1

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

我有重复测量的纵向数据。在很多天里，我对每个主题都进行了反复观察。

目标：

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

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

如果我想做“主题之间”，我将如何编写脚本？

这是我的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

感谢您的帮助！

4

0 回答 0