我正在尝试将我的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模型中的自相关(来自不同的数据集),并且它没有从原始模式调整系数。
在其他地方找不到如何调整phi.start和theta.start参数中的p / q值
任何帮助将不胜感激