使用fitdistplus的fitdist和二项分布

时间:2014-12-23 21:50:26

标签: r distribution

我刚刚发现了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没有)。

有什么想法?

1 个答案:

答案 0 :(得分:5)

您是否总是不知道试验次数(即尺寸参数)?如果是,请尝试

fitBinom=fitdist(data=scorebinom, dist="binom", fix.arg=list(size=8), start=list(prob=0.3))

估计p及其误差。