微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

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

如何解决复制交叉设计中的线性混合模型

model suggested by FDA

我正在努力研究如何使用 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 举报,一经查实,本站将立刻删除。