Dear Augusto,
thank you for your example. Your solutions are the in fact the usual
methods, but these do not apply to my case because I have grouped data.
The good news is, that the solution of Prof. Brian Ripley works
perfectly -- of course :-)
Thank you both for your help.
Thomas P.
[EMAIL PROTECTED] wrote:
> Thomas,
>
> The Gamma distr. can be fitted via ML using:
>
> Library("MASS")
> GF <- fitdistr(given_data,"gamma")
> sh <- GF$estimate[1]
> ra <- GF$estimate[2]
>
> Fitting via Moments, m:
> var <- m[2] - m[1]*m[1]
> sh <- m[1]*m[1]/var
> sc <- m[1]/var
> ra <- 1/sc
>
> G_pdf <- dgamma(breaks,shape=sh,rate=ra,scale=1/ra)
>
>
> Hope it helps,
>
> Augusto
[...] original message deleted
______________________________________________
R-help@stat.math.ethz.ch 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.
