문제

Does anyone know of a package or method in R of carrying out a MANOVA whilst controlling for phylogenetic non-independence?

Thank you!

도움이 되었습니까?

해결책

sos package is your friend:

library('sos')
findFn('phylogenetic MANOVA')

seems like geiger package and more precisely aov.phylo performs phylogenetic ANOVA or MANOVA.

Here an example from the help :

library(geiger)
geo=get(data(geospiza))
dat=geo$dat
d1=dat[,1]
grp<-as.factor(c(rep(0, 7), rep(1, 6)))
names(grp)=rownames(dat)

## MANOVA
x=aov.phylo(dat~grp, geo$phy, nsim=50, test="Wilks")
Multivariate Analysis of Variance Table

Response: dat
          Df   Wilks approx-F num-Df den-Df   Pr(>F) Pr(phy)
group      1 0.27872   3.6229      5      7 0.061584   0.549
Residuals 11      
라이센스 : CC-BY-SA ~와 함께 속성
제휴하지 않습니다 StackOverflow
scroll top