Models

[annaschultze]

This pull request is misleadingly titled models, but is actually a finalisation of descriptives for both asthma and COPD populations.

Requesting review of the Stata code. Specific QC requests are:

• Please check all variable lists as these have changed over time
• Please check the exposure variable for both populations vs. the protocol
• Please comment on logic of code set-up, order of programs etc.
• QC of the egfr calculations

[DarthCTR]

@annaschultze Looks amazing!!!

A couple notes as I’m working through: COPD

• The setup in model.do is great. I know you probably want to reduce the copy/paste across study populations, but I think this is the perfect balance between having the most efficient code and having interpretable code
• I know this is fake data, but what are the expectations set for BMI, smoking, ethnicity, and IMD? They look a bit weird in the test table 1. BMI missing 52%, smoking missing 82%, ethnicity missing 25% and completely missing ethnicities 2 and 4, and IMD completely missing 2,3,and missing.
• In Table 1, meds are listed as no then yes. Once we hit CKD, they are listed as yes then no. Suggest being consistent down the table.
• In KM plot, suggest capitalising Death and Population, or lower case hospital. (I think I prefer the latter).
• Does it get too busy if we do 10 or 15 day ticks on the x axis of KM plot?

Asthma

• Table file doesn’t include “04” like it did for COPD. And I guess you don’t include “05” on KM plots. I think I like labelling the figures with the number of editor its created in.
• Same comments re: table 1 as above

This is great stuff.

[ianjdouglas]

Great work Anna! Just a few thoughts from me on specific dofiles:

01_cr_create_asthma_population Line 29 – should be <18 (it's correct in the code, just the comment needs to match)? For dead before index date are weonly using CPNS deaths, not ONS dead too?

01_cr_create_copd_population For dead before index date are weonly using CPNS deaths, not ONS dead too?

02_cr_create_asthma_exposure SABA only should be not on other asthma meds. Most have been excluded but if anyone is on LABA_LAMA combi product they won't be picked up at the moment. Although not an asthma combi due to LAMA, the LABA is an asthma drug and would mean they're not on SABA only. Also need to keep an eye out for anyone on SABA LAMA in the baseline group as I think that would be weird

Thanks v much and great to see so much progress!