capture log close _all clear * Define some variables; local today=subinstr("`c(current_date)'"," ","",.) local time=subinstr("`c(current_time)'",":","",.) * stata version compatibility version 16.0 * log file number local vnum="1_08" * Syntax name local syntax "Descriptives01" * Projectname local project "C19a" * Open log file; log using "$logfilepath/StataLog_AW_`project'_`syntax'_`vnum'_`today'_`time'", replace text name(Merge`syntax') local logname=r(name) timer clear 1 * Change Log ******************************** * 1.00: Initial * 1.01: Added dispink5 * 1.02: Remove the two people that are missing on dispink3 ( if dispink2018_3!=.) * 1.03: Added more groups Sun3lmic and dispinklmic * 1.04-1.05: Added comments for public extraction & added additional descriptives * 1.06: Changes after reviewer comments, also added mena * 1.07: Added descriptives by age category * 1.08: Added descriptives without sex separation clear use "$datafilepath/AW_`project'_Prepare01" cd "$outputpath" * Replace event variable with only the deaths we care about gen codhold=. replace codhold=0 if _d==0 replace codhold=1 if _d==1 * all cause without Covid replace _d=0 if cod_cov==1 & _d==1 strate agecat, per(365250) output(groupagecatcovall`today'_`vnum', replace) strate sex, per(365250) output(groupSexcovall`today'_`vnum', replace) strate dispink2018_3, per(365250) output(groupdispinkcovall`today'_`vnum', replace) strate marital_status, per(365250) output(groupmaritalstatuscovall`today'_`vnum', replace) strate stolan, per(365250) output(groupstolancovall`today'_`vnum', replace) strate mena, per(365250) output(groupmenacovall`today'_`vnum', replace) strate lmic, per(365250) output(grouplmiccovall`today'_`vnum', replace) strate sun3, per(365250) output(groupsun3covall`today'_`vnum', replace) * Occurrence Exposure tables strate groupSexSun3, per(365250) output(groupSexSun3covall`today'_`vnum', replace) strate groupSexDispink, per(365250) output(groupSexDispink3covall`today'_`vnum', replace) strate groupSexMaritalStatus, per(365250) output(groupSexMaritalStatuscovall`today'_`vnum', replace) strate groupSexlmic, per(365250) output(groupSexlmiccovall`today'_`vnum', replace) strate groupSexmena, per(365250) output(groupSexmenacovall`today'_`vnum', replace) strate groupSexStolan, per(365250) output(groupSexStolancovall`today'_`vnum', replace) strate groupDispink3lmic, per(365250) output(groupDispink3lmiccovall`today'_`vnum', replace) strate groupSun3lmic, per(365250) output(groupSun3lmiccovall`today'_`vnum', replace) * Occurrence and exposure tables by agecat strate groupSexagecatSun3, per(365250) output(groupSexagecatSun3covall`today'_`vnum', replace) strate groupSexagecatDispink, per(365250) output(groupSexagecatDispink3covall`today'_`vnum', replace) strate groupSexagecatMaritalStatus, per(365250) output(groupSexagecatMaritalStatuscovall`today'_`vnum', replace) strate groupSexagecatlmic, per(365250) output(groupSexagecatlmiccovall`today'_`vnum', replace) strate groupSexagecatmena, per(365250) output(groupSexagecatmenacovall`today'_`vnum', replace) strate groupSexagecatStolan, per(365250) output(groupSexagecatStolancovall`today'_`vnum', replace) * Analyses for COVID-19 deaths replace _d=codhold replace _d=0 if cod_cov!=1 & _d==1 strate agecat, per(365250) output(groupagecatcov`today'_`vnum', replace) strate sex, per(365250) output(groupSexcov`today'_`vnum', replace) strate dispink2018_3, per(365250) output(groupdispinkcov`today'_`vnum', replace) strate marital_status, per(365250) output(groupmaritalstatuscov`today'_`vnum', replace) strate stolan, per(365250) output(groupstolancov`today'_`vnum', replace) strate mena, per(365250) output(groupmenacov`today'_`vnum', replace) strate lmic, per(365250) output(grouplmiccov`today'_`vnum', replace) strate sun3, per(365250) output(groupsun3cov`today'_`vnum', replace) * Occurrence Exposure tables strate groupSexSun3, per(365250) output(groupSexSun3cov`today'_`vnum', replace) strate groupSexDispink, per(365250) output(groupSexDispink3cov`today'_`vnum', replace) strate groupSexMaritalStatus, per(365250) output(groupSexMaritalStatuscov`today'_`vnum', replace) strate groupSexlmic, per(365250) output(groupSexlmiccov`today'_`vnum', replace) strate groupSexmena, per(365250) output(groupSexmenacov`today'_`vnum', replace) strate groupSexStolan, per(365250) output(groupSexStolancov`today'_`vnum', replace) strate groupDispink3lmic, per(365250) output(groupDispink3lmiccov`today'_`vnum', replace) strate groupSun3lmic, per(365250) output(groupSun3lmiccov`today'_`vnum', replace) * Occurrence and exposure tables by agecat strate groupSexagecatSun3, per(365250) output(groupSexagecatSun3cov`today'_`vnum', replace) strate groupSexagecatDispink, per(365250) output(groupSexagecatDispink3cov`today'_`vnum', replace) strate groupSexagecatMaritalStatus, per(365250) output(groupSexagecatMaritalStatuscov`today'_`vnum', replace) strate groupSexagecatlmic, per(365250) output(groupSexagecatlmiccov`today'_`vnum', replace) strate groupSexagecatmena, per(365250) output(groupSexagecatmenacov`today'_`vnum', replace) strate groupSexagecatStolan, per(365250) output(groupSexagecatStolancov`today'_`vnum', replace) * Drop first episodes for those who celebrated birthday drop if logistReg==0 * Tabulate relevant variables tab1 sex sun3 dispink2018_3 lmic mena stolan marital_status agecat, m timer list 1 log close `logname'