Dear all,
I'm testing the effect of species and sex in my sample by using the principal
component scores of a PCA analysis.
I have 30 PCs and I tried to see if there is any significant difference from
males to females, given that there is a significant effect of phylogeny (factor
with several species).
I didi it like this:
Y<-PCA$pc.scores[,1:30]
fit <- manova(Y ~ sp*sex)
summary(fit, test="Wilks")
I get a barely significant p-value for the effect of sex and I'd like to know
for which of the species there is a difference between males and females.
I tried TukeyHSD(fit) but I get the following error:
Error in model.tables.aov(x, "means") :
'model.tables' is not implemented for multiple responses
So this has to do with the fact that I have a multivariate independent
variable. Is there an alternative function to this?
Thanks in advance,
Sérgio.
--
Institut des Sciences de l'Evolution
UMR5554, CNRS, IRD, EPHE
Université de Montpellier
Place Eugène Bataillon
34095 Montpellier Cedex 05
France
Email: sergio.ferreira-card...@umontpellier.fr
Tel: +33 (4 ) 67 14 46 52
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
