stringtie_assemble: Assemble transcripts with StringTie

View source: R/05_04_assembly_quantification_stringtie.R

stringtie_assembleR Documentation

Assemble transcripts with StringTie

Description

Assemble transcripts with StringTie

Usage

stringtie_assemble(
  sample_info = NULL,
  qc_table = NULL,
  mappingdir = "results/04_read_mapping",
  gff_path = NULL,
  stringtiedir = "results/05_quantification/stringtie",
  threads = 1
)

Arguments

sample_info

Data frame of sample metadata created with the functions create_sample_info and infer_strandedness. The function infer_strandedness adds a column named "Orientation" with library strandedness information, which is mandatory.

qc_table

Data frame of fastp summary statistics as returned by summary_stats_fastp().

mappingdir

Directory where .bam files are stored.

gff_path

Path to GFF/GTF file with annotations.

stringtiedir

Directory where StringTie output files will be stored.

threads

Number of threads to use. Default: 1.

Value

A 2-column data frame with samples in the first column and status in the second column, with "OK" if transcripts were assembled, and NA otherwise.

Examples

data(sample_info)
qc_table <- summary_stats_fastp(system.file("extdata", package = "bears"))
mappingdir <- system.file("extdata", package="bears")
gff_path <- system.file("extdata", "Homo_sapiens.GRCh37.75_subset.gtf",
                        package="bears")
stringtiedir <- tempdir()
if(stringtie_is_installed()) {
    a <- stringtie_assemble(sample_info, qc_table, mappingdir, 
                            gff_path, stringtiedir)
}

almeidasilvaf/gead documentation built on April 26, 2023, 3:47 p.m.