Thanks that worked perfectly. One thing if I may. Is it possible to make the
center dot red and the lines connecting the dots black?
Thanks,
Chris
On Jul 7, 2011, at 5:10 PM, Abhijit Dasgupta, PhD wrote:
> 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.
