tintoodiel               package:ade4               R Documentation

_T_i_n_t_o _a_n_d _O_d_i_e_l _e_s_t_u_a_r_y _g_e_o_c_h_e_m_i_s_t_r_y

_D_e_s_c_r_i_p_t_i_o_n:

     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).

_U_s_a_g_e:

     data(tintoodiel)

_F_o_r_m_a_t:

     'tintoodiel' is a list containing the following objects : 

     _x_y : a data frame that contains spatial coordinates of the 52
          sites   

     _t_a_b : a data frame with 12 columns (concentration of heavy metals)
          and 52 rows (sites)

     _n_e_i_g : an object of class 'neig'

_S_o_u_r_c_e:

     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.

_E_x_a_m_p_l_e_s:

     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))
     }

