***Analyze Linkage indicator and paradata *** use "Z:\biases-PASS\data\call_record_and_linkage_data_cross.dta", clear * overall mode switches: dy_md_sw_tot * overall interview switches: dy_int_sw_tot * overall contact attempts: nr_contacts_tot * ever refused: ever_refused * AAPOR groups (final): FIN_GRP * Mode of last contact attempt: FINMET * Household interview status: REALSTAT * Linkage categories: cat_linkage gen calls_1_2 = 0 replace calls_1_2 = 1 if inlist(nr_contacts_tot,1,2) gen calls_3_5 = 0 replace calls_3_5 = 1 if inlist(nr_contacts_tot,3,4,5) gen calls_6_10 = 0 replace calls_6_10 = 1 if inlist(nr_contacts_tot,6,7,8,9,10) gen calls_11_more = 0 replace calls_11_more = 1 if nr_contacts_tot >=11 & nr_contacts_tot <= 500 gen calls_grp = 1 if calls_1_2==1 replace calls_grp = 2 if calls_3_5 == 1 replace calls_grp = 3 if calls_6_10 == 1 replace calls_grp = 4 if calls_11_more == 1 gen cati_ct = 0 replace cati_ct = 1 if inlist(FINMET,1) gen linked_conserv = 0 replace linked_conserv = 1 if inrange(new_quality_fin,6,17) gen linked_liberal = 0 replace linked_liberal = 1 if inrange(new_quality_fin,1,17) *Begin bivariate analysis (liberal) svyset psu, strata(strpsu) singleunit(scaled) svy: mean linked_liberal, over(calls_grp) svy: logit linked_liberal i.calls_grp tab calls_grp linked_liberal svy: mean linked_liberal, over(dy_md_sw_tot) svy: logit linked_liberal dy_md_sw_tot tab dy_md_sw_tot linked_liberal svy: mean linked_liberal, over(dy_int_sw_tot) svy: logit linked_liberal dy_int_sw_tot tab dy_int_sw_tot linked_liberal svy: mean linked_liberal, over(nr_contacts_tot) svy: logit linked_liberal nr_contacts_tot svy: mean linked_liberal, over(ever_refused) svy: logit linked_liberal ever_refused tab ever_refused linked_liberal svy: mean linked_liberal, over(cati_ct) svy: logit linked_liberal cati_ct tab cati_ct linked_liberal svy: mean linked_liberal, over(FIN_GRP) svy: logit linked_liberal ib7.FIN_GRP tab FIN_GRP linked_liberal svy: logit linked_liberal i.calls_grp ever_refused cati_ct ib7.FIN_GRP svy: logit linked_liberal ib7.FIN_GRP svy: mean linked_liberal **Bivariate analysis (conservative) svyset psu, strata(strpsu) singleunit(scaled) svy: mean linked_conserv, over(calls_grp) svy: logit linked_conserv i.calls_grp tab calls_grp linked_conserv svy: mean linked_conserv, over(dy_md_sw_tot) svy: logit linked_conserv dy_md_sw_tot tab dy_md_sw_tot linked_conserv svy: mean linked_conserv, over(dy_int_sw_tot) svy: logit linked_conserv dy_int_sw_tot tab dy_int_sw_tot linked_conserv svy: mean linked_conserv, over(nr_contacts_tot) svy: logit linked_conserv nr_contacts_tot svy: mean linked_conserv, over(ever_refused) svy: logit linked_conserv ever_refused tab ever_refused linked_conserv svy: mean linked_conserv, over(cati_ct) svy: logit linked_conserv cati_ct tab cati_ct linked_conserv svy: mean linked_conserv, over(FIN_GRP) svy: logit linked_conserv ib7.FIN_GRP tab FIN_GRP linked_conserv svy: logit linked_conserv i.calls_grp ever_refused cati_ct ib7.FIN_GRP svy: logit linked_conserv ib7.FIN_GRP svy: mean linked_conserv tab linked_conserv ****Household Interview Analysis use "Z:\biases-PASS\data\para_data_linkage_indicator_survey_data_household.dta", clear gen ub2_since09 = 0 if HA0300 == 2 replace ub2_since09 = 1 if HA0300 == 1 gen hhtype_cat = . replace hhtype_cat = 1 if type_HH == 1 replace hhtype_cat = 2 if type_HH == 2 replace hhtype_cat = 3 if inlist(type_HH,3,4,5,6) replace hhtype_cat = 4 if inlist(type_HH,-5,7,8) gen hhchild = 1 if inlist(HEK1700,-2,1,2) replace hhchild = 0 if HEK1700 == -3 gen hhsizecat = 1 if HA0100 == 1 replace hhsizecat = 2 if HA0100 == 2 replace hhsizecat = 3 if HA0100 == 3 replace hhsizecat = 4 if inrange(HA0100,4,50) gen HLS0100a_rc = 0 replace HLS0100a_rc = 1 if HLS0100a == 2 gen HLS0200a_rc = 0 replace HLS0200a_rc = 1 if HLS0200a == 2 gen HLS0300a_rc = 0 replace HLS0300a_rc = 1 if HLS0300a == 2 gen HLS0400a_rc = 0 replace HLS0400a_rc = 1 if HLS0400a == 2 gen HLS0500a_rc = 0 replace HLS0500a_rc = 1 if HLS0500a == 2 gen HLS0600a_rc = 0 replace HLS0600a_rc = 1 if HLS0600a == 2 gen HLS0700a_rc = 0 replace HLS0700a_rc = 1 if HLS0700a == 2 gen HLS0800a_rc = 0 replace HLS0800a_rc = 1 if HLS0800a == 2 gen HLS0900a_rc = 0 replace HLS0900a_rc = 1 if HLS0900a == 2 gen HLS1000a_rc = 0 replace HLS1000a_rc = 1 if HLS1000a == 2 gen HLS1100a_rc = 0 replace HLS1100a_rc = 1 if HLS1100a == 2 gen HLS1200a_rc = 0 replace HLS1200a_rc = 1 if HLS1200a == 2 gen HLS1300a_rc = 0 replace HLS1300a_rc = 1 if HLS1300a == 2 gen depriv_item = HLS0100a_rc + HLS0200a_rc + HLS0300a_rc + HLS0400a_rc + HLS0500a_rc + HLS0600a_rc + HLS0700a_rc + HLS0800a_rc + HLS0900a_rc + HLS1000a_rc + HLS1100a_rc + HLS1200a_rc + HLS1300a_rc gen depriv_item_cat = 1 if inlist(depriv_item,0) replace depriv_item_cat = 2 if inlist(depriv_item,1) replace depriv_item_cat = 3 if inlist(depriv_item,2) replace depriv_item_cat = 4 if inrange(depriv_item,3,50) gen HLS1400a_rc = 0 replace HLS1400a_rc = 1 if HLS1400a == 2 gen HLS1500a_rc = 0 replace HLS1500a_rc = 1 if HLS1500a == 2 gen HLS1600a_rc = 0 replace HLS1600a_rc = 1 if HLS1600a == 2 gen HLS1700a_rc = 0 replace HLS1700a_rc = 1 if HLS1700a == 2 gen HLS1800a_rc = 0 replace HLS1800a_rc = 1 if HLS1800a == 2 gen HLS1900a_rc = 0 replace HLS1900a_rc = 1 if HLS1900a == 2 gen HLS2000a_rc = 0 replace HLS2000a_rc = 1 if HLS2000a == 2 gen HLS2100a_rc = 0 replace HLS2100a_rc = 1 if HLS2100a == 2 gen HLS2200a_rc = 0 replace HLS2200a_rc = 1 if HLS2200a == 2 gen HLS2300a_rc = 0 replace HLS2300a_rc = 1 if HLS2300a == 2 gen HLS2400a_rc = 0 replace HLS2400a_rc = 1 if HLS2400a == 2 gen HLS2500a_rc = 0 replace HLS2500a_rc = 1 if HLS2500a == 2 gen HLS2600a_rc = 0 replace HLS2600a_rc = 1 if HLS2600a == 2 gen depriv_activ = HLS1400a_rc + HLS1500a_rc + HLS1600a_rc + HLS1700a_rc + HLS1800a_rc + HLS1900a_rc + HLS2000a_rc + HLS2100a_rc + HLS2200a_rc + HLS2300a_rc + HLS2400a_rc + HLS2500a_rc + HLS2600a_rc gen depriv_activ_cat = 1 if depriv_activ == 0 replace depriv_activ_cat = 2 if inlist(depriv_activ,1,2) replace depriv_activ_cat = 3 if inlist(depriv_activ,3,4) replace depriv_activ_cat = 4 if inrange(depriv_activ,5,50) gen neighborhoodcat=1 if HW2100 == 1 replace neighborhoodcat = 2 if HW2100 == 2 replace neighborhoodcat = 3 if inlist(HW2100,3,4,5) gen hhinc_cat = 1 if inrange(HEK0600,0,999) replace hhinc_cat = 2 if inrange(HEK0600,1000,1999) replace hhinc_cat = 3 if inrange(HEK0600,2000,2999) replace hhinc_cat = 4 if inrange(HEK0600,3000,10000000) replace hhinc_cat = 4 if HEK0600==-5 replace hhinc_cat = 5 if inlist(HEK0600,-2,-1) gen hhsave_cat = 1 if inlist(HEK1200,0,1) replace hhsave_cat = 2 if inlist(HEK1200,2,3,4) replace hhsave_cat = 3 if inlist(HEK1200,5) replace hhsave_cat = 4 if inlist(HEK1200,6,7) replace hhsave_cat = 5 if inlist(HEK1200,-2,-1) gen ownhh = 1 if HW0300 == 3 replace ownhh = 0 if inlist(HW0300,-3,1,2,4) *Begin bivariate analysis (liberal) svyset psu [pweight=wqmihh], strata(strpsu) singleunit(scaled) gen linked_conserv = 0 replace linked_conserv = 1 if inrange(new_quality_fin,6,17) gen linked_liberal = 0 replace linked_liberal = 1 if inrange(new_quality_fin,1,17) svy: mean linked_liberal, over(hhtype_cat) svy: logit linked_liberal hhtype_cat tab hhtype_cat linked_liberal svy: mean linked_liberal, over(hhchild) svy: logit linked_liberal hhchild tab hhchild linked_liberal svy: mean linked_liberal, over(hhsizecat) svy: logit linked_liberal i.hhsizecat tab hhsizecat linked_liberal svy: mean linked_liberal, over(ownhh) svy: logit linked_liberal i.ownhh tab ownhh linked_liberal svy: mean linked_liberal, over(depriv_item_cat) svy: logit linked_liberal i.depriv_item_cat tab depriv_item_cat linked_liberal svy: mean linked_liberal, over(depriv_activ_cat) svy: logit linked_liberal i.depriv_activ_cat tab depriv_activ_cat linked_liberal svy: mean linked_liberal, over(neighborhoodcat) svy: logit linked_liberal i.neighborhoodcat tab neighborhoodcat linked_liberal svy: mean linked_liberal, over(hhinc_cat) svy: logit linked_liberal i.hhinc_cat tab hhinc_cat linked_liberal svy: mean linked_liberal, over(hhsave_cat) svy: logit linked_liberal i.hhsave_cat tab hhsave_cat linked_liberal svy: mean linked_liberal, over(ub2_since09) svy: logit linked_liberal i.ub2_since09 tab ub2_since09 linked_liberal svy: mean linked_liberal svy: logit linked_liberal i.hhtype_cat i.ownhh i.depriv_item_cat i.hhinc_cat i.hhsave_cat ub2_since09 *Begin bivariate analysis (conservative) svyset psu [pweight=wqmihh], strata(strpsu) singleunit(scaled) svy: mean linked_conserv, over(hhtype_cat) svy: logit linked_conserv hhtype_cat tab hhtype_cat linked_conserv svy: mean linked_conserv, over(ownhh) svy: logit linked_conserv i.ownhh tab ownhh linked_conserv svy: mean linked_conserv, over(depriv_item_cat) svy: logit linked_conserv i.depriv_item_cat tab depriv_item_cat linked_conserv svy: mean linked_conserv, over(hhinc_cat) svy: logit linked_conserv i.hhinc_cat tab hhinc_cat linked_conserv svy: mean linked_conserv, over(hhsave_cat) svy: logit linked_conserv i.hhsave_cat tab hhsave_cat linked_conserv svy: mean linked_conserv, over(ub2_since09) svy: logit linked_conserv i.ub2_since09 tab ub2_since09 linked_conserv svy: mean linked_conserv tab linked_conserv svy: logit linked_conserv i.hhtype_cat i.ownhh i.depriv_item_cat i.hhinc_cat i.hhsave_cat ub2_since09 ****Person Interview Analysis use "Q:\data\para_data_linkage_indicator_survey_data_household_and person_spell.dta", clear gen agecat = 1 if inrange(palter,0,39) replace agecat = 2 if inrange(palter,40,54) replace agecat = 3 if inrange(palter,55,69) replace agecat = 4 if inrange(palter,70,150) gen female = 1 if zpsex == 2 replace female = 0 if zpsex == 1 gen nevmarried = 0 if inlist(PD0500,2,3,4,5,6,7,-3) replace nevmarried = 1 if inlist(PD0500,1) gen educcat = 2 if PB1300a == 1 | PB1300b == 1 | PB1300c == 1 | PB1300d == 1 | PB1300e == 1 replace educcat = 3 if PB1300f == 1 | PB1300g == 1 | PB1300h == 1 | PB1300i == 1 | PB1300j == 1 replace educcat = 1 if educcat == . gen linkconsent = 1 if RegP0100 == 1 replace linkconsent = 0 if RegP0100 == 2 gen empchange = PET4000 if PET4000 ~=-10 & PET4000~=-3 & PET4000 ~=-2 & PET4000 ~=-1 gen empchangecat = 1 if empchange == 0 replace empchangecat = 2 if inlist(empchange,1,2) replace empchangecat = 3 if inlist(empchange,3,4) replace empchangecat = 4 if inlist(empchange,-5) | inrange(empchange,5,100) gen emplength_mo = PET4100a if inrange(PET4100a,1,500) replace emplength_mo = 0 if PET4100a == -5 gen emplength_yr = PET4100b if inrange(PET4100b,1,500) replace emplength_yr = 0 if PET4100b == -5 gen emplength = emplength_mo + emplength_yr*12 gen emplengthcat = 1 if inrange(emplength,0,131) replace emplengthcat = 2 if inrange(emplength,132,251) replace emplengthcat = 3 if inrange(emplength,252,383) replace emplengthcat = 4 if inrange(emplength,384,1000) gen everregunemp = 1 if PA0435 == 1 replace everregunemp = 0 if PA0435 == 2 gen disability = 0 if PG0500 == 2 replace disability = 1 if inlist(PG0500,1,3) gen privinsure = 1 if inlist(PG1300,2,3,4,5,6) replace privinsure = 0 if PG1300 == 1 gen socpension = 1 if PEK1400 == 1 replace socpension = 0 if PEK1400==2 gen foreignborn = 1 if PMI0100 == 2 replace foreignborn = 0 if PMI0100 == 1 gen employed = 0 replace employed = 1 if empl_status==1 gen occup_cat = 4 if employed == 0 replace occup_cat = 1 if occup_status == 1 replace occup_cat = 2 if occup_status == 2 replace occup_cat = 3 if inlist(occup_status,3,4,5,6,7,8) *Begin bivariate analysis (liberal) svyset psu [pweight=wqmip], strata(strpsu) singleunit(scaled) gen linked_conserv = 0 replace linked_conserv = 1 if inrange(new_quality_fin,6,17) gen linked_liberal = 0 replace linked_liberal = 1 if inrange(new_quality_fin,1,17) svy: mean linked_liberal, over(agecat) svy: logit linked_liberal agecat tab agecat linked_liberal svy: mean linked_liberal, over(female) svy: logit linked_liberal female tab female linked_liberal svy: mean linked_liberal, over(nevmarried) svy: logit linked_liberal nevmarried tab nevmarried linked_liberal svy: mean linked_liberal, over(educcat) svy: logit linked_liberal educcat tab educcat linked_liberal svy: mean linked_liberal, over(linkconsent) svy: logit linked_liberal linkconsent tab linkconsent linked_liberal svy: mean linked_liberal, over(employed) svy: logit linked_liberal employed tab employed linked_liberal svy: mean linked_liberal, over(occup_cat) svy: logit linked_liberal occup_cat tab occup_cat linked_liberal svy: mean linked_liberal, over(empchangecat) svy: logit linked_liberal empchangecat tab empchangecat linked_liberal svy: mean linked_liberal, over(emplengthcat) svy: logit linked_liberal emplengthcat tab emplengthcat linked_liberal svy: mean linked_liberal, over(everregunemp) svy: logit linked_liberal everregunemp tab everregunemp linked_liberal svy: mean linked_liberal, over(disability) svy: logit linked_liberal disability tab disability linked_liberal svy: mean linked_liberal, over(privinsure) svy: logit linked_liberal privinsure tab privinsure linked_liberal svy: mean linked_liberal, over(socpension) svy: logit linked_liberal socpension tab socpension linked_liberal svy: mean linked_liberal, over(foreignborn) svy: logit linked_liberal foreignborn tab foreignborn linked_liberal svy: mean linked_liberal, over(occup_cat) svy: logit linked_liberal occup_cat tab occup_cat linked_liberal svy: logit linked_liberal i.agecat female nevmarried i.educcat foreignborn privinsure i.occup_cat *Begin bivariate analysis (conservative) svyset psu [pweight=wqmip], strata(strpsu) singleunit(scaled) gen linked_conserv = 0 replace linked_conserv = 1 if inrange(new_quality_fin,6,17) gen linked_liberal = 0 replace linked_liberal = 1 if inrange(new_quality_fin,1,17) svy: mean linked_conserv, over(agecat) svy: logit linked_conserv agecat tab agecat linked_conserv svy: mean linked_conserv, over(female) svy: logit linked_conserv female tab female linked_conserv svy: mean linked_conserv, over(nevmarried) svy: logit linked_conserv nevmarried tab nevmarried linked_conserv svy: mean linked_conserv, over(educcat) svy: logit linked_conserv educcat tab educcat linked_conserv svy: mean linked_conserv, over(linkconsent) svy: logit linked_conserv linkconsent tab linkconsent linked_conserv svy: mean linked_conserv, over(employed) svy: logit linked_conserv employed tab employed linked_conserv svy: mean linked_conserv, over(occup_cat) svy: logit linked_conserv occup_cat tab occup_cat linked_conserv svy: mean linked_conserv, over(empchangecat) svy: logit linked_conserv empchangecat tab empchangecat linked_conserv svy: mean linked_conserv, over(emplengthcat) svy: logit linked_conserv emplengthcat tab emplengthcat linked_conserv svy: mean linked_conserv, over(everregunemp) svy: logit linked_conserv everregunemp tab everregunemp linked_conserv svy: mean linked_conserv, over(disability) svy: logit linked_conserv disability tab disability linked_conserv svy: mean linked_conserv, over(privinsure) svy: logit linked_conserv privinsure tab privinsure linked_conserv svy: mean linked_conserv, over(socpension) svy: logit linked_conserv socpension tab socpension linked_conserv svy: mean linked_conserv, over(foreignborn) svy: logit linked_conserv foreignborn tab foreignborn linked_conserv svy: mean linked_conserv, over(occup_cat) svy: logit linked_conserv occup_cat tab occup_cat linked_conserv svy: logit linked_conserv i.agecat female nevmarried i.educcat foreignborn privinsure i.occup_cat