如何在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）[autocorrelation 1协方差模型]，因此我不确定如何在R中进行此操作。

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

由于

Answer 1

我不知道你的意思＆＃34;索引＆＃34;对于不同的模型，但要为残差指定AR（1）协方差结构，您可以将corr=corAR1()添加到您的LME调用中。

Answer 2

滞后$ 1 $的相关性是$ r $，其中$ -1＆lt;对于固定的$ AR（1）$型号，r <1 $。滞后$ k \ geq 1 $的相关性是$ r ^ k $。这只需乘以$ X_t $的方差即可得到自协方差矩阵。