View source: R/spca_randtest.R
spca_randtest | R Documentation |
The function spca_randtest
implements Monte-Carlo tests for the
presence of significant spatial structures in a sPCA object. Two tests are
run, for global (positive autocorrelation) and local (negative
autocorrelation) structures, respectively. The test statistics used are the
sum of the absolute values of the corresponding eigenvalues.
spca_randtest(x, nperm = 499)
x |
A |
nperm |
The number of permutations to be used for the test. |
A list with two objects of the class 'randtest' (see
as.randtest
), the first one for 'global' structures
(positivie autocorrelation) and the second for 'local' structures (negative
autocorrelation).
Original code by Valeria Montano adapted by Thibaut Jombart.
## Not run: ## Load data data(sim2pop) ## Make spca spca1 <- spca(sim2pop, type = 1, scannf = FALSE, plot.nb = FALSE) spca1 plot(spca1) ## run tests (use more permutations in practice, e.g. 999) tests <- spca_randtest(spca1, nperm = 49) ## check results tests plot(tests[[1]]) # global structures ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.