多个分组做Procrustes test

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

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

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

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

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