= Solution
Model A has fifteen separate curve parameters plus one error-scale parameter. With such broad priors, one might have expected more of those sixteen parameters to be effectively estimated, rather than $p_D=12$, and might also have expected its greater unpooled flexibility to give better fit than the displayed $\overline D=248.1$. Model B instead has a smaller mean deviance despite sharing information between trees.
This is a diagnostic concern, not a mathematical impossibility. <Effective parameter count in DIC> depends on posterior curvature and parameterization, while mean posterior deviance is not minimized deviance. A weakly identified nonlinear asymptote–slope trade-off, skewness, multiple posterior regions or slow <Markov chain Monte Carlo> mixing could explain the unexpected comparison. For A, the implied deviance at the posterior mean is $\overline D-p_D=236.1$; for B it is $230.3$. Inspect traces, effective sample sizes, independent chains and fitted curves, and extend sampling where needed. Fifty thousand iterations alone do not establish convergence.
Back to article page