1

为了拟合我的数据集的几个概率分布，我尝试了几个使用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函数吗？

谢谢，

4

1 回答 1

2

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

fitBinom <- fitdist(data=Y1, dist="binom", fix.arg=list(size=100), start=list(prob=0.35))
qqp(Y1, "binom", size=100, fitBinom$estimate)
于 2017-09-22T19:52:17.383 回答