use "E:\covid\Covid_Smoking\Meta\Paper 2\BMC\COVID_DataSingleFile.dta", clear /*All 47 studies*/ metan lnodd lnlci lnuci, eform lcols (Study) effect(OR) random xlabel (.5, 1,2,5) favours ("Less severe" # "More severe") title("Smoking and COVID-19 Disease Progression", size (textsizestyle) linegap(1)) textsize (120) /*Exclude high risk studies*/ metan lnodd lnlci lnuci if averagescore <=1, eform lcols (Study) effect(OR) random xlabel (.5, 1,2,5) favours ("Less severe" # "More severe") title("Smoking and COVID-19 Disease Progression" "Exclude 5 high-risk studies", size (textsizestyle) linegap(1)) textsize (120) /*lowess*/ metareg lnodd meanage, wsse(se) metareg lnodd meanage, wsse(se) eform lowess lnodd meanage /*Smoking and COVID Mortality*/ metan mlnodd mlnlci mlnuci, eform lcols (Study) effect(OR) random xlabel (.5, 1,2,5, 10) favours ("Less severe" # "More severe") title("Smoking and COVID-19 Mortality", size (textsizestyle) linegap(1)) textsize (150) /*Current smoking and Covid progression*/ metan current_severe current_nonsevere never_severe never_nonsevere, or lcols (Study) random xlabel (.2,.5, 1,2,5, 10) favours ("Less severe" # "More severe") title("Current Smoker and COVID-19 Disease Progression", size (textsizestyle) linegap(1)) textsize (150) /*Unadjusted ODDs*/ metan unlnodd unlnlci unlnuci, eform lcols (Study) effect(OR) random xlabel (.2,.5, 1,2,5, 10, 20) favours ("Less severe" # "More severe") title("Smoking and COVID-19 Disease Progression" "Unadjusted ORs", size (textsizestyle) linegap(1)) textsize (150) /*Adjusted ODDs*/ metan lnmodd lnmlci lnmcui , eform lcols (Study) effect(OR) random xlabel (.2,.5, 1,2,5, 10, 20) favours ("Less severe" # "More severe") title("Smoking and COVID-19 Disease Progression" "Adjusted ORs", size (textsizestyle) linegap(1)) textsize (150) clear