2010-03-14 82 views
-2

如何在R中指定隨機因子?如何在R中指定隨機因子?

如果我有一個因子x1這應該是隨機的,我可以嘗試這樣的事情嗎?

lm(y ~ x1, data = p) 
+0

對不起,你想達到什麼目的?你想創建一些隨機因素,還是你想對數據框架的兩個因素做一個線性迴歸?,一些示例數據+代碼將真的幫助... – PaulHurleyuk 2010-03-14 20:31:03

+0

因此,如果我嘗試lme4(y〜x1,data = p),將x1視爲一個隨機因子? – 2010-03-15 18:35:12

+1

我懷疑x1將被視爲「錯誤:無法找到函數」lme4「」,因爲沒有名爲lme4的函數。 – 2010-03-15 19:56:53

回答

4

您偶然尋找lme4包這是對線性混合效應(即固定與隨機)造型重點?

+0

如果我寫lm(y〜x1,data = p), x1將被視爲R中的一個固定因子, 我想知道在R中是否有辦法指定x1是一個隨機因子 – TKBell 2010-03-14 21:28:52

+3

是的,在* * lme4 **包。請按照我的答案中的鏈接閱讀lme4軟件包的文檔。 lm()函數用於基本的線性建模,並且不支持隨機效應。 – 2010-03-14 22:33:50

5

最簡單的(使用舊nlme而不是新lme4),

library(nlme) 
lme(fixed = y ~ 1, random = ~1|x1, data = p) 

lme4相對應的是:

library(lme4) 
lmer(y~1+(1|x1), data = p) 
2

Zuur等2011人,混合效應模型和擴展的生態與R,使用nlme軟件包有很好的隨機效果。他們解釋了固定模型,隨機模型和混合模型之間的差異以及如何指定隨機截距和隨機斜率。