I'm trying to figure out how I should be interpreting the output of a procrustes analysis in vegan. Specifically, what is the difference between the sum of squares results in the summary(procrustes_object)
vs protest
functions?
Here is a toy example:
library(vegan)
#some environmental data
env <- read.table(header=TRUE, row.names="site",text="site temp rainfall ph
A 11 600 6
B 13 550 8
C 13 500 6
D 13 450 8
E 14 400 7
F 15 400 7")
#some animal abundances
animals <- read.table(header=TRUE, row.names="site",text="site frogs birds mammals insects spiders
A 54 45 32 88 77
B 64 54 30 85 74
C 50 49 28 81 50
D 30 43 26 84 60
E 13 44 24 86 40
F 14 51 22 50 22")
#some tree abundances
trees <- read.table(header=TRUE, row.names="site",text="site elm oak maple yew willow pine
A 1 55 44 81 34 88
B 3 58 50 78 40 87
C 7 56 40 74 33 75
D 3 54 24 77 22 80
E 1 55 10 79 14 70
F 7 57 11 43 15 61")
#run CCAs for animals and trees against available environmental data
cca_animals <- cca(animals, env)
cca_trees <- cca(trees, env)
pro1 <- procrustes(cca_animals, cca_trees) #run procrustes to compare animal and tree CCAs
summary(pro1) #procrustes sum of squares = 0.786
protest(cca_animals, cca_trees) #procrustes sum of quares = 0.047
#correlation in a symmetric procrustes rotation=0.976
The two sum of squares results I get are very different, but I'm not sure which I should be using. Also, how is this related to the "correlation" value in the protest output? As ever, I'm sure that this question is largely motivated by my ignorance of basic statistics, but as always, remedial schooling is very appreciated.
protest
runs symmetric Procrustes analysis, whereasprocrustes
defaults to non-symmetric analysis.