2017-01-27 7 views
0

для монтажа несколько распределения вероятностей для моего набора данных, я стараюсь несколько дистрибутивов с помощью 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) 

Не работает функция qqpcar. Существуют и другие варианты установки биномиального распределения, например, qqp?

Спасибо,

ответ

1

qqp фактически работает с fitdistrplus. Вам нужно только добавить все параметры для биномиального распространения. Вам также необходимо заменить «Биномный» на «бином»

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