cor_columns: Number of columns which are highly correlated to other...

In eilslabs/epic: Integrative Analysis and Visualization of Epigenomic Sequencing Data

Description

Number of columns which are highly correlated to other columns

Usage

cor_columns(x, abs_cutoff = 0.5, size = 1000, mc = 1, ...)

Arguments

x	a matrix, correlation is calculated by columns
abs_cutoff	cutoff of absolute correlation
size	size of blocks
mc	number of cores
...	pass to cor

Details

For each column, it looks for number of other columns which correlate with absolute correlation coefficient larger tham abs_cutoff. The calculation involves pair-wise correlation of all columns in the matrix. When number of columns is huge in a matrix, it is out of ability of R to store such a long vector. This function solves this problem by splitting the columns into k blocks and looks at each block sequentially or parallel.

The code is partially adapted from https://rmazing.wordpress.com/2013/02/22/bigcor-large-correlation-matrices-in-r/

Value

A vector that represents how many other columns correlate to current column under the correlation cutoff.

Author(s)

Zuguang Gu <[e-mail]>

Examples

## Not run: 
mat = matrix(rnorm(100000 * 10), ncol = 100000, nrow = 20)
cor_columns(mat)

## End(Not run)
NULL

eilslabs/epic documentation built on May 16, 2019, 1:24 a.m.