plot_sensitivity_analysis_results {StrathE2E2}R Documentation

Create a plot of results from the Morris method sensitivity analysis

Description

Reads results generated by the function Sensitivity_analysis_StrathE2E() and plots the Elementary Effect mean (x-axis; magnitude of sensitivity) against Elementary Effect standard deviation (y-axis; strength of interactions between parameters).

Usage

plot_sensitivity_analysis_results(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

Results from the Sensitivity_analysis_StrathE2E() function are stored in the CSV file /results/Modelname/Variantname/sorted_parameter_elementary_effects-*.csv where * represents the model run identifier (model.ident) text embedded in the R-list object created by the read_model() function.

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

Each symbol in the plot represents a single parameter in the model. The parameters are colour-coded to indicate 6 different types - fitted and fixed parameters of the ecology model, fishing fleet model parameters, fishery harvest ratios, environmental drivers, and physical configuration parameters.

The plot also shows a wedge formed by the two dashed lines. These correspond to +/-2 standard errors of the mean, so for points falling outside of the wedge there is a significant expectation that the distribution of elementary effects is non-zero. For points falling within the wedge the distribution of elementary effects is not significantly different from zero

For details of how the Elementary Effect values are derived for each parameter see ?Sensitivity_analysis_StrathE2E

Value

Graphical display in a new graphics window

References

Morris, M.D. (1991). Factorial sampling plans for preliminary computational experiments. Technometrics, 33, 161-174.

See Also

read_model, Sensitivity_analysis_StrathE2E, concatenate_raw_sensitivity_analysis_results , process_sensitivity_analysis_results_offline

Examples

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

# Plot the example data for this model variant provided with the package
plot_sensitivity_analysis_results(model, use.example=TRUE)

# To direct the graph output to a file rather than the screen, wrap the plot_sensitivity_analysis_results() function call in a graphical device call:
pdf("C:/Users/username/Documents/Foldername/plot.pdf",width=6,height=6)      # or jpeg("plot.jpg"), png("plot.png")
	plot_sensitivity_analysis_results(model, use.example=TRUE)
dev.off()


[Package StrathE2E2 version 2.0.0 Index]