It contains my assignment in Bioinformatics course talking about Differential Expression and this semesters’ code
My assignment in bioinformatics course
Codes using different packages producing differential expression files
1 | setwd('61samplesmatrix') |
some plot from DESeq2
Here is the pdf version of the following codes and its output.
If you can’t see pdf, please try this more interactive way by click here
1 | ### DESeq2 |
the codes use DESeq2 and produces files:
05.diffexp/hcc_example.miRNA.homer.DESeq2.rlog.mx
05.diffexp/hcc_example.miRNA.NCvsHCC.DESeq2.tsv
edgeR
1 | # 0.filter expression matrix |
the codes use edgeR and produces files:
05.diffexp/hcc_example.miRNA.NCvsHCC.edgeR.classic.tsv
05.diffexp/hcc_example.miRNA.homer.edgeR.TMM.mx
05.diffexp/hcc_example.miRNA.NCvsHCC.edgeR.tsv
Wilcox/Mann-Whitney-U Test
1. normalize the reads by library size (edgeR)
2. identify differential expressed gene using wilcoxon.test()
1 | cpmMx <- read.table("04.counts/hcc_example.miRNA.homer.rpm.mx",sep="\t",header=T) |
the codes use wilconx test and produces files:
05.diffexp/hccexample.miRNA.NCvsHCC.wilcox.tsv
exploring results
1 | # filter expression matrix |
Heatmap for count matrix
1 | library("pheatmap") |