Learn R Programming

ade4 (version 1.7-8)

tintoodiel: Tinto and Odiel estuary geochemistry

Description

This data set contains informations about geochemical characteristics of heavy metal pollution in surface sediments of the Tinto and Odiel river estuary (south-western Spain).

Usage

data(tintoodiel)

Arguments

Format

tintoodiel is a list containing the following objects :

xy

: a data frame that contains spatial coordinates of the 52 sites

tab

: a data frame with 12 columns (concentration of heavy metals) and 52 rows (sites)

neig

: an object of class neig

Examples

Run this code
# NOT RUN {
data(tintoodiel)
if(!adegraphicsLoaded()) {
  
# }
# NOT RUN {
  if(requireNamespace("pixmap", quietly = TRUE)) {
    estuary.pnm <- pixmap::read.pnm(system.file("pictures/tintoodiel.pnm", package = "ade4"))
    s.label(tintoodiel$xy, pixmap = estuary.pnm, neig = tintoodiel$neig,
      clab = 0, cpoi = 2, cneig = 3, addax = FALSE, cgrid = 0, grid = FALSE)
  }
# }
# NOT RUN {
  
  estuary.pca <- dudi.pca(tintoodiel$tab, scan = FALSE, nf = 4)
  
  if(requireNamespace("maptools", quietly = TRUE) & requireNamespace("spdep", quietly = TRUE)) {
    estuary.listw <- spdep::nb2listw(neig2nb(tintoodiel$neig))
    estuary.pca.ms <- multispati(estuary.pca, estuary.listw, scan = FALSE, nfposi = 3, nfnega = 2)
    summary(estuary.pca.ms)
    par(mfrow = c(1, 2))
    barplot(estuary.pca$eig)
    barplot(estuary.pca.ms$eig)
    par(mfrow = c(1, 1))
}}
# }

Run the code above in your browser using DataLab