R

如何決定將 REML 設置為 True 還是 False?

  • April 1, 2017

我找到了一個網頁告訴 lmer:

如果您的隨機效應是嵌套的,或者您只有一個隨機效應,並且您的數據是平衡的(即每個因子組中的樣本量相似),請將 REML 設置為 FALSE,因為您可以使用最大似然。如果您的隨機效應被交叉,請不要設置 REML 參數,因為無論如何它默認為 TRUE。

我的lmer模型中有 2 個隨機效應。一種是嵌套的:

(1|Random1A/Random1B) + (1|Random2)

我應該設置REMLTRUE(默認情況下)還是FALSE

在我(並非完全不知情)看來,您從網頁和您收到的評論中得到了一些有問題的建議。

  • 您可以隨時使用 REML(或 ML)(無論隨機效應結構如何 - 單個與多個、平衡與不平衡、交叉與嵌套)
  • 在簡單的情況下(平衡/嵌套/等),REML 可以被證明可以提供方差分量的無偏估計(但不能提供例如標準偏差或對數標準偏差的無偏估計)
  • 如果使用 REML 而不是 ML 擬合,則無法比較固定效應不同的模型;這就是為什麼評論者建議REML=FALSE您在嘗試進行模型選擇時使用
  • 但是,我不建議您首先進行模型選擇,當然如果您要依賴條件置信區間和 p 值(即,在不考慮影響的情況下分析改裝的“最小充分”模型)型號選擇)

從我在Fox et al 2015的章節中:

當您對隨機效應方差的大小感興趣時,通常最好使用 REML(如果可用),但當您通過假設檢驗或信息論標準(如 AIC)比較具有不同固定效應的模型時,則永遠不要使用它。

引用自:https://stats.stackexchange.com/questions/272633

comments powered by Disqus