parserGear: Post-processing for the identification result

In wenbostar/PGA: An package for identification of novel peptides by customized database derived from RNA-Seq

Post-processing for the identification result

Description

This function is mainly for q-value calculation, protein inference and novel peptides spectra annotation.

Usage

parserGear(file = NULL, db = NULL, outdir = "parser_outdir",
  prefix = "pga", fdr = 0.01, novelPrefix = "VAR",
  decoyPrefix = "###REV###", alignment = 1, xmx = NULL, thread = 1,
  msfile = NULL, verbose = 1)

Arguments

file	MS/MS search file. Currently, only XML format file of X!Tandem, DAT result of Mascot and mzIdentML result file (MS-GF+, MyriMatch, IPeak, OMSSA, ...) are supported.
db	A FASTA format database file used for MS/MS searching. Usually, it is from the output of the function dbCreator.
outdir	Output directory.
prefix	The prefix of output file.
fdr	FDR for peptide identification. Default is 0.01 at PSM level.
novelPrefix	The prefix of novel protein ID. Default is "VAR". "VAR" is the prefix which used by function dbCreator. This value should be left to the default when your database is constructed by the function getTrinityDB.
decoyPrefix	The prefix of decoy sequences ID. Default is "###REV###". "###REV###" is the prefix which used by function dbCreator.
alignment	0 or 1 to determine if peptide should be alignment or not. Default is 1.
xmx	The maximum Java heap size. The unit is "G".
thread	This parameter is used to specify the number of threads. "0" represents that all of the available threads are used; "1" represents one thread is used; "2" represents two threads are used,and so on. Default is 1.
msfile	The MS/MS data (mgf format) used.

Value

none

Examples

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)

msfile <- system.file("extdata/input", "pga.mgf",package="PGA")

## X!Tandem as the peptide identification software
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(file = idfile, db = dbfile, decoyPrefix="#REV#",xmx=1,thread=8,
           outdir = "parser_outdir")

## Mascot as the peptide identification software
dat_file <- system.file("extdata/input", "mascot.dat",package="PGA")
parserGear(file = dat_file, db = dbfile, decoyPrefix="#REV#",xmx=1,thread=8,
           outdir = "parser_outdir_mascot")

wenbostar/PGA documentation built on March 24, 2022, 6:48 p.m.