允许gls中的相关参数取决于分组因子

Question

下面是我遇到的问题的MWE。我正在使用nlme软件包中的Orthodont数据集，其中包含4个测量数据，包括27个孩子（16个女孩，11个男孩）。为了模拟相关性，我通过指定correlation = corSymm(form = ~1|Subject)来使用非结构化协方差结构。我允许在不同的测量时间内进行非恒定方差，但我也想允许男孩和女孩的方差 - 协方差参数不同（例如，因为男孩的测量值之间的相关性可能高于女孩）。我可以通过指定weights = varIdent(form = ~1|age*Sex)来允许方差参数的异质性，但有人知道如何允许/指定相关参数的异质性吗？

我知道通过在proc mixed声明（http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect019.htm）中指定group选项可以在SAS repeated中实现这一点，但我还没有找到办法达成协议这个在R。

非常感谢提前！

library(nlme)
head(Orthodont)
fit <- gls(distance ~ age * Sex, correlation = corSymm(form = ~1|Subject),
           weights = varIdent(form = ~1|age*Sex), data = Orthodont, na.action = na.exclude)

Answer 1

我并不完全确定这是您正在寻找的内容，但请尝试一下：

fit <- gls(distance ~ age * Sex,
           correlation = corSymm(form = ~1|Subject/Sex),
           weights = varIdent(form = ~1|age*Sex), 
           data = Orthodont, na.action = na.exclude)

它适合男性和女性的完全独立的日常关联。如果你运行

summary(fit)

它非常清楚地表明了这一点。