汽车包装二项分布的拟合

Question

为了为我的数据集拟合几种概率分布，我尝试了使用car和MASS包的几种分布，例如对于泊松分布：

require(car); require (MASS)

#Data set
X1 <- rnorm(100, 175, 7)
Y1 <- rbinom(100,1, 0.35)
##Fitting Poisson:
poisson <- fitdistr(Y1, "Poisson")
qqp(Y1, "pois", poisson$estimate)

但是，如果我的选项是二项式分发，MASS包不支持这种分发，那么我使用fitdistrplus包：

require(fitdistrplus)
fitBinom=fitdist(data=Y1, dist="binom", fix.arg=list(size=100), start=list(prob=0.35))

当我尝试：

qqp(Y1, "Binomial", fitBinom$estimate)

不工作的qqp函数的car包。还有其他方法来拟合二项分布，如qqp函数吗？

谢谢,

Answer 1

qqp实际上适用于fitdistrplus。您只需要为二项分布添加所有参数。您还需要将“二项式”替换为"binom“。

fitBinom <- fitdist(data=Y1, dist="binom", fix.arg=list(size=100), start=list(prob=0.35))
qqp(Y1, "binom", size=100, fitBinom$estimate)