如何解决复制交叉设计中的线性混合模型
我正在努力研究如何使用 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
很明显，我可能对原始模型规范有什么误解或误解。
我可能建议您在 r-sig-mixed-models@r-project.org 邮件列表中尝试后续问题，那里有广泛的读者群和广泛的专业知识（即，在混合模型主题方面的专业知识比 StackOverflow 上的更多）
版权声明：本文内容由互联网用户自发贡献，该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容， 请发送邮件至 dio@foxmail.com 举报，一经查实，本站将立刻删除。