Center acceptancedist around 0.58

OpenThermochronology · Nov 25, 2022 · ee89eab · ee89eab
1 parent 22066ec
commit ee89eab
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/test/testcomplete.jl b/test/testcomplete.jl
@@ -102,14 +102,17 @@ llmean = mean(@view(lldist[model.burnin:end]))
 @info "Mean ll: $llmean"
 
 @test isa(acceptancedist, AbstractVector{Bool})
-@test isapprox(mean(acceptancedist), 0.5, atol=0.3)
+@test isapprox(mean(acceptancedist), 0.58, atol=0.3)
 @info "Mean acceptance rate: $(mean(acceptancedist[model.burnin:end]))"
 
 @test isa(ndist, AbstractVector{Int})
 @test minimum(ndist) >= 0
 @test maximum(ndist) <= model.maxpoints
 @info "Mean npoints: $(mean(ndist[model.burnin:end]))"
 
+@info "Mean σⱼₜ: $(mean(σⱼtdist[model.burnin:end]))"
+@info "Mean σⱼT: $(mean(σⱼTdist[model.burnin:end]))"
+
 ## ---
 detail = DetailInterval(
     agemin = 0, # Youngest end of detail interval
@@ -173,7 +176,7 @@ llmean = mean(@view(lldist[model.burnin:end]))
 @info "Mean ll: $llmean"
 
 @test isa(acceptancedist, AbstractVector{Bool})
-@test isapprox(mean(acceptancedist), 0.5, atol=0.3)
+@test isapprox(mean(acceptancedist), 0.58, atol=0.3)
 @info "Mean acceptance rate: $(mean(acceptancedist[model.burnin:end]))"
 
 @test isa(ndist, AbstractVector{Int})