- 阅读权限
- 255
- 威望
- 1 级
- 论坛币
- 132048 个
- 通用积分
- 9219.7014
- 学术水平
- 925 点
- 热心指数
- 1073 点
- 信用等级
- 703 点
- 经验
- 130943 点
- 帖子
- 3353
- 精华
- 1
- 在线时间
- 3498 小时
- 注册时间
- 2009-7-18
- 最后登录
- 2025-5-18
|
- set.seed(1234)
- x <- rnorm(20)
- df <- data.frame(x = x,
- y = x + rnorm(20))
- plot(y ~ x, data = df)
- # model
- mod <- lm(y ~ x, data = df)
- # predicts + interval
- newx <- seq(min(df$x), max(df$x), length.out=100)
- preds <- predict(mod, newdata = data.frame(x=newx),
- interval = 'confidence')
- # plot
- plot(y ~ x, data = df, type = 'n')
- # add fill
- polygon(c(rev(newx), newx), c(rev(preds[ ,3]), preds[ ,2]), col = 'grey80', border = NA)
- # model
- abline(mod)
- # intervals
- lines(newx, preds[ ,3], lty = 'dashed', col = 'red')
- lines(newx, preds[ ,2], lty = 'dashed', col = 'red')
复制代码
|
-
总评分: 学术水平 + 5
热心指数 + 5
信用等级 + 5
查看全部评分
|