Learn R Programming

ade4 (version 1.6-2)

coinertia: Coinertia Analysis

Description

The coinertia analysis performs a double inertia analysis of two tables.

Usage

coinertia(dudiX, dudiY, scannf = TRUE, nf = 2)
## S3 method for class 'coinertia':
plot(x, xax = 1, yax = 2, \dots) 
## S3 method for class 'coinertia':
print(x, \dots) 
## S3 method for class 'coinertia':
summary(object, \dots)

Arguments

dudiX
a duality diagram providing from one of the functions dudi.coa, dudi.pca, ...
dudiY
a duality diagram providing from one of the functions dudi.coa, dudi.pca, ...
scannf
a logical value indicating whether the eigenvalues bar plot should be displayed
nf
if scannf FALSE, an integer indicating the number of kept axes
x, object
an object of class 'coinertia'
xax, yax
the numbers of the x-axis and the y-axis
...
further arguments passed to or from other methods

Value

  • Returns a list of class 'coinertia', sub-class 'dudi' containing:
  • callcall
  • rankrank
  • nfa numeric value indicating the number of kept axes
  • RVa numeric value, the RV coefficient
  • eiga numeric vector with all the eigenvalues
  • lwa numeric vector with the rows weigths (crossed table)
  • cwa numeric vector with the columns weigths (crossed table)
  • taba crossed table (CT)
  • liCT row scores (cols of dudiY)
  • l1Principal components (loadings for cols of dudiY)
  • coCT col scores (cols of dudiX)
  • c1Principal axes (cols of dudiX)
  • lXRow scores (rows of dudiX)
  • mXNormed row scores (rows of dudiX)
  • lYRow scores (rows of dudiY)
  • mYNormed row scores (rows of dudiY)
  • aXCorrelations between dudiX axes and coinertia axes
  • aYCorrelations between dudiY axes and coinertia axes

encoding

latin1

WARNING

IMPORTANT : dudi1 and dudi2 must have identical row weights.

References

Dol�dec, S. and Chessel, D. (1994) Co-inertia analysis: an alternative method for studying species-environment relationships. Freshwater Biology, 31, 277--294. Dray, S., Chessel, D. and J. Thioulouse (2003) Co-inertia analysis and the linking of the ecological data tables. Ecology, 84, 11, 3078--3089.

Examples

Run this code
data(doubs)
dudi1 <- dudi.pca(doubs$env, scale = TRUE, scan = FALSE, nf = 3)
dudi2 <- dudi.pca(doubs$fish, scale = FALSE, scan = FALSE, nf = 2)
coin1 <- coinertia(dudi1,dudi2, scan = FALSE, nf = 2)
s.arrow(coin1$l1, clab = 0.7)

s.arrow(coin1$c1, clab = 0.7)

par(mfrow = c(1,2))
s.corcircle(coin1$aX)
s.corcircle(coin1$aY)
par(mfrow = c(1,1))
 
coin1
summary(coin1)
plot(coin1)

Run the code above in your browser using DataLab