Normalization of AlignmentsTrack / coverage plot by library size
1
0
Entering edit mode
wiedemak • 0
@wiedemak-14399
Last seen 7.0 years ago

Is it possible to normalize a coverage plot made by AlignmnetsTrack( file.bam) by the library size? I tried to use transform function in plotTracks( ), but it doesn't change the AlignmentsTrack coverage output.

Gviz normalization • 1.2k views
ADD COMMENT
1
Entering edit mode
Robert Ivanek ▴ 750
@robert-ivanek-5892
Last seen 12 months ago
Switzerland

Hi,

You did not provide any code. It is hard to guess where the problem could be. 

Typical mistake is not using full argument name. You need to provide the transformation function as argument called "transformation". The name has to be exact. For example this works for me:

library(Gviz)
afrom <- 2960000
ato <- 3160000
alTrackOrig <- AlignmentsTrack(system.file(package="Gviz", "extdata", "gapped.bam"), isPaired=TRUE, type="coverage")
alTrackTransformed <- AlignmentsTrack(system.file(package="Gviz", "extdata", "gapped.bam"), isPaired=TRUE, 
                           transformation=function(x) {x*10}, type="coverage")
plotTracks(list(alTrackOrig, alTrackTransformed), from=afrom, to=ato, chromosome="chr12")

 

ADD COMMENT

Login before adding your answer.

Traffic: 893 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6