** ANALYSIS CODE FOR PAPER "DOES BENEFIT FRAMING IMPROVE RECORD LINKAGE CONSENT RATES? A SURVEY EXPERIMENT"

clear

***DATA PREPARATION STEP

import excel "N:\Stata_data.xlsx", sheet("Tabelle1") cellrange(A1:Q51) firstrow

destring lfd, replace

merge 1:1 lfd using "N:\experiment_atypik.dta"

drop if _merge == 1

gen monitor = 1 if _merge == 3
replace monitor = 0 if _merge== 2

drop _merge
merge 1:1 lfd using "N:\experiment_atypik_newvars.dta"

svyset int_id [pweight=hoch], fpc(fpc1) || betnr_atypik, strata(risiko) fpc(fpc2) singleunit(scaled) || lfd, fpc(fpc3)
  
gen treatment = 1 if inlist(f31b,1,2,.a,.b,.d)
replace treatment = 2 if inlist(f31a,1,2,.a,.b,.d)
replace treatment = 3 if inlist(f31c,1,2,.a,.b,.d)


gen consent = 1 if inlist(f31b,1) | inlist(f31a,1) | inlist(f31c,1)
replace consent = 0 if inlist(f31b,2,.a,.b,.d) | inlist(f31a,2,.a,.b,.d) | inlist(f31c,2,.a,.b,.d)

gen sbgrp = 1 if inlist(treatment,1,2)
gen sbgrp2 = 2 if inlist(treatment,1,3)

tab consent treatment, col chi
svy, subpop(sbgrp2): tabulate consent treatment, col se ci

*** Respondent characteristics***
gen age = 2013-gebjahr

gen age_cat = 1 if age >=15 & age <=32
replace age_cat = 2 if age >= 33 & age <=45
replace age_cat = 3 if age >= 46 & age <=53
replace age_cat = 4 if age >=54 & age <= 100

gen female = 1 if sex == 1
replace female = 0 if sex == 2

gen educ = 1 if inlist(f96,1,2,96)
replace educ = 2 if inlist(f96,3)
replace educ = 3 if inlist(f96,4,5)

gen mig = 0 if inlist(migration, 1,2,3)
replace mig = 1 if inlist(migration, 2,3)

gen nevermarried = 0 if inlist(mart_status,1,2,3,4)
replace nevermarried = 1 if inlist(mart_status,1)

gen emp = 1 if inlist(employment,1)
replace emp = 2 if inlist(employment,2)
replace emp = 3 if inlist(employment,3,4)

gen more6wobreak = 1 if inlist(no_break,1)
replace more6wobreak = 0 if inlist(no_break,0,2)

gen workhours = 1 if inlist(az_tg,1,2) /**ref */
replace workhours = 2 if inlist(az_tg,3,4)
replace workhours = 3 if inlist(az_tg,5,6)

gen numchild = 1 if kids_14 == 0
replace numchild = 2 if kids_14 == 1
replace numchild = 3 if inlist(kids_14,2,3,4,5)

gen rightswh = 0 if no_right == 0
replace rightswh = 1 if inlist(no_right,1,2,3)

gen desworktime = 1 if des_wt == 1
replace desworktime = 2 if des_wt == 3
replace desworktime = 3 if des_wt == 2


gen int_female = 1 if intsex==2
replace int_female = 0 if intsex==1

gen int_age = 1 if inlist(intaltgr,1)
replace int_age = 2 if inlist(intaltgr,2,3)
replace int_age = 3 if inlist(intaltgr,4,7)

gen int_educ = 1 if inlist(int_sab,1,6,9)
replace int_educ = 2 if inlist(int_sab,2,4)
replace int_educ = 3 if inlist(int_sab,3,5)

gen int_ten = 1 if inlist(int_dau,3,4)
replace int_ten = 0 if inlist(int_dau,1,2)


**ANALYSIS STEP

prop age_cat, over(treatment)
tab age_cat treatment, chi
svy: prop age_cat, over(treatment)
svy: tabulate age_cat treatment, col se

prop female, over(treatment)
tab female treatment, chi
svy: prop female, over(treatment)
svy: tabulate female treatment, col se

prop educ, over(treatment)
tab educ treatment, chi
svy: prop educ, over(treatment)
svy: tabulate educ treatment, col se

prop mig, over(treatment)
tab mig treatment, chi
svy: prop mig, over(treatment)
svy: tabulate mig treatment, col se

prop nevermarried, over(treatment)
tab nevermarried treatment, chi
svy: prop nevermarried, over(treatment)
svy: tabulate nevermarried treatment, col se

prop emp, over(treatment)
tab emp treatment, chi
svy: prop emp, over(treatment)
svy: tabulate emp treatment, col se

svy: prop desworktime, over(treatment)
svy: tabulate desworktime treatment, col se

svy: prop numchild, over(treatment)
svy: tabulate numchild treatment, col se

svy: prop rightswh, over(treatment)
svy: tabulate rightswh treatment, col se

svy: prop workhours, over(treatment)
svy: tabulate workhours treatment, col se

svy: prop more6wobreak, over(treatment)
svy: tabulate more6wobreak treatment, col se

prop consent if treatment == 1
prop consent if treatment == 2
prop consent if treatment == 3
svy: prop consent if inlist(treatment,1,2,3)


svy: tabulate consent treatment, se col

svy: tabulate consent treatment if inlist(treatment, 1,2), se col


svy: logit consent i.treatment if inlist(treatment,1,2)

test [consent]2.treatment == 0


svy: logit consent i.treatment if inlist(treatment,1,3)

test [consent]3.treatment == 0


tab consent age_cat if treatment == 1, col 
tab consent age_cat if treatment == 2, col 
tab consent age_cat if treatment == 3, col 
tab consent age_cat, col 
logit consent i.age_cat##i.treatment
test 2.age_cat#2.treatment=2.age_cat#3.treatment=3.age_cat#2.treatment=3.age_cat#3.treatment=4.age_cat#2.treatment=4.age_cat#3.treatment=0
svy: tabulate consent age_cat if treatment == 1, col se ci
svy: tabulate consent age_cat if treatment == 2, col se ci
svy: tabulate consent age_cat if treatment == 3, col se ci
svy: tabulate consent age_cat, col se
svy: logit consent i.age_cat##i.treatment
test 2.age_cat#2.treatment=2.age_cat#3.treatment=3.age_cat#2.treatment=3.age_cat#3.treatment=4.age_cat#2.treatment=4.age_cat#3.treatment=0


tab consent female if treatment == 1, col 
tab consent female if treatment == 2, col 
tab consent female if treatment == 3, col 
tab consent female, col 
logit consent female##i.treatment
test 1.female#2.treatment=1.female#3.treatment=0
svy: tabulate consent female if treatment == 1, col se
svy: tabulate consent female if treatment == 2, col se
svy: tabulate consent female if treatment == 3, col se
svy: tabulate consent female, col se
svy: logit consent female##i.treatment
test 1.female#2.treatment=1.female#3.treatment=0

tab consent educ if treatment == 1, col 
tab consent educ if treatment == 2, col 
tab consent educ if treatment == 3, col 
tab consent educ, col 
logit consent educ##i.treatment
test 2.educ#2.treatment=2.educ#3.treatment=3.educ#2.treatment=3.educ#3.treatment
svy: tabulate consent educ if treatment == 1, col se
svy: tabulate consent educ if treatment == 2, col se
svy: tabulate consent educ if treatment == 3, col se
svy: tabulate consent educ, col se
svy: logit consent educ##i.treatment
test 2.educ#2.treatment=2.educ#3.treatment=3.educ#2.treatment=3.educ#3.treatment=0

tab consent mig if treatment == 1, col 
tab consent mig if treatment == 2, col 
tab consent mig if treatment == 3, col 
tab consent mig, col 
logit consent mig##i.treatment
test 1.mig#2.treatment=1.mig#3.treatment=0
svy: tabulate consent mig if treatment == 1, col se
svy: tabulate consent mig if treatment == 2, col se
svy: tabulate consent mig if treatment == 3, col se
svy: tabulate consent mig, col se
svy: logit consent mig##i.treatment
test 1.mig#2.treatment=1.mig#3.treatment=0

tab consent nevermarried if treatment == 1, col 
tab consent nevermarried if treatment == 2, col 
tab consent nevermarried if treatment == 3, col 
tab consent nevermarried, col 
logit consent nevermarried##i.treatment
test 1.nevermarried#2.treatment=1.nevermarried#3.treatment=0
svy: tabulate consent nevermarried if treatment == 1, col 
svy: tabulate consent nevermarried if treatment == 2, col 
svy: tabulate consent nevermarried if treatment == 3, col 
svy: tabulate consent nevermarried, col 
svy: logit consent nevermarried##i.treatment
test 1.nevermarried#2.treatment=1.nevermarried#3.treatment=0

tab consent emp if treatment == 1, col 
tab consent emp if treatment == 2, col 
tab consent emp if treatment == 3, col 
tab consent emp, col 
logit consent emp##i.treatment
test 2.emp#2.treatment=2.emp#3.treatment=3.emp#2.treatment=3.emp#3.treatment
svy: tabulate consent emp if treatment == 1, col 
svy: tabulate consent emp if treatment == 2, col 
svy: tabulate consent emp if treatment == 3, col 
svy: tabulate consent emp, col 
svy: logit consent emp##i.treatment
test 2.emp#2.treatment=2.emp#3.treatment=3.emp#2.treatment=3.emp#3.treatment=0


svy: tabulate consent more6wobreak if treatment == 1, col 
svy: tabulate consent more6wobreak if treatment == 2, col 
svy: tabulate consent more6wobreak if treatment == 3, col 
svy: tabulate consent more6wobreak, col 
svy: logit consent i.more6wobreak##i.treatment
test 1.more6wobreak#2.treatment=1.more6wobreak#3.treatment=0

svy: tabulate consent workhours if treatment == 1, col se
svy: tabulate consent workhours if treatment == 2, col se
svy: tabulate consent workhours if treatment == 3, col se
svy: tabulate consent workhours, col se
svy: logit consent i.workhours##i.treatment
test 2.workhours#2.treatment=2.workhours#3.treatment=3.workhours#2.treatment=3.workhours#3.treatment=0

svy: tabulate consent numchild if treatment == 1, col se
svy: tabulate consent numchild if treatment == 2, col se
svy: tabulate consent numchild if treatment == 3, col se
svy: tabulate consent numchild, col se
svy: logit consent numchild##i.treatment
test 2.numchild#2.treatment=2.numchild#3.treatment=3.numchild#2.treatment=3.numchild#3.treatment=0

svy: tabulate consent rightswh if treatment == 1, col se
svy: tabulate consent rightswh if treatment == 2, col se
svy: tabulate consent rightswh if treatment == 3, col se
svy: tabulate consent rightswh, col se
svy: logit consent rightswh##i.treatment
test 1.rightswh#2.treatment=1.rightswh#3.treatment=0

svy: tabulate consent desworktime if treatment == 1, col 
svy: tabulate consent desworktime if treatment == 2, col 
svy: tabulate consent desworktime if treatment == 3, col 
svy: tabulate consent desworktime, col 
svy: logit consent ib3.desworktime##i.treatment
test 1.desworktime#2.treatment=1.desworktime#3.treatment=2.desworktime#2.treatment=2.desworktime#3.treatment=0



prop consent, over(age_cat)
logit consent i.age_cat
svy: prop consent, over(age_cat)
svy: logit consent i.age_cat

prop consent, over(female)
logit consent i.female
svy: prop consent, over(female)
svy: logit consent i.female

prop consent, over(educ)
logit consent i.educ
svy: prop consent, over(educ)
svy: logit consent i.educ

prop consent, over(mig)
logit consent mig
svy: prop consent, over(mig)
svy: logit consent mig

prop consent, over(nevermarried)
logit consent nevermarried
svy: prop consent, over(nevermarried)
svy: logit consent nevermarried

prop consent, over(emp)
logit consent i.emp
svy: prop consent, over(emp)
svy: logit consent i.emp


svy: prop consent, over(desworktime)
svy: logit consent ib3.desworktime

svy: prop consent, over(numchild)
svy: logit consent i.numchild

svy: prop consent, over(rightswh)
svy: logit consent i.rightswh

svy: prop consent, over(workhours)
svy: logit consent i.workhours

svy: prop consent, over(more6wobreak)
svy: logit consent i.more6wobreak

svyset int_id [pweight=hoch], fpc(fpc1) || betnr_atypik, strata(risiko) fpc(fpc2) singleunit(scaled) || lfd, fpc(fpc3)
  
svy: logit consent i.treatment i.age_cat i.female  i.mig i.educ##i.treatment i.workhours##i.treatment i.numchild##i.treatment i.rightswh##i.treatment 

gen timesaving = 1 if treatment == 2
replace timesaving = 0 if inlist(treatment,1,3)

gen studyvalue = 1 if treatment == 3
replace studyvalue = 0 if inlist(treatment,1,2)


prop consent if treatment == 1, over(age_cat)
prop consent if treatment == 2, over(age_cat)
prop consent if treatment == 3, over(age_cat)
svy: prop consent if treatment == 1, over(age_cat)
svy: prop consent if treatment == 2, over(age_cat)
svy: prop consent if treatment == 3, over(age_cat)

prop consent if treatment == 1, over(female)
prop consent if treatment == 2, over(female)
prop consent if treatment == 3, over(female)
svy: prop consent if treatment == 1, over(female)
svy: prop consent if treatment == 2, over(female)
svy: prop consent if treatment == 3, over(female)


prop consent if treatment == 1, over(educ)
prop consent if treatment == 2, over(educ)
prop consent if treatment == 3, over(educ)
svy: prop consent if treatment == 1, over(educ)
svy: prop consent if treatment == 2, over(educ)
svy: prop consent if treatment == 3, over(educ)

prop consent if treatment == 1, over(mig)
prop consent if treatment == 2, over(mig)
prop consent if treatment == 3, over(mig)
svy: prop consent if treatment == 1, over(mig)
svy: prop consent if treatment == 2, over(mig)
svy: prop consent if treatment == 3, over(mig)

prop consent if treatment == 1, over(nevermarried)
prop consent if treatment == 2, over(nevermarried)
prop consent if treatment == 3, over(nevermarried)
svy: prop consent if treatment == 1, over(nevermarried)
svy: prop consent if treatment == 2, over(nevermarried)
svy: prop consent if treatment == 3, over(nevermarried)

prop consent if treatment == 1, over(emp)
prop consent if treatment == 2, over(emp)
prop consent if treatment == 3, over(emp)
svy: prop consent if treatment == 1, over(emp)
svy: prop consent if treatment == 2, over(emp)
svy: prop consent if treatment == 3, over(emp)


svy: prop consent if treatment == 1, over(more6wobreak)
svy: prop consent if treatment == 2, over(more6wobreak)
svy: prop consent if treatment == 3, over(more6wobreak)

svy: prop consent if treatment == 1, over(desworktime)
svy: prop consent if treatment == 2, over(desworktime)
svy: prop consent if treatment == 3, over(desworktime)

svy: prop consent if treatment == 1, over(workhours)
svy: prop consent if treatment == 2, over(workhours)
svy: prop consent if treatment == 3, over(workhours)

svy: prop consent if treatment == 1, over(numchild)
svy: prop consent if treatment == 2, over(numchild)
svy: prop consent if treatment == 3, over(numchild)

svy: prop consent if treatment == 1, over(rightswh)
svy: prop consent if treatment == 2, over(rightswh)
svy: prop consent if treatment == 3, over(rightswh)

**Interviewer observation analysis

tab F5_tgroup

tab F6_wortlaut

*************Treatment Effects (Differences between 

svy: prop consent if age_cat == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if age_cat == 2, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if age_cat == 3, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if age_cat == 4, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if female == 0, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if female == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if mig == 0, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if mig == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if educ == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if educ == 2, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if educ == 3, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if workhours == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if workhours == 2, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if workhours == 3, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if numchild == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if numchild == 2, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if numchild == 3, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if rightswh == 0, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if rightswh == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if desworktime == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if desworktime == 2, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if desworktime == 3, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1

svy: prop consent if more6wobreak == 1, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1
svy: prop consent if more6wobreak == 0, over(treatment)
lincom [_prop_2]2 - [_prop_2]1
lincom [_prop_2]3 - [_prop_2]1