拟合伽马分布可以这么来做——
x.gam<-rgamma(200,rate=0.5,shape=3.5) library(stats4) ## loading package
stats4 ll<-function(lambda,alfa) {n<-200 x<-x.gam -n*alfa*log(lambda)+n*log(gamma(alfa))-(alfa- 1)*sum(log(x))+lambda*sum(x)} ## -log-likelihood function
est<-mle(minuslog=ll, start=list(lambda=2,alfa=1))
summary(est)
对于常用的分布形式,也可以直接用MASS包里的fitdistr函数直接拟合 如
library(MASS) ## loading package MASS
fitdistr(x.gam,"gamma")