Description Usage Arguments Details Value Examples
View source: R/GenerateToyData.R
Generate artificial data for ensemble verification using a signal-plus-noise model
1 2 3 4 5 6 7 8 9 10 11 12 13 14 |
N |
number of forecasts and observations |
mu.y |
expectation value of the observations |
s.s |
standard deviation of the predictable signal |
s.eps |
standard deviation of the unpredictable noise |
mu.x |
expectation value of the ensemble |
beta |
weighting parameter of the signal in the ensemble forecasting system |
s.eta |
average spread of the ensemble |
K |
number of members of the ensemble |
mu.x.ref |
expectation value of the reference ensemble |
beta.ref |
weighting parameter of the signal in the reference ensemble forecasting system |
s.eta.ref |
average spread of the reference ensemble |
K.ref |
number of members of the reference ensemble |
The function simulates data from the latent variable model:
y_t = mu_y + s_t + eps_t
x_t,r = mu_x + beta * s_t + eta_t,r
where y_t is the observation at time t, and x_t,r is the r-th ensemble member at time t. The latent variable s_t is to be understood as the "predictable signal" that generates correlation between observations and ensemble members. If all arguments that end in ".ref" are specified, a reference ensemble is returned to also test comparative verification.
A list with elements:
N-vector of observations
N*K matrix of ensemble members
N*K.ref matrix of reference ensemble members
1 2 | l <- GenerateToyData()
with(l, EnsCrps(ens, obs))
|
[1] 3.076019 6.148100 4.931804 4.928815 2.855213 4.676788 2.256443
[8] 3.037386 1.937655 10.794279 4.636184 10.585075 7.844967 5.713086
[15] 4.917702 6.593254 4.731927 2.397302 11.501353 1.933082
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.