mkatari-bioinformatics-august-2013-deseq
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-deseq [2013/08/23 16:20] – mkatari | mkatari-bioinformatics-august-2013-deseq [2015/06/17 06:14] – mkatari | ||
---|---|---|---|
Line 9: | Line 9: | ||
source(" | source(" | ||
biocLite(" | biocLite(" | ||
- | </ | ||
- | |||
- | Once we have installed the DESeq package in our local R repository we need to let our script know to look in this directory so if anyone wants to run our script, it will know exactly where to look for DESeq. R uses a variable (.libPaths) to store locations where it should look for packages. By default R will look in the global R library and the user's local path (if one exists). We will simply add our R library path to this variable. The other option is to tell the system administrator to add the packages to the global library path. This is done in the following lines of the code: | ||
- | |||
- | < | ||
- | mannyLibPaths = '/ | ||
- | .libPaths(new='/ | ||
</ | </ | ||
Line 95: | Line 88: | ||
< | < | ||
- | pdf(" | + | maFile = paste(pathToOutputDir, |
+ | pdf(maFile) | ||
plotMA(res) | plotMA(res) | ||
dev.off() | dev.off() | ||
</ | </ | ||
- | #To get the genes that have FDR of 10% | + | #To get the genes that have FDR of 10% and save it in the output directory. |
+ | < | ||
+ | resSig = res[ which(res$padj < 0.1), ] | ||
+ | |||
+ | outfile = paste(pathToOutputDir," | ||
write.table(resSig, | write.table(resSig, | ||
- | | + | |
sep=" | sep=" | ||
col.names=T, | col.names=T, | ||
row.names=F, | row.names=F, | ||
quote=F) | quote=F) | ||
+ | </ | ||
+ | |||
mkatari-bioinformatics-august-2013-deseq.txt · Last modified: 2015/08/21 14:13 by mkatari