box_and_whisker_monthly_plots_with_credible_intervals {StrathE2E2}R Documentation

Box-and-whisker plots of observed monthly data on the state of the ecosystem with equivalent properties derived from a Monte Carlo StrathE2E run

Description

Creates a multi-panel plot comparing annual cycles of monthly averaged observational data on nutrinet and plankton concentrations together with the distribution of credible values of equivalent data derived from the final year of a model run generated by the Monte_Carlo_StrathE2E() function

Usage

box_and_whisker_monthly_plots_with_credible_intervals(model,
  use.example = FALSE)

Arguments

model

R-list object defining the model configuration compiled by the read_model() function

use.example

(TRUE or FALSE) Option to use pre-computed example data from the internal North Sea model rather than user-generated data (default=FALSE)

Details

For details of how the distribution of credible output values from StrathE2E are calculated see ?Monte_Carlo_StrathE2E.

The function plots a multi-panel page of box-and-whisker plots showing the medians and variability ranges (quartiles as box-and-whisker) of observational data on nitrate, ammonia, chlorophyll, omnivorous and carnivorous zooplankton, and benthos larvae concentrations (shown in black), alongside comparable box-and-whisker plots (shown in red) of equivalent measures derived from the final years of an ensemble of model runs generted by a Monte Carlo methodology (Monte_Carlo_StrathE2E() function)

The observational data to be plotted are loaded from the folder Modelname/Variantname/Target_data/monthly_target_data-*.csv as part of a read_model() function call and are built into the R-list object generated by read_model().

Optionally the function can read an example data set for one of the two North Sea model variants supplied with the package.

The corresponding measures derived from the final year of a model run generated by the Monte_Carlo_StrathE2E() function call are located in /results/Modelname/Variantname/CredInt_processed_monthly_mass-*.csv, and in the R-list object generated by Monte_Carlo_StrathE2E()

Value

Graphical display in a new graphics window

See Also

read_model, Monte_Carlo_StrathE2E, box_and_whisker_monthly_plots, box_and_whisker_annual_plots_with_credible_intervals

Examples

# Load the 1970-1999 version of the North Sea model supplied with the package:
model <- read_model("North_Sea", "1970-1999")

# Either run the Monte_Carlo_StrathE2E() to generate some results, or use the internal example data

# Plot the final year comparison of observed and modelled ecosystem state on screen:
box_and_whisker_monthly_plots_with_credible_intervals(model, use.example=TRUE)

# Direct the graphics output to a file (Windows OS)... 
pdf("C:/Users/username/Documents/Foldername/plot.pdf",width=8,height=6)      # or jpeg("plot.jpg"), png("plot.png")
	box_and_whisker_monthly_plots_with_credible_intervals(model, use.example=TRUE)
dev.off()

[Package StrathE2E2 version 2.0.0 Index]