ClimProjDiags (version 0.1.3)

Subset: Subset a Data Array

Description

This function allows to subset (i.e. slice, take a chunk of) an array, in a similar way as done in the function take() in the package plyr. There are two main inprovements:

The input array can have dimension names, either in names(dim(x)) or in the attribute 'dimensions', and the dimensions to subset along can be specified via the parameter along either with integer indices or either by their name.

There are additional ways to adjust which dimensions are dropped in the resulting array: either to drop all, to drop none, to drop only the ones that have been sliced or to drop only the ones that have not been sliced.

If an array is provided without dimension names, dimension names taken from the parameter dim_names will be added to the array. This function computes the threshold based on a quantile value for each day of the year of the daily data input.

Usage

Subset(x, along, indices, drop = FALSE)

Value

An array with similar dimensions as the x input, but with trimmed or dropped dimensions.

Arguments

x

A multidimensional array to be sliced. It can have dimension names either in names(dim(x)) or either in the attribute 'dimensions'.

along

Vector with references to the dimensions to take the subset from: either integers or dimension names.

indices

List of indices to take from each dimension specified in 'along'. If a single dimension is specified in 'along' the indices can be directly provided as a single integer or as a vector.

drop

Whether to drop all the dimensions of length 1 in the resulting array, none, only those that are specified in 'along', or only those that are not specified in 'along'. The possible values are, respectively: 'all' or TRUE, 'none' or FALSE, 'selected', and 'non-selected'.

Examples

##Example synthetic data:
data <- 1:(2 * 3 * 372 * 1)
dim(data) <- c(time = 372, lon = 2, lat = 3, model = 1)
data_subset <- Subset(data, c('time', 'model'), 
                     list(1:10, TRUE), drop = 'selected')
dim(data_subset)