Closed KateJohnson closed 5 years ago
Need to validate treatment by checking that exacerbation rate is lower in diagnosed than undiagnosed. Should create a validation function for treatment.
validate_treatment() is added, but need more rigorous tests of treatment effect. Try scenarios in which everyone is diagnosed, versus everyone undiagnosed, and see if you can recover the treatment effect in the rate of exacerbations.
Treatment effects need to be expanded to the full GOLD ABCD classification:
ABCD treatment classification added. Use validate_medication() to look at treatment trends over time.
Treatment is still not assigned probabilistically
Treatments adherence is fully implemented. Use
input <- init_input()
input$values$medication$medication_adherence <- 0.5
To change adherence, which proportionally influences treatment costs and treatment effectiveness (rate of exacerbations)
The coefficients come from the effect of the therapy on the rate of exacerbations.
This is implemented and working, but note that the impact of treatment on the rate of exacerbations has a lag of 1 annual event.