# HiContacts
<!-- badges: start -->
[](https://lifecycle.r-lib.org/articles/stages.html#experimental)
[](https://github.com/js2264/HiContacts/actions/workflows/check-bioc.yml)
[](https://github.com/js2264/HiContacts/actions)
<!-- badges: end -->
HiContacts provides tools to investigate `(m)cool` matrices imported in R by `HiCExperiment`.
It leverages the `HiCExperiment` class of objects, built on pre-existing Bioconductor objects, namely `InteractionSet`, `GenomicInterations` and `ContactMatrix` (`Lun, Perry & Ing-Simmons, F1000Research 2016`), and provides **analytical** and **visualization** tools to investigate contact maps.
## Installation
`HiContacts` is available in Bioconductor. To install the current release, use:
```r
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("HiContacts")
```
To install the most recent version of `HiContacts`, you can use:
```r
install.packages("devtools")
devtools::install_github("js2264/HiContacts")
library(HiContacts)
```
## Citation
If you are using `HiContacts` in your research, please cite:
> Serizay J (2022). _HiContacts: HiContacts: R interface to cool files_.
> R package version 1.1.0
> <https://github.com/js2264/HiContacts>.
## How to use `HiContacts`
`HiContacts` includes a introduction vignette where its usage is
illustrated. To access the vignette, please use:
```r
vignette('HiContacts')
```
## Overview
### Import a .(m)cool file as `HiCExperiment`
```r
mcool_file <- HiContactsData::HiContactsData('yeast_wt', format = 'mcool')
range <- 'I:20000-80000' # range of interest
lsCoolResolutions(mcool_file, verbose = TRUE)
hic <- HiCExperiment::import(mcool_file, format = 'mcool', focus = range, resolution = 1000)
hic
```
### Plotting matrices
```r
plotMatrix(hic, use.scores = 'raw')
plotMatrix(hic, use.scores = 'balanced', limits = c(-4, -1))
```
### Distance law, a.k.a. P(s)
```r
hic <- Contacts(
mcool_file,
pairs = HiContactsData::HiContactsData('yeast_wt', format = 'pairs')
)
ps <- distanceLaw(hic)
plotPs(ps, aes(x = binned_distance, y = norm_p))
```
### Virtual 4C
```r
hic <- Contacts(mcool_file, res = 1000)
v4C <- virtual4C(hic, viewpoint = GRanges('V:150000-170000'))
gg4C(v4C, aes(x = center, y = score, col = chr)) +
facet_wrap(~chr, scales = 'free_x')
```
### Cis-trans ratios
```r
hic <- Contacts(mcool_file, res = 1000)
cisTransRatio(hic)
```
