Search
Question: CHiPpeakAnno GO and/or KEGG results presenting
0
2.3 years ago by
lokapal0
lokapal0 wrote:

Hello, developers and Bioconductor experts! :-)

I'm trying to annotate my genomic features track (bed format) and build for annotated data gene ontology and/or KEGG pathway. Right now I'm using CHiPpeakAnno package. All is working as it should be (as it seems to me). My problem is that I cannot use the results and manual, unfortunately, doesn't cover it. I loaded my dataset to GRanges, annotated it (by the way, what database is the better for H. sapiens annotations?) and built GO with function getEnrichedGO. How can I use the full results? I mean how to build from "over" object/result something like to this figure:


modified 2.3 years ago • written 2.3 years ago by lokapal0
0
2.3 years ago by
Julie Zhu3.8k
United States
Julie Zhu3.8k wrote:

Hello, Julie!

Thanks a lot for prompt answer. If I сould can give you advice concerning CHiPpeakAnno documentation - it would be wise to create section "Working Recipes" and, probably, reorganize "User guide" according to its appearance. It's especially important for annotation and further getEnrichedGO and/or getEnrichedPATH processing because this functions contain some non-obvious parameters (e.g. multiAdjMethod description is as follows: "The multiple testing procedures, for details, see mt.rawp2adjp in multtest package". From mt.rawp2adjp description we can see that "This function computes adjusted p-values for simple multiple testing procedures from a vector of raw (unadjusted) p-values". Right now it is unclear should we use this method or not and in what cases (I can only guess that this option should be used in the case of mix in one GRange object different genomic tracks). I only know that old CHiPpeakAnno description (e.g. from Bioconductor 2.11) contains the other usage example:

enrichedGO = getEnrichedGO (annotatedPeak, orgAnn="org:Hs:eg:db", maxP=0.01, multiAdj=TRUE, minGOterm=10, multiAdjMethod="BH")

So multiAdj option changed its default behaviour. Why? What should we use?

As for total usage it is not obvious from the manual that there are two (or even three if biomart is included) quite different workflows. As a matter of fact single mention in the user guide about it is the following statement: "If you are using TxDb as annotation data, please try to change feature_id_type to “entrez_id”. Moreover it would be good to have some very brief and concise comparison of different databases (as for your module, of course). It would be good to have in the "Recipes" section completely working examples for each workflow.

I do understand that all my questions and problems looks like quite trivial to you, but I assure you that not all CHiPpeakAnno users are hard gene annotators and pathway investigators/creators. We are not so cool in this area of science so we would like to have a bit of help in the manual.

Anyway thanks for very useful package!