computing corrected DIC

[dokato]

Is there a way to access the contribution to the effective number of parameters from observation p_Di:

as suggested in corrected version of DIC from (Plummer, 2008).

I wanted to add that as extra model comparison criterion, but noticed that p_D so far is computed indirectly:

https://github.com/hddm-devs/kabuki/blob/master/kabuki/hierarchical.py#L672

and not broken down by i.