plot_eval: Plot time evolution of errors

In robertsy/assimilr: Local Ensemble Kalman Particle Filters

Description

Plot a time series of bias-rmse-crps for a given field

Usage

plot_eval(ens_eval, field_ind = 2)

Arguments

ens_eval	is a data frame produced by sweq_eval
field_ind	choose which field (typically 1:fluid height, 2: rain, 3:wind)

Examples

set.seed(1)
sweq_run <- sweq_simulate(5*60, 60, 300)
k <- 20
ens0 <- sweq_ens0(k, sweq_run, klag=1000)
l <- 5
block_run <- da_cycle(ens0, sweq_run, block_LEnKPF, l=l, block_size=l/2, 
                      taper=sweq_GC(sweq_run$ndim, l, 0.9), ndim=sweq_run$ndim,
                      get_partition = sweq_partition,
                      manifold_proj = sweq_proj)
block_df <- sweq_eval(sweq_run, block_run)
block_df %>% mutate(method='block-LEnKPF') %>% plot_eval(2)

robertsy/assimilr documentation built on May 27, 2019, 10:33 a.m.