0

cumulativeで family を使用すると、 invglmのmaxit引数vglm.controlが無視されることに気付きました。次のコードを参照してください

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)

最後の行は、フィッティングに使用された反復回数を示しています。maxit = 1に設定しても、出力は4です。誰かがこれを修正する方法を知っていますか?

別のファミリ (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)

誰かが私を助けてくれることを願っています。

4

1 に答える 1

0

知りたい人向け: たとえば、control 引数を使用せずに 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
于 2013-07-16T11:08:38.897 に答える