vectorized function for getting outliers in rows of a matrix

vectorize.r

data <- cbind(rnorm(100),rnorm(100),rnorm(100))

outlierMat <- function(mat) { 
m <- rowMeans(mat) 
devs <- abs(mat - m) 
val <- apply(mat, 1, max) 
pos <- which(mat == val, arr.ind = TRUE) 
out <- cbind(pos,val) 
return(out) 
}

out <- outlierMat(data)
out