Minor tweak to event-dependent exposure simulations

OHDSI · Oct 7, 2024 · fd71ffd · fd71ffd
1 parent 3891ed9
commit fd71ffd
Showing 1 changed file with 9 additions and 2 deletions.
diff --git a/extras/EndOfExposureSimulations.R b/extras/EndOfExposureSimulations.R
@@ -30,7 +30,7 @@ for (trueRr in c(1, 2, 4)) {
                                                    maxBaselineRate = baseLineRate,
                                                    eraIds = 1,
                                                    usageRate = usageRate,
-                                                   usageRateSlope = usageRateSlope,,
+                                                   usageRateSlope = usageRateSlope,
                                                    meanPrescriptionDurations = 30,
                                                    sdPrescriptionDurations = 14,
                                                    simulationRiskWindows = list(rw),
@@ -55,7 +55,9 @@ writeLines(sprintf("Number of simulation scenarios: %d", length(scenarios)))
 # Run simulations ----------------------------------------------------------------------------------
 folder <- "e:/SccsEdeSimulations100"
 
-scenario = scenarios[[1]]
+scenario = scenarios[[20]]
+scenario$censorType
+
 simulateOne <- function(seed, scenario) {
   set.seed(seed)
   sccsData <- simulateSccsData(1000, scenario$settings)
@@ -101,7 +103,12 @@ simulateOne <- function(seed, scenario) {
       arrange(caseId, eraStartDay)
   } else if (scenario$censorType == "Permanent when exposed") {
     probability <- if_else(scenario$censorStrength == "Strong", 0.8, 0.25)
+    # Only outcomes during exposure can lead to exposure censoring:
     outcomeEras <- outcomeEras |>
+      inner_join(sccsData$eras |>
+                   filter(eraType == "rx"),
+                 by = join_by(caseId, between(outcomeDay, eraStartDay, eraEndDay))) |>
+      select(caseId, outcomeDay) |>
       filter(runif() > probability)
 
     filteredExposureEras <- sccsData$eras |>