plot_for_each_guild_catch_by_gear {StrathE2E2} | R Documentation |
Multi-panel bar-plots of the distributions of annual integrated catches of each guild across gears for the final year of a model run generated by the StrathE2E() function
plot_for_each_guild_catch_by_gear(model, results)
model |
R-list object defining the model configuration compiled by the read_model() function |
results |
R-list object containing model results generated by the StrathE2E() function |
The function plots a multi-panel page of bar-plots of annual integrated catches of the guilds in the ecology model which are potentially catchable by the gears in the fishing fleet model: macrophytes, carnivorous zooplankton, suspension/deposit and carnivore/scavenge feeding benthos, planktivorous, demersal and migratory fish, birds, pinnipeds and cetaceans.
Each panel represesents a single guild. In each panel, the x-axis represents the (up to) 12 gear types parameterised in the fishing fleet model, and the y-axis is the annual integrated catch. For each gear there are two bars corresponding to the inshore and offshore zones. Each bar is subdivided into landed and discarded components of the catch. Green and blue sections are offshore and inshore landings respectively; black and grey sections are offshore and inshore discards respectively.
Units for the catch are mMN/m2/year
Graphical display in a new graphics window
read_model
, StrathE2E
, plot_time_series_annual_land_disc
, plot_for_each_gear_catch_by_guild
# Load the 2003-2013 version of the North Sea model supplied with the package: model <- read_model("North_Sea", "2003-2013") #Run the model and generate the results object results <- StrathE2E(model,nyears=10) # Plot the time series of annual landings and discards plot_for_each_guild_catch_by_gear(model, results) # Time series plot of state varaiables in a new window leaving any existing windows open dev.new() plot_for_each_guild_catch_by_gear(model, results) #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") plot_for_each_guild_catch_by_gear(model, results) dev.off() # Read the North Sea/1970-1999 model and set the identifier for output files to "Baseline", run the model and plot the distribution of catches across gears model <- read_model("North_Sea", "1970-1999",model.ident="Baseline") results <- StrathE2E(model,nyears=10) plot_for_each_guild_catch_by_gear(model, results) # Create a new scenario version of the North Sea/1970-1999 model by increasing the activity rate of gear 1 (pelagic trawls and seines) by a factor of 2 model$data$fleet.model$gear_mult[1] <- 2 model$setup$model.ident <- "gear1x2" # Set a new identifier for the outputs results <- StrathE2E(model,nyears=10) dev.new() # Open a new graphics window plot_for_each_guild_catch_by_gear(model, results)