AggrHiCExperiment
S4 class
Source: R/AggrHiCExperiment-class.R
, R/AggrHiCExperiment-methods.R
AggrHiCExperiment.Rd
The AggrHiCExperiment
extends HiCExperiment
class.
Usage
AggrHiCExperiment(
file,
resolution = NULL,
targets,
flankingBins = 50,
metadata = list(),
topologicalFeatures = S4Vectors::SimpleList(),
pairsFile = NULL,
bed = NULL,
maxDistance = NULL,
BPPARAM = BiocParallel::bpparam()
)
# S4 method for class 'AggrHiCExperiment,missing'
slices(x)
# S4 method for class 'AggrHiCExperiment,character'
slices(x, name)
# S4 method for class 'AggrHiCExperiment,numeric'
slices(x, name)
# S4 method for class 'AggrHiCExperiment'
show(object)
Arguments
- file
CoolFile or plain path to a Hi-C contact file
- resolution
Resolution to use with the Hi-C contact file
- targets
Set of chromosome coordinates for which interaction counts are extracted from the Hi-C contact file, provided as a GRanges object (for diagnoal-centered loci) or as a GInteractions object (for off-diagonal coordinates).
- flankingBins
Number of bins on each flank of the bins containing input targets.
- metadata
list of metadata
- topologicalFeatures
topologicalFeatures provided as a named SimpleList
- pairsFile
Path to an associated .pairs file
- bed
Path to regions file generated by HiC-Pro
- maxDistance
Maximum distance to use when compiling distance decay
- BPPARAM
BiocParallel parameters
- x, object
A
AggrHiCExperiment
object.- name
The name/index of slices to extract.
Slots
fileName
Path of Hi-C contact file
resolutions
Resolutions available in the Hi-C contact file.
resolution
Current resolution
interactions
Genomic Interactions extracted from the Hi-C contact file
scores
Available interaction scores.
slices
Available interaction slices.
topologicalFeatures
Topological features associated with the dataset (e.g. loops (\<Pairs\>), borders (\<GRanges\>), viewpoints (\<GRanges\>), etc...)
pairsFile
Path to the .pairs file associated with the Hi-C contact file
metadata
metadata associated with the Hi-C contact file.
Examples
fpath <- HiContactsData::HiContactsData('yeast_wt', 'mcool')
#> see ?HiContactsData and browseVignettes('HiContactsData') for documentation
#> loading from cache
data(centros_yeast)
x <- AggrHiCExperiment(
file = fpath,
resolution = 8000,
targets = centros_yeast[c(4, 7)]
)
#> Going through preflight checklist...
#> Parsing the entire contact matrice as a sparse matrix...
#> Modeling distance decay...
#> Filtering for contacts within provided targets...
x
#> `AggrHiCExperiment` object over 2 targets
#> -------
#> fileName: "/github/home/.cache/R/ExperimentHub/190530f4def5_7752"
#> focus: 2 targets
#> resolutions(5): 1000 2000 4000 8000 16000
#> active resolution: 8000
#> interactions: 10201
#> scores(4): count balanced expected detrended
#> slices(4): count balanced expected detrended
#> topologicalFeatures: targets(2)
#> pairsFile: N/A
#> metadata(0):
slices(x, 'count')[1:10, 1:10, 1]
#> [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
#> [1,] 1364 1054 679 470 571 450 230 132 135 83
#> [2,] 0 832 1430 646 617 432 232 138 142 90
#> [3,] 0 0 1457 1575 864 603 290 227 219 123
#> [4,] 0 0 0 1227 1398 634 356 210 218 109
#> [5,] 0 0 0 0 1285 1544 678 400 361 226
#> [6,] 0 0 0 0 0 1261 1479 591 649 381
#> [7,] 0 0 0 0 0 0 1186 1349 885 433
#> [8,] 0 0 0 0 0 0 0 1252 1561 489
#> [9,] 0 0 0 0 0 0 0 0 1245 1313
#> [10,] 0 0 0 0 0 0 0 0 0 1356