generated from opensafely/research-template
/
measures_test_7.py
67 lines (45 loc) · 1.68 KB
/
measures_test_7.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
###################################################
# This script creates monthly counts/rates of opioid
# prescribing for any opioid prescribing by mode of administration
###################################################
from ehrql import months, INTERVAL, Measures, case, when
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
codes = ["322524004", "322523005", "322526002", "322525003", "322527006"]
dia_in_interval = (
medications.where(medications.dmd_code.is_in(codes))
.where(medications.date.is_on_or_between(index_date, INTERVAL.end_date))
)
any_dia_rx = dia_in_interval.exists_for_patient()
first_dia_rx = dia_in_interval.sort_by(medications.date).last_for_patient().dmd_code
##########
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())
)
#########################
## Overall
measures.define_measure(
name="count",
numerator=any_dia_rx,
denominator=denominator,
group_by={"dmd_code": first_dia_rx}
)