Pablo, Check the qqPlot function in "car":
require(car) qqPlot(x, dist = "gamma", shape = 1.7918012, rate = 0.9458022) Best, Jorge.- On Tue, Apr 2, 2013 at 4:41 AM, pablo.castano <> wrote: > 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.8765314 ) # fitted distribution, > using the loglikelihood estimated parameters plot(ppoints(n=100), > sort(dist)) # create p-p plot > abline(0,1) # diagonal line to check fit > > Is there an implementation of the delta method in R that I could use for > this in order to estimate the variance of the predicted probability ? > > > Regards > > > > -- > View this message in context: > http://r.789695.n4.nabble.com/95-Confidence-Interval-for-a-p-p-plot-tp4662982.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > 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. > [[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.