tweak default max proposal limit and improve max proposal error msg

Author: Red-Portal

src/SliceSampling.jl

@@ -20,6 +20,8 @@ export sample, MCMCThreads, MCMCDistributed, MCMCSerial
 # Interfaces
 abstract type AbstractSliceSampling <: AbstractMCMC.AbstractSampler end
 
+const DEFAULT_MAX_PROPOSALS = 10_000
+
 """
     struct Transition
 
@@ -52,10 +54,9 @@ Return the initial sample for the `model` using the random number generator `rng
 """
 function initial_sample(::Random.AbstractRNG, ::Any)
     error(
-        "`initial_sample` is not implemented but an initialization wasn't provided. " *
+        "`initial_sample` is not implemented but an initialization wasn't provided. ",
         "Consider supplying an initialization to `initial_params`."
     )
-    println("fuck!!!")
 end
 
 # If target is from `LogDensityProblemsAD`, unwrap target before calling `initial_sample`.
@@ -66,10 +67,15 @@ initial_sample(
 ) = initial_sample(rng, parent(wrap))
 
 function exceeded_max_prop(max_prop::Int)
-    error("Exceeded maximum number of proposal $(max_prop).\n", 
-          "Here are possible causes:\n",
-          "- The model might be broken or pathologic.\n",
-          "- There might be a bug in the sampler.")
+    error("Exceeded maximum number of proposal $(max_prop), ",
+          "which indicates an acceptance rate less than $(1/max_prop*100)%. ", 
+          "A quick fix is to increase `max_prop`, ",
+          "but an acceptance rate that is too low often indicates that there is a problem. ",
+          "Here are some possible causes:\n",
+          "- The model might be broken or degenerate. (most likely cause)\n",
+          "- The initialization is pathologic. (try supplying a different `initial_params`)\n",
+          "- There might be a bug in the sampler. (if this is suspected, file an issue to `SliceSampling`)\n"
+          )
 end
 
 ## Univariate Slice Sampling Algorithms

src/multivariate/gibbspolar.jl

@@ -9,14 +9,14 @@ Gibbsian polar slice sampling algorithm by P. Schär, M. Habeck, and D. Rudolf [
 
 # Keyword Arguments
 - `w::Real`: Initial window size for the radius shrinkage procedure
-- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `typemax(Int)`).
+- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `10^6`).
 """
 struct GibbsPolarSlice{W <: Real} <: AbstractMultivariateSliceSampling
     w::W
     max_proposals::Int
 end
 
-GibbsPolarSlice(w::Real; max_proposals::Int = typemax(Int)) = GibbsPolarSlice(w, max_proposals)
+GibbsPolarSlice(w::Real; max_proposals::Int = 10^6) = GibbsPolarSlice(w, max_proposals)
 
 struct GibbsPolarSliceState{T <: Transition, R <: Real, D <: AbstractVector} 
     "Current [`Transition`](@ref)."

src/multivariate/latent.jl

@@ -8,14 +8,14 @@ Latent slice sampling algorithm by Li and Walker[^LW2023].
 - `beta::Real`: Beta parameter of the Gamma distribution of the auxiliary variables.
 
 # Keyword Arguments
-- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `typemax(Int)`).
+- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `$(DEFAULT_MAX_PROPOSALS)`).
 """
 struct LatentSlice{B <: Real} <: AbstractMultivariateSliceSampling
     beta         ::B
     max_proposals::Int
 end
 
-function LatentSlice(beta::Real; max_proposals::Int = typemax(Int))
+function LatentSlice(beta::Real; max_proposals::Int = DEFAULT_MAX_PROPOSALS)
     @assert beta > 0 "Beta must be strictly positive"
     LatentSlice(beta, max_proposals)
 end

src/univariate/doublingout.jl

@@ -9,7 +9,7 @@ Univariate slice sampling by automatically adapting the initial interval through
 
 # Keyword Arguments
 - `max_doubling_out`: Maximum number of "doubling outs" (default: 8).
-- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `typemax(Int)`).
+- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `$(DEFAULT_MAX_PROPOSALS)`).
 """
 struct SliceDoublingOut{W <: Real} <: AbstractUnivariateSliceSampling
     window          ::W
@@ -20,7 +20,7 @@ end
 function SliceDoublingOut(
     window          ::Real;
     max_doubling_out::Int = 8,
-    max_proposals   ::Int = typemax(Int),
+    max_proposals   ::Int = DEFAULT_MAX_PROPOSALS,
 )
     @assert window > 0
     SliceDoublingOut(window, max_doubling_out, max_proposals)

src/univariate/fixedinterval.jl

@@ -8,7 +8,7 @@ Univariate slice sampling with a fixed initial interval (Scheme 2 by Neal[^N2003
 - `window::Real`: Proposal window.
 
 # Keyword Arguments
-- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `typemax(Int)`).
+- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `$(DEFAULT_MAX_PROPOSALS)`).
 """
 struct Slice{W <: Real} <: AbstractUnivariateSliceSampling
     window       ::W
@@ -17,7 +17,7 @@ end
 
 function Slice(
     window       ::Real;
-    max_proposals::Int = typemax(Int), 
+    max_proposals::Int = DEFAULT_MAX_PROPOSALS, 
 )
     @assert window > 0
     Slice(window, max_proposals)

src/univariate/steppingout.jl

@@ -9,7 +9,7 @@ Univariate slice sampling by automatically adapting the initial interval through
 
 # Keyword Arguments
 - `max_stepping_out::Int`: Maximum number of "stepping outs" (default: 32).
-- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `typemax(Int)`).
+- `max_proposals::Int`: Maximum number of proposals allowed until throwing an error (default: `$(DEFAULT_MAX_PROPOSALS)`).
 """
 struct SliceSteppingOut{W <: Real} <: AbstractUnivariateSliceSampling
     window          ::W
@@ -20,7 +20,7 @@ end
 function SliceSteppingOut(
     window          ::Real;
     max_stepping_out::Int = 32,
-    max_proposals   ::Int = typemax(Int),
+    max_proposals   ::Int = DEFAULT_MAX_PROPOSALS,
 )
     @assert window > 0
     SliceSteppingOut(window, max_stepping_out, max_proposals)

test/multivariate.jl

@@ -78,7 +78,7 @@ end
         LatentSlice(5),
 
         # Gibbsian polar slice sampling
-        GibbsPolarSlice(10),
+        GibbsPolarSlice(100),
     ]
         @testset "initial_params" begin
             model  = MultiModel(1.0, 1.0, [0.0])
@@ -87,7 +87,7 @@ end
 
             θ0    = [1.0, 0.1]
             chain = sample(
-                model,
+                model′,
                 sampler,
                 10;
                 initial_params=θ0,