## Synoptic table of Scheden vegetation data
library(cluster)
pam1 <- pam(schedenveg, 4) # PAM clustering with 4 clusters output
## 1) Unordered synoptic percentage frequency table
percfreq <- syntable(schedenveg, pam1$clustering, abund = "percentage",
type = "percfreq")
percfreq # view results
## 2) Differential species analysis
differential <- syntable(schedenveg, pam1$clustering, abund = "percentage",
type = "diffspec")
# show complete table with differential character of species
differential$syntable
# list differential species for second cluster
differential$differentials[2]
## 3) Synoptic table with phi fidelity
phitable <- syntable(schedenveg, pam1$clustering, abund = "percentage",
type = "phi")
phitable
## 4) Synoptic percentage frequency table based on historical classification from 1997
percfreq <- syntable(schedenveg, schedenenv$comm, abund = "percentage",
type = "percfreq")
percfreq
Run the code above in your browser using DataLab