[R] QQ plotting of various distributions...
Petar Milin
pmilin at ff.uns.ac.rs
Thu Sep 17 17:46:10 CEST 2009
Hello!
I am trying with this question again:
I would like to test few distributional assumptions for some behavioral
response data. There are few theories about true distribution of those
data, like: normal, lognormal, gamma, ex-Gaussian
(exponential-Gaussian), Wald (inverse Gaussian) etc. The best way would
be via qq-plot, to show to students differences. First two are trivial:
qqnorm(dat$X)
qqnorm(log(dat$X))
Then, things are getting more "hairy". I am not sure how to make plots
for the rest. I tried gamma with:
qqmath(~ X, data=dat, distribution=function(X)
qgamma(X, shape, scale))
Which should be the same as:
plot(qgamma(ppoints(dat$X), shape, scale), sort(dat$X))
Shape and scale parameters I got via mhsmm package that has gammafit()
for shape and scale parameters estimation.
Am I on right track? Does anyone know how to plot the rest: ex-Gaussian
(exponential-Gaussian), Wald (inverse Gaussian)?
Thanks,
PM
More information about the R-help
mailing list