getTx2gene: Retrieve the data frame linking genes to transcripts.

In statOmics/stageR: stageR: stage-wise analysis of high throughput gene expression data in R

Description

This functions returns a data frame that links the genes with the transcripts being analysed.

Usage

getTx2gene(object, ...)

## S4 method for signature 'stageRTx'
getTx2gene(object, ...)

Arguments

object	an object of the stageRTxClass class.
...	Additional arguments

Value

A matrix linking gene to transcript identifiers.

References

Van den Berge K., Soneson C., Robinson M.D., Clement L. (2017). stageR: a general stage-wise method for controlling the gene-level false discovery rate in differential expression and differential transcript usage. Genome Biology 18:151. https://doi.org/10.1186/s13059-017-1277-0

Examples

#make identifiers linking transcripts to genes
set.seed(1)
genes=paste0("gene",sample(1:200,1000,replace=TRUE))
nGenes=length(table(genes))
transcripts=paste0("tx",1:1000)
tx2gene=data.frame(transcripts,genes)
#gene-wise q-values
pScreen=c(seq(1e-10,1e-2,length.out=nGenes-100),seq(1e-2,.2,length.out=50),seq(50))
names(pScreen)=names(table(genes)) #discards genes that are not simulated
pConfirmation=matrix(runif(1000),nrow=1000,ncol=1)
rownames(pConfirmation)=transcripts
stageRObj <- stageRTx(pScreen=pScreen, pConfirmation=pConfirmation ,pScreenAdjusted=TRUE, tx2gene=tx2gene)
getTx2gene(stageRObj)

statOmics/stageR documentation built on Dec. 15, 2021, 2:58 a.m.