Improve readability #33

# old
cols <- 3:5
x <- colSums(areas * as.matrix(output[, cols])) / sum(areas)
stats::setNames(as.list(x), names(output)[cols])

# new
values <- as.matrix(output[, 3:5])
means <- colSums(areas * values) / sum(areas)
stats::setNames(as.list(means), colnames(values))