1

JAGS allows reporting of DIC and also PED (penalized expected deviance). What are differences between these two statistics? Which is better for what purpose? DIC was the first one, PED was evidently added later, so what was the reason for introducing PED? In what regard and in which situations/type of models is it more appropriate than DIC?

I did a little research and found an answer from Martyn Plummer (JAGS author) stating that the difference is in the penalty for complexity of the model. However, when I looked at my models, I noticed that even if I compare models with the same structure, just different covariate (but always the same dependent variable), I notice that although DIC and PED are correlated many times, there is still some variability and sometimes they both select a different models. And sometimes they are not even correlated:

enter image description here

Corresponding DIC and PED are always from one model run, so MCMC stochasticity cannot explain the variability. Model complexity is the same in all cases, as well as the response variable, the models only differ in the covariate used. Thus, Martin Plummers' answer doesn't explain it for me.

Subquestion: In AIC, if the difference of AIC between two models was greater than 2, we said that the models "differ enough". What would be the equivalent "rule of thumb" for DIC and PED? Is it still 2, or is it 1?

Tomas
  • 5,735
  • 11
  • 52
  • 93
  • 1
    Perhaps this paper (also from Martyn Plummer) will clear up the confusion? http://biostatistics.oxfordjournals.org/content/9/3/523.full – Jacob Socolar Dec 09 '16 at 21:42
  • 1
    A brief summary: Plummer claims (and I assume he's right, but I've not read in detail) that the DIC can be placed on firmer theoretical ground by interpreting it as an approximation to the PED. When this approximation fails, PED should presumably be used as it remains well-justified. – Jacob Socolar Dec 09 '16 at 21:48
  • @user43849, thanks, sounds like being close to an answer, right? You can still win the bounty if you compose an answer! :-) Just few things: 1) I am lost in the terminology a bit. "Penalized loss functions" mentioned in paper - is it what the PED is? 2) does this mean that it is always better and more correct to use PED than DIC? 3) does this mean that PED also works well (unlike DIC) in more complex (hierarchical) models? Thanks a lot! – Tomas Dec 15 '16 at 10:23

0 Answers0