Мы также можем использовать pmin/pmax
из base R
cbind(pmax(mat[,1], mat[,2]), pmin(mat[,1], mat[,2]))
# [,1] [,2]
#[1,] 25 15
#[2,] 55 35
#[3,] 65 45
#[4,] 50 40
данные
mat <- structure(c(25L, 55L, 65L, 50L, 15L, 35L, 45L, 40L), .Dim = c(4L,
2L), .Dimnames = list(NULL, c("co1", "co2")))