Gang

If you want to analyze all variables simultaneously and account for some correlational structure among the different response variables, then the best strategy is to pre-whiten the data and then use lmer. The methods for pre-whitening are described in detail in Pinhiero and Bates in the GLS chapter.  


-----Original Message-----
From: r-help-bounces@r-project.org on behalf of Gang Chen
Sent: Tue 8/5/2008 2:26 PM
To: r-help@stat.math.ethz.ch
Subject: [R] Mixed model with multiple response variables?
 
Hi,

I have a data set collected from 10 measurements (response variables)
on two groups (healthy and patient) of subjects performing 4 different
tasks. In other words there are two fixed factors (group and task),
and 10 response variables. I could analyze the data with aov() or
lme() in package nlme for each response variable separately, but since
most likely there are correlations among the 10 response variables,
would it be more meaningful to run a MANOVA? However manova() in R
seems not to allow an error term in the formula. What else can I try
for this kind of multivariate mixed model?

Also, if I want to find out which response variables (among the 10
measurements) are statistically significant in terms of acting as
indicators for group difference, what kind of statistical analysis
would help me sort them out?

Thanks in advance,
Gang

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


	[[alternative HTML version deleted]]

