Top panel: Distribution of pull statistic (Eq. (23)) computed from 125 f_NL = 12 simulations (purple). A Gaussian fit to the central region (−3 < Pull < 3; dashed line) yields a mean of −0.20 and standard deviation of 1.31. Bottom panel: Mean and standard deviation of pull distribution as a function of s_max, the largest radial bin included in the analysis. The dashed vertical line marks the adopted value s_max = 275 h⁻¹ Mpc used in the test case.