cts <- matrix(data=rnbinom(28,2,.4), ncol=4, nrow=7)
colnames(cts) <- apply(expand.grid(c('a','b'),1:2), 1, paste, collapse='_')
x <- manta(cts, makeSampleDF(cts, group=rep(c('a','b'),2)))
collapseRepliCounts(x, pair=c('a','b'))
Run the code above in your browser using DataLab