Peak plots using R
Can any one tell me how to generate peak plots using R.
I have the chromosome start and end positions. Based on those peak start and end position. How can i generate a plot. Can anyone advice on it.
Thanks
• 5,886 views
•
link
updated
by
Ram
•
written
by
pinky_pinkpinky •
•
written
by
pinky_pinkpinky •
1 answer
You could pretty easily do this with either Gviz or ggbio (they're bioconductor packages). With ggbio, the general idea is to use the autoplot() function with a BED file of peaks that could just be subset (or split) by chromosome and the coordinate sorted and indexed BAM file with the alignments. See section 2.4 of the ggbio vignette for an example of this.
• 0 views
•
link
updated
by
Ram
•
written
by
Devon Ryan
•
written
by
Devon Ryan
Log in to answer this question.
More posts like this
-
How to plot chipseq peak relative to TSS for some genes ?
written by Sarah •
-
Creating a bed file from a de novo fasta file
written by c.e.chong •
-
Identify transcripts from genome positions
written by jomiles26 •
-
Moving average plot of gene expression ratio along the chromosome length
written by Susmita Mandal
-
Forum: Genetic distances Vs Physical positions
written by myprogramming2016 •
-
Color intensity plot peak positions in R from BED file
written by j.mccarter1 •
-
Bed file data Vs HOMER output
written by pinky_pinkpinky •
-
Generation of R Plots
written by pinky_pinkpinky •
-
Chip-Seq GO analysis
written by pinky_pinkpinky •
-
Plotting transposons start and end positions
written by pinky_pinkpinky •
Do you just want to plot the position of the peaks on the genome or their positions relative to some feature (e.g., the transcription start site)?
I want to plot the read counts with respective to the position of the peaks.
Within an individual peak or summed across all peaks? In the latter case, I assume you want the peak width to be normalized (otherwise it'd be difficult to compare things).
It's within an individual peak and plot has to be generated chromosome wise
Yaa I saw the manual. there are lot of information. Thanks.
But I would like to know how to go about with bed file.
What information should I have in bed file.
The bed file should just have the locations of the peaks.
I think I am doing something wrong.Can you please correct me.
f1.bamis just a character vector, trying to autoplot it like that won't make much sense. A better example would be something like:No, this isn't using a BED file to get the regions, but this should be a short enough example to get you started.