Displaying 1 result from an estimated 1 matches for "gf_pois".
2010 Mar 02
1
test the goodness of it for negative binomial type 2
[code]library(MASS)
x=c(rep(0,8096),
rep(1,1629),
rep(2,233),
rep(3,38),
rep(4,4)
)
x.bar=round(mean(x),4)
x.var=round(var(x),4)
p.hat=round(x.bar/x.var,4)
alpha.hat=round(x.bar*p.hat/(1-p.hat),4)
fitdistr(x, "Negative Binomial")
fitdistr(x, "Poisson")[/code]
1- fitdistr(x, "Negative Binomial")
the parameters got here,