[R] glm() function not finding the maximum

Brett Presnell presnell at stat.ufl.edu
Mon Apr 22 21:46:17 CEST 2002

In message <Pine.GSO.4.44.0204221641030.4377-100000 at moeran> you write:
> #Fits data to a gamma distribution using glm()
> gamma1_function(data){
> n_length(data)
> m_summary(glm(data~1, family=Gamma(link=identity)))
> shape_1/as.numeric(m$disp)

You're not using the mle of the gamma shape parameter.  See the
function shape.gamma in library(MASS).

Brett Presnell
Department of Statistics
University of Florida
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch

More information about the R-help mailing list