Question: Visualizing average read coverage across features with varying length
1
gravatar for rxzlmn
5.7 years ago by
rxzlmn10
Singapore
rxzlmn10 wrote:

Hi,

I am trying to figure out a way to plot read densities over features of varying length. I have a set of ~100 genomic features that are between ~3kb to 7kb in length, and I want to plot per-base coverage or raw read coverage across those elements on the same x scale, i.e. from 0 to 100% of feature length.

The only tool so far I found is ngsplot, but it is kind of limited as it is stand-alone and does not have many options.

Is there a way to do that with existing packages or base tools in R?

rna-seq R • 2.0k views
ADD COMMENTlink modified 5.7 years ago by dariober11k • written 5.7 years ago by rxzlmn10
1

For "base tools", are we including just writing your own functions? This shouldn't be too terrible to do yourself (it's just a matter of smoothing things so interpolation makes sense). Have a look at the approx() function.

ADD REPLYlink written 5.7 years ago by Devon Ryan96k

I think that's exactly what I was looking for! Many thanks

ADD REPLYlink written 5.7 years ago by rxzlmn10
2
gravatar for dariober
5.7 years ago by
dariober11k
WCIP | Glasgow | UK
dariober11k wrote:

Have a look at theĀ  deeptools suite in particular computeMatrix.

something like:

computeMatrix scale-regions -R exons.bed -S coverage.bw -out exons.mat.gz --outFileNameData profile.tab -B 100 -A 100

will put in profile.tab the raw numbers you can use in R.

ADD COMMENTlink modified 5.7 years ago • written 5.7 years ago by dariober11k

Thanks, I'll try that!

ADD REPLYlink written 5.7 years ago by rxzlmn10
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: 1578 users visited in the last hour