[R-sig-ME] can mixed effects models handle two dependent variables

Afshartous, David DAfshartous at med.miami.edu
Mon Nov 23 22:08:34 CET 2009

A good starting point is Doran & Lockwood (2006), Journal of Educational and Behavioral Statistics, 31:205-230 (see p.223).   There are also a couple of recent threads on this issue, e.g., see
[R-sig-ME] Multivariate mixed effects model
[R-sig-ME] weighting nlme and multivariate outcomes


On 11/23/09 3:57 PM, "Jude Phillips" <birdlists at gmail.com> wrote:


I am carrying out a stable isotope study, where I have measured stable
isotopes of nitrogen and carbon in plants at different sites.  I have been
able to create mixed models for each isotope separately - eg

fit1.c<-lme(X13c~soil+manure, random=~1|site, weights=varIdent(form = ~ 1 |
soil), data)
fit1.n<-lme(X15n~soil+manure, random=~1|site, weights=varIdent(form = ~ 1 |
soil), data)

What I'm wondering is if I can test both stable isotopes at the same time,
in a MANOVA style analysis.  I can't seem to find any info on this - if
someone could point me in the direction of some info, or tell me if this
isn't possible , I'd much appreciate it.

Thanks, Sarah Adams

        [[alternative HTML version deleted]]

R-sig-mixed-models at r-project.org mailing list

More information about the R-sig-mixed-models mailing list