Matrix Analysis and Normalization by Concordant Information Enhancement
This function removes noise in the main matrix by utilizing information available from the supplementary matrix or summarized supplementary matrix.
The main matrix or data frame. Rows are features (genes/peaks/etc) and cols are samples (conditions/replicates)
The supplementary matrix or data frame.
The higher cutoff. See below for explanation.
The lower cutoff. See below for explanation.
If the supplementary dataset have the same genomic features on rows and samples on columns as
mat_main, it can be directly fed to
mancie. An example is RNA-Seq data of the same cell lines from two labs. If the supplementary dataset has different rows from
mat_main. It need to be first summarized using
summarize_mat to be compatible with
mat_main. An example is RNA-Seq data and DNase-seq data of the same tissue types.
The underlying rationale for using MANCIE is that the variation of genomic features in
mat_supp are concordant with and can be used to remove noise in the variation of genomic features in
(a) If the correlation between row i of
mat_main and row i of
mat_supp is larger than
cutoff1, the new row vector will be the first PC of the matrix formed by these two row vectors.
(b) If the correlation is between
cutoff2, the new row vector will be the weighted average of these two rows. The weight for row i of
mat_main is 1 and the weight for row i of
mat_supp is the correlation between these two row vectors.
(c)If the correlation is smaller than
cutoff2, the new row vector is the original row i of
There should be a reasonable portion of rows that fall into the first and second category. If not, the user should check if the data they would like to try MANCIE on really fits the aforementioned rationale. The user may also vary the default values of
cutoff2 if they see fit. The
mancie function will report percentage of rows falling into each category.
A modified matrix with the same dimensions as the main matrix
1 2 3
Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.