2013-07-02 52 views
0

我注意到,如果我使用家庭cumulativevglm,參數將被忽略。請參閱以下代碼VGAM與家庭累計忽略設置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))) 
[email protected] # 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))) 
[email protected] # output is 1, because I set control=vglm.control(maxit=1) 

我希望有人能夠幫助我。

回答

0

只是對於那些想知道誰: 必須設置在vglm功能的麥克斯特參數不使用的控制參數,例如

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)) 
[email protected] # output is 1