[R] 95% Confidence Interval for a p-p plot

2013-04-01 Thread pablo.castano
Hi, I want to create upper and lower 95% confidence intervals for a p-p plot of an empirical distribution with a theoretical gamma distribution. This is my code: x<-rgamma(100,shape=2, rate=1) # empirical data fitdistr(x,"gamma") # fit a gamma distribution dist<-pgamma(x,shape=1.9884256 ,rate=0

[R] label the mean in bwplot and divide a panel in two plots

2013-03-08 Thread pablo.castano
Hi, I have the following bwplots, and want to add a label with the mean value at each dote. Can you give me suggestions? Also, if I wanted to add another bwplot below the first, can I use the par(mfrow=c(2,1)) function? This is my code: Code #Data x<-c('Small','Large','Large','Medium','Medium