randtest.discrimin {ade4} | R Documentation

## Monte-Carlo Test on a Discriminant Analysis (in C).

### Description

Test of the sum of a discriminant analysis eigenvalues (divided by the rank). Non parametric version of the Pillai's test. It authorizes any weighting.

### Usage

    randtest.discrimin(xtest, nrepet = 999, ...)

### Arguments

`xtest` | an object of class `discrimin`
`nrepet` | the number of permutations
`...` | further arguments passed to or from other methods

### Value

returns a list of class `randtest`

### Examples

    data(meaudret)
    pca1 <- dudi.pca(meaudret$mil, scan = FALSE, nf = 3)
    rand1 <- randtest(discrimin(pca1, meaudret$plan$dat, scan = FALSE), 99)
    rand1
    #Monte-Carlo test
    #Observation: 0.3035
    #Call: as.randtest(sim = sim, obs = obs)
    #Based on 999 replicates
    #Simulated p-value: 0.001
    plot(rand1, main = "Monte-Carlo test")
    summary.manova(manova(as.matrix(meaudret$mil)~meaudret$plan$dat), "Pillai")
    #                   Df Pillai approx F num Df den Df  Pr(>F)
    # meaudret$plan$dat  3   2.73    11.30     27     30 1.6e-09 ***
    # Residuals         16
    # ---
    # Signif. codes:  0 `***' 0.001 `**' 0.01 `*' 0.05 `.' 0.1 ` ' 1
    # 2.731/9 = 0.3034

[Package Contents]