使用lme4爲隨機效應

-5

我有3個隨機變量，X，YZ（所有隨機效應）使用lme4爲隨機效應

x被嵌套在Y，但y爲沿z

我使用下面的函數在lme4劃線，但它不起作用。

<- lmer(A ~ 1 + (1 | x/y) + (1 | y*z) + (1|x/y*z), my data)

有人幫我嗎？非常感謝

'my data'不能是數據框的名稱。並從一個更簡單的模型開始，檢查它是否運行。「不起作用」對獲得良好迴應沒有多大用處。顯示整個工作示例，將其降至最低。 – 2014-11-20 19:10:24

我的數據是在lme4中使用的矩陣。這是正確的，我認爲它與其他模型一起工作。當變量與其他變量嵌套時，此函數不起作用。我的問題是使用的功能。我們如何解決函數中的嵌套變量。 – user3546966 2014-11-20 19:58:48

*請註明*「不起作用」是指... – 2014-11-20 21:56:14

恐怕這還是很不清楚。更多的上下文將是有用的。我的猜測是，你想

A ~ 1 + (1|y)+ (1|z) + (1|y:z) + (1|y:x)

或等價

A ~ 1 + (1|y*z) + (1|y:x)

，但它幾乎是不可能知道。

表達A ~ 1 + (1|y/x) + (1|z/y)應該給你同樣的結果，因爲a/b一般擴展到a + a:b（爲了事項/而不是:），但它不太清楚。

交叉隨機效應一般是通過(1|y) + (1|z)來表示，或由(1|y*z)（其擴展到(1|y) + (1|z) + (1|y:z)）如果如上所討論的有每{Y，Z}組合的多個觀測。

2014-11-20 21:30:39

非常感謝你，本，非常有用。在我的模型中，y在z中交叉。我在函數中沒有看到y * z，請您澄清一下。 – user3546966 2014-11-20 21:50:33

謝謝你Ben，隨機效應應該被編碼，比如說1,2,3,4 ...是否正確？ – user3546966 2014-11-20 22:03:55

分組變量（'x'，'y'，'z'）理想情況下應該是*因素*，但是如果它們不是，lme4會將它們轉換。（如果你不知道什麼是因素，那麼閱讀一些介紹性的R材料可能是一個好主意......） – 2014-11-20 22:49:51

回答