如何解决使用 nlme 指定 2x2 重复测量模型
我有来自 2x2 重复测量实验的数据。四个实验条件由 x1 和 x2 虚拟变量定义。 IE,四个条件是:
x1 = 1 & x2 = 1
x1 = 0 & x2 = 0
x1 = 1 & x2 = 0
x1 = 0 & x2 = 1
此外,还有一个代表性别的gend dummy。连续响应变量是 resp。我正在尝试对此进行建模,并且在我读过的一本书(Field)中展示了 nlme 中的 2x2 重复测量,作者使用 x1 和 x2 变量作为嵌套在人(id)中的级别,例如:
m1 <-
lme(resp ~ x1 + x2 + x1:x2 + gend,random = ~gend | id/x1/x2,data = df,method = "ML")
上述模型包括 x1、x2、交互和性别的固定效应。它还包括性别上的随机斜率,都嵌套在人和 4 个实验条件中。这对我来说似乎很奇怪,因为它在一个人体内嵌套了两个实验条件。
为什么不只是这样定义:
m2 <-
lme(resp ~ x1 + x2 + x1:x2 + gend,random = ~gend | id,method = "ML")
IE,为什么嵌套级别 |id/x1/x2 而不仅仅是 |id 并且包含 x1 和 x2 作为固定效果,因为 x1 和 x2 变量是级别 1(它们因条件而不是每个参与者而异)?
第二个问题:是否有关于指定 lme() 公式的明确参考?我没有找到任何小插曲,包手册似乎也没有指定所有不同的方式来配置 lme() 固定和随机公式。
数据样本:
df <- structure(list(
id = c(1,1,2,3,3),x1 = c(0L,1L,0L,0L),x2 = c(1L,1L),gend = c(0L,resp = c(12L,15L,13L,17L,19L,20L,11L,14L,10L,17L)),class = "data.frame",row.names = c(NA,-12L))
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。