ggplot2 - R - tidy augment confidence interval -


i wondering how can compute confidence interval using broom package.

what trying simple , standard :

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) 

using visreg can plot regression models ci :

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

enter image description here

i interesting in reproducing using broom , ggplot2, far achieved :

 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, aes(x, .fitted, colour = .fitted))  

enter image description here

the augment funciton doesn't compute conf.int. clue how can add smooth confidence invervals ?

 geom_smooth(data=dt, aes(x, y, ymin=lcl, ymax=ucl), size = 1.5,          colour = "red", se = true, stat = "smooth") 

using broom output, can this:

ggplot(data = dt, aes(x, y)) +    geom_ribbon(aes(ymin=.fitted-1.96*.se.fit, ymax=.fitted+1.96*.se.fit), alpha=0.2) +   geom_point(aes(colour = y)) +    geom_line(aes(x, .fitted, colour = .fitted)) +   theme_bw() 

i moved colour=y geom_point() because can't apply colour aesthetic geom_ribbon.

enter image description here


Comments

Popular posts from this blog

aws api gateway - SerializationException in posting new Records via Dynamodb Proxy Service in API -

asp.net - Problems sending emails from forum -