Fortify a

Fortifies an object of class "metaMDS" to produce a data frame of the selected axis scores in long format, suitable for plotting with ggplot2::ggplot().

Usage

# S3 method for class 'metaMDS'
fortify(model, data, axes = 1:2, layers = c("sites", "species"), ...)

# S3 method for class 'metaMDS'
tidy(x, data, layers = c("sites", "species"), ...)

Arguments

model, x: an object of class "metaMDS", the result of a call to vegan::metaMDS().
data: currently ignored.
axes: numeric; which axes to extract scores for.
layers: character; the scores to extract in the fortified object. Passed to display in the respective vegan::scores() method.
...: additional arguments passed to vegan::scores.metaMDS(). Note you can't use display.

Value

A data frame in long format containing the ordination scores. The first two components are the axis scores.

Details

TODO

Author

Gavin L. Simpson

Examples


library("vegan")

data(dune)

ord <- metaMDS(dune)
#> Run 0 stress 0.1192678 
#> Run 1 stress 0.1192678 
#> ... Procrustes: rmse 3.906726e-05  max resid 9.79081e-05 
#> ... Similar to previous best
#> Run 2 stress 0.1183186 
#> ... New best solution
#> ... Procrustes: rmse 0.02026975  max resid 0.06495941 
#> Run 3 stress 0.1939202 
#> Run 4 stress 0.2066373 
#> Run 5 stress 0.1808911 
#> Run 6 stress 0.1183186 
#> ... Procrustes: rmse 2.245981e-06  max resid 7.540827e-06 
#> ... Similar to previous best
#> Run 7 stress 0.1183186 
#> ... Procrustes: rmse 4.917067e-06  max resid 1.600171e-05 
#> ... Similar to previous best
#> Run 8 stress 0.1183186 
#> ... New best solution
#> ... Procrustes: rmse 5.76826e-06  max resid 1.830903e-05 
#> ... Similar to previous best
#> Run 9 stress 0.119268 
#> Run 10 stress 0.1192678 
#> Run 11 stress 0.1183186 
#> ... New best solution
#> ... Procrustes: rmse 1.764602e-06  max resid 5.686181e-06 
#> ... Similar to previous best
#> Run 12 stress 0.1183186 
#> ... Procrustes: rmse 1.111141e-06  max resid 2.854289e-06 
#> ... Similar to previous best
#> Run 13 stress 0.1183186 
#> ... Procrustes: rmse 1.382648e-05  max resid 4.402376e-05 
#> ... Similar to previous best
#> Run 14 stress 0.2192919 
#> Run 15 stress 0.1192678 
#> Run 16 stress 0.1192678 
#> Run 17 stress 0.1192678 
#> Run 18 stress 0.1183186 
#> ... Procrustes: rmse 8.63055e-06  max resid 2.735838e-05 
#> ... Similar to previous best
#> Run 19 stress 0.1192678 
#> Run 20 stress 0.1183186 
#> ... Procrustes: rmse 2.414588e-06  max resid 6.802962e-06 
#> ... Similar to previous best
#> *** Best solution repeated 5 times
head(fortify(ord))
#> # A tibble: 6 × 4
#>   score label   nmds1   nmds2
#>   <fct> <chr>   <dbl>   <dbl>
#> 1 sites 1     -0.841  -0.716 
#> 2 sites 2     -0.505  -0.409 
#> 3 sites 3     -0.0827 -0.437 
#> 4 sites 4     -0.116  -0.522 
#> 5 sites 5     -0.627  -0.0867
#> 6 sites 6     -0.543   0.113