a <- t(matrix(c(4,3,4,1,2,2,3,3,1), 3, 3))
a
b <- t(matrix(c(2,3,5,1,5,9,6,7,2), 3, 3))
b
idx <- do.call(rbind, lapply(1:nrow(a), function(n.row) {
return(which(a[n.row,]==max(a[n.row,])))
}))
idx
out <- do.call(rbind, lapply(1:nrow(idx), function(n.row) {
return(sum(b[n.row...