aaa

Ruby code posted
created at 17 Jul 14:45

Edit | Back 
1
2
3
4
5
6
7
8
9
10
11
12
13
 
require(vegan)
data(dune)
data(dune.env)

# NMDS on fouth-root transformed abundances
un_nmds <- metaMDS(dune^0.25, distance = "bray")

# plot it
cols = rainbow(4)
plot(un_nmds, display = "sites", type = "n")
points(un_nmds, col = cols[dune.env$Management], pch = 16, cex = 1.5)
ordihull(un_nmds, groups = dune.env$Management, lty = "dotted")
ordispider(un_nmds, groups = dune.env$Management, label = TRUE)
420 Bytes in 2 ms with coderay