我有一个包含 12 天日记数据的数据集。我试图用lme参与者的随机截距效应和睡眠质量的随机斜率效应来模拟睡眠质量对压力的影响。我对询问从日记 1 到 12 是否随时间变化并不特别感兴趣,只是考虑了时间变量。
我对lme4多级模型更加熟悉,但考虑到我的日记数据结构,考虑方差中的时间自相关似乎很重要。
我已经尝试从这里遵循模型语法,但我并不完全清楚语法的实际含义。
我的问题:具体的语法是correlation = corAR1(form = ....)什么?corAR1我在以下两个模型中编写的语法有什么区别?
型号 1:
lme(stress_diary ~ sleepquality_diary + time.diary + diaryday + sex,
data = data_diary,
random = ~ 1 + sleepquality_diary| randomid,
correlation = corAR1(form = ~ diaryday | randomid), #What exactly am I specifying?
na.action=na.exclude)
型号 2:
lme(stress_diary ~ sleepquality_diary + time.diary + diaryday + sex,
data = data_diary,
random = ~ 1 + sleepquality_diary| randomid,
correlation = corAR1(), #What exactly am I specifying?
na.action=na.exclude)
变量键:stress_diary = 1 到 5 之间的值 sleepquality_diary = 1 到 5 之间的值 time.diary = 响应时间,从 00:00(午夜)开始 diaryday = 响应日,1:12 性别 = M,F randomid = 参与者
我对时间序列和纵向日记数据很陌生,所以我知道我缺少一些直觉。