Skip to content

Graphing Envfit Output

2 messages · Stephen LeDuc, Jari Oksanen

#
On 9/11/11 00:52 AM, "Stephen LeDuc" <sdleduc at gmail.com> wrote:

            
Steve,

There are two different scaling issues: (1) scaling arrows against each
other, and (2) scaling the bunch of arrows to fit the ordination graph.

(1) The coordinates above are direction cosines or the coordinates of the
heads of unit-length arrows. The envfit plot() function scales these arrows
by their correlations (square root of column r2) so that 'P' will be longer
than 'SON' in your example above. If you want to see the scaled arrows, you
should use function scores() for the envfit() result. If you want to get the
unit length arrows, you cannot use scores() but you must directly access
envfit_result$vectors$arrows.

(2) The bunch of r-adjusted arrows is scaled to fill the current graph. An
internal vegan function ordiArrowMul() that finds the multiplier for the
current plot. Replotting ordination in different size will change the
multiplier. The ordiArrowMul() is unexported in vegan 2.0-* so that you must
call it as vegan:::ordiArrowMul() to see the current multiplier. This
multiplier concerns the r-adjusted arrows which you can see with scores(...,
display="vectors").

The arrow scaling has been asked before, and therefore I thought this is a
FAQ, and therefore I have something in
http://vegan.r-forge.r-project.org/FAQ-vegan.html

This seems to cover only point (2) above. If you think that this FAQ needs
improvement and clarification, please let me know, and we'll fix the FAQ.

Cheers, Jari