「回帰の予測区間のプロット」の編集履歴(バックアップ)一覧はこちら

回帰の予測区間のプロット」(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]]

表示オプション

横に並べて表示:
変化行の前後のみ表示:
記事メニュー
目安箱バナー