granulo {ade4}R Documentation

Granulometric Curves

Description

This data set gives the repartition in diameter classes of deposit samples.

Usage

data(granulo)

Format

granulo is a list of 2 components.

tab
contains the 49 deposit samples, 9 diameter classes, weight of grains by size class
born
contains the boundaries of the diameter classes

Source

Gaschignard-Fossati, O. (1986) Répartition spatiale des macroinvertébrés benthiques d'un bras vif du Rhône. Rôle des crues et dynamique saisonnière. Thèse de doctorat, Université Lyon 1.

Examples

data(granulo)
w <- t(apply(granulo$tab, 1, function (x) x / sum(x)))
w <- data.frame(w)
wtr <- data.frame(t(w))
wmoy <- data.frame(matrix(apply(wtr, 1, mean), 1))
d1 <- dudi.pca(w, scal = FALSE, scan = FALSE)
wmoy <- suprow.default(d1, wmoy)$lisup
s.arrow(d1$c1, clab = 1.5)
s.distri(d1$c1, wtr, cstar = 0.33, cell = 0, 
    axesell = FALSE, add.p = TRUE, clab = 0.75)
s.label(wmoy, cpoi = 5, clab = 0, add.p = TRUE)

Worked out examples


> library(ade4)
> ### Name: granulo
> ### Title: Granulometric Curves
> ### Aliases: granulo
> ### Keywords: datasets
> 
> ### ** Examples
> 
> data(granulo)
> w <- t(apply(granulo$tab, 1, function (x) x / sum(x)))
> w <- data.frame(w)
> wtr <- data.frame(t(w))
> wmoy <- data.frame(matrix(apply(wtr, 1, mean), 1))
> d1 <- dudi.pca(w, scal = FALSE, scan = FALSE)
> wmoy <- suprow.default(d1, wmoy)$lisup
> s.arrow(d1$c1, clab = 1.5)
> s.distri(d1$c1, wtr, cstar = 0.33, cell = 0, 
+     axesell = FALSE, add.p = TRUE, clab = 0.75)
> s.label(wmoy, cpoi = 5, clab = 0, add.p = TRUE)
> 
> 
> 
> 

[Package ade4 Index]