getPearsonMap: Pearson correlation map

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/pca.R

Description

Generate pearson correlation map, usually used to call chromosomal compartments

Usage

1

Arguments

x

object that inherits from class HTCexp

normPerExpected

normalized by expected interaction using the loess calculation of distance dependency. see normPerExpected

center

default=true. center the observed/expected map before calculating the Pearson correlation

...

additional parameters passed to normPerExpected

Details

The function returns an HTCexp object with Pearson correlation map. This is usually the first step of the Principal Component Analysis (see pca.hic). Centering the rows of the observed/expected matrix allows to avoid bias to due to ranges of interaction counts. If true, the correlation of small values should be as valuable as correlation of large values

Value

A HTCexp object

Author(s)

N. Servant, B. Lajoie, R. McCord

See Also

normPerExpected

Examples

1
2
3
4
5
6
7
8
9
## Get Lieberman-Aiden Hi-C data
exDir <- system.file("extdata", package="HiTC")
l <- sapply(list.files(exDir, pattern=paste("HIC_gm06690_"), full.names=TRUE),
import.my5C)
hiC <- HTClist(l)

## get Pearson correlation map
pm <- getPearsonMap(hiC$chr14chr14)
mapC(HTClist(pm), maxrange=1, col.pos=c("black","red"), col.neg=c("black","blue"))

HiTC documentation built on Nov. 8, 2020, 8:27 p.m.