mkatari-bioinformatics-august-2013-loadingphytozome
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionNext revisionBoth sides next revision | ||
mkatari-bioinformatics-august-2013-loadingphytozome [2014/01/06 14:37] – mkatari | mkatari-bioinformatics-august-2013-loadingphytozome [2015/06/18 07:24] – 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 25: | Line 25: | ||
===== Run Fisher Exact Test ===== | ===== Run Fisher Exact Test ===== | ||
- | If you do not have a sample gene list to work with, simply create a random one. | + | Let's use the significant genes identified from Deseq as our example. [[https:// |
< | < | ||
- | gene.list< | + | |
+ | gene.list< | ||
</ | </ | ||
- | Start by sourcing the [[https:// | + | Start by sourcing the [[https:// |
To run the function: | To run the function: | ||
Line 38: | Line 40: | ||
goterm_annot< | goterm_annot< | ||
row.names=1, | row.names=1, | ||
- | goresults=doFisherTest(gene.list, | + | goresults=doFisherTest(as.character(gene.list[,1]), goAnnot=goterm_annot) |
write.table(goresults, | write.table(goresults, | ||
</ | </ |
mkatari-bioinformatics-august-2013-loadingphytozome.txt · Last modified: 2015/06/18 09:05 by mkatari