Procrustes Rotation of Two Configurations and PROTEST
procrustes_analyse(b_res1, b_res2, nperm = 999, ...)
pro_res
data(otutab, package = "pcutils")
b_analyse(otutab, method = "pca") -> b_res1
#> four dataframes in a list, 1 is eig, 2 is sample_site, 3 is var, 4 is var contribution
b_analyse(otutab * abs(rnorm(10)), method = "pca") -> b_res2
#> four dataframes in a list, 1 is eig, 2 is sample_site, 3 is var, 4 is var contribution
pro_res <- procrustes_analyse(b_res1, b_res2)
#> b_res1 use the PCA1
#> b_res2 use the PCA1
plot(pro_res, "Group", metadata)