I am trying to estimate the point prevalence of TB per 100,000 population for all as well as childhood TB from PR file of NFHS-5 (2019-21). The variable I have used is sh29aa. I am using STATA 14 for the analysis with the following command:

gen tb_mtreated =0

replace tb_mtreated =1 if sh29aa>=1 & sh29aa<=3

proportion tb_mtreated if hv102==1 [iw=hv005/1000000]

proportion tb_mtreated if hv102==1 & hv024==23 [iw=hv005/1000000], over (shdist)

The point prevalence by age groups, rural/urban and total is matching with the national report at India level. However, when I am doing the same analysis with same command at state level, I found huge differences in point prevalence by age groups and rural/urban at state level. Also, I would like to mention that the total number of sample is matching with the state reports but not the estimates of point prevalence.

I humbly request to solve this problem on priority basis. ]]>

The problem may be that you are using hv005 and iweights (iw). Please repeat, using the state-specific weight for within-state estimates. Please let us know if you still get a discrepancy.]]>

I have already applied the state weight and IW. The number of usual residents for states which is "N" is matching with the state report while point prevalence by age-group and gender is not. For instance, point prevalence among <15 years for UP was estimated as 34/100,000 from the data but in the report, it is 356/100,000.

Here is the full command of STATA, which I have employed for the state of Uttar Pradesh (if hv024==9):

gen tb_mtreated =0

replace tb_mtreated =1 if sh29aa>=1 & sh29aa<=3

proportion tb_mtreated if hv102==1 [iw=hv005/1000000]

recode hv105 (0/14=1 "0-14 years") (15/59=2 "15-59 years") (else=3 ">=60 years"), gen (age)

proportion tb_mtreated if hv102==1 & hv024==9 [iw=shweight/1000000], over (age)

proportion tb_mtreated if hv102==1 & hv024==9 [iw=shweight/1000000], over (hv104)

proportion tb_mtreated if hv102==1 & hv024==9 [iw=shweight/1000000]

//Results from STATA//

proportion tb_mtreated if hv102==1 & hv024==9 [iw=shweight/1000000], over (age)

Proportion estimation Number of obs = 364,194

_prop_1: tb_mtreated = 0

_prop_2: tb_mtreated = 1

_subpop_1: age = <15

_subpop_2: age = 15-59

_subpop_3: age = >=60

Over Proportion Std. Err. [95% Conf. Interval]

_prop_1

_subpop_1 .9996582 .0000553 .9995307 .9997511

_subpop_2 .9975473 .0001065 .9973295 .9977473

_subpop_3 .994669 .0003812 .9938673 .9953665

_prop_2

_subpop_1 .0003418 .0000553 .0002489 .0004693

_subpop_2 .0024527 .0001065 .0022527 .0026705

_subpop_3 .005331 .0003812 .0046335 .0061327

Request you to kindly look at the discrepancies. ]]>