dbCreator: Create customized protein database from RNA-Seq data

View source: R/dbCreator.R

dbCreatorR Documentation

Create customized protein database from RNA-Seq data

Description

The main function to create customized protein database from RNA-Seq data

Usage

dbCreator(gtfFile = NULL, vcfFile = NULL, bedFile = NULL,
  tabFile = NULL, annotation_path = NULL, outdir, outfile_name,
  lablersid = FALSE, COSMIC = FALSE, bool_get_longest = TRUE,
  organism = "Homo sapiens", make_decoy = TRUE, genome = NULL,
  debug = FALSE, ...)

Arguments

gtfFile

A GTF format file containing novel transcripts information

vcfFile

A VCF format file containing SNV and INDEL information

bedFile

A BED format file from Tophat2 containing juction information.

tabFile

A Tab format file from Hisat2 containing juction information (e.g. splicesites.tab). As HISAT2 is a successor to both HISAT and TopHat2, we recommend that users switch to this option.

annotation_path

This directory contains numerous pieces of genome annotation information which can be downloaded by PrepareAnnotationEnsembl2 or PrepareAnnotationRefseq2.

outdir

Output directory

outfile_name

Output file name

lablersid

A logical indicating whether to do the SNV annotation(dbSNP)

COSMIC

A logical indicating whether to do the SNV annotation(COSMIC)

bool_get_longest

When it's set as TRUE, the longest sequences will be retained after the DNA sequences are six-frame translated into protein sequences. Otherwise, the protein sequences more than 30 aa are retained.

organism

What is the Genus and species of this organism.Please use proper scientific nomenclature for example: "Homo sapiens" and not "human", default is "Homo sapiens".

make_decoy

A logical indicating whether to add the decoy sequences

genome

Genome information. This is a BSgenome object(e.g. Hsapiens). Default is NULL.

debug

Output the debug info (Default is <False>).

...

Additional arguments

Value

The database file

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)
dbCreator(gtfFile=gtffile,vcfFile=vcffile,bedFile=bedfile,
          annotation_path=annotation,outfile_name=outfile_name,
          genome=Hsapiens,outdir=outfile_path)

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