plot-mmcfits.Rd
Method to plot fitting statistics of multiple fits, useful to compare fits.
# S4 method for class 'a4aFitMCMCs,missing'
plot(x, y = missing, ...)
a plot
with fitting statistics
data(ple4)
data(ple4.index)
qmods <- list(list(~s(age, k=6)))
fmods = list()
for(i in 1:6) {
fmods[[paste0(i)]] <-
as.formula(
paste0("~te(age, year, k = c(6,", i+14,"), bs = 'tp') + s(age, k = 6)")
)
}
myFits <-
scas(
FLStocks(ple4), list(FLIndices(ple4.index)),
fmodel = fmods, qmodel=qmods, fit="MCMC", mcmc=SCAMCMC()
)
plot(myFits)