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
.
.
.
.
.

f <- function(x) {
  c(Mean = mean(x),
    Min = min(x),
    Max = max(x))
}

out <- t(sapply(dat, f))

out
#    Mean Min Max
#X1  47.0  21  91
#X2  64.4   7  95
#X3  44.8  18  77
#X4  67.8  39 100
#X5  44.4  13  94
#X6  40.4   2  87
#X7  52.2  19  83
#X8  54.6  11  80
#X9  67.2  53  83
#X10 54.8   3  79
#X11 42.4   8  87
#X12 40.4  10  67
#X13 53.6  30  92

set.seed(1)
dat <- data.frame(matrix(sample(1:100, 13*5, TRUE), ncol = 13))