clear
set mem 300m
cd e:\otherresearch\BrazilNikkei\ParanaData\

* FIRST GET HOUSEHOLD CHARACTERISTICS

* Round 1
use "e:\otherresearch\BrazilNikkei\ParanaData\bddomiciliarparana_l.dta", clear

rename v08a v08_1
rename v08b v08_2
rename v08c v08_3
rename v08d v08_4
rename v08e v08_5
rename v08f v08_6
rename v08g v08_7
rename v08h v08_8
rename v08i v08_9

* Household characteristics
gen hhsize=1 if v10a1<=10
for num 2/9: replace hhsize=hhsize+1 if v10aX<=10
gen size18to60=1 if v10a1<=10 & v08_1>=18 & v08_1<=60
for num 2/9: replace size18to60=size18to60+1 if v10aX<=10 & v08_X>=18 & v08_X<=60
replace size18to60=0 if hhsize~=. & size18to60==.

gen readspaper=v78a==1
gen listenradio=v78b==1
gen watchtv=v78c==1
gen readbook=v78d==1
gen readbunk=v78e==1
gen checkweb=v78f==1

* current migrant
gen currentin=1 if v09b1==2
for num 2/9: replace currentin=1 if v09bX==2
replace currentin=0 if currentin==.

* return migrant
gen hhinjapan=1 if v09b1==3 
for num 2/9: replace hhinjapan=1 if v09bX==3
replace hhinjapan=0 if hhinjapan==. 

* Remittance receipts
gen receiveremit=1 if v66<=8|v66==10|v66==12
replace receiveremit=0 if receiveremit==.
gen refuseremit=v66==.
gen amountreceive=v62a
replace amountreceive=0 if v62a==.
replace amountreceive=amountreceive+v62b if v62b~=.
replace amountreceive=amountreceive+v62c if v62c~=.
gen refuseamount=1 if receiveremit==1 & amountreceive==0
replace refuseamount=0 if receiveremit==1 & amountreceive>0 & amountreceive~=.

#delimit ;
keep quest va vb hhsize size18to60 readspaper listenradio watchtv readbook readbunk checkweb
hhinjapan currentin receiveremit refuseremit refuseamount amountreceive
;
#delimit cr
gen round1p=1
save poolhh1p.dta, replace


* Round 2 - long
use "e:\otherresearch\BrazilNikkei\ParanaData\domiciliarparana2.dta", clear

rename v08a v08_1
rename v08b v08_2
rename v08c v08_3
rename v08d v08_4
rename v08e v08_5
rename v08f v08_6
rename v08g v08_7
rename v08h v08_8
rename v08i v08_9

* Household characteristics
gen hhsize=1 if v10a1<=10
for num 2/9: replace hhsize=hhsize+1 if v10aX<=10
gen size18to60=1 if v10a1<=10 & v08_1>=18 & v08_1<=60
for num 2/9: replace size18to60=size18to60+1 if v10aX<=10 & v08_X>=18 & v08_X<=60
replace size18to60=0 if hhsize~=. & size18to60==.

gen readspaper=v66a==1
gen listenradio=v66b==1
gen watchtv=v66c==1
gen readbook=v66d==1
gen readbunk=v66e==1
gen checkweb=v66f==1

* current migrant
gen currentin=1 if v09b1==2
for num 2/9: replace currentin=1 if v09bX==2
replace currentin=0 if currentin==.

* return migrant
gen hhinjapan=1 if v09b1==3 
for num 2/9: replace hhinjapan=1 if v09bX==3
replace hhinjapan=0 if hhinjapan==. 

* Remittance receipts
gen receiveremit=1 if  v54<=8|v54==10|v54==12
replace receiveremit=0 if receiveremit==.
gen refuseremit=v54==.
gen amountreceive=v50a
replace amountreceive=0 if v50a==.
replace amountreceive=amountreceive+v50b if v50b~=.
replace amountreceive=amountreceive+v50c if v50c~=.
gen refuseamount=1 if receiveremit==1 & amountreceive==0
replace refuseamount=0 if receiveremit==1 & amountreceive>0 & amountreceive~=.

#delimit ;
keep quest va vb hhsize size18to60 readspaper listenradio watchtv readbook readbunk checkweb
hhinjapan currentin receiveremit refuseremit refuseamount amountreceive;
#delimit cr
gen round2p=1
save poolhh2p.dta, replace


* Round 2 short/aka round 3
use "e:\otherresearch\BrazilNikkei\ParanaData\bddomiciliarparana3_nol.dta", clear
drop if va==.

gen hhsize=v08
replace hhsize=. if hhsize==99
gen size18to60=v09
replace size18to60=. if v09==99|v09==98

gen readspaper=v16==1
replace readspaper=. if v16>2
gen listenradio=v17==1
replace listenradio=. if v17>2
gen watchtv=v18==1
replace watchtv=. if v18>2
gen readbook=v19==1
replace readbook=. if v19>2
gen readbunk=v20==1
replace readbunk=. if v20>2
gen checkweb=v21==1
replace checkweb=. if v21>2

* Migration
gen workinjapan=v13==1
replace workinjapan=. if v13>2
gen hhinjapan=workinjapan
replace hhinjapan=1 if v14==1
gen currentin=v15==1
replace currentin=. if v15>2

*Remittances
gen receiveremit=v25<=60
gen refuseremit=v25==80|v25==90|v25==.
rename v2b v26b
gen amountreceive=v26b if v26b~= 999997 & v26b~=999999
replace amountreceive=0 if amountreceive==.
replace amountreceive=amountreceive+v27b if v27b~= 999997 & v27b~=999999
replace amountreceive=amountreceive+v28b if v28b~= 999997 & v28b~=999999
replace amountreceive=amountreceive+v29b if v29b~= 999997 & v29b~=999999
replace amountreceive=amountreceive+v30b if v30b~= 999997 & v30b~=999999
replace amountreceive=amountreceive+v31b if v31b~= 999997 & v31b~=999999
replace amountreceive=amountreceive+v32b if v32b~= 999997 & v32b~=999999
replace amountreceive=amountreceive+v33b if v33b~= 999997 & v33b~=999999
replace amountreceive=amountreceive+v34b if v34b~= 999997 & v34b~=999999
gen refuseamount=1 if receiveremit==1 & amountreceive==0
replace refuseamount=0 if refuseamount==.



#delimit ;
keep quest va vb hhsize size18to60 readspaper listenradio watchtv readbook readbunk checkweb
hhinjapan currentin receiveremit refuseremit refuseamount amountreceive;
#delimit cr
gen round3p=1
save poolhh3p.dta, replace
append using poolhh1p.dta
append using poolhh2p.dta
gen round=1 if round1p==1
replace round=2 if round2p==1
replace round=3 if round3p==1
save poolhhp.dta, replace

*****************
*****************
* Merge Weights and Matching file from listing
use bd123parana.dta, clear
gen double ap=v503
sort ap
merge ap using weightp.dta
tab _merge
drop if _merge==2
drop _merge
sort quest
save mastermatchp.dta, replace

use poolhhp.dta, clear
sort quest
merge quest using mastermatchp.dta
tab _merge
rename _merge mergemaster
save poolhhp.dta, replace


rename weight sampleweight
#delimit ;
keep quest va vb ap hhsize size18to60 readspaper listenradio watchtv readbook readbunk checkweb
hhinjapan currentin receiveremit refuseremit refuseamount amountreceive sampleweight round;
#delimit cr
gen method=1
gen state="Parana"
save poolhhp.dta, replace

use poolhhp.dta, clear
svyset ap [pweight=sampleweight]

svy: mean hhsize 
svy: mean size18to60 
svy: mean readspaper
svy: mean listenradio
svy: mean watchtv 
svy: mean readbook 
svy: mean readbunk 
svy: mean checkweb
svy: mean age 
svy: mean female 
svy: mean married 
svy: mean currentin 
svy: mean hhinjapan 
svy: mean receiveremit 
svy: mean refuseremit 
svy: mean refuseamount if receiveremit==1
svy: mean amountreceive if receiveremit==1 & amountreceive>0
table refuseamount [pw=sampleweight] if amountreceive>0 & receiveremit==1, c(median amountreceive)


save "E:\otherresearch\BrazilNikkei\PooledSamples\paranastratpool.dta"
append using E:\otherresearch\BrazilNikkei\PooledSamples\paranainterceptpool.dta
drop if vb==. & method==1
replace readbook=readbooks if method==2
replace readbunk=readbunkyo if method==2
replace currentin=currentinjapan if method==2
replace receiveremit=receiveremittances if method==2 
replace refuseremit=refuseremitt if method==2
svyset [pweight=sampleweight]
svy: mean hhsize, over(method)
lincom [hhsize]1-[hhsize]2
svy: mean size18to60, over(method)
lincom [size18to60]1-[size18to60]2
svy: mean readspaper, over(method)
lincom [readspaper]1-[readspaper]2
svy: mean listenradio, over(method)
lincom [listenradio]1-[listenradio]2
svy: mean watchtv, over(method)
lincom [watchtv]1-[watchtv]2
svy: mean readbook, over(method)
lincom [readbook]1-[readbook]2
svy: mean readbunk, over(method)
lincom [readbunk]1-[readbunk]2
svy: mean checkweb, over(method)
lincom [checkweb]1-[checkweb]2

svy: mean currentin, over(method) 
lincom [currentin]1-[currentin]2
svy: mean hhinjapan, over(method) 
lincom [hhinjapan]1-[hhinjapan]2
svy: mean receiveremit, over(method) 
lincom [receiveremit]1-[receiveremit]2
svy: mean refuseremit, over(method) 
lincom [refuseremit]1-[refuseremit]2
svy: mean refuseamount if receiveremit==1, over(method)
lincom [refuseamount]1-[refuseamount]2
svy: mean amountreceive if receiveremit==1 & amountreceive>0, over(method)
lincom [amountreceive]1-[amountreceive]2


gen unweightweight=sampleweight if method==1
replace unweightweight=1 if method==2
svyset [pweight=unweightweight]
svy: mean hhsize, over(method)
lincom [hhsize]1-[hhsize]2
svy: mean size18to60, over(method)
lincom [size18to60]1-[size18to60]2
svy: mean readspaper, over(method)
lincom [readspaper]1-[readspaper]2
svy: mean listenradio, over(method)
lincom [listenradio]1-[listenradio]2
svy: mean watchtv, over(method)
lincom [watchtv]1-[watchtv]2
svy: mean readbook, over(method)
lincom [readbook]1-[readbook]2
svy: mean readbunk, over(method)
lincom [readbunk]1-[readbunk]2
svy: mean checkweb, over(method)
lincom [checkweb]1-[checkweb]2

svy: mean currentin, over(method) 
lincom [currentin]1-[currentin]2
svy: mean hhinjapan, over(method) 
lincom [hhinjapan]1-[hhinjapan]2
svy: mean receiveremit, over(method) 
lincom [receiveremit]1-[receiveremit]2
svy: mean refuseremit, over(method) 
lincom [refuseremit]1-[refuseremit]2
svy: mean refuseamount if receiveremit==1, over(method)
lincom [refuseamount]1-[refuseamount]2
svy: mean amountreceive if receiveremit==1 & amountreceive>0, over(method)
lincom [amountreceive]1-[amountreceive]2


* Look at original sample vs those in round 2/3
use "E:\otherresearch\BrazilNikkei\PooledSamples\paranastratpool.dta", clear
gen roundsum=round==1
drop if vb==.
svyset [pweight=sampleweight]
svy: mean hhsize, over(roundsum)
lincom [hhsize]1-[hhsize]0
svy: mean readspaper, over(roundsum)
lincom [readspaper]1-[readspaper]0
svy: mean listenradio, over(roundsum)
lincom [listenradio]1-[listenradio]0
svy: mean watchtv, over(roundsum)
lincom [watchtv]1-[watchtv]0
svy: mean readbook, over(roundsum)
lincom [readbook]1-[readbook]0
svy: mean readbunk, over(roundsum)
lincom [readbunk]1-[readbunk]0
svy: mean checkweb, over(roundsum)
lincom [checkweb]1-[checkweb]0

svy: mean currentin, over(roundsum)
lincom [currentin]1-[currentin]0
svy: mean hhinjapan, over(roundsum)
lincom [hhinjapan]1-[hhinjapan]0
svy: mean receiveremit, over(roundsum)
lincom [receiveremit]1-[receiveremit]0
svy: mean refuseremit, over(roundsum)
lincom [refuseremit]1-[refuseremit]0
svy: mean refuseamount if receiveremit==1, over(roundsum)
lincom [refuseamount]1-[refuseamount]0
svy: mean amountreceive if receiveremit==1 & amountreceive>0, over(roundsum)
lincom [amountreceive]1-[amountreceive]0


* PERSON LEVEL DATA


gen issei1=(v16b1==1 & v08_1>=18 & v08_1<99)
gen nissei1=(v16b1==2 & v08_1>=18 & v08_1<99)
gen sansei1=(v16b1==3 & v08_1>=18 & v08_1<99)
gen yonsei1=(v16b1==4 & v08_1>=18 & v08_1<99)

for num 2/9: gen isseiX=(v16bX==1 & v08_X>=18 & v08_X<99)
for num 2/9: gen nisseiX=(v16bX==2 & v08_X>=18 & v08_X<99)
for num 2/9: gen sanseiX=(v16bX==3 & v08_X>=18 & v08_X<99)
for num 2/9: gen yonseiX=(v16bX==4 & v08_X>=18 & v08_X<99)
for num 1/9: gen nikkeiadultX=(v16bX<=4 & v08_X>=18 & v08_X<99)

#delimit ;
gen hhadultnikkei=nikkeiadult1+nikkeiadult2+nikkeiadult3+nikkeiadult4+nikkeiadult5
+nikkeiadult6+nikkeiadult7+nikkeiadult8+nikkeiadult9;
gen numissei=issei1+issei2+issei3+issei4+issei5+issei6+issei7+issei8+issei9;
gen numnissei=nissei1+nissei2+nissei3+nissei4+nissei5+nissei6+nissei7+nissei8+nissei9;
gen numsansei=sansei1+sansei2+sansei3+sansei4+sansei5+sansei6+sansei7+sansei8+sansei9;
gen numyonsei=yonsei1+yonsei2+yonsei3+yonsei4+yonsei5+yonsei6+yonsei7+yonsei8+yonsei9;
#delimit cr
sum hhadultnikkei numissei numnissei numsansei numyonsei if interviewyes==1