importC: Import high-htroughput 'C' data

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

View source: R/import.R

Description

Import 5C or Hi-C data from list file

Usage

1
2
importC(con, xgi.bed, ygi.bed=NULL, allPairwise=FALSE,
rm.trans=FALSE, lazyload=FALSE)

Arguments

con

input csv file. See details

xgi.bed

BED file describing the 'x' Intervals (i.e. column names) of the contact map

ygi.bed

BED file describing the 'y' intervals (i.e. row names) of the contact map

allPairwise

logical; generate all pairwise chromosomal contact maps, i.e chr1-chr2, chr2-chr1

rm.trans

if true, returns only intra-chromosomal maps

lazyload

logical; see details

Details

This function import high-throughput data from a tab list file.
The standard format for 5C/Hi-C data is the following :
** One list file (tab delimited) bin1 bin2 x12 bin1 bin3 x13 ... ** The BED file(s) describing the intervals ('xgi.bed' and 'ygi.bed' are usually the same for Hi-C but can be different for 5C data) chr1 1 1000000 bin1 chr1 1000001 2000000 bin2 ...

Note that this format is particularly interesting for sparse data as only non null values are stored.

The lazyload option allow to reduce the memory size of imported object. Therefore, only half of inter-chromosomal maps are stored. And intra-chromosomal maps are stored as sparse triangular matrix. Note that even if the contact maps are stored as triangular matrix, the indata method always returns a symmetrical map.

Value

A HTClist object

Author(s)

N. Servant

See Also

exportC,import.my5C, HTCexp-class

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
## Not run: 
data(Nora_5C)

## Data binning
E14.bin<-binningC(E14$chrXchrX)

## Export the new intervals definition
exportC(E14.bin, file="E14")

##Import
importC(con="E14.count", xgi="E14_xgi.bed", ygi="E14_ygi.bed")

## End(Not run)

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