Question: Reproducibility Exon imbalance
gravatar for int11ap1
5 weeks ago by
int11ap1320 wrote:

Dear all,

I am trying to reproduce the exon imbalance analysis from this paper (specifically, Figure 1C):

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 • 71 views
ADD COMMENTlink modified 5 weeks ago • written 5 weeks ago by int11ap1320
Please log in to add an answer.


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