An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-sig-ecology/attachments/20111108/25625482/attachment.pl>
Graphing Envfit Output
2 messages · Stephen LeDuc, Jari Oksanen
On 9/11/11 00:52 AM, "Stephen LeDuc" <sdleduc at gmail.com> wrote:
Dear R Listserve,
I'd like to take my results from the envfit function in the vegan package,
such as those below, and make a vector plot in a graphing package like
Sigma Plot. When graphed in R, it is my understanding that the results
are scaled depending upon the plot axes. Does R just multiply each
dimensional value (Dim1, Dim2) by a scaler to get X and Y coordinates for
each arrow endpoint? If so, is there a way to view the scaler for each
given axis? Any help on this would be greatly appreciated. Best, Steve
***VECTORS
Dim1 Dim2 r2 Pr(>r)
Mineral.N 0.125899 0.992043 0.5715 0.3127
AA_N 0.998794 -0.049102 0.8172 0.1283
SON 0.984886 0.173203 0.7901 0.1235
FF_C -0.222764 0.974872 0.5052 0.3463
P -0.034731 -0.999397 0.1737 0.7537
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
Jari Oksanen, Dept Biology, Univ Oulu, 90014 Finland jari.oksanen at oulu.fi, Ph. +358 400 408593, http://cc.oulu.fi/~jarioksa