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 at 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.