You can easily do this by:
qplot(x=as.factor(sch),y=est, geom='point', colour='red') +
geom_pointrange(aes(x=as.factor(sch), y=est, ymin=lower.95ci, ymax=upper.95ci))+
xlab('School') + ylab("Value-added")+theme_bw()
On 07/07/2011 05:55 PM, Christopher Desjardins wrote:
Hi,
I have the following data:
est
sch190 sch107 sch290 sch256 sch287 sch130
sch139
4.16656026 2.64306071 4.22579866 6.12024789 4.49624748 11.12799127
1.17353917
sch140 sch282 sch161 sch193 sch156 sch288
sch352
3.48197696 -0.29659410 -1.99194986 10.23489859 7.77342138 6.77624539
9.66795001
sch368 sch225 sch301 sch105 sch353 sch291
sch179
7.20229569 4.41989204 5.61586860 5.99460203 -2.65019242 -9.42614560
-0.25874193
sch134 sch135 sch324 sch360 bb1
3.26432479 10.52555091 -0.09637968 2.49668858 -3.24173545
se
sch190 sch107 sch290 sch256 sch287 sch130 sch139 sch140
3.165127 3.710750 4.680911 6.335386 3.896302 4.907679 4.426284 4.266303
sch282 sch161 sch193 sch156 sch288 sch352 sch368 sch225
3.303747 4.550193 3.995261 5.787374 5.017278 7.820763 7.253183 4.483988
sch301 sch105 sch353 sch291 sch179 sch134 sch135 sch324
4.076570 7.564359 10.456522 5.705474 4.247927 5.671536 10.567093 4.138356
sch360 bb1
4.943779 1.935142
sch
[1] "190" "107" "290" "256" "287" "130" "139" "140" "282" "161" "193" "156"
"288"
[14] "352" "368" "225" "301" "105" "353" "291" "179" "134" "135" "324" "360"
"BB"
From this data I have created 95% confidence intervals assuming a normal
distribution.
lower.95ci<- est - se*qnorm(.975)
upper.95ci<- est + se*qnorm(.975)
What I'd like to do is plot the estimate (est) and have lines attach to the
points located in lower.95ci and upper.95ci. Presently I am doing the
following:
qplot(x=as.factor(sch),y=lower.95ci) + geom_point(aes(x=as.factor(sch),y=upper.95ci),colour="black") +
geom_point(aes(x=as.factor(sch), y=est),colour="red") + ylab("Value-Added") +
xlab("School") + theme_bw()
Which creates this graph --->
http://dl.dropbox.com/u/1501309/value_added_test.pdf
That's fine except that it doesn't connect the points vertically. Does anyone
know how I could make the 'black' points connect to the 'red' point, i.e. show
confidence bands?
Thanks,
Chris
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.