splitIntoGroupsPseudo {edgeR}R Documentation

Split pseudodata according to group

Description

Given a pair of groups, split pseudodata for these groups

Usage

 
splitIntoGroupsPseudo(pseudo, group, pair)

Arguments

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

Value

list in which each element is a matrix of count data for an individual group

Author(s)

Davis McCarthy

Examples

# 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))

[Package edgeR version 1.2.4 Index]