Skip to contents

NMDS plot unsing grindr layers

Source: R/gr-NMDS.R
plot_NMDS.Rd

Quickly vizualise MDS and NMDS objects and build customs plots using the layers. See examples.

Usage

plot_NMDS(
  x,
  f = NULL,
  axes = c(1, 2),
  points = TRUE,
  points_transp = 1/4,
  chull = TRUE,
  chullfilled = FALSE,
  labelgroups = FALSE,
  legend = TRUE,
  title = "",
  box = TRUE,
  axesnames = TRUE,
  palette = pal_qual
)

plot_MDS(
  x,
  f = NULL,
  axes = c(1, 2),
  points = TRUE,
  points_transp = 1/4,
  chull = TRUE,
  chullfilled = FALSE,
  labelgroups = FALSE,
  legend = TRUE,
  title = "",
  box = TRUE,
  axesnames = TRUE,
  palette = pal_qual
)

Arguments

x

the result of MDS or NMDS

f

factor specification to feed fac_dispatcher

axes

numeric of length two to select PCs to use (c(1, 2) by default)

points

logical whether to draw this with layer_points

points_transp

numeric to feed layer_points (default:0.25)

chull

logical whether to draw this with layer_chull

chullfilled

logical whether to draw this with layer_chullfilled

labelgroups

logical whether to draw this with layer_labelgroups

legend

logical whether to draw this with layer_legend

title

character if specified, fee layer_title (default to "")

box

logical whether to draw this using layer_box

axesnames

logical whether to draw this using layer_axesnames

palette

color palette to use col_summer by default

Value

a plot

See also

Other grindr: drawers, layers, layers_morphospace, mosaic_engine(), papers, pile(), plot_LDA(), plot_PCA()

Examples

### First prepare an NMDS object
x <- bot %>% efourier %>% NMDS
#> 'norm=TRUE' is used and this may be troublesome. See ?efourier #Details
#> 'nb.h' set to 10 (99% harmonic power)
#> 'comm' has negative data: 'autotransform', 'noshare' and 'wascores' set to FALSE
#> Warning: results may be meaningless because data have negative entries
#>                  in method “bray”
#> Run 0 stress 0.07227125 
#> Run 1 stress 0.07227125 
#> ... Procrustes: rmse 8.460091e-06  max resid 3.090637e-05 
#> ... Similar to previous best
#> Run 2 stress 0.07227125 
#> ... New best solution
#> ... Procrustes: rmse 5.108024e-06  max resid 1.882959e-05 
#> ... Similar to previous best
#> Run 3 stress 0.1558547 
#> Run 4 stress 0.1653725 
#> Run 5 stress 0.07227125 
#> ... Procrustes: rmse 1.464749e-06  max resid 5.407828e-06 
#> ... Similar to previous best
#> Run 6 stress 0.1566038 
#> Run 7 stress 0.07227125 
#> ... Procrustes: rmse 1.490112e-06  max resid 3.743155e-06 
#> ... Similar to previous best
#> Run 8 stress 0.07227125 
#> ... Procrustes: rmse 5.69994e-06  max resid 2.264085e-05 
#> ... Similar to previous best
#> Run 9 stress 0.07227125 
#> ... Procrustes: rmse 3.375566e-06  max resid 1.412437e-05 
#> ... Similar to previous best
#> Run 10 stress 0.07227125 
#> ... Procrustes: rmse 1.702584e-06  max resid 7.024297e-06 
#> ... Similar to previous best
#> Run 11 stress 0.07227125 
#> ... Procrustes: rmse 4.198101e-06  max resid 1.568591e-05 
#> ... Similar to previous best
#> Run 12 stress 0.07227125 
#> ... Procrustes: rmse 2.782198e-06  max resid 9.31499e-06 
#> ... Similar to previous best
#> Run 13 stress 0.07227125 
#> ... Procrustes: rmse 6.140461e-06  max resid 2.422899e-05 
#> ... Similar to previous best
#> Run 14 stress 0.07227125 
#> ... Procrustes: rmse 2.371183e-06  max resid 6.726446e-06 
#> ... Similar to previous best
#> Run 15 stress 0.07227125 
#> ... Procrustes: rmse 4.036543e-06  max resid 1.59661e-05 
#> ... Similar to previous best
#> Run 16 stress 0.07227125 
#> ... Procrustes: rmse 5.818498e-06  max resid 2.334124e-05 
#> ... Similar to previous best
#> Run 17 stress 0.07227125 
#> ... Procrustes: rmse 1.884869e-06  max resid 7.993225e-06 
#> ... Similar to previous best
#> Run 18 stress 0.07227125 
#> ... Procrustes: rmse 9.897688e-06  max resid 3.87602e-05 
#> ... Similar to previous best
#> Run 19 stress 0.1595332 
#> Run 20 stress 0.07227125 
#> ... Procrustes: rmse 6.944369e-06  max resid 2.757848e-05 
#> ... Similar to previous best
#> *** Best solution repeated 15 times

plot_NMDS(x)
#> Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...): plot.new has not been called yet
plot_NMDS(x, ~type) %>% layer_stars() %>% layer_labelpoints()
#> Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...): plot.new has not been called yet

### Same on MDS object
x <- bot %>% efourier %>% MDS
#> 'norm=TRUE' is used and this may be troublesome. See ?efourier #Details
#> 'nb.h' set to 10 (99% harmonic power)

plot_MDS(x)
#> Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...): plot.new has not been called yet
plot_MDS(x, ~type) %>% layer_stars() %>% layer_labelpoints()
#> Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...): plot.new has not been called yet