R 如何在lme中分别拟合两个随机效应?
我正在用nlme软件包中的REML进行线性混合效果模型拟合。这些代码对我来说很有用:R 如何在lme中分别拟合两个随机效应?,r,random,nlme,R,Random,Nlme,我正在用nlme软件包中的REML进行线性混合效果模型拟合。这些代码对我来说很有用: # Linear mixed-effects model fit by REML (intercept and not slope) x <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker) summary(x) # Linear mixed-effects model fit by REML (slope and no i
# Linear mixed-effects model fit by REML (intercept and not slope)
x <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker)
summary(x)
# Linear mixed-effects model fit by REML (slope and no intercept)
x1 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3-1|speaker)
summary(x1)
# Linear mixed-effects model fit by REML (slope and intercept)
x2 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~IV3|speaker)
summary(x2)
#nested random effect
x5 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~1|speaker/item)
summary(x5)
你知道这是什么意思吗?如何将扬声器和项目分别作为随机效果进行匹配?我认为可以通过以下代码使用
lme()
分别包含两个随机效果(一个用于扬声器,一个用于时间):
x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~ speaker + item -1 | id),
x4使用包lme4。它为您定义模型的随机部分提供了更大的灵活性,但是我不能直接从该包中获取p值,因为我使用的是新版本的R,对吗?我以前尝试过该软件包,但无法使用“pval”功能这个软件包很方便,因为它直接给我p值我尝试过以下公式:x如果您将lme4
与lmerTest
软件包结合使用,您可以指定模型并获得所需的p值(通过Satterthwaite或Kenward-Roger近似)。无论怎样,lme
都很有可能不会给出正确的p值……谢谢Ben Bolker。现在我似乎不必再使用lmerTest了,因为我在使用lsmeans的LMM之后进行两两比较,这直接给了我p值。
Warning message:
In Ops.factor(speaker, 1) : + not meaningful for factors
x4 <- lme (DV ~ IV1 + IV2 + IV1*IV2, data=a.frame, random=~ speaker + item -1 | id),