fastCorrDiff (version 0.5)

SS.boot: Bootstraping function for spectral screening

Description

The function initialize a bootstraping procedure for spectral score calculation, which can be used to determine a threshold for SS selection.

Usage

SS.boot(X1, X2, K, B, sv = FALSE, spearman = FALSE)

Arguments

X1

Matrix of the first group of observations. Each row is variable and each column is an observation. Note that this is different from the usual way of statistical data matrices, because usually number of variables is much larger than the number of observations in this case.

X2

Matrix of the second group of observations. Each row is variable and each column is an observation. Note that this is different from the usual way of statistical data matrices, because usually number of variables is much larger than the number of observations in this case.

K

Rank K in spectral screening

B

Number of bootstrap samples

sv

logical value. If TRUE, the square roots of singular values will be used in the score calculation.

spearman

logical value. If TRUE, Spearman's correlation is used. Otherwise (default), the calculation is based on Pearson's correlation.

Value

A matrix with B columns. Each column is a bootstraped score for the p variables.

Details

See the details in the reference paper.

References

Tianxi Li, Xiwei Tang, and Ajay Chatrath. Compressed spectral screening for large-scale differential correlation analysis with application in selecting Glioblastoma gene modules. arXiv preprint arXiv:2111.03721, 2021.

See Also

SS

Examples

# NOT RUN {
X1 <- matrix(rnorm(200),10,20)
X2 <- matrix(rnorm(200),10,20)


bt <- SS.boot(X1,X2,3, 10)

# }