*final analysis (combined version) cd C:\Users\jcanner1\Documents\TMJ // Figure 1 Flowchart (Joe) use MUDS_PCO_data, clear keep if TMJ==1 count *511 total patients with TMJ count if (age==. | sex=="") *31 patients didn't report age and/or sex gen demogx=(age==. | sex=="") count if painranking==. *7 patients didn't report pain ranking gen painx=(painranking==.) tab demogx painx, m count if (age==. | sex=="" | painranking==.) *38 patients didn't report age, sex, and/or painranking (excluded from analysis) gen demogpainx=(age==. | sex=="" | painranking==.) count *473 patients with age, sex, and pain count if (decisionrank_1==. | decisionrank_2==. | decisionrank_3==. | decisionrank_4==. | /// decisionrank_5==. | decisionrank_6==. | decisionrank_7==. ) *101 all ranks missing, excluded from analysis of potential benefits gen potbenx=(decisionrank_1==. | decisionrank_2==. | decisionrank_3==. | decisionrank_4==. | /// decisionrank_5==. | decisionrank_6==. | decisionrank_7==. ) count if icludeAll!=1 *110 missing some SE ranks, excluded from analysis of potential harms gen potharmx=(icludeAll!=1) tab demogpainx potbenx,m tab demogpainx potharmx,m count if TMJ==1 & icludeAll==1 & age!=. & sex!="" & painranking!=. *363 patients included in analysis of potential harms count if TMJ==1 & (decisionrank_1!=. & decisionrank_2!=. & decisionrank_3!=. & decisionrank_4!=. & /// decisionrank_5!=. & decisionrank_6!=. & decisionrank_7!=.) & age!=. & sex!="" & painranking!=. *372 patients included in analysis of potential benefit count if TMJ==1 & age!=. & sex!="" & painranking!=. & /// ((decisionrank_1!=. & decisionrank_2!=. & decisionrank_3!=. & decisionrank_4!=. & /// decisionrank_5!=. & decisionrank_6!=. & decisionrank_7!=.) | icludeAll==1) *388 patients included in either analysis of potential benefits or harms //table 2 (Asieh) sum agedx age if TMJ==1 & age!=. & sex!="" & painranking!=. & /// ((decisionrank_1!=. & decisionrank_2!=. & decisionrank_3!=. & decisionrank_4!=. & /// decisionrank_5!=. & decisionrank_6!=. & decisionrank_7!=.) | icludeAll==1),d tab sex if TMJ==1 & age!=. & sex!="" & painranking!=. & /// ((decisionrank_1!=. & decisionrank_2!=. & decisionrank_3!=. & decisionrank_4!=. & /// decisionrank_5!=. & decisionrank_6!=. & decisionrank_7!=.) | icludeAll==1) sum painranking disorders Medications_c Medications_p if TMJ==1 & age!=. & sex!="" & painranking!=. & /// ((decisionrank_1!=. & decisionrank_2!=. & decisionrank_3!=. & decisionrank_4!=. & /// decisionrank_5!=. & decisionrank_6!=. & decisionrank_7!=.) | icludeAll==1),d tab preuseinfo if TMJ==1 & age!=. & sex!="" & painranking!=. & /// ((decisionrank_1!=. & decisionrank_2!=. & decisionrank_3!=. & decisionrank_4!=. & /// decisionrank_5!=. & decisionrank_6!=. & decisionrank_7!=.) | icludeAll==1), m *barplot graph hbox decisionrank_1 decisionrank_2 decisionrank_4 decisionrank_3 decisionrank_7 decisionrank_5 decisionrank_6 if TMJ==1 & icludeAll==1 & age!=. & sex!="" & painranking!=., /// ascategory ylabel(#7, labels valuelabel) yscale(range(1 7)) legend(off) plotregion(fcolor(white) lcolor(white) lpattern(solid)) /// scheme(s2color) graphregion(fcolor(white) lcolor(white) lpattern(solid) ifcolor(none) ilpattern(solid)) // Figure 3 (Joe) use MUDS_PCO_data, clear keep if TMJ==1 & icludeAll==1 & age!=. & sex!="" & painranking!=. foreach var of varlist SErank?_* SErank??_* { // Denominator egen `var'_N=count(`var') // Count number of 1s and 7s egen `var'_tot1=total(`var'==1) egen `var'_tot7=total(`var'==7) // Proportions gen `var'_p1=`var'_tot1/`var'_N gen `var'_p7=`var'_tot7/`var'_N // Difference in proportions gen `var'_d17=`var'_p1-`var'_p7 set more off list `var'_* in 1 } // All records should now have the same values, so only need to keep 1 list *_d17 in 1 keep in 1 // Take the mean over all birthmonths for each SE forvalues x=1/21 { egen mean`x'=rowmean(SErank*_`x'_d17) local varlabel`x' : var lab SErank_`x' } list mean* // Reshape so that each SE is a separate row keep mean* gen id=1 reshape long mean, i(id) j(SE) // Create a variable for the SE description gen varlabel="" forvalues x=1/21 { replace varlabel="`varlabel`x''" in `x' } replace varlabel=strltrim(varlabel) format varlabel %-80s sort mean list export excel using MUDS_PCO_SErank.xlsx, first(var) replace // Supplement 7 (Asieh, reproduced by Joe) use MUDS_PCO_data, clear keep if TMJ==1 & age!=. & sex!="" & painranking!=. & /// ((decisionrank_1!=. & decisionrank_2!=. & decisionrank_3!=. & decisionrank_4!=. & /// decisionrank_5!=. & decisionrank_6!=. & decisionrank_7!=.) | icludeAll==1) tab birthmonthNum,m tab birthmonthNum sex, row foreach var of varlist age agedx painranking conditions disorders Medications_c Medications_p { preserve collapse (p50) med_`var'=`var' (p25) p25_`var'=`var' (p75) p75_`var'=`var', by(birthmonthNum) di "`var'" list restore } // Footnote 2 (missing age at diagnosis by month) tab birthmonthNum if mi(agedx) // Second to last paragraph of Results tab preuseinfo, m