如何在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 $，其中對於靜態$ AR（1）$模型，$-1 <r <1 $。 滯后$ k \\ geq 1 $的相關值為$ r ^ k $。 通過乘以$ X_t $的方差，可以得到自協方差矩陣。