GetConfidentCalls: Obtain confident gene calls for classifying genes into...

Description Usage Arguments Details Value Note Author(s) Examples

View source: R/GetConfidentCalls.R

Description

Obtain confident gene calls for classifying genes into expression paths

Usage

1
GetConfidentCalls(EBSeqHMMOut, FDR=.05, cutoff=0.5, OnlyDynamic=TRUE,Paths=NULL)

Arguments

EBSeqHMMOut

output from EBSeqHMMTest function

FDR

Target FDR, default is 0.05.

cutoff

cutoff to use for defining a confident call. Genes with PP_path greater or equal to cutoff will be called as a confident call. Default is 0.5.

OnlyDynamic

if specifies as T, only dynamic paths will be shown

Paths

paths that are of interest. Default is NULL. If it is not specified, all possible paths will be considered.

Details

Function GetConfidentCalls() can be used to obtain a list of DE genes/isoforms with user specific cutoffs. To obtain a list of DE genes/isoforms with a target FDR alpha, the user may specify FDR=alpha. To further choose genes/isoforms with high posterior probability of being its most likely path, the user may specify the option cutoff (default is 0.5). Then genes or isoforms with PP(most likely path ) > = 0.5 will be selected

Value

Overall: a list of genes/isoforms that are identified as DE under the target FDR, shown are their names and PPs; EachPath: a list object, each sublist contains confident calls (genes/isoforms) that have PP(path)>=cutoff for a particular expression path, shown are their names and PPs; NumEach: length of each sublist in EachPath. EachPathName: gene/isoform names in each of the sublists in EachPath

Note

Output: output a list of genes that are classified to a expression path as a confident assignment.

Author(s)

Ning Leng

Examples

1
2
3
4
5
6
7
8
data(GeneExampleData)
CondVector <- rep(paste("t",1:5,sep=""),each=3)
Conditions <- factor(CondVector, levels=c("t1","t2","t3","t4","t5"))
Sizes <- MedianNorm(GeneExampleData)
EBSeqHMMGeneOut <- EBSeqHMMTest(Data=GeneExampleData, sizeFactors=Sizes, Conditions=Conditions,
          UpdateRd=2)
GeneDECalls <- GetDECalls(EBSeqHMMGeneOut, FDR=.05)
GeneConfCalls <- GetConfidentCalls(EBSeqHMMGeneOut, FDR=.05,cutoff=.5, OnlyDynamic=TRUE)

EBSeqHMM documentation built on Nov. 1, 2018, 2:35 a.m.