ade4 (version 1.4-16)

withincoinertia: Within-class coinertia analysis

Description

Performs a within-class analysis after a coinertia analysis

Usage

withincoinertia(obj, fac, scannf = TRUE, nf = 2)
## S3 method for class 'witcoi':
plot(x, xax = 1, yax = 2, ...)
## S3 method for class 'witcoi':
print(x, ...)

Arguments

obj

a coinertia analysis (object of class coinertia) obtained by the function coinertia

fac

a factor partitioning the rows in classes

scannf

a logical value indicating whether the eigenvalues barplot should be displayed

nf

if scannf FALSE, an integer indicating the number of kept axes

xax

the axis number to be plotted for the x-axis

yax

the axis number to be plotted for the y-axis

x

an object of class witcoi

...

further arguments passed to or from other methods

Value

• An object of the class witcoi. Outputs are described by the print function

Details

This analysis is equivalent to do a within-class analysis on each initial dudi, and a coinertia analysis on the two within analyses. This function returns additional outputs for the interpretation.

References

Franquet E., Doledec S., and Chessel D. (1995) Using multivariate analyses for separating spatial and temporal effects within species-environment relationships. Hydrobiologia, 300, 425--431.

See Also

coinertia, within

Examples

data(meaudret)
pca1 <- dudi.pca(meaudret$mil, scan = FALSE, nf = 4)
pca2 <- dudi.pca(meaudret$fau, scal = FALSE, scan = FALSE, nf = 4)
   
wit1 <- within(pca1, meaudret$plan$sta, scan = FALSE, nf = 2)
wit2 <- within(pca2, meaudret$plan$sta, scan = FALSE, nf = 2)
coiw <- coinertia(wit1, wit2, scannf = FALSE)

coi <- coinertia(pca1,pca2,scannf = FALSE,nf=3)
coi.w <- withincoinertia(coi,meaudret$plan$sta, scannf = FALSE)
## coiw and coi.w are equivalent

plot(coi.w)