Performs a Monte-Carlo test on on the percentage of explained (i.e. constrained) inertia. The statistic is the ratio of the inertia (sum of eigenvalues) of the constrained analysis divided by the inertia of the unconstrained analysis.
# S3 method for pcaiv
randtest(xtest, nrepet = 99, ...)
# S3 method for pcaivortho
randtest(xtest, nrepet = 99, ...)a list of the class randtest
an object of class pcaiv, pcaivortho or caiv
the number of permutations
further arguments passed to or from other methods
Stéphane Dray stephane.dray@univ-lyon1.fr, original code by Raphaël Pélissier
data(rpjdl)
millog <- log(rpjdl$mil + 1)
coa1 <- dudi.coa(rpjdl$fau, scann = FALSE)
caiv1 <- pcaiv(coa1, millog, scan = FALSE)
randtest(caiv1)
Run the code above in your browser using DataLab