0

我使用 R glm 或Countr (renewalCount) 得到的可能性与我使用 R 估计参数手动计算可能性不同。我试图用 R 中的一些函数来拟合数据集,但似然值与我从其他语言中得到的不一致。另一个奇怪的事情是,对于不同的数据集,我使用 glm(...,family=poisson) 得到了相同的结果!

D = matrix(c(1:9))
Numb = matrix(c(10, 4, 5, 3, 1, 2, 1, 1, 1))
Number = data.frame(Numb)
agg_poiss <- glm(formula = D ~ 1, family = poisson, data = Number)
agg_wei <- renewalCount(formula = D ~ 1, data = Number, dist = "weibull", 
computeHessian = FALSE, control = renewal.control(trace = 0))

对于泊松情况,我得到:(截距)1.6094,这意味着 exp(1.6094) 的泊松参数,但 t 给了我 R 中的可能性 -21.53659,当我手动计算时它会变得不同。

4

0 回答 0