信頼区間で予測をプロットするには少し助けが必要です。私は他の全ての変数は(/モードを意味する)一定に保持されているときvs
(0と1)の効果をプロットしたい次の例線形モデル予測の95%信頼区間をプロットするときにプロットが正しくない
library(Hmisc)
data("mtcars")
mfit = lm(mpg ~ vs + disp + cyl, data = mtcars)
#disp and cyl at their mean
newcar = data.frame(vs = c(0,1), disp = 230, cyl = 6.188)
pmodel <- predict(mfit, newcar, se.fit=TRUE)
を考えてみましょう。私はここで間違ってやっている何
plot(1:2, pmodel$fit[1:2], ylim=c(0,1), pch=19, xlim=c(.5,2.5), xlab="X",
ylab = "Predicted values", xaxt = "n", main = "Figure1")
arrows(1:2, (pmodel$fit[1:2] - 1.96 * pmodel$fit[1:2]),
1:2, (pmodel$fit[1,1] + 1.96 * pmodel$fit[1:2]),
length=0.05, angle=90, code=3)
axis(1, at=c(1,2), labels=c("0","1"))
:
は、私は以下のこのコードを実行するこれを行うには?ありがとう!
: 'ライブラリ(visreg)。 visreg(mfit、xvar = "vs") ' – agenis