fdr: Compute the fdr (q-value) from posterior probabilities

In RGLab/MIMOSA: Mixture Models for Single-Cell Assays

Description

Given the z's from a MIMOSA model, calculates the q-values for each observation.

Usage

fdr(z)

## S3 method for class 'matrix'
fdr(z)

## S3 method for class 'MIMOSAResult'
fdr(z)

## S3 method for class 'MIMOSAResultList'
fdr(z)

Arguments

z	matrix of posterior probabilties, or a MIMOSAResult, or MIMOSAResultList

Value

a vector of q-values or a list of vectors of q-values.

Examples

data(ICS)
E<-ConstructMIMOSAExpressionSet(ICS,
                              reference=ANTIGEN%in%'negctrl',measure.columns=c('CYTNUM','NSUB'),
                               other.annotations=c('CYTOKINE','TCELLSUBSET','ANTIGEN','UID'),
                               default.cast.formula=component~UID+ANTIGEN+CYTOKINE+TCELLSUBSET,
                               .variables=.(TCELLSUBSET,CYTOKINE,UID),
                               featureCols=1,ref.append.replace='_REF')
result<-MIMOSA(NSUB+CYTNUM~UID+TCELLSUBSET+CYTOKINE|ANTIGEN,
   data=E, method='EM',
   subset=RefTreat%in%'Treatment'&ANTIGEN%in%'ENV',
   ref=ANTIGEN%in%'ENV'&RefTreat%in%'Reference')
qvalues<-fdr(result)

RGLab/MIMOSA documentation built on Nov. 13, 2020, 5:04 a.m.