1

FDA推荐的型号

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

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

1 回答 1

0

我会说公式应该是

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 上的更多）

于 2020-12-21T16:58:52.680 回答