# (1) simulating expression data: 100 genes, 12 samples
eset <- make.example.data(what="eset")
# plot heatmap
exprs.heatmap(expr=exprs(eset), grp=as.factor(pData(eset)$GROUP))
# (2) DE analysis
eset <- de.ana(eset)
pdistr(fData(eset)$ADJ.PVAL)
volcano(fc=fData(eset)$FC, p=fData(eset)$ADJ.PVAL)
Run the code above in your browser using DataLab