R中的重复测量Anova（混合模型？）

Question

我有一个数据集，我试图在JMP和R中分析，包含4列和468行。

    ID     Treatment    Time     Response
H1                 T      10             28
H1                 T      20             0
H1                 T      30             0
H1                 T      40             0
H1                 Y      10             12
H1                 Y      20             5
H1                 Y      30             24
H1                 Y      40             28
H2                 T      10             4
H2                 T      20             12
H2                 T      30             4
H2                 T      40             0
H2                 Y      10             44
H2                 Y      20             12
H2                 Y      30             23
H2                 Y      40             11

我的完整数据集包含13个受试者（由ID表示），每个受试者接受两次治疗并且被测量18次。我试图找出对这些数据进行重复测量分析的最佳方法。

我正在使用公式am1 <- aov(Response ~ Treatment + (ID/Time), data = df)，但我的同事出现了我的自由度，所以我请求你帮助解决这个问题的最佳方法。我将我的数据设置为单变量方法我相信，但任何帮助都表示赞赏。感谢

Answer 1

假设您的数据是平衡的，您可以在公式中指定Error项，以下是来自？aov的示例：

 npk.aovE <- aov(yield ~  N*P*K + Error(block), npk)