相关疑难解决方法(0)

我想知道如何使用broom包计算置信区间.

我想做的是简单和标准:

set.seed(1)
x <- runif(50)
y <- 2.5 + (3 * x) + rnorm(50, mean = 2.5, sd = 2)
dat <- data.frame(x = x, y = y)
mod <- lm(y ~ x, data = dat)

使用visreg我可以CI非常简单地绘制回归模型:

library(visreg)
visreg(mod, 'x',  overlay=TRUE) 

我很有兴趣使用broom和重现这个ggplot2,到目前为止我只实现了这个:

 library(broom) 

 dt = lm(y ~ x, data = dat) %>% augment(conf.int = TRUE)  
 ggplot(data = dt, aes(x, y, colour = y)) + 
  geom_point() + geom_line(data = dt, …