Calculates factor-stratified sums for each column
colSumByFac(sY, rowSel)
sparse matrix (dgCmatrix)
integer factor. Note that the 0-th column will return sums for any NA values; 0 or negative values will be omitted
Matrix