Question: Reproducibility Exon imbalance
0
gravatar for int11ap1
7 months ago by
int11ap1360
Barcelona
int11ap1360 wrote:

Dear all,

I am trying to reproduce the exon imbalance analysis from this paper (specifically, Figure 1C): https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4184927/

Authors normalize raw coverages as described below (from supplementary materials): enter image description here

Therefore, considering a matrix of coverages where rows are exons of gene X and columns are samples, the R code to perform the equations A and B should be:

#equation A, norm across exons
coverage <- read.table("coverage.tsv", header=T, row.names=1, sep="\t", check.names=F)
coverage <- t(t(coverage*nrow(coverage))/colSums(coverage))
#equation B, norm across samples
coverage <- (coverage*ncol(coverage))/rowSums(coverage)
coverage <- log10(coverage)

Is that right?

The supplementary table 4 provides raw coverage for gene NOTCH1 and NOTCH2 and the normalized values they use to perform the exon imbalance. However, using the code written above, I do not obtain the normalized values they have.

Do you have any clue what I am doing wrong?

Thank you,

imbalance exon • 170 views
ADD COMMENTlink modified 7 months ago • written 7 months ago by int11ap1360
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 806 users visited in the last hour