powered by
matrix_variability(mat_lst, mean_mat = apply(abind::abind(mat_lst, along = 3), c(1, 2), mean))