[R-sig-eco] "partial" ordination

[Duarte Viana]

Hello all,

I'm interested in applying an ordination method (e.g. PCA or NMDS) to
a species matrix. However, I want to remove the potential effect of
the environment, which I'm not interested in modelling, from the
ordination (i.e. partialling out the environmental effect).

I've thought in the following approach (with the vegan package):

# spe: species matrix
# env: environment matrix
partial.pca<-rda(spe~1+Condition(env))
plot(partial.pca)

Is this a reasonable approach?

Some feedback would be great.
Thanks,

Duarte