generated from opensafely/research-template
/
measures_test.py
98 lines (71 loc) · 2.47 KB
/
measures_test.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
##################################################
# This script creates monthly counts/rates of opioid
# prescribing for any opioid prescribing by mode of administration
###################################################
from ehrql import Dataset, months, INTERVAL, Measures
from ehrql.tables.beta.tpp import (
patients,
medications,
practice_registrations)
import codelists
##########
from argparse import ArgumentParser
parser = ArgumentParser()
parser.add_argument("--start-date", type=str)
parser.add_argument("--intervals", type=int)
args = parser.parse_args()
start_date = args.start_date
intervals = args.intervals
##########
index_date = INTERVAL.start_date
dataset = Dataset()
number_diamorph_rx = medications.where(
medications.dmd_code.is_in(codelists.diamorph_opioid_codes)
).where(
medications.date.is_on_or_between(index_date, INTERVAL.end_date)
).count_for_patient()
number_par_rx = medications.where(
medications.dmd_code.is_in(codelists.par_opioid_codes)
).where(
medications.date.is_on_or_between(index_date, INTERVAL.end_date)
).count_for_patient()
diamorph_opioid_any = medications.where(
medications.dmd_code.is_in(codelists.diamorph_opioid_codes)
).where(
medications.date.is_on_or_between(index_date, INTERVAL.end_date)
).exists_for_patient()
par_opioid_any = medications.where(
medications.dmd_code.is_in(codelists.par_opioid_codes)
).where(
medications.date.is_on_or_between(index_date, INTERVAL.end_date)
).exists_for_patient()
measures = Measures()
measures.define_defaults(intervals=months(intervals).starting_on(start_date))
denominator = (
(patients.age_on(index_date) >= 18)
& (patients.age_on(index_date) < 110)
& ((patients.sex == "male") | (patients.sex == "female"))
& (patients.date_of_death.is_after(index_date) | patients.date_of_death.is_null())
& (practice_registrations.for_patient_on(index_date).exists_for_patient())
)
#########################
measures.define_measure(
name="par_opioid",
numerator=number_par_rx,
denominator=denominator
)
measures.define_measure(
name="diamorph_opioid",
numerator=number_diamorph_rx,
denominator=denominator
)
measures.define_measure(
name="par_opioid_any",
numerator=par_opioid_any,
denominator=denominator
)
measures.define_measure(
name="diamorph_opioid_any",
numerator=diamorph_opioid_any,
denominator=denominator
)