我刚刚发现了fitdistrplus
包,我把它和Poisson发行版一起运行等等。但是在尝试使用二项式时我遇到了问题:
set.seed(20)
#Binomial distributed, mean score of 2
scorebinom <- rbinom(n=40,size=8,prob=.25)
fitBinom=fitdist(data=scorebinom, dist="binom", start=list(size=8, prob=mean(scorebinom)/8))
我收到错误:
Error in fitdist(data = scorebinom, dist = "binom", start = list(size = 8, :
the function mle failed to estimate the parameters,
with the error code 100
In addition: There were 50 or more warnings (use warnings() to see the first 50)
我从这个软件包中看到了很多关于负二项分布的文档,但关于二项式的分析并不多。此功能似乎支持此分发(尽管fitdistr
中的MASS
没有)。
有什么想法?
答案 0 :(得分:5)
您是否总是不知道试验次数(即尺寸参数)?如果是,请尝试
fitBinom=fitdist(data=scorebinom, dist="binom", fix.arg=list(size=8), start=list(prob=0.3))
估计p及其误差。