Question

为了为我的数据集拟合几个概率分布，我尝试使用car和MASS包进行多次分发，例如Poisson分布：

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。您只需要为二项分布添加所有参数。您还需要替换＆＃34;二项式＆＃34;用＆＃34; binom＆＃34;

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

用汽车包装拟合二项分布

1 个答案: