复制交叉设计中的线性混合模型

问题描述

我正在努力研究如何使用 REML 方法为复制交叉设计拟合模型。 FDA建议的模型如上所述，有人可以帮助将其编码为R编码吗？这是我的编码，不知道是对还是错？

samplePK.lmer <- lmer(ykir2~1+Treatment:Sequence:Replication+
    (1|Subject:Sequence:Treatment),data=samplePK,REML=TRUE)

解决方法

我会说公式应该是

response ~ trt + trt:seq:rep + (trt|subj:seq) 

与您的规范的主要区别在于，(trt|subj:seq) 拟合了诸如随机区组或随机斜率模型之类的东西，我们允许跨主题/序列组合的治疗效果发生变化。

在尝试将此模型拟合到模拟数据时，我遇到/注意到了一些问题：

• 如果我们使用“现代”混合模型方法对此进行拟合，则在处理效果 (trt) 和 trt:seq:rep 项之间会有一些参数别名。在瞬间方法中，这并不重要（因为您从未明确估计参数），但它会导致对排名不足模型的抱怨（如果您知道自己在做什么，则可以忽略不计...... ).
• 随机效应 delta_{ij} 的平均值为 {mu_R,mu_T} 似乎是错误的；这是多余的固定效应 mu_k

很明显，我可能对原始模型规范有什么误解或误解。

我可能建议您在 [e-mail] 邮件列表中尝试后续问题，那里有广泛的读者群和广泛的专业知识（即，在混合模型主题方面的专业知识比 StackOverflow 上的更多）