Merge pull request #248 from stan-dev/fix_min_neff

show NA for n_eff/ESS if k > k_threshold
stan-dev · Feb 21, 2024 · 8108ac4 · 8108ac4
2 parents 68ee019 + ef4d236
commit 8108ac4
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/R/diagnostics.R b/R/diagnostics.R
@@ -130,6 +130,7 @@ pareto_k_table <- function(x) {
   S <- dim(x)[1]
   k_threshold <- ps_khat_threshold(S)
   kcut <- k_cut(k, k_threshold)
+  n_eff[k>k_threshold] <- NA
   min_n_eff <- min_n_eff_by_k(n_eff, kcut)
   count <- table(kcut)
   out <- cbind(