benchmarkWrapper: A wrapper that performs the whole benchmark analysis

In plger/RNAontheBENCH: RNAontheBENCH - Benchmark of RNAseq quantification and DEA

Description

Performs the whole series of benchmark analysis (see compareWithNanostring, analyzeSpikein, compareWithPCR, and compareSimulated). The function expects quantification files with the right column headers ("AJ80" and so on, or "s1","s2" and so on for the simulated data) to be in the rpath folder, and to bear some kind of recognizable name. However, to avoid confusion, we suggest using the filename 'transcripts.quant' for transcript-level quantification, 'genes.quant' for gene-level quantification (optional), and 'simulated.quant' for transcript-level quantification of the simulated dataset (optional). If you are benchmarking both the core (12-samples) dataset and the validation (6-samples) dataset, prefix the files with, respectively, 'w12.' and 'w6.' (see example below). If you instead want to specify files manually, use the individual underlying functions.

Usage

benchmarkWrapper(rpath, ANALYSIS_NAME, qt)

Arguments

rpath	The path were the quantification files are stored, and where the output files will be saved.
ANALYSIS_NAME	The name of the analysis pipeline
qt	A string indicating the unit of the expression matrix (either "FPKM", "TPM" or "COUNTS").

Value

Nothing, but saves a bunch of files in 'rpath' and opens an html page to browse the results.

Examples

# first we create a directory and put the example quantification file in it:
data(exampledata)
dir.create("example")
write.table(exampleTranscriptLevel,"example/w12.transcripts.quant",sep="\t",quote=FALSE)
write.table(exampleGeneLevel,"example/w12.genes.quant",sep="\t",quote=FALSE)
# run the wrapper, specifying that folder:
benchmarkWrapper("example", "tophat.featureCount", qt="COUNTS")

plger/RNAontheBENCH documentation built on May 25, 2019, 8:22 a.m.