/* 

easySHARE 

Zakladni kod pro Stata 

Vice informaci na www.share-project.org
http://share.cerge-ei.cz

Kontakt:
Radim Bohacek
SHARE CZ
NHU AV CR
Politickych veznu 7
111 21 Praha 1
radim.bohacek@gmail.com
tel: 775 379 336

Nize je kod pro prohlizeni dat, kodovani,
regresi panelovou a pro jednu vlnu, a probit.

Priklad pouziti dat neni nijak slavny.

*/

clear
clear matrix
version 12 
set mem 400m
set more off
set matsize 800

/* zde upravte dle ulozeni souboru */
use D:\work\SHARE\easySHARE\data\stata\easyshare_rel1-0-0.dta, clear

/* Priprava dat pro analyzu (BMI) panel */
sum bmi_mod,d  /* vyskytuji se missing variables */
tab bmi_mod if bmi_mod<0
mvdecode bmi_mod, mv(-1/-15)  /* zakodovat missing variables */
bysort country: sum bmi_mod  /* zobrazit podle zemi */

/* dalsi promenne */
mvdecode eduyears_mod, mv(-1/-15)  /* zakodovat missing variables */
mvdecode income_pct_w4, mv(-1/-15)  /* zakodovat missing variables */

mvdecode maxgrip, mv(-1/-15)  /* zakodovat missing variables */
mvdecode age, mv(-1/-15)  /* zakodovat missing variables */
mvdecode ep005_, mv(-1/-15)  /* zakodovat missing variables */
mvdecode smoking, mv(-1/-15)  /* zakodovat missing variables */
mvdecode br010_mod, mv(-1/-15)  /* zakodovat missing variables */
mvdecode income_pct_w4, mv(-1/-15)  /* zakodovat missing variables */
mvdecode euro1, mv(-1/-15)  /* zakodovat missing variables */

/* zobrazit nektere */
tab eduyears_mod
tab income_pct_w4
tab female

/* prekodovat: prijem v top 20 percentilech */
tab  income_pct_w4
gen top_income = .
replace top_income = 1 if income_pct_w4 > 8
replace top_income = 0 if income_pct_w4 < 9

/* prekodovat: alkohol */ 
tab br010_mod
gen alcohol = .
replace alcohol = 1 if br010_mod > 4
replace alcohol = 0 if br010_mod < 5 

/* prekodovat: retired */
tab ep005_
gen retired = .
replace retired = 1 if ep005_ == 1
replace retired = 0 if ep005_ > 1

/* ulozit data pro analyzu */
save D:\work\SHARE\easySHARE\data\stata\bmi_data.dta, replace

/* ************************************************************ */
/* Priklad: body-mass index (BMI) panel */
/* upravit data pro panel */
use D:\work\SHARE\easySHARE\data\stata\bmi_data.dta, clear
keep if wavepart==234 /* CR ma pouze vlnu 2, 3 a 4 */

bysort mergeid: gen mergeid_n = _n ==1  /* kodovat cislo respondenta pro panel (viz ReleaseGuide) */
replace mergeid_n = sum(mergeid_n)

xtset mergeid_n wave  /* priprava dat pro panelovou regresi */

xtreg bmi_mod  eduyears_mod top_income female  retired maxgrip euro1 alcohol smoking if country~=28, mle
xtreg bmi_mod  eduyears_mod top_income female  retired maxgrip euro1 alcohol smoking if country==28, mle

/* ************************************************************ */
/* Priklad: body-mass index (BMI) pouze pro respondenty 4. vlny */
use D:\work\SHARE\easySHARE\data\stata\bmi_data.dta, clear
keep if wave==4

regress bmi_mod  eduyears_mod top_income female  retired maxgrip euro1 alcohol smoking if country~=28
regress bmi_mod  eduyears_mod top_income female  retired maxgrip euro1 alcohol smoking if country==28 /* pouze pro CR */

/* obesita pro BMI>30 */
gen obese = .
replace obese = 1 if bmi_mod >= 30
replace obese = 0 if bmi_mod < 30

tab obese if country==28

probit obese  eduyears_mod top_income female  retired maxgrip euro1 alcohol smoking if country~=28
probit obese  eduyears_mod top_income female  retired maxgrip euro1 alcohol smoking if country==28 /* pouze pro CR */