kcponds {ade4} | R Documentation |
This data set contains informations about 33 ponds in De Maten reserve (Genk, Belgium).
data(kcponds)
: a data frame with 15 environmental variables(columns) on 33 ponds(rows)
: an object of class area
: a data frame with the coordinates of ponds
: an object of class neig
Variables of kcponds$tab
are the following ones : depth, area, O2 (oxygen concentration),
cond (conductivity), pH, Fe (Fe concentration), secchi (Secchi disk depth), N (NNO concentration),
TP (total phosphorus concentration), chla (chlorophyll-a concentration), EM (emergent macrophyte cover),
FM (floating macrophyte cover), SM (submerged macrophyte cover), denMI (total density of macroinvertebrates),
divMI (diversity macroinvertebrates)
Cottenie, K. (2002) Local and regional processes in a zooplankton metacommunity. PhD, Katholieke Universiteit Leuven, Leuven, Belgium.
http://www.kuleuven.ac.be/bio/eco/phdkarlcottenie.pdf
data(kcponds) w <- as.numeric(scalewt(kcponds$tab$N)) if(adegraphicsLoaded()) { if(requireNamespace("sp", quiet = TRUE)) { g1 <- s.label(kcponds$xy, Sp = kcponds$Spatial, pSp.col = "white", nb = kcponds$nb, plab.cex = 0, paxes.asp = "fill", plot = F) g2 <- s.label(kcponds$xy, Sp = kcponds$Spatial, pSp.col = "white", plabels.cex = 0.8, paxes.asp = "fill", plot = F) g3 <- s.value(kcponds$xy, w, psub.text = "Nitrogen concentration", paxe.asp = "fill", plot = F) G <- rbindADEg(g1, g2, g3, plot = T) } } else { par(mfrow=c(3, 1)) area.plot(kcponds$area) s.label(kcponds$xy, add.p = TRUE, cpoi = 2, clab = 0) s.label(kcponds$xy, add.p = TRUE, cpoi = 3, clab = 0) s.label(kcponds$xy, add.p = TRUE, cpoi = 0, clab = 0, neig = kcponds$neig, cneig = 1) area.plot(kcponds$area) s.label(kcponds$xy, add.p = TRUE, clab = 1.5) s.value(kcponds$xy, w, cleg = 2, sub = "Nitrogen concentration", csub = 4, possub = "topright", include = FALSE) par(mfrow = c(1, 1)) } ## Not run: par(mfrow = c(3, 1)) pca1 <- dudi.pca(kcponds$tab, scan = FALSE, nf = 4) if(requireNamespace("maptools", quiet = TRUE) & requireNamespace("spdep", quiet = TRUE)) { multi1 <- multispati(pca1, nb2listw(neig2nb(kcponds$neig)), scan = F, nfposi = 2, nfnega = 1) summary(multi1) } par(mfrow = c(1, 1)) ## End(Not run)