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

2013-04-01 Thread Jorge I Velez
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 di

[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