>100% CMA_sliding_window

[dffyfe]

Hi Dan,

I've been comparing CMA7 to Ambers manual Adherence work and it seems pretty much bang on (I think any differences can be put down to rounding or human error). One thing that would be good as an addition is the option to show >100% adherence on the sliding window with all the overuse going into the last window. i.e. in the example below show 100%, 100%, 100%, 364%. Obviously if it was 1 window we'd want to show overuse in that.

GT100 <- tribble(
  ~PATIENT_ID,  ~DATE,  ~PERDAY,    ~CATEGORY,  ~DURATION,
  1,    "2020-06-01",   2,  "medA", 40,
  1,    "2020-06-29",   2,  "medA", 40,
  1,    "2020-07-27",   2,  "medA", 40,
  1,    "2020-08-24",   2,  "medA", 40,
  1,    "2020-09-21",   2,  "medA", 40,
  1,    "2020-10-19",   2,  "medA", 40,
  1,    "2020-11-16",   2,  "medA", 40,
  1,    "2020-12-14",   2,  "medA", 40,
  1,    "2021-01-11",   2,  "medA", 40,
  1,    "2021-02-08",   2,  "medA", 40,
  1,    "2021-03-08",   2,  "medA", 40,
  1,    "2021-04-05",   2,  "medA", 40,
  1,    "2021-05-03",   2,  "medA", 40,
  1,    "2021-05-31",   2,  "medA", 40,
  1,    "2021-06-28",   2,  "medA", 40,
  1,    "2021-07-26",   2,  "medA", 40,
  1,    "2021-08-23",   2,  "medA", 40,
  1,    "2021-09-20",   2,  "medA", 40,
  1,    "2021-10-18",   2,  "medA", 40,
  1,    "2021-11-15",   2,  "medA", 40,
  1,    "2021-12-13",   2,  "medA", 40
)

GT100$DATE <- as.Date(GT100$DATE)

cma7 <- CMA_sliding_window(CMA.to.apply = "CMA7",
                           data=GT100,
                           ID.colname = "PATIENT_ID",
                           event.date.colname = "DATE",
                           event.duration.colname = "DURATION",
                           medication.class.colname = "CATEGORY",
                           carry.only.for.same.medication = FALSE,
                           followup.window.start = as.Date("2020-06-01"),
                           followup.window.start.unit = "days",
                           followup.window.duration = 600,
                           followup.window.duration.unit = "days",
                           observation.window.start = as.Date("2021-01-01"), 
                           observation.window.start.unit = "days",
                           observation.window.duration = 365,
                           observation.window.duration.unit = "days",
                           sliding.window.start = 0,
                           sliding.window.start.unit = "days",
                           sliding.window.duration = 91,
                           sliding.window.duration.unit = "days",
                           #sliding.window.step.duration = 91,
                           #sliding.window.step.unit = "days",
                           sliding.window.no.steps = 4,
                           date.format = "%Y-%m-%d")

plot(cma7, show.period = "dates", 
     period.in.days = 28,
     show.legend = FALSE,
     ylab = c(withCMA = "Section - Drug Name"),
     cex.axis = 1,
     title = "CHIno",
     CMA.cex = 0.8,
     plot.CMA = T,
     plot.CMA.as.histogram = F,
     CMA.plot.ratio = 0.000000000000000000000001,
     min.plot.size.in.characters.horiz=0,
     min.plot.size.in.characters.vert=0)

Best,

David

[ddediu]

Actually seems to be a more general question about a new (?) CMA that can allow (controlled) overuse...