METADATA RECORD
Authors: S. Volzke & R. Dunlop, The University of Queensland, Brisbane, QLD, Australia 
Contact: r.dunlop@uq.edu.au
In collaboration with C. Booth (SMRU Consulting), 
With funding provided by the National Environmental Science Program (NESP) Marine and Coastal Hub 2024-2025
(NESP MaC Hub Project 4.9)
https://www.nespmarinecoastal.edu.au/project/4-9/


Files: 
(1) summary_statistics_SV2.R provides code to reproduce graphs and visualise results, supported by 
(2) SUMMARY_FUNCTIONS_SV.R which includes the functions used in the above file.

The above access the following output data files:
(3) iPCoD_EE_S1_BW_Output.rdata
(4) iPCoD_EE_S1_SRW_Output.rdata
(5) iPCoD_EE_S2_BW_Output.rdata
(6) iPCoD_EE_S2_SRW_Output.rdata

The .rdata files contain the outputs from iPCoD in a file called "dat.out".
dat.out is a matrix called where the results of each of 1,000 paired population trajectories are stored. 
The matrix dat.out contains 5 columns and 26 rows for each of the 1,000 simulations.
The columns are:
• dat.out [,1,] "Impacted" is the total size of the impacted population (numeric)
• dat.out [,2,] "Unimpacted" is the total size of the baseline population (numeric)
• dat.out [,3,] "Unimpacted-Impacted" is the difference in total population size between the two paired populations (baseline minus impacted) (numeric)
• dat.out [,4,] "average pa% change Impacted" is the average % annual change in the size of the impacted population from time.point 1. (decimal between 0-1)
• dat.out [,5,] "average pa% change Unimpacted" is the average % annual change in the size of the baseline population from time.point 1. (decimal between 0-1)
The rows represent the population information for each year in the simulation. When running this code it is important to understand that for a 25 year simulation, the dat.out file will have 26 time points for each of the 1,000 paired simulations:
time.point=1 the start of year 1 (ie: before any time has passed or any impact has occurred)
time.point=2 is the start of year 2 (ie: after 1 year of impact)
time.point=3 is the start of year 3 (ie: after 2 years of impact has occurred)
… etc …
time.point=26 is the start of year 26 (ie: after the full 25 simulated years).