Home » Topics » Nutrition and Anthropometry » nutritional status
nutritional status [message #16186] |
Mon, 19 November 2018 11:25 |
radhs
Messages: 4 Registered: November 2018
|
Member |
|
|
HI!!
I am trying to relate the nutritional status haz whz waz mean z score by age in months in kr file.Can you please improvise my code, It would be really helpful.
gen weight = v005/1000000
gen psu =v021
gen strata = v022
svyset psu [pw = weight], strata(strata) vce(linearized)
recode hw1 (0/5 = 1 "<6") (6/8 = 2 "6-8") (9/11 = 3 "9-11") (12/17 = 4 "12-17") (18/23 = 5 "18-23") (24/35 = 6 "24/35") (36/47 = 7 "36-47") (48/59 = 8 "48-59"), gen(child_age)
gen whz=hw72
replace whz=. if whz>=9996
histogram whz
gen wasted=.
replace wasted=0 if whz ~=.
replace wasted=1 if whz <-200
gen waz=hw71
replace waz=. if waz>=9996
histogram waz
gen under=.
replace under=0 if waz ~=.
replace under=1 if waz <-200
gen HAZ=hw70
replace HAZ=. if HAZ>=9996
histogram HAZ
gen stunted=.
replace stunted=0 if HAZ ~=.
replace stunted=1 if HAZ <-200
tab child_age stunted [iw=v005/1000000],miss
qui regr stunted under wasted if stunted !=. & under !=. & wasted !=. [pw=v005/1000000]
drop if e(sample)!=1
tab child_age stunted [iw=v005/1000000],miss
<6 | 14,098.17 3,500.687 | 17,598.86
6-8 | 9,147.855 2,325.8678 | 11,473.72
9-11 | 7,896.852 2,771.066 | 10,667.92
12-17 | 13,420.65 8,316.5437 | 21,737.2
18-23 | 11,639.71 10,335.97 | 21,975.68
24/35 | 24,758.8 18,443.34 | 43,202.14
36-47 |25,827.762 19,671.92 | 45,499.69
48-59 | 26,023.61 17,310.86 | 43,334.46
-----------+----------------------+----------
Total | 132,813.4 82,676.25 | 215,489.7
IT should be equal to 219796.
17644 for <6.
svy:mean HAZ if child_age==1
[Updated on: Mon, 19 November 2018 14:16] Report message to a moderator
|
|
|
Goto Forum:
Current Time: Mon Dec 23 14:24:01 Coordinated Universal Time 2024
|