In [1]:
library(ggplot2)
In [7]:
dmod <- lm(price ~ cut, data=diamonds)
cuts <- data.frame(cut=unique(diamonds$cut), predict(dmod, data.frame(cut = unique(diamonds$cut)), se=TRUE)[c("fit","se.fit")])
se <- ggplot(cuts, aes(cut, fit,
ymin = fit - se.fit, ymax=fit + se.fit, colour = cut))
se + geom_linerange() + geom_pointrange() + geom_errorbar(width=0.5) + geom_crossbar(width=0.5)
Out[7]:
In [ ]: