bnbc: Normalize Contact Matrices with BNBC
In bnbc: Bandwise normalization and batch correction of Hi-C data

Description Usage Arguments Details Value References See Also Examples

View source: R/bnbc.R

Applies BNBC method to normalize contact matrices.

1 2	bnbc(cg, batch, threshold = NULL, step = NULL, qn = TRUE, nbands = NULL, mod = NULL, mean.only = FALSE, tol = 5, bstart = 2, verbose = TRUE)

`cg`	A ContactGroup object.
`batch`	A single batch indicator variable.
`threshold`	The maximum distance interacting loci are allowed to be separated by.
`step`	The step size, or the number of bases a contact matrix cell represents.
`qn`	Whether to apply quantile normalization on each band matrix. Defaults to TRUE.
`bstart`	The first band to normalize. Defaults to 2.
`nbands`	The last band to normalize. Defaults to `nrow(cg) - 1`.
`mod`	A model matrix specifying which sample information is to be preserved by ComBat. Optional.
`mean.only`	Whether ComBat should not correct for batch effect in the variances of band matrix rows. Defaults to FALSE, which means variances are corrected. Set to TRUE if there is only one observation per batch.
`tol`	The number of significant digits for which the mean value of a band matrix must be greater than 0 to be processed by ComBat.
`verbose`	Should the function print progress?

Normalization and batch correction is performed in a band-wise manner, correcting all samples' observations of one matrix off-diagonal (which we refer to as a matrix "band") at a time. For each matrix band, we collect all samples' observations into a single matrix. We then apply quantile normalization to ensure distributional similarity across samples. Finally, we perform batch effect correction using ComBat on this matrix. Each samples' matrix band is then replaced with its corrected version. We refer to this process of Band-Wise Normalization and Batch Correction as BNBC.

This function applies BNBC to the set of contact matrices and returns a ContactGroup object with matrix bands bstart:nbands corrected. For those rows in the matrix bands which cannot be corrected we set all elements to 0.

Very high bands contain little data in Hi-C experiments, and we don't recommend to analyze those or apply this function to high bands, see the nbands argument to the function.

We recommend performing bnbc on contact matrices which have been converted to log-CPM and smoothed, see the example.

A ContactGroup object for which matrix bands bstart:nbands have had BNBC applied.

Johnson, W.E., Li, C. and Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 2007, 8:118-127. doi: 10.1093/biostatistics/kxj037

Fletez-Brant et al. Distance-dependent between-sample normalization for Hi-C experiments. In preparation.

ContactGroup, logCPM, boxSmoother, band

data(cgEx)
batches <- colData(cgEx)$Batch
cgEx.cpm <- logCPM(cgEx)
cgEx.smooth <- boxSmoother(cgEx, 5, mc.cores=1)
cgEx.bnbc <- bnbc(cgEx.smooth, batches, 1e7, 4e4, bstart=2, nbands=4)