nrow <- length(colnames(df13))
ncol = 4
output <- matrix(1,mrows, mcol)

for (i in colnames(df13)){
  min <- min$df13
  max <- max$df13
  mean <- mean$df13
}

Feature name | Mean | Min | Max
Var 1        | 2    |  4  |  5
Var 2        | 6    |  5  |  9
.
.
.
.
.