有没有人知道在控制系统发育非独立性的同时执行MANOVA的R包或方法?
谢谢!
发布于 2013-05-31 20:11:49
sos包是您的朋友:
library('sos')
findFn('phylogenetic MANOVA')看起来像geiger包,更准确地说,aov.phylo执行系统发生方差或MANOVA。
下面是帮助中的一个示例:
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 https://stackoverflow.com/questions/16856408
复制相似问题