「回帰の予測区間のプロット」の編集履歴(バックアップ)一覧はこちら
「回帰の予測区間のプロット」(2008/03/22 (土) 16:19:55) の最新版変更点
追加された行は緑色になります。
削除された行は赤色になります。
#--- サンプルデータ作成 ---#
set.seed(2801)
data <- data.frame(x=rep(seq(5), rep(6,5)))
data$y <- data$x * 10 + 5 + rnorm(nrow(data), sd=10)
#--- 予測したい x の範囲 ---#
x.pred <- data.frame(x=seq(1, 6, length=101))
#--- 直線回帰 ---#
fm1 <- lm(y ~ x, data)
#--- 信頼区間 ---#
fm1.conf <- predict(fm1, x.pred, interval="confidence") # 平均値
fm1.pred <- predict(fm1, x.pred, interval="prediction") # 個別
#--- プロット ---#
matplot(x.pred, fm1.pred, lty=c(1,2,2), type="l", col="black",
xlab="x", ylab="y")
points(data$x, data$y)
lines(x.pred$x, fm1.conf[,"lwr"], lty=2, col="blue")
lines(x.pred$x, fm1.conf[,"upr"], lty=2, col="blue")
*参照
[[R]]
表示オプション
横に並べて表示:
変化行の前後のみ表示: