多个分组做Procrustes test

2020-06-16 15:40:49 浏览数 (1)

前文:Procrustes test 介绍了其在vegan中的基本用法。但是存在一个局限性,只能输入两个矩阵且没办法分组。 正好昨天有读者问我相关问题,我顺手查了一下发现还有其他一些包可以做,并可以实现更一般化的分析Generalised Procrustes Analysis。

代码语言:javascript复制
library(FactoMineR)
?GPA
##GPA,Generalised Procrustes Analysis数据中可包含缺失值,并可以添加多个分组
data(wine)
res.gpa <- GPA(wine[,-(1:2)], group=c(5,3,10,9,2),
               name.group=c("olf","vis","olfag","gust","ens"))
X11()
plot (res.gpa)

另外还发现三个包shapes,Evomorph和smacof也可以做GPA,但是不能分组。

END

一个环境工程专业却做生信分析的深井冰博士,深受拖延症的困扰。想给自己一点压力,争取能够不定期分享学到的生信小技能,亦或看文献过程中的一些笔记与小收获,记录生活中的杂七杂八。

目前能力有限,尚不能创造知识,只是知识的搬运工。

欢迎分享,转载请联系我。

0 人点赞