如何拟合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)

由于模型是利用AR （1）SAS中的[自相关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 $的方差给出了自协方差矩阵。