mkatari-bioinformatics-august-2013-deseq
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revisionLast revisionBoth sides next revision | ||
mkatari-bioinformatics-august-2013-deseq [2015/06/17 08:19] – mkatari | mkatari-bioinformatics-august-2013-deseq [2015/08/21 13:53] – mkatari | ||
---|---|---|---|
Line 26: | Line 26: | ||
Then we will load the experimental design. An example is provided [[https:// | Then we will load the experimental design. An example is provided [[https:// | ||
< | < | ||
- | expdesign = read.table(pathToExpDesign) | + | expdesign = read.table(" |
</ | </ | ||
Line 53: | Line 53: | ||
An important part of DESeq is to estimate dispersion. This is simply a form of variance for the genes. | An important part of DESeq is to estimate dispersion. This is simply a form of variance for the genes. | ||
< | < | ||
+ | # if you have replicates do the following: | ||
cds = estimateDispersions( cds ) | cds = estimateDispersions( cds ) | ||
+ | ### HOWEVER, If you have NO replicates, then try this | ||
+ | cds = estimateDispersions( cds, method=" | ||
</ | </ | ||
Line 85: | Line 88: | ||
#To get the genes that have FDR of 10% and save it in the output directory. | #To get the genes that have FDR of 10% and save it in the output directory. | ||
< | < | ||
- | resSig | + | resSigind |
+ | resSigrep = res[ which(res$padj < 0.1 & res$log2FoldChange < -1), ] | ||
- | outfile | + | indoutfile |
+ | repoutfile = " | ||
- | write.table(resSig, | + | write.table(resSigind, |
- | | + | |
+ | sep=" | ||
+ | col.names=T, | ||
+ | row.names=F, | ||
+ | quote=F) | ||
+ | |||
+ | write.table(resSigrep, | ||
+ | repoutfile, | ||
sep=" | sep=" | ||
col.names=T, | col.names=T, |
mkatari-bioinformatics-august-2013-deseq.txt · Last modified: 2015/08/21 14:13 by mkatari