使用順序包裝中的clmm2難以獲得擬合和置信區間

Question

我在R中使用序數包中的clmm函數，以使累積混合模型適合我的數據。 在我嘗試獲得預測的概率之前，它一直運行良好。 通過指定se.fit = TRUE和interval = TRUE，我無法獲得SE或置信區間。 看起來像這樣：

mod1<-clmm2(response~X0+X1+X2+X3+X4+X5+X7+X0*X2*X3+X2*X3*X4+X0:X4, random=X6, 
                      data=df,link ="logistic", threshold ="flexible", 
                      Hess=TRUE, nAGQ=7)

如您所見，那里有很多互動（都很重要）。 我試圖為我的問題創建一個可重復的虛擬數據集，但clmm無法與更簡單的數據集實現收斂。 我以序數形式獲取了包裝中包含的葡萄酒數據集，並對公式進行了一些更改以模仿我自己的（雖然我認為這沒有任何意義）：

library(ordinal)
data(wine)
fm1 <- clmm2(rating ~ temp + contact+bottle+temp:contact:bottle+temp:contact+ temp:bottle+bottle:contact,random=judge, data=wine,link ="logistic", threshold ="flexible", 
         Hess=TRUE, nAGQ=7)

head(do.call("cbind", predict(fm1, se.fit=TRUE, interval=TRUE)))

然后我得到這個錯誤：

Error in head(do.call("cbind", predict(fm1, se.fit = TRUE, interval = TRUE))) : 
 error in evaluating the argument 'x' in selecting a method for function 'head' : Erreur dans do.call("cbind", predict(fm1, se.fit = TRUE, interval = TRUE)) : second argument must be a list

我的猜測是，在這種情況下，預測甚至不會計算SE和IC。 有人知道為什么嗎？ 無論如何，要獲得這些價值？

非常感謝！

Answer 1

clmm2對象的預測方法不提供標准錯誤。 請參閱其幫助頁面。 這與R包作者處理混合效果模型時的慣常做法保持一致。