小编Aur*_*lie的帖子

我一直在尝试将重复测量模型从 SAS 转换为 R，因为合作者将进行分析但没有 SAS。我们处理 4 组，每组 8 到 10 只动物，然后每只动物有 5 个时间点。模拟数据文件可在此处获得https://drive.google.com/file/d/0B-WfycVUQyhaVGU2MUpuQkg4Mk0/edit?usp=sharing作为 Rdata 文件和此处https://drive.google.com/file/d/ 0B-WfycVUQyhaR0JtZ0V4VjRkTk0/edit?usp=作为excel文件共享：

原始 SAS 代码 (1) 是：

proc mixed data=essai.data_test method=reml;
    class group time mice;
    model param = group time group*time / ddfm=kr;
    repeated time / type=un subject=mice group=group;
run;

这使 ：

    Type 3 Tests des effets fixes
               DDL     DDL     Valeur
Effet         Num.    Res.          F    Pr > F
group            3    15.8       1.58    0.2344
time             4    25.2      10.11    <.0001
group*time      12 …