Displaying 1 result from an estimated 1 matches for "lmc2".
Did you mean:
lmc
2011 Jun 02
0
allowing individual level correlations to differ by cluster in lme in R
...s and effects
differ by cluster. I have worked out the code to account for the former but not both (see code below).
lme(yl~
-1+cons.1+treat.1+cons.2+treat.2,
random=~-1+cons.1+cons.2|cidl,
weights=varIdent(form=~1|resp.cidl),
corr=corCompSymm(form=~1|cidl/indl),
data=datalong,
control=lmc2) # SANN: 20k iterations
where resp.cidl = response-specific cluster id, lmc2 = a lmcControl object specifying SANN as optimisation method among other settings
I would like to find out if there is a way to allow for the correlations to differ by cluster via the corr option. Any suggestion or p...