Я пытаюсь подогнать мою модель временных рядов Poission в garmaFit (пакет gamlss), но мои параметры не работают.
model6.0 <-garmaFit(aces ~ time+ smokban + offset(log(stdpop)), data=data, order=c(0,0),phi.start =11 ,theta.start = 11, family=PO)
Это запускает модель, но не регулирует ее на основе коэффициентов, которые я могу сгенерировать из
model1 <- glm(aces ~ offset(log(stdpop)) + smokban + time, family=poisson, data)
Я даже попытался использовать известный метод для контроля автокорреляции в модели ols ( из другого набора данных), и он не корректировал коэффициенты из исходного режима.
Не могу найти в другом месте, как настроить значения p / q в параметрах phi.start и theta.start
Любая помощь будет принята с благодарностью