col_stats function

Apply a statistic function to all column vectors