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)

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

Source

Borrego, J., Morales, J.A., de la Torre, M.L. and Grande, J.A. (2002) Geochemical characteristics of heavy metal pollution in surface sediments of the Tinto and Odiel river estuary (south-western Spain). Environmental Geology, 41, 785–796.

Examples

data(tintoodiel)

## Not run: 
if (require(pixmap, quiet = TRUE)){
estuary.pnm <- 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)
}

## End(Not run)

estuary.pca <- dudi.pca(tintoodiel$tab, scan = FALSE, nf = 4)
if (require(maptools, quiet = TRUE) & require(spdep, quiet = TRUE)) {
estuary.listw <- 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))
}

Worked out examples

> library(ade4)
> ### Name: tintoodiel
> ### Title: Tinto and Odiel estuary geochemistry
> ### Aliases: tintoodiel
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(tintoodiel)
> 
> if (require(pixmap, quiet = TRUE)){
+ estuary.pnm <- 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)
+ }

> 
> 
> estuary.pca <- dudi.pca(tintoodiel$tab, scan = FALSE, nf = 4)
> if (require(maptools, quiet = TRUE) & require(spdep, quiet = TRUE)) {
+ estuary.listw <- 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))
+ }
deldir 0.0-12 

     Please note: The process for determining duplicated points
     has changed from that used in version 0.0-9 (and previously).


Multivariate Spatial Analysis
Call: multispati(dudi = estuary.pca, listw = estuary.listw, scannf = FALSE, 
    nfposi = 3, nfnega = 2)

Scores from the initial duality diagramm:
         var       cum     ratio        moran
RS1 4.401792  4.401792 0.2751120  0.222781944
RS2 3.133961  7.535753 0.4709846  0.217365304
RS3 1.818129  9.353883 0.5846177  0.004424936
RS4 1.379856 10.733739 0.6708587 -0.209208617

Multispati eigenvalues decomposition:
            eig      var      moran
CS1   1.5292058 3.707268  0.4124886
CS2   0.6486201 2.415632  0.2685095
CS3   0.4969485 1.126228  0.4412505
CS15 -0.2964507 1.648786 -0.1797994
CS16 -0.4623624 1.365391 -0.3386300
> 
> 
> 
>