Update initial proposal ll in hiatus case

brenhinkeller · Jul 16, 2024 · 3780a5c · 3780a5c · brenhinkeller · Jul 16, 2024
1 parent 971dd92
commit 3780a5c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "Chron"
 uuid = "68885b1f-77b5-52a7-b2e7-6a8014c56b98"
 authors = ["C. Brenhin Keller <[email protected]>"]
-version = "0.6.2"
+version = "0.6.3"
 
 [deps]
 DelimitedFiles = "8bb1440f-4735-579b-a4ab-409b98df4dab"

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