*Rodriguez
*Epi 265, Spring 2017

*HW #8 4 scenarios; 1st S2 = 0.9; 2nd S2 = 0.5:; 3rd S1 = 0.9; 4th S1 = 0.5 
set obs 1000
set seed 11717
gen error = rnormal()
gen u_rural = runiform()>0.7 + error
gen S2 = runiform()>0.2 + error if u_rural ==1
replace S2 = runiform()>0.05 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0 
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 15*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.05 + error if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.1 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.15 + error if S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S2==1
margins, dydx(severe_HFI) 
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S2==1


*2nd scenario  S2 = ~0.57:
set obs 1000
set seed 11718
gen error = rnormal()
gen u_rural = runiform()>0.7 + error
gen S2 = runiform()>0.6 + error if u_rural ==1
replace S2 = runiform()>0.4 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 15*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.05 + error if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.1 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.15 + error if  S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S2==1
margins, dydx(severe_HFI)
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S2==1

*3rd scenario S1 = 0.9: 
set obs 1000
set seed 11719
gen error = rnormal()
gen u_rural = runiform()>0.7 + error 
gen S2 = runiform()>0.2 + error if u_rural ==1
replace S2 = runiform()>0.05 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 15*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.02 + error if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.95 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.6 + error if S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S1==1
margins, dydx(severe_HFI)
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S1==1
 
*4th scenario S1 = 0.5: 
set obs 1000
set seed 11720
gen error = rnormal()
gen u_rural = runiform()>0.7 + error
gen S2 = runiform()>0.2 + error if u_rural ==1
replace S2 = runiform()>0.05 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 15*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.2 + error if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.6 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.7 + error if S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S1==1
margins, dydx(severe_HFI)
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S1==1

****Now assume the null under each of the 4 scenarios above:
set obs 1000
gen error = rnormal()
gen u_rural = runiform()>0.7 + error
gen S2 = runiform()>0.2 + error if u_rural ==1
replace S2 = runiform()>0.05 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0 
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 0*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.05 + error if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.1 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.15 + error if S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S2==1
margins, dydx(severe_HFI) 
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S2==1


*2nd scenario  S2 = ~0.57:
set obs 1000
set seed 11718
gen error = rnormal()
gen u_rural = runiform()>0.7 + error
gen S2 = runiform()>0.5 + error if u_rural ==1
replace S2 = runiform()>0.4 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 0*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.05 + error if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.1 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.15 + error if  S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S2==1
margins, dydx(severe_HFI)
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S2==1

*3rd scenario S1 = 0.9: 
set obs 1000
set seed 11719
gen error = rnormal()
gen u_rural = runiform()>0.7 + error 
gen S2 = runiform()>0.2 + error if u_rural ==1
replace S2 = runiform()>0.05 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 0*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.1 if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.35 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.2 + error if S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S1==1
margins, dydx(severe_HFI)
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S1==1
 
*4th scenario S1 = 0.5: 
set obs 1000
set seed 11720
gen error = rnormal()
gen u_rural = runiform()>0.7 + error
gen S2 = runiform()>0.2 + error if u_rural ==1
replace S2 = runiform()>0.05 + error if u_rural ==0
gen severe_HFI = runiform()>0.8 + error if u_rural ==0
replace severe_HFI = runiform()>0.7 + error if u_rural ==1
gen diet_qual = 66.8 - 0*severe_HFI + 9*u_rural + 3*error
gen S1 = runiform()>0.2 + error if diet_qual>68 & severe_HFI ==0 
replace S1 = runiform()>0.6 + error if diet_qual>68 & severe_HFI ==1 
replace S1 = runiform()>0.7 + error if S1==.
gen ID = _n
order ID diet_qual severe_HFI u_rural S1 S2 error
regress diet_qual severe_HFI if S1==1
margins, dydx(severe_HFI)
bootstrap, reps(1000) size(500): regress diet_qual severe_HFI if S1==1