Last chance! 50% off unlimited learning
Sale ends in
Applies the wavelet-based weather simulation algorithm to multiple sites (single site possible as well)
prsim.weather(data_p, data_t, station_id_p="Precip",
station_id_t="Temp", number_sim=1, win_h_length=15,
n_wave=100,verbose=TRUE,t_margin='sep',p_margin='egpd',...)A list with elements temperature and precipitation of
data frames with time information, observations, and
number_sim columns containing the simulated data.
list of precipitation data frames. One list entry, i.e. data frame, corresponds to one station/grid cell. Each data frame contains the time indications and precipitation of one station. See ‘Details’.
list of temperature data frames. One list entry, i.e. data frame, corresponds to one station/grid cell. Each data frame contains the time indications and temperature of one station. See ‘Details’.
identifies the precipitation variable name in case several time series are present in data_p. See ‘Details’.
identifies the temperature variable name in case several time series are present in data_t. See ‘Details’.
number of simulations to be carried out.
(half-)length of moving window size.
marginal distribution to be used for the backtransformation of temperature. Can be either "sep" or any type of CDF (see ‘Details’). "sep" uses the four-parameter skewed-exponential power for backtransformation. CDF allows for specifying any distribution ‘Examples’.
marginal distribution to be used for the backtransformation of precipitation. Can be either "egpd" or any type of CDF (see ‘Details’). "egpd" uses the extended GPD for backtransformation. CDF allows for specifying any distribution ‘Examples’.
logical. Should progress be reported?
number of scales to be considered in the continuous wavelet transform.
any other argument passed to the sub-function specifying the cdf for fitting. See ‘Details’ and ‘Examples’.
Manuela Brunner
Time can be given with three columns named "YYYY", "MM", "DD", or as in POSIXct format YYYY-MM-DD.
All leap days (Feb 29th) will be omitted from the analysis, but no missing observations are allowed.
Stations are identified by list index.
Alternative distributions can be specified by providing three functions: (1) a function fitting the parameters of a distributions and providing a vector of these parameters as output (CDF_fit), (2) a function simulating random numbers from this distribution (rCDF), and (3) a function specifying the distribution (pCDF). See ‘Examples’ for the generalized beta for the second kind and for the Generalized Extreme Values (GEV) distribution.
Brunner, M. I., and E. Gilleland (2021). Spatial compound hot-dry events in the United States: assessment using a multi-site multi-variable weather generator, in preparation.