matcontrib <- matcontrib0[is_normal, ] # wipe out aberrant data
matcontrib <- t(apply(matcontrib, 1, function(x) x / sum(x)))
matcontrib <- matcontrib0[is_normal, ] # wipe out aberrant data
matcontrib <- t(apply(matcontrib, 1, function(x) x / sum(x)))