r - lme 中的协方差结构 - AR(1)

Question

我的响应变量是 Yijk 对应的recovery时间

• patient我（我=1，...，我）
• 与treatmentj (j=1,...,J)
• 并在timek (k=1,...,K)处测量

我想拟合以下模型：模型方程，其中：

• μ 是全局固定截距
• αj 是治疗的固定效应
• bik 是具有以下协方差结构的随机效应。表示患者 i 的 K 维效应向量，则其方差-协方差矩阵将具有以下 AR(1) 结构。 方差协方差矩阵
• uijk 是方差为 σ² 的通常误差项

考虑以下命令：

lme(recovery ~ treatment, method="REML", random=~1|patient, correlation=corAR1,form=~time|patient,data=data)

几个问题：

1. 这个correlation论点对应的是什么？协方差的结构是什么？那是我定义为 R 的 var-cov 矩阵吗？
2. 这条线真的做我想做的事吗？
3. 如果不是，它有什么作用？
4. 如果没有，有没有办法做我想做的事？

先感谢您！

Answer 1

首先，你有一个命令lme，我假设这是nlme因为 a)lme在我知道的任何包中或者 R 可以找到的包中都不是 R 命令，并且 b) 相关性不是lme4

其次，在nlme他们的文档中：

描述组内相关结构的可选 corStruct 对象。有关可用的 corStruct 类的描述，请参阅 corClasses 的文档。默认为 NULL，对应于没有组内相关性。

里面corClasses说

corAR1 1 阶自回归过程。

因此，前两个问题的答案似乎是“是”。