evolqg (version 0.1-7)

BootstrapRep: Bootstrap analysis via resampling

Description

Calculates the repeatability of the covariance matrix of the suplied data via bootstrap resampling

Usage

BootstrapRep(ind.data, ComparisonFunc, ..., iterations = 1000,
  correlation = FALSE, parallel = FALSE)

Arguments

ind.data

Matrix of residuals or indiviual measurments

ComparisonFunc

comparison function

...

Aditional arguments passed to ComparisonFunc

iterations

Number of resamples to take

correlation

If TRUE, correlation matrix is used, else covariance matrix.

parallel

if TRUE computations are done in parallel. Some foreach backend must be registered, like doParallel or doMC.

Value

• returns the mean repeatability, that is, the mean value of comparisons from samples to original statistic.

Details

Samples with replacement are taken from the full population, a statistic calculated and compared to the full population statistic.

See Also

MonteCarloStat, AlphaRep

Examples

BootstrapRep(iris[,1:4], MantelCor, iterations = 5, correlation = TRUE)

BootstrapRep(iris[,1:4], RandomSkewers, iterations = 50)

BootstrapRep(iris[,1:4], KrzCor, iterations = 50, correlation = TRUE)

BootstrapRep(iris[,1:4], PCAsimilarity, iterations = 50)

#Multiple threads can be used with some foreach backend library, like doMC or doParallel
#library(doParallel)
##Windows:
#cl <- makeCluster(2)
#registerDoParallel(cl)
##Mac and Linux:
#registerDoParallel(cores = 2)
#BootstrapRep(iris[,1:4], PCAsimilarity,
#             iterations = 5,
#             parallel = TRUE)