R/gammaMME.R

Defines functions gammaMME

Documented in gammaMME

## the canonical gamma MME -- but the cMLE is lower variance and just as fast!
gammaMME <- function(x) { 
 return(c(shape=(mean(as.matrix(x),na.rm=T)/apply(as.matrix(x),2,sd,na.rm=T))^2,
          scale=apply(as.matrix(x),2,var,na.rm=T)/mean(as.matrix(x),na.rm=T)))
} 

Try the rGammaGamma package in your browser

Any scripts or data that you put into this service are public.

rGammaGamma documentation built on May 29, 2017, 8:01 p.m.