Equivalent of covariance pattern model with a heterogeneous compound symmetric
Hi Pardis, Four say 4 occasions you could use: HCS.gls <- gls (read ~ 1+occ2+occ3+occ4, data=DA68, method="REML", correlation = corCompSymm(form = ~ occasion|id ), weights = varIdent (form = ~1 | occasion) ) Best regards, Ben.
On 30-5-2019 4:57, Pardis Miri wrote:
Hi everyone,
I am trying to figure out how to replicate the results I've gotten from SPSS when using the Repeated Covariance Type : Compound Symmetry: Heterogeneous. Can anyone help me figure this out?
So far I know the following does't produce the same results.
comsym <- gls(SC_slope~Bodysite+Pattern+Strength,data = df.na.removed.,
correlation=corCompSymm(form = ~ 1 |id), method="REML", na.action = na.omit)
summary(comsym)
[[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models