alpha.approxeb {edgeR}R Documentation

Estimate the prior weight, alpha

Description

Estimate the prior weight, using an approximate empirical Bayes rule

Usage

 
alpha.approxeb(object, verbose=TRUE) 

Arguments

object DGEList object containing the raw data with elements data (table of counts), group (vector indicating group) and lib.size (vector of library sizes)
verbose whether to write comments, default true

Value

EBList object with elements sigma2.0.est (numeric scale sigma_0^2 estimate), alpha (estimate for the prior weight, alpha), scores (likelihood scores), infos (Fisher expected information), quantileAdjusted (list from output of quantileAdjust)

Author(s)

Mark Robinson, Davis McCarthy

Examples

y<-matrix(rnbinom(20,size=1,mu=10),nrow=5)
d<-DGEList(data=y,group=rep(1:2,each=2),lib.size=rep(c(1000:1001),2))
alpha<-alpha.approxeb(d)

[Package edgeR version 1.2.4 Index]