sleuth_gene_table: Create a gene table from a sleuth object

In pachterlab/sleuth: Tools for investigating RNA-Seq

Create a gene table from a sleuth object

Description

Take in a sleuth object with added genes, and return a data table in which genes list the most significant transcript mapping to themselves

Usage

sleuth_gene_table(obj, test, test_type = "lrt", which_model = "full",
  which_group = "ens_gene")

Arguments

obj	a sleuth object
test	a character string denoting which beta to use
test_type	either 'wt' for Wald test or 'lrt' for likelihood ratio test
which_model	a character string denoting which model to use
which_group	a character string denoting which gene group to use

Value

a data.frame with the following columns

gene name; if ext_gene name specified, it will be legible gene name. If ens_gene name, it will be an Ensemble gene (assuming user followed the vignette)

most_sig_trancript: Most significant transcript for the given gene

pval: p-value for the test chosen

qval: False discovery rate normalized p-value (Benjamini-Hochberg, see: p.adjust)

num_transcripts: Total number of transcripts for the given gene

list_of_transcripts: All transcripts associated with this gene

Examples

sleuth_genes <- sleuth_gene_table(so, 'conditionIP', test_type ='wt',
                                  which_group = 'ext_gene')
head(sleuth_genes) # show info for first 5 genes
sleuth_genes[1:5, 6] # show transcripts for first 5 genes

pachterlab/sleuth documentation built on Nov. 17, 2022, 4:51 p.m.