getTx2gene: Retrieve the data frame linking genes to transcripts.

Description Usage Arguments Value References Examples

Description

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

Usage

1
2
3
4
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
#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.