Распределение Пуассона в пакете R mixdist

Question

Я думаю, что мои данные могут быть описаны суммой распределений Пуассона, и я узнал о пакете mixdist для R. Мне удалось подогнать распределения гаммы и норм, но я не могу понять, как использовать Пуассона. Сначала я попробовал данные примера:

    library(mixdist) 
    data(poisdat)
    data(poispar)
    plot.mixdata(poisdat) #this works
    fitp<-mix(poisdat, coef(poispar), "pois") #this doesn't

, но я получаю сообщение об ошибке Error in if (usecondit & ncol(mixdat) - 2 != k) stop("Conditional data are not consistent with mixpar.") : argument is of length zero

Как получить рабочий пример mixdist и Poisson?

Answer 1

Вы ищете это:

mix(poisdat, poispar, "pois", constr = mixconstr(consigma = "POIS"))

Просто mix(poisdat, poispar, "pois") выдает ошибку:

Ошибка в testconstr (mixdat, mixpar, dist, constr): нужды Пуассонаconsigma = POIS.

, поэтому этот constr параметр.