连续变量是随机效应吗？(R中的lme函数)

Question

荷尔蒙水平被样本质量夸大，即使在通过样本质量校正荷尔蒙水平之后(这是内分泌学家的一个常见问题)。

我在试着确定治疗是否会影响荷尔蒙水平，“修正”样本质量。lme(hormone levels ~ treatment, random= list(~1|INDIVIDUAL, ~1|sample mass), na.action="na.omit", method = "ML",dados)

然而，评论者说我不能使用连续变量作为随机效果。另一种选择是什么？

Answer 1

欢迎使用堆栈溢出。这个问题可能更适合交叉验证，因为它更多的是关于统计而不是编码。不管怎样，我还是要回答。

评论者是正确的，你不能有一个连续的预测者作为一个随机的效果。在这里看一些关于这方面的讨论：https://stats.stackexchange.com/questions/105698/how-do-i-enter-a-continuous-variable-as-a-random-effect-in-a-linear-mixed-effect

为了直接回答您的问题，另一种方法是将预测sample mass作为固定效果添加到模型中，其中它将是协变量。这意味着该模型将考虑荷尔蒙水平如何随大小和治疗方法的变化。这就是user63230的建议，我认为这是前进的正确方式。如果您对sample mass如何改变treatment有一个特定的预测，您可以包括一个交互。如果您希望治疗对不同的个人产生不同的影响，您可以为个人的治疗设置一个随机斜率。