Enriched GO questions
0
0
Entering edit mode
Julie Zhu ★ 4.3k
@julie-zhu-3596
Last seen 5 months ago
United States
Dear Christopher, For the dataset you sent to me, you could use read.table in R and then use the 5th column (x[,5]) to parse the EnsemblID with substr or strsplit (https://stat.ethz.ch/pipermail/r-help/2008-October/178335.html). If you follow the following steps to obtain annotated overlapping peaks, you will have a RangedData to be directly plugged into getEnrichedGO function at step 3. 1. Apply findOverlappingPeaks to find overlapping peaks 2. Apply annotatePeakInBatch to annotate the overlapping peaks 3. Apply getEnrichedGO to the annotated peaks from step 2. For detailed information, please refer to http://bioconductor.org/pack ages/2.6/bioc/vignettes/ChIPpeakAnno/inst/doc/ChIPpeakAnno.pdf and htt p://bioconductor.org/packages/2.6/bioc/vignettes/ChIPpeakAnno/inst/doc /ChIPpeakAnno.pdf Best regards, Julie On 2/16/11 12:29 PM, "Christopher Ricupero" <ricupero@eden.rutgers.edu> wrote: Dear Julie, I have a quick question that I think will be very easy for you to answer. I have been making great progress using ChIPPeak Anno and the overlapping functions. Through this process I have generated many Overlapping and Unique Peak datasets. However as you can see in the attached excel file example, I have unique peaks that do contain the Ensemble IDs needed for Go Annotation, but they are not contained separately. Through my comparison, the names contain the EnsemblID within the peak IDs (Ex –“ Conrol K9 Dy1-chr1_peak12171 ENSRNOG00000030145-Conrol K9 Dy3-chr1_peak22483 ENSRNOG00000030145” However, before I run the go annotation, I will need to parse out the Ensemble ID. I can do this using the right function in Excel ‘=RIGHT(E2,18)’, since all Ensemble Ids are 18 characters. However, I would rather just complete this in R, by accessing the ranged data file and the column names that contains the Ensembl ID. Can you tip me off on how I can just access the ‘names’ column and then just parse out the EnsemblID? I need this so I can then run the function: getEnrichedGo. Thank you Christopher Ricupero Rutgers University [[alternative HTML version deleted]]
Annotation GO annotate PROcess Annotation GO annotate PROcess • 754 views
ADD COMMENT

Login before adding your answer.

Traffic: 655 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