ktab.list.dudi {ade4}R Documentation

Creation of a K-tables from a list of duality diagrams

Description

creates a list of class ktab from a list of duality diagrams.

Usage

ktab.list.dudi(obj, rownames = NULL, colnames = NULL, tabnames = NULL)

Arguments

obj a list of objects of class 'dudi'. Each element of the list must have the same row names for $tab and even for $lw
rownames the row names of the K-tables (otherwise the row names of the $tab)
colnames the column names of the K-tables (otherwise the column names of the $tab)
tabnames the names of the arrays of the K-tables (otherwise the names of the obj if they exist, or else "Ana1", "Ana2", ...)

Value

returns a list of class ktab. See ktab

Author(s)

Daniel Chessel
Anne B Dufour dufour@biomserv.univ-lyon1.fr

Examples

data(euro123)
pca1 <- dudi.pca(euro123$in78, scale = FALSE, scann = FALSE)
pca2 <- dudi.pca(euro123$in86, scale = FALSE, scann = FALSE)
pca3 <- dudi.pca(euro123$in97, scale = FALSE, scann = FALSE)
ktabeuro <- ktab.list.dudi(list(pca1, pca2, pca3), 
    tabnames = c("1978","1986","1997"))
kplot(sepan(ktabeuro), mfr = c(2,2), clab.c = 1.5)

data(meaudret)
w1 <- split(meaudret$mil,meaudret$plan$dat)
ll <- lapply(w1, dudi.pca, scann = FALSE)
kta <- ktab.list.dudi(ll, rownames <- paste("Station", 1:5, sep=""))
kplot(sepan(kta), clab.r = 1.5, clab.c = 0.75)

data(jv73)
w <- split(jv73$poi, jv73$fac.riv)
wjv73poi <- lapply(w, dudi.pca, scal = FALSE, scan = FALSE)
wjv73poi <- lapply(wjv73poi, t)
wjv73poi <- ktab.list.dudi(wjv73poi)
kplot(sepan(wjv73poi), permut = TRUE, traj = TRUE)

Worked out examples


> library(ade4)
> ### Name: ktab.list.dudi
> ### Title: Creation of a K-tables from a list of duality diagrams
> ### Aliases: ktab.list.dudi
> ### Keywords: multivariate
> 
> ### ** Examples
> 
> data(euro123)
> pca1 <- dudi.pca(euro123$in78, scale = FALSE, scann = FALSE)
> pca2 <- dudi.pca(euro123$in86, scale = FALSE, scann = FALSE)
> pca3 <- dudi.pca(euro123$in97, scale = FALSE, scann = FALSE)
> ktabeuro <- ktab.list.dudi(list(pca1, pca2, pca3), 
+     tabnames = c("1978","1986","1997"))
> kplot(sepan(ktabeuro), mfr = c(2,2), clab.c = 1.5)
> 
> data(meaudret)
> w1 <- split(meaudret$mil,meaudret$plan$dat)
> ll <- lapply(w1, dudi.pca, scann = FALSE)
> kta <- ktab.list.dudi(ll, rownames <- paste("Station", 1:5, sep=""))
> kplot(sepan(kta), clab.r = 1.5, clab.c = 0.75)
> 
> data(jv73)
> w <- split(jv73$poi, jv73$fac.riv)
> wjv73poi <- lapply(w, dudi.pca, scal = FALSE, scan = FALSE)
> wjv73poi <- lapply(wjv73poi, t)
> wjv73poi <- ktab.list.dudi(wjv73poi)
> kplot(sepan(wjv73poi), permut = TRUE, traj = TRUE)
> 
> 
> 
> 

[Package ade4 Index]