Get bootstrapped estimates of SparCC correlation coefficients. To get empirical p-values, pass this output to pval.sparccboot
.
sparccboot(data, sparcc.params = list(), statisticboot = function(data,
indices) triu(do.call("sparcc", c(list(data[indices, , drop = FALSE]),
sparcc.params))$Cor), statisticperm = function(data, indices)
triu(do.call("sparcc", c(list(apply(data[indices, ], 2, sample)),
sparcc.params))$Cor), R, ncpus = 1, ...)
Community count data
named list of parameters to pass to sparcc
function which takes data and bootstrap sample indices and results the upper triangle of the bootstapped correlation matrix
function which takes data and permutated sample indices and results the upper triangle of the null correlation matrix
number of bootstraps
number of cores to use for parallelization
additional arguments that are passed to boot::boot