Update initial proposal ll in hiatus case

brenhinkeller · Jul 16, 2024 · 296d931 · 296d931
1 parent 971dd92
commit 296d931
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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