首页 文章

如何用lmer混合模型表示配对观察

提问于
浏览 次
0

我是lmer的新手，我正在尝试生成一个混合模型，我在其中测试以下效果：

  • 组（固定）：2级

  • 治疗（固定）：2级（未刺激和刺激）

  • 组*治疗

关于因变量“结果”，考虑“主体”的随机效应 .

在该实验中，两组中的每个受试者具有一个受刺激的臂和一个未刺激的受试者 .

到目前为止，我提出的模型是

lmer(Outcome ~ Group + Treatment + Group*Treatment + (1|Subject), REML=FALSE, data= data)

但是，我不确定如何指定每个受试者的一只手臂未受刺激而一只受刺激 .

有人可以帮忙吗？

r statistics lme4 mixed-models bioconductor

1 回答

  • 0

    如果您的问题更多地是针对您的案例的适当模型规范，我会说这取决于您的学习和您的目标 . 你所描述的是符合你的公式，它是有意义的 . 您已经使用（1 |学生）计算学生效果，而治疗指定了治疗组和未治疗组 . 我建议你看看this post which discusses fixed and mixed effects

    关于使用公式在lmer中指定模型的方法，我的第一个评论是以下3个是等价的：

    • 结果〜团体治疗组*治疗

    • 结果〜团体治疗组：治疗

    • 结果〜组*治疗

    第三种是第二种的紧凑形式，第一种是多余的 . 然后我建议您尝试以下有效的替代方案，以便您更熟悉公式表示法

    model2 <- lmer(Outcome ~ Treatment +(1+Treatment|Group)+(1|Subject), REML=FALSE, data= data);coef(summary(model2));ranef(model2)
model3 <- lmer(Outcome ~ Treatment +(0+Treatment|Group)+(1|Subject), REML=FALSE, data= data);coef(summary(model3));ranef(model3)
model4 <- lmer(Outcome ~ Treatment +(1|Group)+(1|Subject), REML=FALSE, data= data);coef(summary(model4));ranef(model4)
    回复于

相关问题