From 296d9314b2cd2314baee0e25a560b6c92012db7c Mon Sep 17 00:00:00 2001 From: "C. Brenhin Keller" Date: Mon, 15 Jul 2024 21:32:54 -0400 Subject: [PATCH] Update initial proposal ll in hiatus case --- src/StratMetropolis.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/StratMetropolis.jl b/src/StratMetropolis.jl index 3c2bec3..e670ce2 100644 --- a/src/StratMetropolis.jl +++ b/src/StratMetropolis.jl @@ -783,8 +783,8 @@ end # Add log likelihood for hiatus duration - duration = @. min(model_ages[closest_hiatus_unique - 1] - model_ages[closest_hiatus_unique], Hiatus_duration) - ll += normpdf_ll(Hiatus_duration, Hiatus_duration_sigma, duration) + duration = @. model_ages[closest_hiatus_unique - 1] - model_ages[closest_hiatus_unique] + ll += normcdf_ll(Hiatus_duration, Hiatus_duration_sigma, duration) # Preallocate variables for MCMC proposals llₚ=ll