alpha.approxeb {edgeR} | R Documentation |
Estimate the prior weight, using an approximate empirical Bayes rule
alpha.approxeb(object, verbose=TRUE)
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 |
EBList
object with elements p
(overall proportion), p1
(estimates for first group), p2
(estimates for second group)
Mark Robinson
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)