View source: R/matching_function.R
prune_event_durations | R Documentation |
Flags or removes leftover supply durations after dosage changes, the end of a special period,
or treatment interruption.
The function accepts the raw list output of compute_event_durations
and additional arguments
to specify event durations that need to be removed.
prune_event_durations( data, include = c("special periods", "treatment interruptions", "dosage changes"), medication.class.colnames = data$medication.class.colnames, days.within.out.date.1, days.within.out.date.2, keep.all = TRUE, suppress.warnings = FALSE, return.data.table = FALSE, ... )
data |
A |
include |
A |
medication.class.colnames |
A |
days.within.out.date.1 |
event durations from before the dosage change, special period, or treatment interruptions are removed if there is a new dispensing event within the number of days specified as integer after the dosage change or end of the special period/treatment interruption. |
days.within.out.date.2 |
event durations from before dosage change, special period, or treatment interruption are removed if there is NO new dispensing event within the number of days specified as integer after the dosage change or end of the special period/treatment interruption. |
keep.all |
Logical, should events be kept and marked for removal?
If |
suppress.warnings |
Logical, if |
return.data.table |
Logical, if |
... |
other possible parameters. |
Dosage changes, special periods, and treatment interruptions may lead to overestimation of implementation, e.g. if patients get a refill after discharge from hospital and don't continue to use their previous supply. Likewise, it may also lead to overestimation of persistence, e.g. when patients discontinue treatments after the end of a special period or treatment interruption.
A data.frame
or data.table
, the pruned event_durations.
## Not run: # select medication class of interest and compute event durations disp_data <- durcomp.dispensing[ID == 3 & grepl("J01EE01", ATC.CODE)] presc_data <- durcomp.prescribing[ID == 3 & grepl("J01EE01", ATC.CODE)] # compute event durations event_durations_list <- compute_event_durations(disp.data = disp_data, presc.data = presc_data, special.periods.data = durcomp.hospitalisation, ID.colname = "ID", presc.date.colname = "DATE.PRESC", disp.date.colname = "DATE.DISP", date.format = "%Y-%m-%d", medication.class.colnames = c("ATC.CODE", "UNIT", "FORM"), total.dose.colname = "TOTAL.DOSE", presc.daily.dose.colname = "DAILY.DOSE", presc.duration.colname = "PRESC.DURATION", visit.colname = "VISIT", force.init.presc = TRUE, force.presc.renew = TRUE, split.on.dosage.change = TRUE, trt.interruption = "carryover", special.periods.method = "carryover", suppress.warnings = FALSE, return.data.table = TRUE, progress.bar = FALSE) # prune event durations event_durations <- prune_event_durations(event_durations_list, include = c("special periods"), medication.class.colnames = "ATC.CODE", days.within.out.date.1 = 7, days.within.out.date.2 = 30, keep.all = FALSE) ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.