8

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

我知道通过在语句中proc mixed指定group选项（ http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect019.htm），这在 SAS 中是可能的，但我还没有找到在 R 中处理这个问题的方法。repeated

提前谢谢了！

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)
4

1 回答 1

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)

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

于 2015-03-31T17:12:07.480 回答