[R] glm() function not finding the maximum

[Brett Presnell]

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
http://www.stat.ufl.edu/~presnell/
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._