Вероятность того, что я получу использование 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 дает мне вероятность -21.53659 в R, которая становится другой, когда я вычисляю вручную.