R mixdist包中的泊松分布

Question

我认为我的数据可能由Poisson分布的总和描述，我发现了R的mixdist包。我设法拟合了gamma和lnorm分布，但我不知道如何使用Poisson 。我首先尝试了示例数据：

    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参数。