spict v1.3.8 (2023-12-08)

New features

• Added function calc.process.resid() that calculates the process residuals for
  the biomass and fishing mortality processes. The residuals are added to the
  fitted spict object in the data frame called rep$process.resid. The
  residuals can be summarised with the function sumspict.diagnostics() and
  visualised with the function plotspict.diagnostic.process(). Note, that the
  process residuals have mainly been tested for scenarios without seasonal fishing patterns.

• Added function plotspict.compare that shows estimates from multiple
  fits, of B, F, catch, F/Fmsy, B/Bmsy and the production curve on the same plot. The
  fits can be given separately, like plotspict.compare(fit1, fit2) or in a list, like
  plotspict.compare(list(fit1, fit2)).

Bug fixes

• make.man.inp now checks all fractiles to be between 0 and 0.5 (closes #167)
• add.manlines correctly reports shorter management periods than 1 year (closes #168)
• add.manlines checks for errflag for multiple management scenarios
• annual function was rewritten to accommodate cases where intime does not include full years. This happens when the first data point is not on a full year (closes #170, thanks to Paz Sampedro)
• man.select checks for NAs for multiple management scenarios
• summary now respects the CI argument. Added CI argument for sumspict.* functions

Minor changes

• Retro runs have names: Baseline, -1, -2, -3, ... (closes #172)