add DIC

Author: itsdfish

ParetoSmootherImportanceSampling/Project.toml

@@ -0,0 +1,2 @@
+[deps]
+Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"

examples/arsenic/compare.png

@@ -32,17 +32,17 @@ arXiv preprint arXiv:1507.02646.
 
 module ParetoSmoothedImportanceSampling
 
-using StatsFuns
+using StatsFuns, Statistics
 
 psis_path = @__DIR__
 
 include("psisloo.jl")
 include("psislw.jl")
 include("gpdfitnew.jl")
 include("gpinv.jl")
-#include("logsumexp.jl")
 include("waic.jl")
 include("pk_utilities.jl")
+include("dic.jl")
 
 export
     psis_path

src/ParetoSmoothedImportanceSampling.jl

@@ -0,0 +1,18 @@
+
+"""
+    dic(loglike::AbstractVector{<:Real})
+
+Computes Deviance Information Criterion (DIC).
+
+# Arguments
+* `log_lik::AbstractArray`: Array of size n x m containing n posterior samples of the log likelihood terms p(y_i|\theta^s).
+
+# Returns
+* `dic::Real`: DIC value
+"""
+function dic(loglike::AbstractVector{<:Real})
+    D = deviance.(loglike)
+    return mean(D) + 0.5 * var(D)
+end
+
+deviance(loglikelihood::Real) = -2 * loglikelihood