Is there a way to perform a variable reduction for a partial canonical ordination (either redundancy analysis or correspondence analysis) with the function ordistep from the package vegan ? I checked Borcard et al. (2011) Numerical Ecology with R and I could not find an answer. I tried the following syntax and it returned the next error :
Error in formula.default(prda1) : invalid formula
prda1 <- rda(spp.h, env.partial, avu.rda)
rda.ordistep <- ordistep(rda(spp ~ 1, data = env),
scope = formula(partial.rda),
direction= "forward",
pstep=10000,
trace = FALSE)
?