In R: I'm doing a partial RDA, using a forward selection procedure to identify the explanatory variables which are most important for a coleoptera community.
But how do I get not only adjusted R2 values, but also p-values for all tested variables when using function ordiR2step?
And why is function ordiR2step so much faster than ordistep in doing so??
I first answer to the last question: ordiR2step
is so much faster than ordistep
because it does not estimate all p-values.
The p-values are found with permutation tests, and running those tests to all variables takes time.
There is no way of getting all so-called p-values in ordiR2step
because the method does not use those, as is explained in the source paper cited (Blanchet et al. 2008 in the help page) and the help page. That is the reason why the method is called ordi
R2step
: it uses primarily R2 instead of p-values. If you want to get all p-values, you must use a method that was written to evaluate those, that is ordistep
.