randtest.between {ade4}R Documentation

Monte-Carlo Test on the between-groups inertia percentage (in C).

Description

Performs a Monte-Carlo test on the between-groups inertia percentage.

Usage

## S3 method for class 'between':
randtest(xtest, nrepet = 999, ...)

Arguments

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

Value

a list of the class randtest

Author(s)

Jean Thioulouse ade4-jt@biomserv.univ-lyon1.fr

References

Romesburg, H. C. (1985) Exploring, confirming and randomization tests. Computers and Geosciences, 11, 19–37.

Examples

data(meaudret)
pca1 <- dudi.pca(meaudret$mil, scan = FALSE, nf = 3)
rand1 <- randtest(between(pca1, meaudret$plan$dat, scan = FALSE), 99)
rand1
plot(rand1, main = "Monte-Carlo test")

Worked out examples


> library(ade4)
> ### Name: randtest.between
> ### Title: Monte-Carlo Test on the between-groups inertia percentage (in
> ###   C).
> ### Aliases: randtest.between
> ### Keywords: multivariate nonparametric
> 
> ### ** Examples
> 
> data(meaudret)
> pca1 <- dudi.pca(meaudret$mil, scan = FALSE, nf = 3)
> rand1 <- randtest(between(pca1, meaudret$plan$dat, scan = FALSE), 99)
> rand1
Monte-Carlo test
Call: randtest.between(xtest = between(pca1, meaudret$plan$dat, scan = FALSE), 
    nrepet = 99)

Observation: 0.3722686 

Based on 99 replicates
Simulated p-value: 0.01 
Alternative hypothesis: greater 

    Std.Obs Expectation    Variance 
3.302810875 0.158975973 0.004170464 
> plot(rand1, main = "Monte-Carlo test")

> 
> 
> 
>
[Package ade4 Index]