From 69314005184a36ab7e6591a7b0448a6ce9103412 Mon Sep 17 00:00:00 2001
From: alexmccreight <57416850+alexmccreight@users.noreply.github.com>
Date: Sun, 9 Nov 2025 11:17:23 -0600
Subject: [PATCH] add point mass at zero

---
 R/sufficient_stats_methods.R | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/R/sufficient_stats_methods.R b/R/sufficient_stats_methods.R
index 06c095f9..c5078193 100644
--- a/R/sufficient_stats_methods.R
+++ b/R/sufficient_stats_methods.R
@@ -357,7 +357,7 @@ update_variance_components.ss <- function(data, params, model, ...) {
     }
 
     # Region 1: Very dense near zero (tau2/100 to tau2/10)
-    small_grid <- exp(seq(log(mom_result$tau2/100), log(mom_result$tau2/10), length.out = 11))
+    small_grid <- exp(seq(log(mom_result$tau2/100), log(mom_result$tau2/10), length.out = 10))
 
     # Region 2: Moderate density around tau2 (tau2/10 to tau2*3)
     medium_grid <- exp(seq(log(mom_result$tau2/10), log(mom_result$tau2*3), length.out = 6))
@@ -365,7 +365,7 @@ update_variance_components.ss <- function(data, params, model, ...) {
     # Region 3: Sparse coverage towards minimum sparse effect
     large_grid <- exp(seq(log(mom_result$tau2*3), log(min_sparse_var), length.out = 3))
 
-    est_sa2 <- c(unique(c(small_grid, medium_grid, large_grid)))
+    est_sa2 <- c(0, unique(c(small_grid, medium_grid, large_grid)))
 
     # Call mr.ash with residuals
     mrash_output <- mr.ash.alpha::mr.ash(