mkatari-bioinformatics-august-2013-loadingphytozome
Differences
This shows you the differences between two versions of the page.
Next revision | Previous revision | ||
mkatari-bioinformatics-august-2013-loadingphytozome [2013/08/21 09:26] – created mkatari | mkatari-bioinformatics-august-2013-loadingphytozome [2015/06/18 09:05] (current) – [Run Fisher Exact Test] mkatari | ||
---|---|---|---|
Line 5: | Line 5: | ||
===== Read Phytozome annotation file. ===== | ===== Read Phytozome annotation file. ===== | ||
- | This step was tricky because there were characters such as quotes and hash marks. We need to reset the default values to allow the file to load without a problem. | + | This step was tricky because there were characters such as quotes and hash marks. We need to reset the default values to allow the file to load without a problem. |
< | < | ||
- | anno< | + | anno< |
| | ||
| | ||
Line 22: | Line 22: | ||
source(" | source(" | ||
</ | </ | ||
+ | |||
+ | ===== Run Fisher Exact Test ===== | ||
+ | |||
+ | Let's use the significant genes identified from Deseq as our example. [[https:// | ||
+ | |||
+ | < | ||
+ | |||
+ | gene.list< | ||
+ | |||
+ | </ | ||
+ | |||
+ | Start by sourcing the [[https:// | ||
+ | |||
+ | To run the function: | ||
+ | < | ||
+ | source(" | ||
+ | goterm_annot< | ||
+ | row.names=1, | ||
+ | goresults=doFisherTest(as.character(gene.list[, | ||
+ | write.table(goresults, | ||
+ | </ | ||
+ | |||
+ | Visualizing your results in the GO-term graph. Use the RamiGO package. Simply provide the GO terms that are enriched and the corresponding p-value for the go-terms. By default it will save a file to your working directory called RamiGO.png | ||
+ | |||
+ | < | ||
+ | source(" | ||
+ | biocLite(" | ||
+ | library(RamiGO) | ||
+ | getAmigoTree(rownames(gotermresults)[1: | ||
+ | | ||
+ | </ | ||
+ |
mkatari-bioinformatics-august-2013-loadingphytozome.1377077181.txt.gz · Last modified: 2013/08/21 09:26 by mkatari