# Example 1: Model predictions only ---------------------------------------
a_model <- dmc_dm()
deltas <- calc_stats(
a_model,
type = "delta_funs",
minuends = "incomp",
subtrahends = "comp"
)
plot(deltas)
plot(deltas, col = "black", lty = 2, xlim = c(0.2, 0.65))
# Example 2: Observed and predicted data ----------------------------------
obs_data(a_model) <- dmc_synth_data
deltas <- calc_stats(
a_model,
type = "delta_funs",
minuends = "incomp",
subtrahends = "comp"
)
plot(deltas)
# Example 3: Observed data only -------------------------------------------
deltas <- calc_stats(
dmc_synth_data,
type = "delta_funs",
minuends = "incomp",
subtrahends = "comp"
)
plot(deltas)
# Example 4: Observed data with intervals ---------------------------------
deltas <- calc_stats(
dmc_synth_data,
type = "delta_funs",
minuends = "incomp",
subtrahends = "comp",
resample = TRUE
)
plot(deltas)
Run the code above in your browser using DataLab