data(SEQC.Example)
exDat <- initDat(datType="count", isNorm=FALSE, exTable=MET.CTL.countDat,
filenameRoot="testRun", sample1Name="MET",
sample2Name="CTL", erccmix="RatioPair",
erccdilution=1/100, spikeVol=1, totalRNAmass=0.500,
choseFDR=0.1)
exDat <- est_r_m(exDat)
exDat <- dynRangePlot(exDat)
exDat <- geneExprTest(exDat)
exDat <- erccROC(exDat)
exDat <- estLODR(exDat, kind="ERCC", prob=0.9)
exDat <- annotLODR(exDat)
#to print 4 main plots to a single page pdf file
saveERCCPlots(exDat, plotsPerPg = "manuscript",saveas = "pdf")
#to print 4 plots to a jpeg file
saveERCCPlots(exDat, plotsPerPg = "manuscript",saveas = "jpeg")
# or to create a multiple page pdf of all plots produced
saveERCCPlots(exDat, plotsPerPg = "single", plotlist = exDat$Figures)
# or to create a multiple page pdf of just 2 plots
saveERCCPlots(exDat, plotsPerPg = "single",
plotlist = list(exDat$Figures$lodrPlot, exDat$Figures$maPlot))
Run the code above in your browser using DataLab