*1. randomization evaluation (Appendix Table A3) *---------------------------------------------------- *prepare SP level data for randomization evaluate *-------------------------------------------------- use "[Location]\NHATS_Round_11_SP_File.dta", clear keep spid r11d2intvrage r5dgender rl5dracehisp w11anfinwgt0 w11varstrat w11varunit merge 1:1 spid using "[Location]\NHATS_Round_11_SP_Sen_Dem_File.dta" keep spid r11d2intvrage r5dgender rl5dracehisp r11dintvwrage r11ddeathage w11anfinwgt0 w11varstrat w11varunit merge 1:1 spid using "[Location]\NHATS_Round_5_SP_File_V2.dta" keep if _merge == 3 keep spid r11d2intvrage r5dgender rl5dracehisp r11dintvwrage el5higstschl r11ddeathage w11anfinwgt0 w11varstrat w11varunit merge 1:1 spid using "[Location]\NHATS_Round_1_SP_File.dta" keep if _merge==1 | _merge==3 gen sp_edush = . replace sp_edush = el5higstschl if el5higstschl != -1 replace sp_edush = el1higstschl if el5higstschl == -1 keep spid r11d2intvrage r5dgender rl5dracehisp r11dintvwrage el5higstschl sp_edush r11ddeathage w11anfinwgt0 w11varstrat w11varunit describe el5higstschl label copy el5higstschl sp_edush label values sp_edush sp_edush save spvar, replace *merge with NSOC data - NSOC OP file *-------------------------------------------------- cd "[Location]\Data" import delimited "[Location]\Data\nsoc_r11_sample_info_20230117.csv", clear save "[Location]\Data\nsoc_r11_sample_info_20230117.dta", replace tostring opid, replace format(%03.0f) save contact, replace cd "[Location]\Data\Data for paper" use "[Location]\NHATS_Round_11_OP_File.dta", clear gen op11tkplhlp = -1 replace op11tkplhlp = 1 if op11tkplhlp1 == 1 | op11tkplhlp2 == 1 foreach var in op11outhlp op11insdhlp op11bedhlp op11eathlp op11bathhlp op11toilhlp op11dreshlp /// op11launhlp op11shophlp op11mealhlp op11bankhlp op11medshlp /// op11tkplhlp op11moneyhlp op11dochlp op11insurhlp { recode `var' (-1=0) /*(.=0) */ } gen actsumr11 = op11outhlp + op11insdhlp + op11bedhlp + op11eathlp + op11bathhlp + op11toilhlp +op11dreshlp + op11launhlp +op11shophlp +op11mealhlp +op11bankhlp +op11medshlp + op11tkplhlp + op11moneyhlp + op11dochlp + op11insurhlp keep spid opid op11relatnshp r11dresid op11dgender actsumr11 merge m:1 spid using spvar drop _merge *merge with NSOC tracker data merge 1:1 spid opid using "[Location]\NSOC_R11_Crss_OP_Tracker_File.dta" keep if _merge == 3 drop _merge merge 1:1 spid opid using "[Location]\Data\contact.dta" *keep only eligible caregivers keep if op11dnsocelig == 1 *recode variables gen CGrelation = . replace CGrelation = 1 if op11relatnshp == 2 replace CGrelation = 2 if op11relatnshp >=3 & op11relatnshp <= 8 replace CGrelation = 3 if (op11relatnshp >=9 & op11relatnshp <= 29) | op11relatnshp == 91 replace CGrelation = 4 if (op11relatnshp >=30 & op11relatnshp <= 41) | op11relatnshp == 92 label define rela 1 "1 Spouse / partners" 2 "2 Children" 3 "3 Other relatives" 4 "4 Non-relatives" label values CGrelation rela gen sp_educ = . replace sp_educ = 1 if sp_edush == 1 | sp_edush == 2 | sp_edush == 3 replace sp_educ = 2 if sp_edush == 4 replace sp_educ = 3 if sp_edush >= 5 & sp_edush <= 9 label define edu 1 "1 Below high school" 2 "2 High school" 3 "3 Above high school" label values sp_educ edu gen sp_race = . replace sp_race = rl5dracehisp if rl5dracehisp == 1 | rl5dracehisp == 2 replace sp_race = 3 if rl5dracehisp == 4 replace sp_race = 4 if rl5dracehisp == 3 | rl5dracehisp == 5 label define race 1 "1 White, non-Hispanic" 2 "2 Black, non-Hispanic" 3 "3 Hispanic" 4 "4 Other" label values sp_race race gen sp_age = r11dintvwrage replace sp_age = r11ddeathage if r11dintvwrage == -1 *weighted svyset w11varunit [pweight=w11anfinwgt0], strata(w11varstrat) singleunit(scaled) svy: tab op11dgender op11dmixmode, col svy: tab CGrelation op11dmixmode, col phi op11dgender op11dmixmode, V cchi2 phi CGrelation op11dmixmode, V cchi2 svy: mean sp_age if op11dmixmode == 1 svy: mean sp_age if op11dmixmode == 2 svy: regress sp_age i.op11dmixmode svy: mean sp_age , over(op11dmixmode) estat sd esize two sp_age, by(op11dmixmode) svy: tab r5dgender op11dmixmode, col svy: tab sp_educ op11dmixmode, col svy: tab sp_race op11dmixmode, col phi r5dgender op11dmixmode, V cchi2 phi sp_educ op11dmixmode, V cchi2 phi sp_race op11dmixmode, V cchi2 svy: tab hasphone op11dmixmode, col svy: tab hasaddress op11dmixmode, col svy: tab hasemail op11dmixmode, col phi hasphone op11dmixmode, V cchi2 phi hasaddress op11dmixmode, V cchi2 phi hasemail op11dmixmode, V cchi2 ssc install stddiff help stddiff /************************************************************************/ *Prepare data for propensity score analysis /************************************************************************/ cd "[Location]\Data" import delimited "[Location]\Data\nsoc_r11_sample_info_20230117.csv", clear save "[Location]\Data\nsoc_r11_sample_info_20230117.dta", replace tostring opid, replace format(%03.0f) save contact, replace *this data has information on HasAddress HasPhone HasEmail use "[Location]\NHATS_Round_11_OP_File.dta", clear keep spid opid op11proxy op11dgender op11prsninhh op11outhlp op11insdhlp op11bedhlp op11tkplhlp1 op11tkplhlp2 op11launhlp op11shophlp op11mealhlp op11bankhlp op11moneyhlp op11eathlp op11bathhlp op11toilhlp op11dreshlp op11medshlp op11dochlp op11dochlpmst op11insurhlp op11spcaredfr op11chhlpfin op11hlpchfin op11ishelper op11helpsched op11numdayswk op11numdaysmn op11numhrsday op11paidhelpr op11sppayshlp op11govpayhlp op11inspayhlp op11othpayhlp op11payunit op11hourlypay op11weeklypay op11monthlypy op11progmpaid op11dhrsmth op11lmlhelper op11lmlpaidhelpr op11contctflg op11dflmiss op11dmissadd op11dprobdup op11ddupeid op11ddiffop op11dnotfield save nhatsop_R11var, replace use "[Location]\NHATS_Round_11_SP_File.dta", clear keep spid r11dresid r11dresidlml r11dpt1mode r11dintdays fl11pt2miss is11resptype r5dgender rl5dracehisp r11d2intvrage r11dmetnonmet te11cellphone te11othrphone te11computer te11compoth te11tablet te11emailtext te11oftnemail te11shoponli1 te11shoponli2 te11shoponli3 te11socialnet te11shoponli5 te11intrntmd1 te11intrntmd2 te11intrntmd3 te11intrntmd4 ha11shop ha11shoponlipers mc11regdoclyr mc11telehealth mc11tpersevr5 save nhatsSP_R11var, replace use "[Location]\NSOC_R11_Crss_OP_Tracker_File.dta", clear sort spid opid merge 1:1 spid opid using "[Location]\NSOC_R11_Crss_File.dta" keep if op11dnsoccomplete==1 drop _merge *destring opid, replace merge 1:1 spid opid using contact keep if _merge == 1 | _merge == 3 drop _merge merge 1:1 spid opid using nhatsop_R11var keep if _merge == 3 drop _merge merge m:1 spid using nhatsSP_R11var keep if _merge == 3 drop _merge *create weights for both caregivers to living and LML SPs gen w11cgfinwgt=0 replace w11cgfinwgt=w11cgfinwgt0 replace w11cgfinwgt=w11cglmlfinwgt0 if fl11spdied==1 /*control for age and relationship*/ /*recode age */ gen agegrp=0 replace agegrp=1 if chd11dage>=18 & chd11dage<=40 replace agegrp=2 if chd11dage>40 & chd11dage<=60 replace agegrp=3 if chd11dage>60 & chd11dage<=80 replace agegrp=4 if chd11dage>80 replace agegrp=6 if agegrp==0 /*missing age*/ tab agegrp label define age 1 "1 18-40" 2 "2 41-60" 3 "3 61-80" 4 "4 >80" 6 "6 Missing" label value agegrp age gen relgrp=0 replace relgrp=1 if c11relatnshp==2 /*spouse*/ replace relgrp=2 if c11relatnshp>=3 & c11relatnshp<=8 /*all children, including inlaw, step*/ replace relgrp=3 if (c11relatnshp>=19 & c11relatnshp<=22) | c11relatnshp==26 |c11relatnshp==28 | c11relatnshp==29 /*grandchild, nices nephews*/ replace relgrp=4 if c11relatnshp==25 | c11relatnshp==11 | c11relatnshp==10 | c11relatnshp==9 | c11relatnshp==15 | c11relatnshp==91 replace relgrp=5 if relgrp==0 tab relgrp label define rela 1 "1 spouse" 2 "2 children" 3 "3 grandchild/nices/nephews" 4 "4 other relatives" 5 "5 other nonrelatives" label value relgrp rela gen edu4cat = . replace edu4cat = 1 if chd11educ == 1 | chd11educ == 2 | chd11educ == 3 replace edu4cat = 2 if chd11educ == 4 replace edu4cat = 3 if chd11educ == 5 | chd11educ == 6 | chd11educ == 7 | chd11educ == 8 replace edu4cat = 4 if chd11educ == 9 replace edu4cat = 9 if edu4cat == . label define edu4cat 1 "1 HS, <=College" 4 "4 >College" label value edu4cat edu4cat tab edu4cat replace chd11dmartstat = 9 if chd11dmartstat < 0 gen childunder18 = 0 replace childunder18 = 1 if chd11numchu18 > 0 replace childunder18 = 9 if chd11chldlvng < 0 | chd11numchu18 == -6 replace chd11dnuminhh = . if chd11dnuminhh < 0 replace op11prsninhh = 2 if op11prsninhh == -1 replace te11tablet = . if te11tablet < 0 /*deceased SP not answer TE questions*/ replace te11othrphone = . if te11othrphone < 0 /*deceased SP not answer TE questions*/ replace te11computer = . if te11computer < 0 /*deceased SP not answer TE questions*/ replace te11cellphone = . if te11cellphone < 0 /*deceased SP not answer TE questions*/ replace ha11shop = . if ha11shop < 0 /*deceased SP not asked*/ gen resid = r11dresid replace resid = r11dresidlml if r11dresid == 6 gen shoponline = ha11shoponlipers replace shoponline = 4 if ha11shop == 3 | ha11shop == 5 replace shoponline = . if shoponline < 0 label define shoponline 1 "1 ONLY IN PERSON" 2 "2 ONLY ONLINE" 3 "3 IT VARIED / BOTH WAYS" 4 "4 someone else/not done last month" label value shoponline shoponline gen telhealth = mc11telehealth replace telhealth = 4 if mc11regdoclyr == 2 replace telhealth = . if telhealth < 0 label define telhealth 1 "1 ONLY IN PERSON" 2 "2 ONLY ONLINE" 3 "3 IT VARIED / BOTH WAYS" 4 "4 did not see doc last year" label value telhealth telhealth replace op11bankhlp = 0 if op11bankhlp == -1 replace r11d2intvrage = . if r11d2intvrage == -1 gen mode=0 replace mode=1 if op11dmixmode==1 & c11dintmode==2 replace mode=2 if op11dmixmode==1 & c11dintmode==1 replace mode=3 if op11dmixmode==2 replace mode=4 if op11dmixmode==1 & c11dintmode>2 /*16 cases used both modes*/ gen inimode_web = . replace inimode_web = 0 if c11dintmode == 1 | c11dintmode == 3 replace inimode_web = 1 if c11dintmode == 2 | c11dintmode == 4 replace cdc11hlphrsdy = . if cdc11hlphrsdy < 0 replace fl11spdied = 0 if fl11spdied == -1 foreach var of varlist che11lifemng che11gaveup che11feelconf che11lifelone che11lifegood che11rcvrqick che11easychng { replace `var' = . if `var' < 0 gen `var'_2cat = . replace `var'_2cat = 1 if `var' == 1 | `var' == 2 replace `var'_2cat = 0 if `var' == 3 | `var' == 4 } *wellbeing score foreach x of varlist che11lifemng che11feelconf che11lifegood che11easychng che11rcvrqick { replace `x' = . if `x' < 0 gen `x'new = 5 - `x' tab `x'new `x' } replace che11gaveup = . if che11gaveup < 0 replace che11lifelone = . if che11lifelone < 0 gen wbscore = che11lifemngnew + che11feelconfnew + che11gaveup + che11lifegoodnew + che11lifelone + che11easychngnew + che11rcvrqicknew tab wbscore gen mode_3grp = . replace mode_3grp = 1 if inimode_web == 1 replace mode_3grp = 2 if inimode_web == 0 & op11dmixmode == 1 replace mode_3grp = 3 if op11dmixmode == 2 label define mode_3grp 1 "1 Initial web" 2 "2 Initial phone" 3 "3 phone-only" label value mode_3grp mode_3grp save "[Location]\Data\Data for paper\data_4_p_score.dta", replace /************************************************************************/ *propensity score analysis /************************************************************************/ cd "[Location]\Data\Data for paper" use data_4_p_score, clear /*all*/ svyset c11varunit [pweight=w11cgfinwgt], strata(c11varstrat) singleunit(scaled) *subpop size: 958 svy, subpop(if op11dmixmode==1):logit inimode_web i.agegrp i.relgrp i.fl11spdied i.c11gender i.edu4cat i.chd11dmartstat i.childunder18 chd11dnuminhh i.crl11dcgracehisp i.hasaddress i.hasphone i.hasemail i.r11dmetnonmet i.op11prsninhh predict p1, pr somersd inimode_web p1 if op11dmixmode==1, tr(c) twoway histogram p1, color(*.5) || kdensity p1 ||, by(mode_3grp) save data_with_p_score, replace cd "[Location]\Data\Data for paper" use data_with_p_score, clear drop if mode_3grp == 2 *egen id = concat(spid opid), format(%10.0f) punct(:) *destring id, replace order spid opid gen id = _n save tomatch_grp1and3, replace //Bivariate test before matching use tomatch_grp1and3, clear replace chd11dage = . if chd11dage < 0 tab relgrp inimode_web if p1 !=., chi2 row tab fl11spdied inimode_web if p1 !=., chi2 row tab c11gender inimode_web if p1 !=., chi2 row tab edu4cat inimode_web if p1 !=., chi2 row tab chd11dmartstat inimode_web if p1 !=., chi2 row tab childunder18 inimode_web if p1 !=., chi2 row tab crl11dcgracehisp inimode_web if p1 !=., chi2 row tab hasaddress inimode_web if p1 !=., chi2 row tab hasphone inimode_web if p1 !=., chi2 row tab hasemail inimode_web if p1 !=., chi2 row tab r11dmetnonmet inimode_web if p1 !=., chi2 row tab op11prsninhh inimode_web if p1 !=., chi2 row stddiff chd11dage if agegrp < 6 & mode_3grp != 2, by(inimode_web) stddiff i.relgrp fl11spdied c11gender i.edu4cat i.chd11dmartstat childunder18 i.crl11dcgracehisp hasaddress hasphone hasemail r11dmetnonmet op11prsninhh chd11dnuminhh if mode_3grp != 2, by(inimode_web) mean chd11dage if agegrp2 < 6 & mode_3grp != 2 ttest chd11dage if agegrp2 < 6, by(inimode_web) ttest chd11dnuminhh, by(inimode_web) cd "[Location]\Data\Data for paper revision" use "[Location]\Data\Data for paper\data_4_p_score", clear replace chd11dage = . if chd11dage < 0 *appendix table 1: phone vs. phone *SMD before matching stddiff chd11dage i.relgrp fl11spdied c11gender i.edu4cat i.chd11dmartstat childunder18 i.crl11dcgracehisp hasaddress hasphone hasemail r11dmetnonmet op11prsninhh chd11dnuminhh if mode_3grp == 2 | mode_3grp == 3, by(mode_3grp) ttest chd11dage if mode_3grp == 2 | mode_3grp == 3, by(mode_3grp) ttest chd11dnuminhh if mode_3grp == 2 | mode_3grp == 3, by(mode_3grp) *. Propensity Score Analysis *---------------------------------------------------- cd "[Location]\Data\Data for paper" use tomatch_grp1and3, clear //PSM based on Logistic regression 1 //Scheme 1 - Nearest neighbor within caliper .25*SD use tomatch_grp1and3, clear drop if p1==. gen logit1=log((1-p1)/p1) sum logit1 display .25*1.503159 set seed 13455 generate x=uniform() sort x psmatch2 inimode_web , pscore(logit1) caliper(.37578975) noreplacement descending list spid opid inimode_web mode_3grp _pscore _treated _support _weight _id _n1 _nn _pdif in 1/5 sort _id g match=id[_n1] g treat=id if _nn==1 drop if treat== . sum treat keep treat match save scheme1, replace //scheme1 clear use scheme1, replace keep treat gen tx=1 rename treat id save t1,replace use scheme1, replace keep match gen tx=0 rename match id append using t1 sort id save t2, replace clear use tomatch_grp1and3, replace sort id merge 1:1 id using t2 tab tx drop if tx==. save post1, replace //bivarate analysis after matching cd "[Location]\Data\Data for paper" use post1, clear tab relgrp inimode_web if p1 !=., chi2 row tab fl11spdied inimode_web if p1 !=., chi2 row tab c11gender inimode_web if p1 !=., chi2 row tab edu4cat inimode_web if p1 !=., chi2 row tab chd11dmartstat inimode_web if p1 !=., chi2 row tab childunder18 inimode_web if p1 !=., chi2 row tab crl11dcgracehisp inimode_web if p1 !=., chi2 row tab hasaddress inimode_web if p1 !=., chi2 row tab hasphone inimode_web if p1 !=., chi2 row tab hasemail inimode_web if p1 !=., chi2 row tab r11dmetnonmet inimode_web if p1 !=., chi2 row tab op11prsninhh inimode_web if p1 !=., chi2 row mean chd11dnuminhh if p1 !=. & inimode_web == 1 mean chd11dnuminhh if p1 !=. & inimode_web == 0 replace chd11dage = . if chd11dage < 0 mean chd11dage ttest chd11dage if agegrp < 6, by(inimode_web) ttest chd11dnuminhh, by(inimode_web) *SMD stddiff chd11dage if agegrp < 6, by(inimode_web) stddiff i.relgrp fl11spdied c11gender i.edu4cat i.chd11dmartstat childunder18 i.crl11dcgracehisp hasaddress hasphone hasemail r11dmetnonmet op11prsninhh chd11dnuminhh, by(inimode_web) *phone in mixed-mode vs more likely to choose phone in phone-only cd "[Location]\Data\Data for paper" use data_with_p_score, clear drop if mode_3grp == 1 order spid opid gen id = _n save tomatch_grp2and3, replace //PSM based on Logistic regression 1 //Scheme 1 - Nearest neighbor within caliper .25*SD use tomatch_grp2and3, clear drop if p1==. gen logit1=log((1-p1)/p1) sum logit1 display .25*1.608234 set seed 13455 generate x=uniform() sort x gen group2_iniphone = 1 replace group2_iniphone = 0 if mode_3grp == 3 psmatch2 group2_iniphone , pscore(logit1) caliper(.4020585) noreplacement descending list spid opid inimode_web mode_3grp _pscore _treated _support _weight _id _n1 _nn _pdif in 1/5 sort _id g match=id[_n1] g treat=id if _nn==1 drop if treat== . sum treat keep treat match save scheme1_grp2and3, replace //scheme1 clear use scheme1_grp2and3, replace keep treat gen tx=1 rename treat id save t1,replace use scheme1_grp2and3, replace keep match gen tx=0 rename match id append using t1 sort id save t2, replace clear use tomatch_grp2and3, replace sort id merge 1:1 id using t2 tab tx drop if tx==. save post1_grp2and3, replace //bivarate analysis after matching - phone vs. phone cd "[Location]\Data\Data for paper" use post1_grp2and3, clear replace chd11dage = . if chd11dage < 0 ttest chd11dage if mode_3grp == 2 | mode_3grp == 3, by(mode_3grp) mean chd11dage ttest chd11dage if mode_3grp == 2 | mode_3grp == 3, by(mode_3grp) ttest chd11dnuminhh if mode_3grp == 2 | mode_3grp == 3, by(mode_3grp) stddiff chd11dage i.relgrp fl11spdied c11gender i.edu4cat i.chd11dmartstat childunder18 i.crl11dcgracehisp hasaddress hasphone hasemail r11dmetnonmet op11prsninhh chd11dnuminhh, by(mode_3grp) /************************************************************************/ *variable level analysis for experimantal data /************************************************************************/ cd "[Location]\Data" use "[Location]\NSOC_R11_Crss_OP_Tracker_File.dta", clear sort spid opid merge 1:1 spid opid using "[Location]\NSOC_R11_Crss_File.dta" keep if op11dnsoccomplete==1 gen w11cgfinwgt=0 replace w11cgfinwgt=w11cgfinwgt0 replace w11cgfinwgt=w11cglmlfinwgt0 if fl11spdied==1 drop _merge /*control for age and relationship*/ /*recode age */ gen agegrp=0 replace agegrp=1 if chd11dage>=18 & chd11dage<=40 replace agegrp=2 if chd11dage>40 & chd11dage<=60 replace agegrp=3 if chd11dage>60 & chd11dage<=80 replace agegrp=4 if chd11dage>80 replace agegrp=6 if agegrp==0 /*missing age*/ tab agegrp gen relgrp=0 replace relgrp=1 if c11relatnshp==2 /*spouse*/ replace relgrp=2 if c11relatnshp>=3 & c11relatnshp<=8 /*all children, including inlaw, step*/ replace relgrp=3 if (c11relatnshp>=19 & c11relatnshp<=22) | c11relatnshp==26 |c11relatnshp==28 | c11relatnshp==29 /*grandchild, nices nephews*/ replace relgrp=4 if c11relatnshp==25 | c11relatnshp==11 | c11relatnshp==10 | c11relatnshp==9 | c11relatnshp==15 | c11relatnshp==91 replace relgrp=5 if relgrp==0 tab relgrp /*all*/ svyset c11varunit [pweight=w11cgfinwgt], strata(c11varstrat) singleunit(scaled) svy: tab crl11dcgracehisp op11dmixmode, col svy: tab crl11dcgracehisp op11dmixmode, row log using df.log foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11busntype cec11dmswkunit cec11doccpdscr cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11occuptn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno crl11yourrace1 crl11yourrace2 crl11yourrace3 crl11yourrace4 crl11yourrace5 cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit chi11dgifttospcat crl11dcgracehisp cca11hwoftshp cdc11dhlpyrs chd11dmartstat chd11dnuminhh chi11dfromspcat chi11dpaycat chi11dgifttospcat chi11incimf crl11dcgracehisp { recode `var' (-6=-8) (-7=-8) svy: tab `var' op11dmixmode if `var'~=-1 & `var'~=-2 & `var'~=-4 & `var'~=-3 & `var'~=-5, col } log close *create variable level data in excel import excel "[Location]\Coding of question level data\var_level_data_based_on_technical reportv3.xlsx", sheet("Sheet1") firstrow clear save "[Location]\Data\Data for paper\var_level_data_based_on_tech_reportv3.dta", replace use "[Location]\Coding of question level data\grid summary_effect size_12.14.2022.dta", clear drop if Variable == "cca11hwoftshp" & Item == "Derived/imputed" keep Variable Chisquare Sample_size CramersVV df p_new ttest t_p_value Sample_size_con cohen_d cohendcat q_type_2 Type q_type_3 replace p_new = "" if p_new == "NA" | p_new== "na" destring p_new, replace replace Variable = "cca11hlpteeth" if Variable == "cca11hlpteeth " merge 1:1 Variable using "[Location]\Data\Data for paper\var_level_data_based_on_tech_reportv3.dta" list Variable _merge if _merge == 1 | _merge == 2 keep if _merge == 3 /*keep only variables included in the technical paper*/ drop _merge keep Variable Variablelabel Chisquare n df CramersVV p_new p_sig_cat Effectsize_Cat ttest t_p_value cohen_d Effectsize_con p_sig_con gen cat_var = 0 replace cat_var = 1 if p_sig_cat != . order Variable Variablelabel Type q_type_2 q_type_3 n Chisquare df p_new p_sig_cat CramersVV Effectsize_Cat ttest t_p_value p_sig_con cohen_d Effectsize_con cd "[Location]\Data\Data for paper" save var_level_sig_tests, replace *merge with coding of the objective and subjective questions, which as originally coded in excel and transfered in Stata "\coding_a_m.dta" use "[Location]\Coding of question level data\coding_a_m.dta", clear replace variable = "cca11hlpteeth" if variable == "cca11hlpteeth " save "[Location]\Coding of question level data\coding_a_m2.dta" cd "[Location]\Data\Data for paper" use var_level_sig_tests, clear rename Variable variable merge 1:1 variable using "[Location]\Coding of question level data\coding_a_m2.dta" keep if _merge == 1 | _merge == 3 *correct coding errors in coding_a_m2.dta gen objective = objective_m *code for those variables not being coded by M and A: replace objective = 1 if variable == "cdc11dhlpyrs" | variable == "chd11dage" | variable == "chd11dmartstat" | variable == "chd11dnuminhh" | variable == "chd11dspouage" | variable == "chi11dfromspcat" | variable == "chi11dgifttospcat" | variable == "chi11dpaycat" | variable == "chi11incim1" | variable == "chi11incim2" | variable == "chi11incim3" | variable == "chi11incim4" | variable == "chi11incim5" | variable == "chi11incimf" replace objective = 0 if variable == "cpp11hlpkptca" | variable == "cpp11hlpkptgo" | variable == "cpp11hlpkptgr" | variable == "cpp11hlpkptrl" | variable == "cpp11hlpkptvo" | variable == "cpp11hlpkptvs" | variable == "cpp11hlpkptwk" | variable == "cpp11impvst" | variable == "crl11primarace" | variable == "cse11frfamtlk" | variable == "crl11yourrace1" | variable == "crl11yourrace2" | variable == "crl11yourrace3" | variable == "crl11yourrace4" | variable == "crl11yourrace5" | variable == "crl11dcgracehisp" gen demo_exclude = 0 replace demo_exclude = 1 if variable == "chd11dage" | variable == "chd11dspouage" | variable == "crl11dcgracehisp" | variable == "crl11yourrace1" | variable == "crl11yourrace2" | variable == "crl11yourrace3" | variable == "crl11yourrace4" | variable == "crl11yourrace5" *correct df replace df = "" if df == "na" destring df, gen(newdf) recast int newdf, force replace newdf = 1 if variable == "cel11flexhrs" | variable == "cel11mswkhelp" replace newdf = 2 if variable == "cec11wrkmulti" | variable == "cpp11hlpkptca" | variable == "cpp11volun" | variable == "chc11regdoc" replace newdf = 3 if variable == "cca11hlpordmd" | variable == "cel11misswork" | variable == "cll11training" replace newdf = 4 if variable == "cec11look4wrk" | variable == "che11brethlmt" replace newdf = 5 if variable == "cac11diffemlv" | variable == "cac11vaccinedif" replace newdf = 6 if variable == "ch11incomeest" encode Type, gen(newtype) encode q_type_2, gen(newq_type_2) encode q_type_3, gen(newq_type_3) gen type_num = . replace type_num = 1 if q_type_3 == "Dummy" replace type_num = 2 if q_type_3 == "Ordinal" replace type_num = 3 if q_type_3 == "Nominal" replace type_num = 4 if q_type_3 == "Open text" | q_type_3 == "Open-ended Year 2cat" | q_type_3 == "choose all" label define type_num 1 "1 Dummy" 2 "2 Ordinal" 3 "3 Nominal" 4 "4 Other" label values type_num type_num gen n_centered = (n-1388) gen num_option = newdf + 1 save var_level_sig_tests2, replace *******create the table****** *the BH adjustment of significance was done in excel cd "[Location]\Data\Data for paper" use var_level_sig_tests2, clear drop _merge merge 1:1 variable using table4count_n_w_chi2 // this data is created in 0.1 code drop _merge merge 1:1 variable using cnt_options drop _merge gen type_num3cat = type_num replace type_num3cat = 3 if type_num == 4 label define type_num3cat 1 "1 Binary" 2 "2 Ordinal" 3 "3 Nominal & Others" label values type_num3cat type_num3cat replace type_num3cat = 1 if variable == "cpp11hlpkptca" | variable == "cec11wrkmulti" | variable == "cpp11volun" replace type_num3cat = 3 if variable == "cca11hlpordmd" | variable == "cel11misswork" | variable == "cll11training" replace countofoptions = 3 if variable == "cca11hlpordmd" gen num_options_cat = countofoptions replace num_options_cat = 5 if countofoptions > 4 drop if demo_exclude == 1 | cat_var == 0 | newtype == 9 mean table4_esize if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 tabulate objective if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, summarize(table4_esize) anova table4_esize objective if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 tabulate type_num3cat if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, summarize(table4_esize) anova table4_esize type_num3cat if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 tabulate num_options_cat if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, summarize(table4_esize) anova table4_esize num_options_cat if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 mean table4_count_n if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 tab table4_sig_adjusted_designgrp if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 tab objective table4_sig_adjusted_designgrp if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, chi2 row tab type_num3cat table4_sig_adjusted_designgrp if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, chi2 row tab num_options_cat table4_sig_adjusted_designgrp if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, chi2 row logit table4_sig_adjusted_designgrp ib1.objective ib1.type_num3cat countofoptions table4_count_n /* ----------------------------------------------------------------------------------------------- table4_sig_adjusted_designgrp | Coefficient Std. err. z P>|z| [95% conf. interval] ------------------------------+---------------------------------------------------------------- 0.objective | 1.188469 1.141253 1.04 0.298 -1.048346 3.425284 | type_num3cat | 2 Ordinal | 1.060117 1.539546 0.69 0.491 -1.957338 4.077573 3 Nominal & Others | 2.032326 1.486323 1.37 0.172 -.8808136 4.945465 | countofoptions | -.2693375 .4950409 -0.54 0.586 -1.2396 .7009249 table4_count_n | .0003103 .0007075 0.44 0.661 -.0010763 .001697 _cons | -4.230681 1.712197 -2.47 0.013 -7.586526 -.874837 ----------------------------------------------------------------------------------------------- */ betareg table4_esize ib1.type_num3cat ib1.objective countofoptions table4_count_n if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 /* ------------------------------------------------------------------------------------- table4_esize | Coefficient Std. err. z P>|z| [95% conf. interval] --------------------+---------------------------------------------------------------- table4_esize | type_num3cat | 2 Ordinal | .5188256 .1401256 3.70 0.000 .2441845 .7934667 3 Nominal & Others | .2224315 .1652769 1.35 0.178 -.1015052 .5463682 | 0.objective | .087016 .1191863 0.73 0.465 -.1465848 .3206168 countofoptions | .0891391 .0260622 3.42 0.001 .0380581 .1402202 table4_count_n | -.0005336 .0000697 -7.66 0.000 -.0006702 -.000397 _cons | -2.587952 .1259541 -20.55 0.000 -2.834817 -2.341086 --------------------+---------------------------------------------------------------- scale | _cons | 3.449443 .0995036 34.67 0.000 3.254419 3.644466 ------------------------------------------------------------------------------------- */ margins, dydx(type_num3cat) margins, dydx(objective) margins, dydx(countofoptions) margins, dydx(table4_count_n) /************************************************************************/ *variable level analysis for matched data /************************************************************************/ *variable level analysis cd "[Location]\Data\Data for paper" use post1, clear phi cpp11hlpkptgr inimode_web if cpp11hlpkptgr >0 /* Pearson chi2(1) = 1.9307 Pr = 0.165 phi = Cohen's w = fourfold point correlation = 0.0469 */ foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11busntype cec11dmswkunit cec11doccpdscr cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11occuptn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno crl11yourrace1 crl11yourrace2 crl11yourrace3 crl11yourrace4 crl11yourrace5 cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit chi11dgifttospcat crl11dcgracehisp cca11hwoftshp cdc11dhlpyrs chd11dmartstat chd11dnuminhh chi11dfromspcat chi11dpaycat chi11dgifttospcat chi11incimf crl11dcgracehisp { recode `var' (-6=-8) (-7=-8) tab `var' inimode_web if `var'~=-1 & `var'~=-2 & `var'~=-4 & `var'~=-3 & `var'~=-5 & p1 != . & `var'~=-8, chi2 col matrix test = (r(p)) matrix rownames test = `var' matrix colnames test = "p-value" mat2txt, matrix(test) saving(chisqtest_without_missing) append } import excel "[Location]\Data\Data for paper\postpsm_p_value_no_missing.xlsx", sheet("Sheet1") firstrow clear save postpsm_p_value_no_missing, replace **************************************************** clear input count_n variable 999 1 end tostring variable, replace save pscore_count_n, replace use post1, clear //*sample size in the model foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11busntype cec11dmswkunit cec11doccpdscr cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11occuptn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno crl11yourrace1 crl11yourrace2 crl11yourrace3 crl11yourrace4 crl11yourrace5 cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit chi11dgifttospcat crl11dcgracehisp cca11hwoftshp cdc11dhlpyrs chd11dmartstat chd11dnuminhh chi11dfromspcat chi11dpaycat chi11dgifttospcat chi11incimf crl11dcgracehisp { use post1, clear recode `var' (-6=-8) (-7=-8) replace `var' = . if `var'==-1 | `var'==-2 | `var'==-4 |`var'==-3 | `var'==-5 | `var'==-8 collapse (count) count_n=`var' gen variable = "`var'" save temp1, replace use pscore_count_n, clear append using temp1 save pscore_count_n, replace } use pscore_count_n, clear drop if variable == "1" save pscore_count_n2, replace use pscore_count_n2, clear duplicates tag variable, generate(duple) sort variable quietly by variable: gen dup = cond(_N==1,0,_n) drop if dup>1 drop dup duple save pscore_count_n3, replace log using "[Location]\Data\Data for paper\effect_size_without_missing.log", replace use post1, clear foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11busntype cec11dmswkunit cec11doccpdscr cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11occuptn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno crl11yourrace1 crl11yourrace2 crl11yourrace3 crl11yourrace4 crl11yourrace5 cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit chi11dgifttospcat crl11dcgracehisp cca11hwoftshp cdc11dhlpyrs chd11dmartstat chd11dnuminhh chi11dfromspcat chi11dpaycat chi11dgifttospcat chi11incimf crl11dcgracehisp { replace `var' = . if `var'==-1 | `var'==-2 | `var'==-4 |`var'==-3 | `var'==-5 | `var'==-6 | `var'==-7 | `var'==-8 phi `var' inimode_web } log close // save effect size data import excel "[Location]\Data\Data for paper\variable list for Mason_v2_first2filled.xlsx", sheet("without missing") firstrow clear cd "[Location]\Data\Data for paper" save post_phi_without_missing, replace cd "[Location]\Data\Data for paper" use var_level_sig_tests2, clear drop _merge merge 1:1 variable using postpsm_p_value_no_missing drop _merge keep if cat_var == 1 & demo_exclude == 0 count keep variable p_without_missing post_psm_sig_no_misisng export excel using "[Location]\Data\Data for paper\post_psm_Benjamini_Hochberg.xlsx" *note: BH adjustments are done using excel; three variables were removed due to countofoptions > 15: cec11busntype cec11doccpdscr cec11occuptn import excel using "[Location]\Data\Data for paper\post_psm_Benjamini_Hochberg2.xlsx" , sheet("Sheet1") firstrow clear save postpsm_p_BHadjusted, replace cd "[Location]\Data\Data for paper" use var_level_sig_tests2, clear drop _merge merge 1:1 variable using postpsm_p_BHadjusted drop _merge merge 1:1 variable using pscore_count_n3 drop _merge merge 1:1 variable using post_phi_without_missing drop _merge merge 1:1 variable using cnt_options drop _merge list variable Variablelabel objective type_num countofoptions count_n if cat_var == 1 & countofoptions > 6 *during checking, it is found that in the excel file, the phi for cpp11hlpkptgr is accidentally copied as the chisquare, make corrections below replace phi_binary = 0.0469 if phi_binary == 1.9307 gen type_num3cat = type_num replace type_num3cat = 3 if type_num == 4 label define type_num3cat 1 "1 Dummy" 2 "2 Ordinal" 3 "3 Nominal & Others" label values type_num3cat type_num3cat *correct coding error below replace type_num3cat = 1 if variable == "cpp11hlpkptca" | variable == "cec11wrkmulti" | variable == "cpp11volun" replace type_num3cat = 3 if variable == "cca11hlpordmd" | variable == "cel11misswork" | variable == "cll11training" replace countofoptions = 3 if variable == "cca11hlpordmd" gen num_options_cat = countofoptions replace num_options_cat = 5 if countofoptions > 4 // BH adjusted: post1_in_sig_adjusted drop if demo_exclude == 1 | cat_var == 0 | newtype == 9 gen catesize = phi_binary replace catesize = Cohens_w if phi_binary == . & Cohens_w != . mean catesize tabulate objective if cat_var == 1 , summarize(catesize) anova catesize objective tabulate type_num3cat if cat_var == 1 , summarize(catesize) anova catesize type_num3cat tabulate num_options_cat if cat_var == 1 , summarize(catesize) anova catesize num_options_cat * (Estimates) *---------------------------------- *post_psm_sig_no_misisng is before adjustment, we used the BH adjusted significance: sig_after_adjustment tab objective if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 tab sig_after_adjustment if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 tabulate objective if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, summarize(sig_after_adjustment) anova sig_after_adjustment objective tabulate type_num3cat if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, summarize(sig_after_adjustment) anova sig_after_adjustment type_num3cat tabulate num_options_cat if cat_var == 1 & demo_exclude == 0 & countofoptions < 15, summarize(sig_after_adjustment) anova sig_after_adjustment num_options_cat *------------------------------- mean countofoptions if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 mean count_n if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 gen type_2cat = type_num3cat replace type_2cat = 2 if type_num3cat == 3 gen interaction = 0 if type_2cat != . & objective != . replace interaction =1 if type_2cat == 1 & objective == 0 tab type_num3cat objective if cat_var == 1 & demo_exclude == 0 //p-value adjusted using BH method logit sig_after_adjustment ib1.objective ib1.type_num3cat countofoptions count_n if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 /* -------------------------------------------------------------------------------------- sig_after_adjustment | Coefficient Std. err. z P>|z| [95% conf. interval] ---------------------+---------------------------------------------------------------- 0.objective | .1703912 .6361014 0.27 0.789 -1.076345 1.417127 | type_num3cat | 2 Ordinal | 2.160663 .872166 2.48 0.013 .4512492 3.870077 3 Nominal & Others | 2.750239 .871221 3.16 0.002 1.042677 4.457801 | countofoptions | -.0986515 .1801142 -0.55 0.584 -.4516689 .2543658 count_n | .0014917 .0007891 1.89 0.059 -.0000549 .0030383 _cons | -4.524062 .9586741 -4.72 0.000 -6.403029 -2.645095 -------------------------------------------------------------------------------------- */ mean count_n betareg catesize ib1.type_num3cat ib1.objective countofoptions sample_size if cat_var == 1 & demo_exclude == 0 & countofoptions < 15 /* ------------------------------------------------------------------------------------- catesize | Coefficient Std. err. z P>|z| [95% conf. interval] --------------------+---------------------------------------------------------------- catesize | type_num3cat | 2 Ordinal | .914355 .1280602 7.14 0.000 .6633617 1.165348 3 Nominal & Others | .7499406 .1477555 5.08 0.000 .4603451 1.039536 | 0.objective | .0972362 .1072323 0.91 0.365 -.1129354 .3074077 countofoptions | .0510011 .0244624 2.08 0.037 .0030558 .0989465 sample_size | -.0006902 .0001129 -6.11 0.000 -.0009115 -.0004688 _cons | -2.642023 .1180185 -22.39 0.000 -2.873336 -2.410711 --------------------+---------------------------------------------------------------- scale | _cons | 3.484509 .0978863 35.60 0.000 3.292655 3.676362 ------------------------------------------------------------------------------------- */ margins, dydx(type_num3cat) margins, dydx(objective) margins, dydx(countofoptions) margins, dydx(sample_size) /******************************************************************************************/ *Illustration of % of differences betweeen mixed-mode and phone by varying % of web intake /******************************************************************************************/ cd "[Location]\Data\Data for paper" use post1, clear drop _merge gen group4cat = . replace group4cat = 1 if mode_3grp == 3 replace group4cat = 2 if mode_3grp == 1 label define group4cat 1 "1 telephone-only matched to web respondents in mixed mode" 2 "2 web respondents in mixed-mode design" label values group4cat group4cat save post1_2_merge_post23, replace use post1_grp2and3, clear drop _merge gen group4cat = . replace group4cat = 3 if mode_3grp == 3 replace group4cat = 4 if mode_3grp == 2 merge 1:1 spid opid group4cat using post1_2_merge_post23 drop _merge label define group4cat2 /// 1 "1 telephone-only matched to web respondents in mixed mode" /// 2 "2 web respondents in mixed-mode design" /// 3 "3 telephone-only matched to phone respondents in mixed mode" /// 4 "4 phone respondents in mixed-mode design" label values group4cat group4cat2 tab1 group4cat tab mode_3grp op11dmixmode *w11cgfinwgt is final weight including LML weights egen sum_weights = total( w11cgfinwgt), by(group4cat) tab sum_weights group /* sum_weight | group4cat s | 1 telepho 2 web res 3 telepho 4 phone r | Total -----------+--------------------------------------------+---------- 4571788 | 0 0 0 418 | 418 4618325 | 0 0 418 0 | 418 5629064 | 528 0 0 0 | 528 6023749 | 0 528 0 0 | 528 -----------+--------------------------------------------+---------- Total | 528 528 418 418 | 1,892 */ *sum of 1 & 3 weights is 10247389 (telephone only) *sum of 2 & 4 weights is 10595537 (mixed-mode) *sum of group 1: 5629064 *sum of group 2: 6023749 *sum of group 3: 4618325 *sum of group 4: 4571788 /* We could multiply by adjustment factors weights for groups 1, 2 so they are 10% and 3,4 so they are 90%: group 1 is 10% of 1+3; group 2 is 10% of 2+4 Groups 1 & 3 Groups 2 & 4 10247389 0 0 10595537 0 0 10247389 0.1 1024738.9 10595537 0.1 1059553.7 10247389 0.2 2049477.8 10595537 0.2 2119107.4 10247389 0.3 3074216.7 10595537 0.3 3178661.1 10247389 0.4 4098955.6 10595537 0.4 4238214.8 10247389 0.5 5123694.5 10595537 0.5 5297768.5 10247389 0.6 6148433.4 10595537 0.6 6357322.2 10247389 0.7 7173172.3 10595537 0.7 7416875.9 10247389 0.8 8197911.2 10595537 0.8 8476429.6 10247389 0.9 9222650.1 10595537 0.9 9535983.3 10247389 1 10247389 10595537 1 10595537 */ gen wgt_10per_web = w11cgfinwgt replace wgt_10per_web = w11cgfinwgt * (1024738.9 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_10per_web = w11cgfinwgt * (1059553.7 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_10per_web = w11cgfinwgt * (9222650.1 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_10per_web = w11cgfinwgt * (9535983.3 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_10per = total( wgt_10per_web), by(group4cat) tab sum_weights_10per group4cat gen wgt_20per_web = w11cgfinwgt replace wgt_20per_web = w11cgfinwgt * (2049477.8 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_20per_web = w11cgfinwgt * (2119107.4 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_20per_web = w11cgfinwgt * (8197911.2 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_20per_web = w11cgfinwgt * ( 8476429.6 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_20per = total( wgt_20per_web), by(group4cat) tab sum_weights_20per group4cat gen wgt_30per_web = w11cgfinwgt replace wgt_30per_web = w11cgfinwgt * (3074216.7 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_30per_web = w11cgfinwgt * (3178661.1 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_30per_web = w11cgfinwgt * (7173172.3 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_30per_web = w11cgfinwgt * (7416875.9 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_30per = total( wgt_30per_web), by(group4cat) tab sum_weights_30per group4cat gen wgt_40per_web = w11cgfinwgt replace wgt_40per_web = w11cgfinwgt * (4098955.6 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_40per_web = w11cgfinwgt * (4238214.8 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_40per_web = w11cgfinwgt * (6148433.4 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_40per_web = w11cgfinwgt * (6357322.2 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_40per = total( wgt_40per_web), by(group4cat) tab sum_weights_40per group4cat gen wgt_50per_web = w11cgfinwgt replace wgt_50per_web = w11cgfinwgt * (5123694.5 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_50per_web = w11cgfinwgt * (5297768.5 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_50per_web = w11cgfinwgt * (5123694.5 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_50per_web = w11cgfinwgt * (5297768.5 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_50per = total( wgt_50per_web), by(group4cat) tab sum_weights_50per group4cat gen wgt_60per_web = w11cgfinwgt replace wgt_60per_web = w11cgfinwgt * (6148433.4 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_60per_web = w11cgfinwgt * (6357322.2 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_60per_web = w11cgfinwgt * (4098955.6 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_60per_web = w11cgfinwgt * (4238214.8 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_60per = total( wgt_60per_web), by(group4cat) tab sum_weights_60per group4cat gen wgt_70per_web = w11cgfinwgt replace wgt_70per_web = w11cgfinwgt * (7173172.3 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_70per_web = w11cgfinwgt * (7416875.9 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_70per_web = w11cgfinwgt * (3074216.7 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_70per_web = w11cgfinwgt * (3178661.1 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_70per = total( wgt_70per_web), by(group4cat) tab sum_weights_70per group4cat gen wgt_80per_web = w11cgfinwgt replace wgt_80per_web = w11cgfinwgt * (8197911.2 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_80per_web = w11cgfinwgt * (8476429.6 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_80per_web = w11cgfinwgt * (2049477.8 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_80per_web = w11cgfinwgt * (2119107.4 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_80per = total( wgt_80per_web), by(group4cat) tab sum_weights_80per group4cat gen wgt_90per_web = w11cgfinwgt replace wgt_90per_web = w11cgfinwgt * (9222650.1 / 5629064) if group4cat == 1 // x% of 1+3 replace wgt_90per_web = w11cgfinwgt * (9535983.3 / 6023749) if group4cat == 2 // x% of 2+4 replace wgt_90per_web = w11cgfinwgt * (1024738.9 / 4618325) if group4cat == 3 // 1-x% of 1+3 replace wgt_90per_web = w11cgfinwgt * (1059553.7 / 4571788) if group4cat == 4 // 1-x% of 2+4 egen sum_weights_90per = total( wgt_90per_web), by(group4cat) tab sum_weights_90per group4cat foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { replace `var' = . if `var' < 0 } replace cll11informed = . if cll11informed == 5 replace cll11respect = . if cll11respect == 5 save data_4group_4_percent_web_wgts, replace //get count for sample size excluding missing clear input table5_count_n variable 999 1 end tostring variable, replace save table5_count_n_06302023, replace use data_4group_4_percent_web_wgts, clear foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { use data_4group_4_percent_web_wgts, clear collapse (count) table5_count_n=`var' gen variable = "`var'" save temp1, replace use table5_count_n_06302023, clear append using temp1 save table5_count_n_06302023, replace } drop if variable == "1" save table5_count_n_per_web_06302023, replace * 10% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_10percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_10per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close *note: all BH adjustment were completed using excel. import excel "[Location]\Data\Data for paper\web_10percent_sig_BH_code6.xlsx", sheet("web10") firstrow clear keep variable chi2_web10 sigweb10 sig_web10per_adjusted save sig_web10per_adjusted_code6, replace use sig_web10per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_10per_web = sqrt(chi2_web10/table5_count_n) mean esize_10per_web tab sig_web10per_adjusted save web10per_final_code6, replace * 20% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_20percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_20per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_20percent_sig_BH_code6.xlsx", sheet("web20") firstrow clear keep variable chi2_web20 sigweb20 sig_web20per_adjusted save sig_web20per_adjusted_code6, replace use sig_web20per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_20per_web = sqrt(chi2_web20/table5_count_n) mean esize_20per_web tab sig_web20per_adjusted save web20per_final_code6, replace * 30% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_30percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_30per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_30percent_sig_BH_code6.xlsx", sheet("web30") firstrow clear keep variable chi2_web30 sigweb30 sig_web30per_adjusted save sig_web30per_adjusted_code6, replace use sig_web30per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_30per_web = sqrt(chi2_web30/table5_count_n) mean esize_30per_web tab sig_web30per_adjusted save web30per_final_code6, replace * 40% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_40percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_40per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_40percent_sig_BH_code6.xlsx", sheet("web40") firstrow clear keep variable chi2_web40 sigweb40 sig_web40per_adjusted save sig_web40per_adjusted_code6, replace use sig_web40per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_40per_web = sqrt(chi2_web40/table5_count_n) mean esize_40per_web tab sig_web40per_adjusted save web40per_final_code6, replace * 50% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_50percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_50per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_50percent_sig_BH_code6.xlsx", sheet("web50") firstrow clear keep variable chi2_web50 sigweb50 sig_web50per_adjusted save sig_web50per_adjusted_code6, replace use sig_web50per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_50per_web = sqrt(chi2_web50/table5_count_n) mean esize_50per_web tab sig_web50per_adjusted save web50per_final_code6, replace * 60% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_60percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_60per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_60percent_sig_BH_code6.xlsx", sheet("web60") firstrow clear keep variable chi2_web60 sigweb60 sig_web60per_adjusted save sig_web60per_adjusted_code6, replace use sig_web60per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_60per_web = sqrt(chi2_web60/table5_count_n) mean esize_60per_web tab sig_web60per_adjusted save web60per_final_code6, replace * 70% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_70percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_70per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_70percent_sig_BH_code6.xlsx", sheet("web70") firstrow clear keep variable chi2_web70 sigweb70 sig_web70per_adjusted save sig_web70per_adjusted_code6, replace use sig_web70per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_70per_web = sqrt(chi2_web70/table5_count_n) mean esize_70per_web tab sig_web70per_adjusted save web70per_final_code6, replace * 80% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_80percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_80per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_80percent_sig_BH_code6.xlsx", sheet("web80") firstrow clear keep variable chi2_web80 sigweb80 sig_web80per_adjusted save sig_web80per_adjusted_code6, replace use sig_web80per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_80per_web = sqrt(chi2_web80/table5_count_n) mean esize_80per_web tab sig_web80per_adjusted save web80per_final_code6, replace * 90% web *-------------------------------------- log using "[Location]\Data\Data for paper\web_90percent_chi_code6.log", replace use data_4group_4_percent_web_wgts, clear svyset c11varunit [pweight=wgt_90per_web], strata(c11varstrat) singleunit(scaled) foreach var of varlist cac11arguelv cac11careamt cac11closr2sp cac11coorddif cac11dealbetr cac11decsdif cac11diffemlv cac11diffemo cac11diffinc cac11diffinlv cac11diffphlv cac11diffphy cac11exhaustd cac11fmlydisa cac11hospdif cac11insurdif cac11joylevel cac11medsdif cac11moreconf cac11moresat cac11nerveslv cac11notime cac11spapprlv cac11toomuch cac11uroutchg cac11vaccinedif cca11cmpbnkng cca11cmpgrcry cca11cmpordrx cca11hlpbnkng cca11hlpdiet cca11hlpexrcs cca11hlpfeet cca11hlpmdtk cca11hlpmed cca11hlpordmd cca11hlpshot cca11hlpskin cca11hlpteeth cca11hwoftchs cca11hwoftdrv cca11hwofthom cca11hwoftott cca11hwoftpc cca11hwoftshp cca11training ccv11covid ccv11effects ccv11limitsp ccv11nohelp ccv11quarantine ccv11vaccine ccv11vaccinemth ccv11vaccinetry ccv11vaccineyr cdc11hlphrmvf cdc11hlphrpce cdc11hlpmthst cdc11hlpsched cdc11hlpyrpls cdi11hwget2sp cdi11tm2spunt cec11dmswkunit cec11dworktype cec11flexhrs cec11hlpafwk1 cec11look4wrk cec11misswork cec11msswrkmt cec11mswkfmly cec11mswkhelp cec11mswkpd cec11mswkpers cec11mswksick cec11mswkvctn cec11ownbusns cec11unpdwrk cec11wrk4pay cec11wrkmulti cec11wrksched cel11flexhrs cel11misswork cel11mswkhelp cel11mswkpd cel11unpdwrk cel11wrk4pay cel11wrksched ch11incomeest chc11coordnat chc11hlpdramt chc11hlpinsrn chc11hlpmanage chc11hlpmdapt chc11hlpothin chc11hlpspkdr chc11hosphlp chc11hospset chc11hosptrain chc11hvaccappt chc11listen chc11online chc11regdoc chc11talkoft chc11understood chd11cgbrthmt chd11chldlvng chd11dmartstat chd11marendmt chd11martstat chd11spbrthmt che11armlmt che11armstr che11arthrits che11brethlmt che11brethprb che11cancer che11cncrbrst che11cncroth che11cncrpros che11cncrskin che11diabetes che11easychng che11enrgylmt che11feelconf che11fltdown che11fltltlin che11fltnervs che11fltworry che11gaveup che11health che11hearing che11highbld che11hrtattck che11leglmt che11legstr che11lifegood che11lifelone che11lifemng che11lost10lb che11lowenrgy che11lungdis che11moodbord che11moodcher che11moodfull che11moodlone che11moodpcfl che11moodupst che11osteoprs che11othheart che11pain che11painlmt che11rcvrqick che11seeing che11sleepint che11sleeptrb che11try2lose chi11certdpst chi11chckacct chi11dfromspcat chi11dgifttospcat chi11dpaycat chi11fromsp chi11gifttosp chi11incimf chi11medicaid chi11medicare chi11medigap chi11mutlfnds chi11ownhome chi11payasdvc chi11payhmhlp chi11payinsur chi11paymbdvc chi11paymeds chi11paysfty chi11privinsr chi11retrplan chi11svgsacct chi11tricare chi11uninsurd cll11brehlp cll11caredecis cll11carenowan cll11informed cll11makedecis cll11painhlp cll11respect cll11sadhlp cll11talkmed cll11training cpp11careothr cpp11go4fun cpp11gorelgs cpp11hlpkptca cpp11hlpkptgo cpp11hlpkptgr cpp11hlpkptrl cpp11hlpkptvo cpp11hlpkptvs cpp11hlpkptwk cpp11impgo cpp11impgroup cpp11imprel cpp11impvst cpp11prtgroup cpp11volun cpp11vstfrfam cpp11wrk4pay crl11hisplatno cse11fndfnhlp cse11frfamact cse11frfamhlp cse11frfamtlk cse11hlpgtdev cse11hlppdhlp cse11homsafer cse11spptgrp cse11srv2hlp cse11training cvs11email cvs11phone cvs11video cvs11visit { display "`var'" svy: tab `var' op11dmixmode, notab } log close import excel "[Location]\Data\Data for paper\web_90percent_sig_BH_code6.xlsx", sheet("web90") firstrow clear keep variable chi2_web90 sigweb90 sig_web90per_adjusted save sig_web90per_adjusted_code6, replace use sig_web90per_adjusted_code6, clear merge 1:1 variable using table5_count_n_per_web_06302023 drop _merge gen esize_90per_web = sqrt(chi2_web90/table5_count_n) mean esize_90per_web tab sig_web90per_adjusted save web90per_final_code6, replace ********** merge all data together ********** use web10per_final_code6, clear merge 1:1 variable using web20per_final_code6 drop _merge merge 1:1 variable using web30per_final_code6 drop _merge merge 1:1 variable using web40per_final_code6 drop _merge merge 1:1 variable using web50per_final_code6 drop _merge merge 1:1 variable using web60per_final_code6 drop _merge merge 1:1 variable using web70per_final_code6 drop _merge merge 1:1 variable using web80per_final_code6 drop _merge merge 1:1 variable using web90per_final_code6 drop _merge mean sig_web10per_adjusted sig_web20per_adjusted sig_web30per_adjusted sig_web40per_adjusted sig_web50per_adjusted sig_web60per_adjusted sig_web70per_adjusted sig_web80per_adjusted sig_web90per_adjusted /* ----------------------------------------------------------------------- | Mean Std. err. [95% conf. interval] ----------------------+------------------------------------------------ sig_web10per_adjusted | 0 (omitted) sig_web20per_adjusted | 0 (omitted) sig_web30per_adjusted | 0 (omitted) sig_web40per_adjusted | .0089686 .0063275 -.003501 .0214382 sig_web50per_adjusted | .0134529 .007732 -.0017846 .0286904 sig_web60per_adjusted | .0179372 .0089078 .0003825 .0354919 sig_web70per_adjusted | .0313901 .0117029 .0083271 .0544532 sig_web80per_adjusted | .0493274 .0145339 .0206852 .0779695 sig_web90per_adjusted | .0493274 .0145339 .0206852 .0779695 ----------------------------------------------------------------------- */ mean esize_10per_web esize_20per_web esize_30per_web esize_40per_web esize_50per_web esize_60per_web esize_70per_web esize_80per_web esize_90per_web