mvgauss_to_Xarray: Computes an ensemble of realisations from a multivariate...

View source: R/tools.R

mvgauss_to_XarrayR Documentation

Computes an ensemble of realisations from a multivariate gaussian distribution parameters

Description

For multiple types of external forcings (natural, all in the current version), mvgauss_to_Xarray generates a sample of time series corresponding to the response to each forcing from the gaussian parameters.

Usage

mvgauss_to_Xarray(mu, Sigma, Nres)

Arguments

mu

a vector of time series corresponding to the mean of the multivariate gaussian distribution. The names of each element must correspond to the folowing pattern : year_forcing, eg 1850_nat for the mean natural response in 1850. The current allowed forcings are nat and all. mu is tipically a vector of concatenated time series relative to the response to the natural and all external forcings for a given period.

Sigma

a matrix of size [length(mu),length(mu)] corresponding to the variance-covariance matrix of the multivariate gaussian distribution

Nres

the size of the wished sample

Value

returns a 3-D array of dimension [Nres,length(year),length(forcing)] where year is a time series of years and forcing is the type of forcing


saidqasmi/KCC documentation built on July 8, 2022, 6:02 a.m.