## ----style, eval=TRUE, echo=FALSE, results='asis'------------------------
BiocStyle::latex()

## ----env, echo=FALSE,warning=FALSE,message=FALSE-------------------------
suppressPackageStartupMessages(library("PGA"))
#suppressPackageStartupMessages(library("R.utils"))

## ----bdb, eval=TRUE, warning=FALSE, message=FALSE------------------------
vcffile <- system.file("extdata/input", "PGA.vcf",package="PGA")
bedfile <- system.file("extdata/input", "junctions.bed",package="PGA")
gtffile <- system.file("extdata/input", "transcripts.gtf",package="PGA")
annotation <- system.file("extdata", "annotation",package="PGA")
outfile_path<-"db/"
outfile_name<-"test"
library(BSgenome.Hsapiens.UCSC.hg19)
dbfile <- dbCreator(gtfFile=gtffile,vcfFile=vcffile,bedFile=bedfile,
                    annotation_path=annotation,outfile_name=outfile_name,
                    genome=Hsapiens,outdir=outfile_path)

## ----databasesearching, echo=TRUE, cache=FALSE, tidy=FALSE,eval=TRUE, warning=FALSE, message=FALSE----
msfile <- system.file("extdata/input", "pga.mgf",package="PGA")
idfile <- runTandem(spectra = msfile, fasta = dbfile, outdir = "./", cpu = 6,
                    enzyme = "[KR]|[X]", varmod = "15.994915@M",itol = 0.05,
                    fixmod = "57.021464@C", tol = 10, tolu = "ppm",
                    itolu = "Daltons", miss = 2, maxCharge = 8, ti = FALSE)


## ----parserGear, echo=TRUE, cache=FALSE, tidy=FALSE, eval=TRUE, warning=FALSE, message=FALSE----
parserGear(file = idfile, db = dbfile, decoyPrefix="#REV#",xmx=1,thread=8,
           outdir = "parser_outdir")


## ----mascotParser, eval=FALSE, echo=TRUE, cache=FALSE, tidy=FALSE, warning=FALSE, message=FALSE----
## dat_file<-"mascot_raw.dat"
## parserGear(file = dat_file, db = dbfile, decoyPrefix="#REV#",xmx=1,thread=8,
##            outdir = "parser_outdir")

## ----reportg, echo=TRUE, cache=FALSE, tidy=FALSE, eval=TRUE, warning=FALSE, message=FALSE----
reportGear(parser_dir = "parser_outdir", tab_dir = outfile_path,
           report_dir = "report")

## ----auto, echo=TRUE, cache=FALSE, tidy=FALSE, eval=TRUE, warning=FALSE, message=FALSE----
vcffile <- system.file("extdata/input", "PGA.vcf",package="PGA")
bedfile <- system.file("extdata/input", "junctions.bed",package="PGA")
gtffile <- system.file("extdata/input", "transcripts.gtf",package="PGA")
annotation <- system.file("extdata", "annotation",package="PGA")
library(BSgenome.Hsapiens.UCSC.hg19)
msfile <- system.file("extdata/input", "pga.mgf",package="PGA")
easyRun(gtfFile=gtffile,vcfFile=vcffile,bedFile=bedfile,spectra=msfile,
        annotation_path=annotation,genome=Hsapiens,cpu = 6,
        enzyme = "[KR]|[X]", varmod = "15.994915@M",itol = 0.05,
        fixmod = "57.021464@C", tol = 10, tolu = "ppm", itolu = "Daltons",
        miss = 2, maxCharge = 8, ti = FALSE,xmx=1)

## ----sessioninfo, results='asis', echo=FALSE-----------------------------
toLatex(sessionInfo())