以下は、私が遭遇している問題の MWE です。Orthodont27人の子供（女の子16人、男の子11人）の4つの測定値で構成されるnlmeパッケージのデータセットを使用しています。相関をモデル化するために、 を指定して構造化されていない共分散構造を使用していcorrelation = corSymm(form = ~1|Subject)ます。さまざまな測定機会にわたる非一定の分散を許可していますが、分散共分散パラメーターが男の子と女の子で異なることも許可したいと思います (たとえば、測定間の相関は女の子よりも男の子の方が高い/低いため) ）。を指定することで、分散パラメーターのこの不均一性をweights = varIdent(form = ~1|age*Sex)許可できますが、相関パラメーターのこの不均一性を許可または指定する方法を知っている人はいますか?
ステートメントでオプションをproc mixed指定することにより、SASでこれが可能であることは知っています（ http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_mixed_sect019.htm）が、私はRでこれに対処する方法を見つけていません。grouprepeated
よろしくお願いします！
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)