splitIntoGroupsPseudo {edgeR} | R Documentation |
Given a pair of groups, split pseudodata for these groups
splitIntoGroupsPseudo(pseudo, group, pair)
pseudo |
data matrix to be split (e.g. quantile adjusted pseudodata) |
group |
group indicator, must be same length as ncol(pseudo) |
pair |
pair of groups to be split from the data |
list in which each element is a matrix of count data for an individual group
Davis McCarthy
# generate raw data from NB, create list object y<-matrix(rnbinom(80,size=1,mu=10),nrow=20) d<-DGEList(data=y,group=rep(1:2,each=2),lib.size=rep(c(1000:1001),2)) rownames(d$data)<-paste("tagno",1:nrow(d$data),sep=".") z<-splitIntoGroupsPseudo(d$data,d$group,pair=c(1,2))