具有系列累积的VGAM忽略设置最大值

Question

我注意到，如果我在vglm中使用cumulative家族，那么vglm.control中的maxit参数将被忽略。请参阅以下代码

library(VGAM)
pneumo <- transform(pneumo, let = log(exposure.time))
(fit <- vglm(cbind(normal, mild, severe) ~ let,
             family=cumulative(parallel = TRUE, reverse = TRUE), 
             data=pneumo, control=vglm.control(maxit=1)))
fit@iter # output is 4, even though I set control=vglm.control(maxit=1)

最后一行给出了用于拟合的迭代次数。这里的输出是4，尽管我设置了maxit=1。有人知道如何修复这个问题吗？

如果我使用不同的族(例如poissonff族)，maxit参数将在拟合过程中考虑，请参见：

(fit2 <- vglm(normal ~ let, family=poissonff, data=pneumo, 
             control=vglm.control(maxit=1)))
fit2@iter # output is 1, because I set control=vglm.control(maxit=1)

我希望有人能帮助我。

Answer 1

仅供那些想知道的人使用:您必须在vglm-function中设置maxit参数，而不使用控制参数，例如

library(VGAM)
pneumo <- transform(pneumo, let = log(exposure.time))
(fit <- vglm(cbind(normal, mild, severe) ~ let,
             family=cumulative(parallel = TRUE, reverse = TRUE), 
             data=pneumo, maxit=1))
fit@iter # output is 1