* bmi_setup.do
infile using "F:\flash\academica\bmi.DCT", clear
codebook, compact
renvars _all, lowcase
gen male = r0536300
replace male = 0 if r0536300 == 2
tab1 r0538600 r0538700
gen race_eth = r0538700
replace race_eth = 3 if r0538600 == 1
replace race_eth = 5 if r0538700 == 3
tab race_eth, m
gen white = 1 if race_eth ==1
replace white = 0 if race_eth ~= 1 & race_eth < .
gen black = 1 if race_eth ==2
replace black = 0 if race_eth ~= 2 & race_eth < .
gen hispanic = 1 if race_eth ==3
replace hispanic = 0 if race_eth ~= 3 & race_eth < .
gen asian = 1 if race_eth ==4
replace asian = 0 if race_eth ~= 4 & race_eth < .
gen other = 1 if race_eth ==5
replace other = 0 if race_eth ~= 5 & race_eth < .
tab1 race_eth white black hispanic asian other race_eth, m


mvdecode _all, mv(-1=.a\-2=.b\-3=.c\-4=.d\-5=.e)
label define gender 1 "male" 0 "female" .a "refuse" .b "don't know" ///
  .c "invalid skip" .d "valid skip" .e "non-interview"
label define race 1 "white" 2 "black" 3 "hispanic" 4 "asian" 5 "other" ///
  .a "refuse" .b "don't know" .c "invalid skip" .d "valid skip" ///
  .e "non-interview"
label define other .a "refuse" .b "don't know" .c "invalid skip" ///
  .d "valid skip" .e "non-interview"
label values male gender
label values race_eth race
label values r0322500 other
label values r0322600 other
label values r0322700 other
label values r0536300 other
label values r0536700 other
label values r0538600 other
label values r0538700 other
label values r1487300 other
label values r1487400 other
label values r1487500 other
label values r1487600 other
label values r2164100 other
label values r2164200 other
label values r2164300 other
label values r3482000 other
label values r3482100 other
label values r3482200 other
label values r4880200 other
label values r4880300 other
label values r4880400 other
label values r6497600 other
label values r6497700 other
label values r6497800 other
tab1  r0536300 r0536700 r0538600 r0538700 r1487300 ///
  r1487400 r1487500 r1487600 male race_eth
gen hght97 = (r0322500 * 12) + r0322600
gen hght98 = (r2164100 * 12) + r2164200
gen hght99 = (r3482000 * 12) + r3482100
gen hght00 = (r4880200 * 12) + r4880300
gen hght01 = (r6497600 * 12) + r6497700
gen hght02 = (s0905500 * 12) + s0905600
gen hght03 = (s2978200 * 12) + s2978300
list r0322500 r0322600 hght97 in 3000/3010, nolabel
list r2164100 r2164200 hght98 in 3000/3010, nolabel
list r3482100 r3482200 hght99 in 3000/3010, nolabel
list r4880200 r4880300 hght00 in 3000/3010, nolabel
list r6497600 r6497700 hght01 in 3000/3010, nolabel
list s0905500 s0905600 hght02 in 3000/3010, nolabel
list s2978200 s2978300 hght03 in 3000/3010, nolabel
gen bmi97 = (r0322700/hght97^2) * 703
gen bmi98 = (r2164300/hght98^2) * 703
gen bmi99 = (r3482200/hght99^2) * 703
gen bmi00 = (r4880400/hght00^2) * 703
gen bmi01 = (r6497800/hght01^2) * 703
gen bmi02 = (s0905700/hght02^2) * 703
gen bmi03 = (s2978400/hght03^2) * 703
list bmi97 bmi98 bmi99 bmi00 bmi01 bmi02 bmi03 in 3000/3020
rename r0536700 age
sum bmi97 bmi98 bmi99 bmi00 bmi01 bmi02 bmi03 if age == 12 

save "F:\flash\academica\bmi_working.dta"
keep if age==12
by male, sort : summarize bmi97 bmi98 bmi99 bmi00 bmi01 bmi02 bmi03
mvpatt bmi97-bmi03
by race_eth, sort : summarize bmi97 bmi98 bmi99 bmi00 bmi01 bmi02 bmi03
mvpatt bmi97-bmi03
rename r0000100 id
rename r1487300 grlprb_y
rename r1487400 boyprb_y
rename r1487500 grlprb_p
rename r1487600 boyprb_p
keep id grlprb_y boyprb_y grlprb_p boyprb_p male race_eth ///
  white black hispanic asian other bmi97 bmi98 bmi99 bmi00 ///
  bmi01 bmi02 bmi03
save "F:\flash\academica\bmi_stata.dta"
stata2mplus using "F:\flash\academica\bmi_stata"