doctest failure in ~/.julia/packages/PSIS/fU76x/src/core.jl:212-233 ```jldoctest psis; setup = :(using Random; Random.seed!(42)) julia> using Distributions julia> proposal, target = Normal(), TDist(7); julia> x = rand(proposal, 1_000, 1, 30); # (ndraws, nchains, nparams) julia> log_ratios = @. logpdf(target, x) - logpdf(proposal, x); julia> result = psis(log_ratios) ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— ``` Subexpression: result = psis(log_ratios) Evaluated output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— Expected output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— diff = Warning: Diff output requires color. ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) ——

doctest failure in ~/.julia/packages/PSIS/fU76x/src/core.jl:238-255 ```jldoctest psis julia> using MCMCDiagnosticTools julia> reff = ess(log_ratios; kind=:basic, split_chains=1, relative=true); julia> result = psis(log_ratios, reff) ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— ``` Subexpression: result = psis(log_ratios, reff) Evaluated output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— Expected output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— diff = Warning: Diff output requires color. ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) ——

doctest failure in ~/.julia/packages/PSIS/fU76x/src/core.jl:212-233 ```jldoctest psis; setup = :(using Random; Random.seed!(42)) julia> using Distributions julia> proposal, target = Normal(), TDist(7); julia> x = rand(proposal, 1_000, 1, 30); # (ndraws, nchains, nparams) julia> log_ratios = @. logpdf(target, x) - logpdf(proposal, x); julia> result = psis(log_ratios) ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— ``` Subexpression: result = psis(log_ratios) Evaluated output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— Expected output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— diff = Warning: Diff output requires color. ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 7 (23.3%) 959 (0.5, 0.7] okay 13 (43.3%) 938 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) ——

doctest failure in ~/.julia/packages/PSIS/fU76x/src/core.jl:238-255 ```jldoctest psis julia> using MCMCDiagnosticTools julia> reff = ess(log_ratios; kind=:basic, split_chains=1, relative=true); julia> result = psis(log_ratios, reff) ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— ``` Subexpression: result = psis(log_ratios, reff) Evaluated output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— Expected output: ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) —— diff = Warning: Diff output requires color. ┌ Warning: 9 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:373 ┌ Warning: 1 parameters had Pareto shape values k > 1. Corresponding importance sampling estimates are likely to be unstable and are unlikely to converge with additional samples. └ @ PSIS ~/.julia/packages/PSIS/... ~/.julia/packages/PSIS/fU76x/src/core.jl:376 PSISResult with 1000 draws, 1 chains, and 30 parameters Pareto shape (k) diagnostic values: Count Min. ESS (-Inf, 0.5] good 9 (30.0%) 806 (0.5, 0.7] okay 11 (36.7%) 842 (0.7, 1] bad 9 (30.0%) —— (1, Inf) very bad 1 (3.3%) ——

doctest failure in ~/.julia/packages/InferenceObjects/ATAPR/src/inference_data.jl:289-316 ```jldoctest cat julia> coords = (; a_dim=["x", "y", "z"]); julia> dims = dims=(; a=[:a_dim]); julia> data = Dict(:a => randn(100, 4, 3), :b => randn(100, 4)); julia> idata = from_dict(data; coords=coords, dims=dims) InferenceData with groups: > posterior julia> idata_cat1 = cat(idata, idata; dims=:chain) InferenceData with groups: > posterior julia> idata_cat1.posterior ┌ 100×8×3 Dataset ┐ ├─────────────────┴──────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered ├──────────────────────────────────────────────────── layers ┤ :a eltype: Float64 dims: draw, chain, a_dim size: 100×8×3 :b eltype: Float64 dims: draw, chain size: 100×8 ├────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704" ``` Subexpression: idata_cat1.posterior Evaluated output: ╭─────────────────╮ │ 100×8×3 Dataset │ ├─────────────────┴──────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered ├──────────────────────────────────────────────────── layers ┤ :a eltype: Float64 dims: draw, chain, a_dim size: 100×8×3 :b eltype: Float64 dims: draw, chain size: 100×8 ├────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-12T01:28:14.766" Expected output: ┌ 100×8×3 Dataset ┐ ├─────────────────┴──────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered ├──────────────────────────────────────────────────── layers ┤ :a eltype: Float64 dims: draw, chain, a_dim size: 100×8×3 :b eltype: Float64 dims: draw, chain size: 100×8 ├────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704" diff = Warning: Diff output requires color. ┌ ╭─────────────────╮ │ 100×8×3 Dataset ┐ │ ├─────────────────┴──────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered ├──────────────────────────────────────────────────── layers ┤ :a eltype: Float64 dims: draw, chain, a_dim size: 100×8×3 :b eltype: Float64 dims: draw, chain size: 100×8 ├────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704""2025-02-12T01:28:14.766"

doctest failure in ~/.julia/packages/InferenceObjects/ATAPR/src/inference_data.jl:320-338 ```jldoctest cat julia> idata_cat2 = cat(idata, idata; dims=:run) InferenceData with groups: > posterior julia> idata_cat2.posterior ┌ 100×4×3×2 Dataset ┐ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704" ``` Subexpression: idata_cat2.posterior Evaluated output: ╭───────────────────╮ │ 100×4×3×2 Dataset │ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-12T01:28:14.766" Expected output: ┌ 100×4×3×2 Dataset ┐ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704" diff = Warning: Diff output requires color. ┌ ╭───────────────────╮ │ 100×4×3×2 Dataset ┐ │ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704""2025-02-12T01:28:14.766"

doctest failure in ~/.julia/packages/InferenceObjects/ATAPR/src/inference_data.jl:344-380 ```jldoctest cat julia> observed_data = Dict(:y => randn(10)); julia> idata2 = from_dict(data; observed_data=observed_data, coords=coords, dims=dims) InferenceData with groups: > posterior > observed_data julia> idata_cat3 = cat(idata, idata2; groups=(:posterior,), dims=:run) InferenceData with groups: > posterior > observed_data julia> idata_cat3.posterior ┌ 100×4×3×2 Dataset ┐ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704" julia> idata_cat3.observed_data ┌ 10-element Dataset ┐ ├────────────── dims ┤ ↓ y_dim_1 ├────────────────────┴─────────────── layers ┐ :y eltype: Float64 dims: y_dim_1 size: 10 ├────────────────────────────────────────────┴ metadata ┐ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.71" ``` Subexpression: idata_cat3.posterior Evaluated output: ╭───────────────────╮ │ 100×4×3×2 Dataset │ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-12T01:28:14.766" Expected output: ┌ 100×4×3×2 Dataset ┐ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704" diff = Warning: Diff output requires color. ┌ ╭───────────────────╮ │ 100×4×3×2 Dataset ┐ │ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704""2025-02-12T01:28:14.766"

doctest failure in ~/.julia/packages/InferenceObjects/ATAPR/src/inference_data.jl:344-380 ```jldoctest cat julia> observed_data = Dict(:y => randn(10)); julia> idata2 = from_dict(data; observed_data=observed_data, coords=coords, dims=dims) InferenceData with groups: > posterior > observed_data julia> idata_cat3 = cat(idata, idata2; groups=(:posterior,), dims=:run) InferenceData with groups: > posterior > observed_data julia> idata_cat3.posterior ┌ 100×4×3×2 Dataset ┐ ├───────────────────┴─────────────────────────────────── dims ┐ ↓ draw , → chain, ↗ a_dim Categorical{String} ["x", "y", "z"] ForwardOrdered, ⬔ run ├─────────────────────────────────────────────────────────────┴ layers ┐ :a eltype: Float64 dims: draw, chain, a_dim, run size: 100×4×3×2 :b eltype: Float64 dims: draw, chain, run size: 100×4×2 ├──────────────────────────────────────────────────────────── metadata ┤ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.704" julia> idata_cat3.observed_data ┌ 10-element Dataset ┐ ├────────────── dims ┤ ↓ y_dim_1 ├────────────────────┴─────────────── layers ┐ :y eltype: Float64 dims: y_dim_1 size: 10 ├────────────────────────────────────────────┴ metadata ┐ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.71" ``` Subexpression: idata_cat3.observed_data Evaluated output: ╭────────────────────╮ │ 10-element Dataset │ ├────────────── dims ┤ ↓ y_dim_1 ├────────────────────┴─────────────── layers ┐ :y eltype: Float64 dims: y_dim_1 size: 10 ├────────────────────────────────────────────┴ metadata ┐ Dict{String, Any} with 1 entry: "created_at" => "2025-02-12T01:28:19.299" Expected output: ┌ 10-element Dataset ┐ ├────────────── dims ┤ ↓ y_dim_1 ├────────────────────┴─────────────── layers ┐ :y eltype: Float64 dims: y_dim_1 size: 10 ├────────────────────────────────────────────┴ metadata ┐ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.71" diff = Warning: Diff output requires color. ┌ ╭────────────────────╮ │ 10-element Dataset ┐ │ ├────────────── dims ┤ ↓ y_dim_1 ├────────────────────┴─────────────── layers ┐ :y eltype: Float64 dims: y_dim_1 size: 10 ├────────────────────────────────────────────┴ metadata ┐ Dict{String, Any} with 1 entry: "created_at" => "2025-02-04T10:49:16.71""2025-02-12T01:28:19.299"

doctest failure in ~/.julia/packages/PosteriorStats/cl4WO/src/model_weights.jl:38-57 ```jldoctest model_weights; filter = [r"└.*"] julia> using ArviZExampleData julia> models = ( centered=load_example_data("centered_eight"), non_centered=load_example_data("non_centered_eight"), ); julia> elpd_results = map(models) do idata log_like = PermutedDimsArray(idata.log_likelihood.obs, (2, 3, 1)) return loo(log_like) end; ┌ Warning: 1 parameters had Pareto shape values 0.7 < k ≤ 1. Resulting importance sampling estimates are likely to be unstable. └ @ PSIS ~/.julia/packages/PSIS/... julia> model_weights(elpd_results; method=Stacking()) |> pairs pairs(::NamedTuple) with 2 entries: :centered => 5.34175e-19 :non_centered => 1.0 ``` Subexpression: model_weights(elpd_results; method=Stacking()) |> pairs Evaluated output: pairs(::NamedTuple) with 2 entries: :centered => 3.50538e-31 :non_centered => 1.0 Expected output: pairs(::NamedTuple) with 2 entries: :centered => 5.34175e-19 :non_centered => 1.0 diff = Warning: Diff output requires color. pairs(::NamedTuple) with 2 entries: :centered => 5.34175e-19 3.50538e-31 :non_centered => 1.0