ClimProjDiags (version 0.1.3)

CombineIndices: Combine weighted indices of n-dimensional arrays

Description

Function to combine climate indices for multiple models through addition, subtraction, division or averaging, optionally applying weights to each index.

Usage

CombineIndices(indices, weights = NULL, operation = "mean")

Value

An array of the same dimensions as one of the elements in the parameter indices.

Arguments

indices

List of n-dimensional arrays with equal dimensions to be combined.

weights

Vector of weights for the indices, whose length is the same as the list of parameter indices. If not provided, a weight of 1 is assigned to each index. If operation = 'mean' the weights are normalized to sum 1 all together.

operation

The operation for combining the indices, either "mean" (default), "add", "subtract" or "divide".

Examples

a <- matrix(rnorm(6), 2, 3)
b <- matrix(rnorm(6), 2, 3)

comb_ind <- CombineIndices(indices = list(a, b), weights = c(2, 1), operation = "add")
print(comb_ind)

a <- rnorm(24)
dim(a) <- c(lon = 2, lat = 3, mod = 4)
b <- rnorm(24)
dim(b) <- c(lon = 2, lat = 3, mod = 4)
comb_ind <- CombineIndices(indices = list(a, b), weights = c(2, 1), operation = "add")
print(comb_ind)