Question: what do negative values from bigWigAverageOverBed output mean?
0
gravatar for ApoorvaB
2.8 years ago by
ApoorvaB210
United States
ApoorvaB210 wrote:

Hi,

This might seem basic, but I have a bed file with promoter regions I am interested in and Bigwig files from my ChIP-Seq experiment after subtracting input from treatment files using deeptools bamcompare. I want to know the coverage of my bw files. I came across the bigWigAverageOverBed utility from UCSC that I used but I am not sure I understand the results. Can anyone help me understand how this utility works and what negative values of the score mean?

Thanks

ADD COMMENTlink modified 2.8 years ago by Devon Ryan90k • written 2.8 years ago by ApoorvaB210
1
gravatar for Devon Ryan
2.8 years ago by
Devon Ryan90k
Freiburg, Germany
Devon Ryan90k wrote:

If the normalized input level is, on average, higher than that of the ChIP (at least in the regions of interest) then the average value will be negative.

ADD COMMENTlink written 2.8 years ago by Devon Ryan90k

Thank you Devon. I have data that looks like this

chr chr_start chr_end ENSEMBL Diff_cell undiff_cell

chr2 76542036 76825604 ENSMUSG00000051747 -3.09314 7.20316

-3.093 and 7.2 are average coverage values of the differentiated cell sample and undifferentiated cell sample. Now i need to know if this gene and other genes with similar values are enriched in any GO category. Would it be right to log transform the average coverage (form bigWigAverageOverBed), calculate the fold change and use the FC to create a ranked list of genes for GSEA ? In that case, how do i log transform large negative values ?

ADD REPLYlink modified 2.8 years ago • written 2.8 years ago by ApoorvaB210
1

The only way for fold-changes to make sense is if you put a floor of 0 on the values and add a small value to everything.

ADD REPLYlink written 2.8 years ago by Devon Ryan90k
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: 1053 users visited in the last hour