def runSailfishAddModels(infiles, outfiles): ''' Computes read counts across transcripts and genes based on a fastq file and an indexed transcriptome using Sailfish. Runs the sailfish "quant" function across transcripts with the specified options. Read counts across genes are counted as the total in all transcripts of that gene (based on the getTranscript2GeneMap table) ''' infiles, transcript2geneMap = infiles index, fastqfile = infiles transcript_outfile, gene_outfile = outfiles Quantifier = PipelineRnaseq.SailfishQuantifier( infile=fastqfile, transcript_outfile=transcript_outfile, gene_outfile=gene_outfile, annotations=index, job_threads=PARAMS["alignment_free_threads"], job_memory=PARAMS["sailfish_memory"], options=PARAMS["sailfish_options"], bootstrap=PARAMS["alignment_free_bootstrap"], libtype=PARAMS['sailfish_libtype'], transcript2geneMap=transcript2geneMap) Quantifier.run_all()
def runSailfish(infiles, outfiles): ''' Computes read counts across transcripts and genes based on a fastq file and an indexed transcriptome using Sailfish. Runs the sailfish "quant" function across transcripts with the specified options. Read counts across genes are counted as the total in all transcripts of that gene (based on the getTranscript2GeneMap table) Parameters ---------- infiles: list list with three components 0 - list of strings - paths to fastq files to merge then quantify across using sailfish 1 - string - path to sailfish index file 2 - string - path to table mapping transcripts to genes sailfish_threads: int :term: `PARAMS` the number of threads for sailfish sailfish_memory: str :term: `PARAMS` the job memory for sailfish sailfish_options: str :term: `PARAMS` string to append to the sailfish quant command to provide specific options, see http://sailfish.readthedocs.io/en/master/index.html sailfish_bootstrap: int :term: `PARAMS` number of bootstrap samples to run. Note, you need to bootstrap for differential expression with sleuth if there are no technical replicates. If you only need point estimates, set to 1. sailfish_libtype: str :term: `PARAMS` sailfish library type http://sailfish.readthedocs.io/en/master/library_type.html#fraglibtype outfiles: list paths to output files for transcripts and genes ''' fastqfile, index, transcript2geneMap = infiles transcript_outfile, gene_outfile = outfiles Quantifier = PipelineRnaseq.SailfishQuantifier( infile=fastqfile, transcript_outfile=transcript_outfile, gene_outfile=gene_outfile, annotations=index, job_threads=PARAMS["alignment_free_threads"], job_memory=PARAMS["sailfish_memory"], options=PARAMS["sailfish_options"], bootstrap=PARAMS["alignment_free_bootstrap"], libtype=PARAMS['sailfish_libtype'], transcript2geneMap=transcript2geneMap) Quantifier.run_all()