{求助} 线性回归 ± 置信区间（95%），且置信区间是两条直线形式

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助大神，用R 进行 截距为0的线性回归，如何能作出图中的两条直线形式的置信区间？万分感谢~~~

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：置信区间 线性回归 求助大神 万分感谢 如何 置信区间

1. x <- runif(50, 0, 25)
   
2. y <- 40*x + 200*rnorm(x)
   
3. 
   
4. fit <- lm(y ~ 0 + x)
   
5. newx <- seq(min(x), max(x), length = 200)
   
6. pred <- predict(fit, newdata = data.frame(x = newx), interval = c("confidence"), level = 0.95)
   
7. plot(x, y, type = "n")
   
8. polygon(c(rev(newx), newx), c(rev(pred[ ,3]), pred[ ,2]), col = 'grey', border = NA)
   
9. points(x, y)
   
10. abline(fit, lwd = 2)
    
11. lines(newx, pred[,2], col = "red", lty = 2)
    
12. lines(newx, pred[,3], col = "red", lty = 2)

复制代码

qoiqpwqr 发表于 2016-4-25 10:06

很接近了！这两条线应该不是截距为0的，我这个地方郁闷了很久，搞不清楚是咋出来的~~