如何使用 nlme 包为多级建模编写 2 级数据的代码

问题描述

我正在为如何在 R 中的多级模型中描述级别 2 数据而苦苦挣扎。 我正在使用 nlme 包。

我有重复测量的纵向数据。我对每个科目进行了很多天的重复观察。

目标：

• 级别 1 将是主题 ID 内的个人观察
• 级别 2 将是总体平均值之间的差异主题 ID 之间（集群）。

我正在尝试确定考试成绩是否受学习时间的显着影响，并查看科目内和科目之间是否存在显着差异。

如果我想做“主题之间”，我该怎么写剧本？

这是我的 Level 1 模型的脚本

model1 <- lme(fixed = Testscore~Studytime,random =~1|SubjectID,data=dataframe,na.action=na.omit)

下面是我的示例数据框

      `Subject ID` Observations Testscore Studytime
 1            1            1        50       600
 2            1            2        72       900
 3            1            3        82       627
 4            1            4        90      1000
 5            1            5        81       300
 6            1            6        37       333
 7            2            1        93       900
 8            2            2        97      1000
 9            2            3        99      1200
10            2            4        85       600
11            3            1        92       800
12            3            2        73       900
13            3            3        81      1000
14            3            4        96       980
15            3            5        99      1300
16            4            1        47       600
17            4            2        77       900
18            4            3        85       950

感谢您的帮助！

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）