box_and_whisker_monthly_plots_with_credible_intervals {StrathE2E2} | R Documentation |
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
box_and_whisker_monthly_plots_with_credible_intervals(model, use.example = FALSE)
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) |
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()
Graphical display in a new graphics window
read_model
, Monte_Carlo_StrathE2E
, box_and_whisker_monthly_plots
, box_and_whisker_annual_plots_with_credible_intervals
# 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()