CompatHelper: bump compat for AbstractMCMC to 4, (keep existing compat)

Project.toml

@@ -11,7 +11,7 @@ Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 Statistics = "10745b16-79ce-11e8-11f9-7d13ad32a3b2"
 
 [compat]
-AbstractMCMC = "3.2"
+AbstractMCMC = "3.2, 4"
 ArrayInterface = "2, 3, 4, 5"
 Distributions = "0.22, 0.23, 0.24, 0.25"
 julia = "1"

test/simple.jl

@@ -33,6 +33,11 @@
             @test all(length(x) == 2_000 for x in samples)
             @test mean(mean, samples) ≈ μ atol = 0.05
             @test mean(var, samples) ≈ σ² atol = 0.05
+
+            # initial parameter
+            init_x = randn(5)
+            samples = sample(ESSModel(prior, ℓ), ESS(), alg, 10, 5; progress=false, init_params=init_x)
+            @test map(first, samples) == init_x
         end
 
         # initial parameter
@@ -66,6 +71,11 @@
             @test all(length(x) == 2_000 for x in samples)
             @test mean(mean, samples) ≈ μ atol = 0.05
             @test mean(var, samples) ≈ σ² atol = 0.05
+
+            # initial parameter
+            init_x = randn(5)
+            samples = sample(ESSModel(prior, ℓ), ESS(), alg, 10, 5; progress=false, init_params=init_x)
+            @test map(first, samples) == init_x
         end
 
         # initial parameter
@@ -100,6 +110,11 @@
             @test all(length(x) == 2_000 for x in samples)
             @test mean(mean, samples) ≈ μ atol = 0.05
             @test mean(var, samples) ≈ σ² atol = 0.05
+
+            # initial parameter
+            init_x = [randn(1) for _ in 1:5]
+            samples = sample(ESSModel(prior, ℓvec), ESS(), alg, 10, 5; progress=false, init_params=init_x)
+            @test map(first, samples) == init_x
         end
 
         # initial parameter
@@ -136,6 +151,11 @@
             @test all(length(x) == 2_000 for x in samples)
             @test mean(mean, samples) ≈ μ atol = 0.05
             @test mean(var, samples) ≈ σ² atol = 0.05
+
+            # initial parameter
+            init_x = [randn(1) for _ in 1:5]
+            samples = sample(ESSModel(prior, ℓvec), ESS(), alg, 10, 5; progress=false, init_params=init_x)
+            @test map(first, samples) == init_x
         end
 
         # initial parameter