用汽车包拟合二项分布

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)

无法使用 car 包的 qqp 功能。还有其他选项可以拟合二项分布，例如 qqp function?

谢谢，

Answer 1

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)

用汽车包拟合二项分布

Fitting Binomial distribution with car package

r

r-car

probability-distribution