Hi David
thanks for your suggestions. Will try the different approaches/packages
to make sure I am on the right way.
 
Cheers
David
 
______________________________________________ 

David Costantini, PhD

http://www.davidcostantini.it
NERC Postdoctoral research associate

Institute of Biodiversity, Animal Health and Comparative Medicine
School of Life Sciences
College of Medical, Veterinary and Life Sciences
University of Glasgow
Graham Kerr Building, room 511
Glasgow G12 8QQ, UK

See also my association Ornis italica
http://www.ornisitalica.com
http://www.birdcam.it
____________________________________________________


________________________________

From: David Duffy [mailto:David.Duffy@qimr.edu.au]
Sent: Mon 21/05/2012 1.16
To: David Costantini
Cc: r-sig-mixed-models@r-project.org
Subject: Re: [R-sig-ME] help with MANOVE with mixed effects



On Sat, 19 May 2012, David Costantini wrote:

> Dear All I wonder if anyone can help with this issue? I have two
> fixed factors (experimental group and sex) and a random factor (brood).
> I need the random factor to control for pseudoreplication because each
> brood provides a certain number of siblings, which are
> pseudoreplicates because share genes/environment. Then I have 8
> variables. so I would like to to compare the covariance structures among
> groups while controlling for the random factor. I was told that R might
> do that, I mean a manova with mixed effects (manova(lme(....))

Not as easily as that ;) One can represent a multivariate model as a lmer
or lme univariate model with appropriate random effects representing each
item (there is recent query about this on the list).

MCMCglmm will fit as a multivariate model.

Or you can use structural equation modelling packages eg OpenMx.

Or you can go outside of R to other packages eg WOMBAT
(http://didgeridoo.une.edu.au/km/wombat.php), MENDEL, ASReml...

--
| David Duffy (MBBS PhD)                                         ,-_|\
| email: davidD@qimr.edu.au  ph: INT+61+7+3362-0217 fax: -0101  /     *
| Epidemiology Unit, Queensland Institute of Medical Research   \_,-._/
| 300 Herston Rd, Brisbane, Queensland 4029, Australia  GPG 4D0B994A v



	[[alternative HTML version deleted]]

