Fix diagnostic if model won't fit

OHDSI · Oct 7, 2024 · 17cec00 · 17cec00
1 parent 65ba03f
commit 17cec00
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/R/Diagnostics.R b/R/Diagnostics.R
@@ -529,6 +529,9 @@ computeExposureChangeP <- function(sccsData, studyPopulation, exposureEraId = NU
                                                modelType = "cpr",
                                                quiet = TRUE)
   fit <- Cyclops::fitCyclopsModel(cyclopsData)
+  if (fit$return_flag != "SUCCESS") {
+    return(NA)
+  }
   fit$log_likelihood
   logRr <- coef(fit)
   if (logRr >= bounds[1] && logRr <= bounds[2]) {