R/designCount.R

In sdef: Synthesizing List of Differentially Expressed Features

designCount <-
function(array,design) {

sum1intersects <- function(c1,c2) return(all(c1 == c2))
res <- vector(mode="numeric",length=(nrow(design)-1))

for(i in 2:nrow(design)){
res[i-1] <- sum(apply(array,1,sum1intersects,design[i,]))
}
return(res)
}