I have a BED file from ChIPseq as a Granges object with associated scores values as metadata. In order to plot this region (1Mb) with Gviz, I need firstly to normalize the scores by a normalization factor calculated from each enriched peak present in the 1Mb region (ratio before and after normalization with MAnorm). This is required because I need to compare acetylation marks in two conditions. My question is: how could I multiply every score present in a given interval of the Granges object by the corresponding normalization factor for that interval? Then I should repeat this operation for each interval and associated normalization factor until I have covered the entire 1Mb region.
1 answer
I got it.
Firstly I had to subset by interval, then simply multiply the scores to the corresponding normalization factor. Finally create a grangesList comprising the new objects and turn it again to granges for visualization.
Log in to answer this question.
-
NA values in conumee detail
written by sativus •
-
Seeking feedback on ChIP-seq normalization method: Calculating scaling factors by dividing input by…
written by kalavattam •
-
Making use of relative density of genomic features between gene sets
written by Charles Murtaugh •
-
Differential Peak Analysis for a specific gene using only bigWig + BED files?
written by Maycon •
-
help for dba.count settings for RNAPII-ChIP using GRanges peakset
written by bertb •
-
Merge several bigWig files (with different intervals) into single GRanges object
written by Roman Hillje •
-
how to calculate with metadata of GRanges Objects
written by Assa Yeroslaviz
-
Plotting confidence intervals with NA values
written by viniciushs.z •
-
Filtering de novo transcriptomes - finding transcripts present in N samples
written by m.fletcher •
-
Which statistical test should be used to test for enrichment of gene lists?
written by Ashutosh •
