search

samueldnj / hierProdAnalysis

A simulation-estimation model for that compares the performance of hierarchical multispecies and single species surplus production stock assessment models.
0 stars 0 forks source link

4d arrays to hold MCMC output #11

Closed samueldnj closed 7 years ago

samueldnj commented 8 years ago

Working with MCMC output would be simplified by using arrays to hold MCMC values. This would allow apply functions to be used to compute quantiles.

dim = c(nReps, nS, mcTrials, nPars)

Because of unequal numbers of MC trials that have pospen=0, I'll focus on comparing MLEs for now

samueldnj commented 7 years ago

Updated the blob so MCMC output is held in 4d arrays, with dimensions (nReps,nS,nMCMC,nMCpars). For mcPar, nMCpars is controlled by the control file, while for mcBio, nMCpars is just the number of time steps (mcBio is the biomass posterior).