generated from opensafely/research-template
/
plot.do
69 lines (64 loc) · 1.78 KB
/
plot.do
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
* Open a log file
cap log close
log using ./logs/plot, replace t
clear
use ./output/main.dta
*follow-up time and events*
stset end_date , origin(start_date) failure(failure==1)
keep if _st==1
tab _t,m
tab _t drug,m col
by drug, sort: sum _t ,de
tab _t drug if failure==1,m col
tab failure drug,m col
*K-M curve*
sort drug _t
by drug, sort: gen _dn=sum(_d)
gen _tn=1 if drug==0&_dn<=5
replace _tn=1 if drug==1&_dn<=5
replace _tn=4 if drug==0&_dn>5&_dn<=10
replace _tn=2 if drug==1&_dn>5&_dn<=10
replace _tn=8 if drug==0&_dn>10&_dn<=15
replace _tn=4 if drug==1&_dn>10&_dn<=15
replace _tn=12 if drug==0&_dn>15&_dn<=20
replace _tn=7 if drug==1&_dn>15&_dn<=20
replace _tn=28 if drug==0&_dn>20
replace _tn=14 if drug==1&_dn>20&_dn<=25
replace _tn=20 if drug==1&_dn>25&_dn<=30
replace _tn=28 if drug==1&_dn>30
replace _t=_tn
sts graph, by(drug) ylabel(.98(.01)1)
graph export ./output/kmcurve.svg, as(svg) replace
*monthly count*
gen month=month(start_date)
tab month drug, col
clear
use ./output/main_mol.dta
*follow-up time and events*
stset end_date , origin(start_date) failure(failure==1)
keep if _st==1
tab _t,m
tab _t drug,m col
by drug, sort: sum _t ,de
tab _t drug if failure==1,m col
tab failure drug,m col
*K-M curve*
sort drug _t
by drug, sort: gen _dn=sum(_d)
gen _tn=4 if drug==0&_dn<=5
replace _tn=1 if drug==1&_dn<=5
replace _tn=8 if drug==0&_dn>5&_dn<=10
replace _tn=2 if drug==1&_dn>5&_dn<=10
replace _tn=28 if drug==0&_dn>10
replace _tn=4 if drug==1&_dn>10&_dn<=15
replace _tn=7 if drug==1&_dn>15&_dn<=20
replace _tn=14 if drug==1&_dn>20&_dn<=25
replace _tn=20 if drug==1&_dn>25&_dn<=30
replace _tn=28 if drug==1&_dn>30
replace _t=_tn
sts graph, by(drug) ylabel(.98(.01)1)
graph export ./output/kmcurve_mol.svg, as(svg) replace
*monthly count*
gen month=month(start_date)
tab month drug, col
log close