generated from opensafely/research-template
/
an_covidvsnoncovid_full.do
93 lines (78 loc) · 2.97 KB
/
an_covidvsnoncovid_full.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
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
cap log close
log using ./analysis/output/an_covidvsnoncovid_full, replace t
use ./analysis/cr_create_analysis_dataset_STSET_ONSCSDEATH.dta, clear
gen coviddeath = onsdeath==1 if _d==1
replace coviddeath= 0 if _d==0
gen noncoviddeath = onsdeath>1 if _d==1
replace noncoviddeath= 0 if _d==0
**MULTIVARIATE
cap prog drop baselogistic
prog define baselogistic
syntax , outcome(string) age(string) bp(string) [ethnicity(real 0) if(string)]
if `ethnicity'==1 local ethnicity "i.ethnicity"
else local ethnicity
cap logistic `outcome' `age' ///
i.male ///
i.obese4cat ///
i.smoke_nomiss ///
`ethnicity' ///
i.imd ///
`bp' ///
i.chronic_respiratory_disease ///
i.asthmacat ///
i.chronic_cardiac_disease ///
i.diabcat ///
i.cancer_exhaem_cat ///
i.cancer_haem_cat ///
i.chronic_liver_disease ///
i.stroke_dementia ///
i.other_neuro ///
i.reduced_kidney_function_cat ///
i.organ_transplant ///
i.spleen ///
i.ra_sle_psoriasis ///
i.other_immunosuppression
end
baselogistic, outcome(coviddeath) age("age1 age2 age3") bp("i.htdiag_or_highbp") ethnicity(0)
if _rc==0{
estimates
estimates save ./analysis/output/models/an_covidvsnoncovid_full_covid_agespline_bmicat_noeth, replace
}
else di "WARNING AGE SPLINE MODEL DID NOT FIT (OUTCOME `outcome')"
*Age group model (not adj ethnicity)
baselogistic, outcome(coviddeath) age("ib3.agegroup") bp("i.htdiag_or_highbp") ethnicity(0)
if _rc==0{
estimates
estimates save ./analysis/output/models/an_covidvsnoncovid_full_covid_agegroup_bmicat_noeth, replace
}
else di "WARNING GROUP MODEL DID NOT FIT (OUTCOME `outcome')"
*Complete case ethnicity model
baselogistic, outcome(coviddeath) age("age1 age2 age3") bp("i.htdiag_or_highbp") ethnicity(1)
if _rc==0{
estimates
estimates save ./analysis/output/models/an_covidvsnoncovid_full_covid_agespline_bmicat_CCeth, replace
}
else di "WARNING CC ETHNICITY MODEL WITH AGESPLINE DID NOT FIT (OUTCOME `outcome')"
***
baselogistic, outcome(noncoviddeath) age("age1 age2 age3") bp("i.htdiag_or_highbp") ethnicity(0)
if _rc==0{
estimates
estimates save ./analysis/output/models/an_covidvsnoncovid_full_noncovid_agespline_bmicat_noeth, replace
}
else di "WARNING AGE SPLINE MODEL DID NOT FIT (OUTCOME `outcome')"
*Age group model (not adj ethnicity)
baselogistic, outcome(noncoviddeath) age("ib3.agegroup") bp("i.htdiag_or_highbp") ethnicity(0)
if _rc==0{
estimates
estimates save ./analysis/output/models/an_covidvsnoncovid_full_noncovid_agegroup_bmicat_noeth, replace
}
else di "WARNING GROUP MODEL DID NOT FIT (OUTCOME `outcome')"
*Complete case ethnicity model
baselogistic, outcome(noncoviddeath) age("age1 age2 age3") bp("i.htdiag_or_highbp") ethnicity(1)
if _rc==0{
estimates
estimates save ./analysis/output/models/an_covidvsnoncovid_full_noncovid_agespline_bmicat_CCeth, replace
}
else di "WARNING CC ETHNICITY MODEL WITH AGESPLINE DID NOT FIT (OUTCOME `outcome')"
***
log close