r - 您如何在 R 中拟合具有 AR(1) 随机效应相关结构的线性混合模型？

Question

我正在尝试使用 R 重新运行其他人的项目，因此我们需要在 R 中使用一些宏。

这里有一个非常基本的问题：

m1.nlme = lme(log.bp.dia ~ M25.9to9.ma5iqr + temp.c.9to9.ma4iqr + o3.ma5iqr + sea_spring + sea_summer + sea_fall + BMI + male + age_ini, data=barbara.1.clean, random = ~ 1|study_id)

由于该模型在 SAS 中使用 AR(1) [自相关 1 协方差模型] 进行人内方差，因此我不确定如何在 R 中执行此操作。

我在哪里可以看到不同模型的索引，比如非结构化？

谢谢

Answer 1

我不知道不同模型的“索引”是什么意思，但是要为残差指定 AR(1) 协方差结构，您可以添加corr=corAR1()到 lme 调用中。

Answer 2

滞后 $1$ 处的相关性为 $r$，其中 $-1< r <1$ 对于固定 $AR(1)$ 模型。滞后 $k \geq 1$ 处的相关性为 $r^k$。这只需乘以 $X_t$ 的方差即可得到自协方差矩阵。