From model building 3 using the collapsed MJ exposure variable, MJ2.
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): mlogit BP_cat i.MJ2, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 6, 104) = 2.22
Prob > F = 0.0467
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | 1.12882 .0853281 1.60 0.112 .9717611 1.311263
_IMJ2_2 | 1.367861 .1567383 2.73 0.007 1.089961 1.716616
_cons | .3045146 .0101963 -35.51 0.000 .2849618 .325409
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | .9606186 .0796187 -0.48 0.629 .8150965 1.132121
_IMJ2_2 | .9212273 .1022908 -0.74 0.462 .7392485 1.148003
_cons | .4037964 .0129189 -28.34 0.000 .3789865 .4302306
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | .9290549 .0868315 -0.79 0.433 .7719571 1.118123
_IMJ2_2 | .9196586 .1445331 -0.53 0.595 .6735199 1.255749
_cons | .2133771 .009085 -36.28 0.000 .1961097 .232165
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): mlogit BP_cat i.MJ2 gndr ridageyr, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 12, 98) = 135.85
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | 1.174563 .0941209 2.01 0.047 1.002078 1.376737
_IMJ2_2 | 1.347399 .168559 2.38 0.019 1.051516 1.72654
gndr | 2.222371 .1048229 16.93 0.000 2.024031 2.440147
ridageyr | 1.026825 .0027857 9.76 0.000 1.021319 1.032361
_cons | .0735952 .0079714 -24.09 0.000 .0593769 .0912182
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.100645 .093988 1.12 0.264 .929275 1.303618
_IMJ2_2 | .991614 .1125176 -0.07 0.941 .7919053 1.241687
gndr | 2.362098 .1160804 17.49 0.000 2.142879 2.603742
ridageyr | 1.048138 .0025136 19.60 0.000 1.043168 1.053131
_cons | .0398168 .0040783 -31.47 0.000 .0325014 .0487787
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | 1.218877 .1194473 2.02 0.046 1.003708 1.480172
_IMJ2_2 | 1.122065 .1956506 0.66 0.510 .7942012 1.585278
gndr | 2.668387 .1793587 14.60 0.000 2.335565 3.048636
ridageyr | 1.085236 .0027698 32.05 0.000 1.07976 1.090739
_cons | .0040827 .0005002 -44.90 0.000 .0032025 .0052047
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): mlogit BP_cat i.MJ2 gndr ridageyr i.AL_cat bmxbmi, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.AL_cat _IAL_cat_0-2 (naturally coded; _IAL_cat_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 67,991
Number of PSUs = 214 Population size = 292,944,878
Subpop. no. obs = 18,821
Subpop. size = 127,924,286
Design df = 109
F( 21, 89) = 83.66
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | 1.18907 .0958878 2.15 0.034 1.013433 1.395147
_IMJ2_2 | 1.30479 .1737977 2.00 0.048 1.002047 1.699
gndr | 2.433773 .1273676 17.00 0.000 2.193986 2.699768
ridageyr | 1.027306 .0030228 9.16 0.000 1.021333 1.033315
_IAL_cat_1 | 1.08711 .0796698 1.14 0.257 .9401391 1.257057
_IAL_cat_2 | 1.314678 .1107666 3.25 0.002 1.112493 1.553608
bmxbmi | 1.062911 .0051213 12.66 0.000 1.052809 1.07311
_cons | .0112267 .0020332 -24.79 0.000 .007841 .0160744
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.188509 .1008407 2.04 0.044 1.004547 1.40616
_IMJ2_2 | 1.082407 .1410556 0.61 0.545 .8360256 1.401398
gndr | 2.668435 .1435724 18.24 0.000 2.398526 2.968717
ridageyr | 1.048978 .0028223 17.77 0.000 1.043399 1.054586
_IAL_cat_1 | 1.147197 .0779982 2.02 0.046 1.00257 1.312686
_IAL_cat_2 | 1.311935 .1022692 3.48 0.001 1.124122 1.531125
bmxbmi | 1.078146 .0044717 18.14 0.000 1.069319 1.087045
_cons | .0036868 .0006747 -30.61 0.000 .0025651 .0052988
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | 1.269826 .1374891 2.21 0.029 1.024582 1.573772
_IMJ2_2 | 1.192516 .2176642 0.96 0.337 .8305279 1.712277
gndr | 3.181072 .2311443 15.93 0.000 2.754412 3.673823
ridageyr | 1.088467 .0032648 28.26 0.000 1.082016 1.094957
_IAL_cat_1 | 1.348174 .1035372 3.89 0.000 1.157821 1.569822
_IAL_cat_2 | 1.691413 .1749037 5.08 0.000 1.377975 2.076146
bmxbmi | 1.103586 .0058826 18.49 0.000 1.091988 1.115306
_cons | .0001444 .0000344 -37.15 0.000 .0000901 .0002314
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1 & gndr==1): mlogit BP_cat i.MJ2, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 10,341
Subpop. size = 71,193,951.1
Design df = 109
F( 6, 104) = 2.43
Prob > F = 0.0310
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | 1.168872 .1019609 1.79 0.076 .9832934 1.389476
_IMJ2_2 | 1.303918 .1905679 1.82 0.072 .9760012 1.742008
_cons | .4424707 .0181896 -19.83 0.000 .407849 .4800313
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | .9193642 .0990083 -0.78 0.437 .7426613 1.13811
_IMJ2_2 | .792077 .1263921 -1.46 0.147 .5773191 1.086723
_cons | .6050349 .0253812 -11.98 0.000 .5567646 .6574901
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | .9055911 .1102886 -0.81 0.417 .7113832 1.152818
_IMJ2_2 | .7835272 .1528135 -1.25 0.214 .5323266 1.153268
_cons | .3292496 .0180238 -20.29 0.000 .2953966 .3669822
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1 & gndr==0): mlogit BP_cat i.MJ2, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 10,679
Subpop. size = 70,330,509.5
Design df = 109
F( 6, 104) = 1.77
Prob > F = 0.1125
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | .8097371 .1159149 -1.47 0.143 .6097131 1.075381
_IMJ2_2 | .9541053 .2362193 -0.19 0.850 .5841014 1.558491
_cons | .2161699 .0098768 -33.52 0.000 .1974545 .2366592
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | .8153086 .101456 -1.64 0.104 .6371041 1.043359
_IMJ2_2 | .8416749 .1789191 -0.81 0.419 .5522921 1.282685
_cons | .2749268 .0115081 -30.85 0.000 .2530386 .2987083
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | .7260548 .1032662 -2.25 0.026 .5477023 .9624856
_IMJ2_2 | .8284301 .2217825 -0.70 0.484 .4873252 1.408293
_cons | .1391744 .0077181 -35.56 0.000 .1246881 .1553437
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): mlogit BP_cat i.MJ2 gndr i.MJ2*gndr, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.MJ2*gndr _IMJ2Xgndr_# (coded as above)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 15, 95) = 35.73
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | .8097371 .1159149 -1.47 0.143 .6097131 1.075381
_IMJ2_2 | .9541053 .2362193 -0.19 0.850 .5841013 1.558491
gndr | 2.046866 .1064922 13.77 0.000 1.846319 2.269196
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.443521 .2268949 2.34 0.021 1.057129 1.971144
_IMJ2Xgndr_2 | 1.366641 .4060856 1.05 0.295 .7583814 2.462755
_cons | .2161699 .0098768 -33.52 0.000 .1974545 .2366592
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | .8153086 .101456 -1.64 0.104 .6371041 1.043359
_IMJ2_2 | .8416749 .1789191 -0.81 0.419 .5522921 1.282685
gndr | 2.200713 .1119424 15.51 0.000 1.989664 2.434148
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.127627 .1721399 0.79 0.433 .8332304 1.52604
_IMJ2Xgndr_2 | .9410722 .2729556 -0.21 0.835 .5296176 1.672182
_cons | .2749268 .0115081 -30.85 0.000 .2530386 .2987083
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | .7260556 .1032662 -2.25 0.026 .5477031 .9624865
_IMJ2_2 | .8284303 .2217825 -0.70 0.484 .4873253 1.408293
gndr | 2.365734 .1645127 12.38 0.000 2.061147 2.715331
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.247275 .2372418 1.16 0.248 .8555379 1.818383
_IMJ2Xgndr_2 | .9457973 .3206978 -0.16 0.870 .4829868 1.852085
_cons | .1391744 .0077181 -35.56 0.000 .1246881 .1553437
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1 & gndr==1): mlogit BP_cat i.MJ2 ridageyr, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 10,341
Subpop. size = 71,193,951.1
Design df = 109
F( 9, 101) = 43.35
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | 1.192354 .1076517 1.95 0.054 .9969918 1.425997
_IMJ2_2 | 1.329292 .1994886 1.90 0.061 .9872915 1.789763
ridageyr | 1.004148 .0033446 1.24 0.217 .9975408 1.010799
_cons | .3774005 .0521229 -7.06 0.000 .2870271 .4962289
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.05689 .1143265 0.51 0.610 .8529415 1.309604
_IMJ2_2 | .9080252 .1436316 -0.61 0.543 .6636561 1.242375
ridageyr | 1.029751 .0028605 10.55 0.000 1.024097 1.035436
_cons | .1871393 .0223335 -14.04 0.000 .1477205 .2370768
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | 1.190368 .1461572 1.42 0.159 .9332426 1.518336
_IMJ2_2 | 1.028267 .2194621 0.13 0.896 .6735895 1.5697
ridageyr | 1.06313 .0037561 17.33 0.000 1.055711 1.070601
_cons | .0252106 .0042179 -22.00 0.000 .0180956 .0351231
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1 & gndr==0): mlogit BP_cat i.MJ2 ridageyr, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 10,679
Subpop. size = 70,330,509.5
Design df = 109
F( 9, 101) = 107.17
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | 1.073181 .1508434 0.50 0.616 .812245 1.417942
_IMJ2_2 | 1.208898 .3014948 0.76 0.448 .7374284 1.981799
ridageyr | 1.053534 .003991 13.77 0.000 1.045654 1.061474
_cons | .0263858 .0042208 -22.72 0.000 .0192168 .0362293
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.157565 .1506857 1.12 0.263 .8943271 1.498284
_IMJ2_2 | 1.126787 .2449948 0.55 0.584 .7323004 1.733782
ridageyr | 1.068428 .0042994 16.45 0.000 1.059941 1.076984
_cons | .0180336 .0033011 -21.94 0.000 .0125464 .0259206
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | 1.223032 .1930029 1.28 0.205 .8945499 1.672135
_IMJ2_2 | 1.261376 .3773209 0.78 0.439 .6972083 2.282058
ridageyr | 1.115458 .004349 28.02 0.000 1.106871 1.124111
_cons | .0012028 .0002358 -34.30 0.000 .0008156 .0017738
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): mlogit BP_cat i.MJ2 gndr i.MJ2*gndr ridageyr, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.MJ2*gndr _IMJ2Xgndr_# (coded as above)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 18, 92) = 92.78
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | .9347764 .1320261 -0.48 0.634 .7065386 1.236743
_IMJ2_2 | 1.075654 .2648861 0.30 0.768 .6602461 1.752424
gndr | 2.118695 .1098636 14.48 0.000 1.911764 2.348023
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.423994 .2191215 2.30 0.024 1.049679 1.93179
_IMJ2Xgndr_2 | 1.374698 .40371 1.08 0.281 .7681161 2.460296
ridageyr | 1.026806 .0027888 9.74 0.000 1.021294 1.032348
_cons | .0755296 .0082775 -23.57 0.000 .0607834 .0938534
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.046803 .1279568 0.37 0.709 .8215793 1.333768
_IMJ2_2 | 1.033806 .2154751 0.16 0.874 .683962 1.562593
gndr | 2.340298 .1219319 16.32 0.000 2.110692 2.594881
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.103751 .1708689 0.64 0.525 .8121185 1.500109
_IMJ2Xgndr_2 | .9576629 .2742653 -0.15 0.880 .5428753 1.689372
ridageyr | 1.048124 .0025142 19.59 0.000 1.043152 1.053118
_cons | .0399635 .0040878 -31.48 0.000 .0326301 .048945
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | 1.086063 .1604401 0.56 0.577 .8103985 1.455498
_IMJ2_2 | 1.143104 .3203457 0.48 0.634 .6559421 1.992077
gndr | 2.622011 .1893981 13.34 0.000 2.272263 3.025592
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.209534 .2343342 0.98 0.328 .8238646 1.775745
_IMJ2Xgndr_2 | .9906064 .366021 -0.03 0.980 .4762755 2.060364
ridageyr | 1.08522 .0027752 31.98 0.000 1.079734 1.090734
_cons | .0041179 .0004957 -45.62 0.000 .0032438 .0052276
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.One interaction term was significant, the gender-light-moderate-use interaction term for elevated blood pressure. Other interaction terms were not significant. Women had RR of elevated BP given light-moderate MJ use of 1.07 (0.81, 1.42), whereas men had RR of elevated BP given light-moderate MJ use of 1.19 (1.00, 1.43).
use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1 & gndr==1): mlogit BP_cat i.MJ2 ridageyr i.AL_cat bmxbmi SMK_cat, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.AL_cat _IAL_cat_0-2 (naturally coded; _IAL_cat_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 69,001
Number of PSUs = 214 Population size = 299,084,506
Subpop. no. obs = 9,152
Subpop. size = 63,733,404.8
Design df = 109
F( 21, 89) = 28.17
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | 1.250992 .1221354 2.29 0.024 1.030904 1.518068
_IMJ2_2 | 1.33964 .2373208 1.65 0.102 .9429807 1.903153
ridageyr | 1.005501 .0037389 1.48 0.143 .9981181 1.012939
_IAL_cat_1 | 1.149784 .1129504 1.42 0.158 .946365 1.396927
_IAL_cat_2 | 1.320086 .1484308 2.47 0.015 1.056376 1.649628
bmxbmi | 1.063518 .0079253 8.26 0.000 1.047926 1.079342
SMK_cat | 1.029392 .0435686 0.68 0.495 .9465629 1.119469
_cons | .0574753 .0153199 -10.72 0.000 .0338882 .0974795
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.214722 .1439525 1.64 0.104 .9604429 1.536322
_IMJ2_2 | 1.06719 .1940419 0.36 0.721 .7442767 1.530203
ridageyr | 1.029906 .0031876 9.52 0.000 1.023607 1.036243
_IAL_cat_1 | 1.156478 .0991721 1.70 0.093 .9757182 1.370725
_IAL_cat_2 | 1.349532 .1306974 3.10 0.002 1.113838 1.635101
bmxbmi | 1.089682 .006738 13.89 0.000 1.076409 1.103118
SMK_cat | .9983936 .0351118 -0.05 0.964 .9311731 1.070467
_cons | .0145284 .0031696 -19.40 0.000 .0094282 .0223876
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | 1.370609 .2045816 2.11 0.037 1.01961 1.842439
_IMJ2_2 | 1.146622 .2748104 0.57 0.569 .7130544 1.843816
ridageyr | 1.066729 .0045068 15.29 0.000 1.057834 1.075698
_IAL_cat_1 | 1.452782 .1595815 3.40 0.001 1.168558 1.806138
_IAL_cat_2 | 1.642396 .196858 4.14 0.000 1.295111 2.082804
bmxbmi | 1.126564 .0089798 14.95 0.000 1.108906 1.144503
SMK_cat | 1.016077 .0401589 0.40 0.687 .9395206 1.098871
_cons | .000538 .0001772 -22.85 0.000 .00028 .0010335
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1 & gndr==0): mlogit BP_cat i.MJ2 ridageyr i.AL_cat bmxbmi SMK_cat, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.AL_cat _IAL_cat_0-2 (naturally coded; _IAL_cat_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 69,130
Number of PSUs = 214 Population size = 300,167,662
Subpop. no. obs = 9,619
Subpop. size = 63,953,118.8
Design df = 109
F( 21, 89) = 40.74
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | .9417161 .1476767 -0.38 0.703 .6901428 1.284994
_IMJ2_2 | .8640943 .2136508 -0.59 0.556 .5293408 1.410545
ridageyr | 1.051196 .0043949 11.94 0.000 1.042521 1.059943
_IAL_cat_1 | 1.003735 .0986039 0.04 0.970 .8261539 1.219487
_IAL_cat_2 | 1.173217 .1748829 1.07 0.286 .8731151 1.576468
bmxbmi | 1.069307 .0057842 12.39 0.000 1.057904 1.080833
SMK_cat | 1.093613 .0417663 2.34 0.021 1.013889 1.179606
_cons | .003741 .0009293 -22.50 0.000 .0022865 .0061208
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.158355 .1614896 1.05 0.294 .8787008 1.527012
_IMJ2_2 | 1.022408 .2400318 0.09 0.925 .6420127 1.628191
ridageyr | 1.06938 .0051039 14.05 0.000 1.059312 1.079543
_IAL_cat_1 | 1.116312 .1105469 1.11 0.269 .9173732 1.358393
_IAL_cat_2 | 1.148223 .1756233 0.90 0.368 .8479524 1.554824
bmxbmi | 1.074681 .0054192 14.28 0.000 1.063994 1.085475
SMK_cat | 1.050515 .0387598 1.34 0.184 .976436 1.130214
_cons | .0018363 .0005389 -21.46 0.000 .0010264 .0032852
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | 1.030042 .1694673 0.18 0.858 .7434285 1.427154
_IMJ2_2 | 1.057236 .3383727 0.17 0.862 .5606389 1.993705
ridageyr | 1.116043 .0052849 23.18 0.000 1.105618 1.126567
_IAL_cat_1 | 1.190287 .1504747 1.38 0.171 .926479 1.529212
_IAL_cat_2 | 1.634714 .2800388 2.87 0.005 1.164092 2.295598
bmxbmi | 1.090867 .0078134 12.14 0.000 1.075491 1.106464
SMK_cat | 1.121027 .0593363 2.16 0.033 1.009383 1.24502
_cons | .0000689 .0000242 -27.31 0.000 .0000344 .0001381
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): mlogit BP_cat i.MJ2 gndr i.MJ2*gndr ridageyr i.AL_cat bmxbmi SMK_cat, rrr pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.MJ2*gndr _IMJ2Xgndr_# (coded as above)
i.AL_cat _IAL_cat_0-2 (naturally coded; _IAL_cat_0 omitted)
(running mlogit on estimation sample)
Survey: Multinomial logistic regression
Number of strata = 105 Number of obs = 67,941
Number of PSUs = 214 Population size = 292,707,116
Subpop. no. obs = 18,771
Subpop. size = 127,686,524
Design df = 109
F( 30, 80) = 72.65
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_cat | RRR Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal | (base outcome)
-------------+----------------------------------------------------------------
Elevated |
_IMJ2_1 | .8243411 .1254131 -1.27 0.207 .6097525 1.114449
_IMJ2_2 | .7927396 .1971952 -0.93 0.353 .48419 1.297912
gndr | 2.248691 .1314721 13.86 0.000 2.002648 2.524962
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.670421 .2780228 3.08 0.003 1.201054 2.323215
_IMJ2Xgndr_2 | 1.867233 .5680013 2.05 0.042 1.021788 3.412211
ridageyr | 1.026536 .0030809 8.73 0.000 1.020448 1.03266
_IAL_cat_1 | 1.075728 .0813741 0.96 0.337 .9259551 1.249727
_IAL_cat_2 | 1.263572 .1154815 2.56 0.012 1.054224 1.514493
bmxbmi | 1.063828 .0050965 12.92 0.000 1.053775 1.073977
SMK_cat | 1.054214 .0300462 1.85 0.067 .9963139 1.115478
_cons | .0114571 .0020667 -24.78 0.000 .0080133 .0163811
-------------+----------------------------------------------------------------
Stage_1_HTN |
_IMJ2_1 | 1.044973 .1320084 0.35 0.728 .8135193 1.342276
_IMJ2_2 | .9559892 .2205987 -0.20 0.846 .6051031 1.510346
gndr | 2.590116 .1447862 17.03 0.000 2.31848 2.893578
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.239244 .2044134 1.30 0.196 .8936652 1.718458
_IMJ2Xgndr_2 | 1.199446 .3703414 0.59 0.557 .6504484 2.211813
ridageyr | 1.048655 .0028953 17.21 0.000 1.042932 1.054409
_IAL_cat_1 | 1.14289 .0792043 1.93 0.057 .9962139 1.311163
_IAL_cat_2 | 1.295854 .1039035 3.23 0.002 1.105451 1.519053
bmxbmi | 1.078423 .0044827 18.16 0.000 1.069575 1.087345
SMK_cat | 1.016224 .0256554 0.64 0.525 .9666266 1.068366
_cons | .0037255 .0006781 -30.72 0.000 .0025972 .0053439
-------------+----------------------------------------------------------------
Stage_2_HTN |
_IMJ2_1 | .9551634 .1424034 -0.31 0.759 .7108028 1.283531
_IMJ2_2 | 1.048607 .3054359 0.16 0.871 .5886991 1.867808
gndr | 3.053203 .2385314 14.29 0.000 2.615225 3.564531
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.48713 .3040776 1.94 0.055 .9916246 2.230234
_IMJ2Xgndr_2 | 1.146982 .4518924 0.35 0.728 .5253294 2.50427
ridageyr | 1.087718 .0033221 27.53 0.000 1.081154 1.094323
_IAL_cat_1 | 1.338681 .104944 3.72 0.000 1.146038 1.563705
_IAL_cat_2 | 1.634896 .1695781 4.74 0.000 1.331095 2.008036
bmxbmi | 1.104627 .0058779 18.70 0.000 1.093038 1.116338
SMK_cat | 1.050075 .0332736 1.54 0.126 .9861556 1.118137
_cons | .0001434 .0000334 -37.99 0.000 .0000904 .0002275
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): logit BP_abn i.MJ2, or pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running logit on estimation sample)
Survey: Logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 2, 108) = 0.31
Prob > F = 0.7360
------------------------------------------------------------------------------
| Linearized
BP_abn | Odds Ratio Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_IMJ2_1 | 1.008883 .0533953 0.17 0.868 .9084166 1.120461
_IMJ2_2 | 1.068427 .0898673 0.79 0.433 .9043671 1.262248
_cons | .9216882 .024326 -3.09 0.003 .8747143 .9711848
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): logit BP_abn i.MJ2 gndr ridageyr, or pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
(running logit on estimation sample)
Survey: Logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 4, 106) = 316.01
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_abn | Odds Ratio Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_IMJ2_1 | 1.154821 .0673943 2.47 0.015 1.028683 1.296426
_IMJ2_2 | 1.154975 .110962 1.50 0.137 .9547223 1.397231
gndr | 2.368786 .0882831 23.14 0.000 2.200118 2.550384
ridageyr | 1.047687 .0020572 23.73 0.000 1.043618 1.051772
_cons | .0922189 .0074624 -29.46 0.000 .0785538 .1082612
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): logit BP_abn i.MJ2 gndr ridageyr i.AL_cat bmxbmi, or pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.AL_cat _IAL_cat_0-2 (naturally coded; _IAL_cat_0 omitted)
(running logit on estimation sample)
Survey: Logistic regression
Number of strata = 105 Number of obs = 67,991
Number of PSUs = 214 Population size = 292,944,878
Subpop. no. obs = 18,821
Subpop. size = 127,924,286
Design df = 109
F( 7, 103) = 221.90
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_abn | Odds Ratio Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_IMJ2_1 | 1.20759 .0747217 3.05 0.003 1.068215 1.365149
_IMJ2_2 | 1.194591 .1253127 1.69 0.093 .9703442 1.470662
gndr | 2.662121 .1092002 23.87 0.000 2.454254 2.887593
ridageyr | 1.047906 .0022563 21.73 0.000 1.043444 1.052388
_IAL_cat_1 | 1.15799 .0614299 2.77 0.007 1.04242 1.286373
_IAL_cat_2 | 1.376358 .0929352 4.73 0.000 1.203957 1.573446
bmxbmi | 1.077091 .0039284 20.36 0.000 1.069333 1.084905
_cons | .0089558 .0012704 -33.24 0.000 .0067609 .0118633
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): logit BP_abn i.MJ2 gndr i.MJ2*gndr, or pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.MJ2*gndr _IMJ2Xgndr_# (coded as above)
(running logit on estimation sample)
note: _IMJ2_1 omitted because of collinearity
note: _IMJ2_2 omitted because of collinearity
note: gndr omitted because of collinearity
Survey: Logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 5, 105) = 102.76
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_abn | Odds Ratio Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_IMJ2_1 | .7936889 .0653237 -2.81 0.006 .6742277 .9343166
_IMJ2_2 | .8773115 .146366 -0.78 0.434 .6303034 1.221119
gndr | 2.184386 .0838953 20.34 0.000 2.024279 2.357156
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.255226 .1403023 2.03 0.044 1.005799 1.566508
_IMJ2Xgndr_2 | 1.088019 .2238085 0.41 0.683 .7237291 1.635673
_cons | .630271 .0197932 -14.70 0.000 .5922374 .6707472
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): logit BP_abn i.MJ2 gndr i.MJ2*gndr ridageyr, or pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.MJ2*gndr _IMJ2Xgndr_# (coded as above)
(running logit on estimation sample)
note: _IMJ2_1 omitted because of collinearity
note: _IMJ2_2 omitted because of collinearity
note: gndr omitted because of collinearity
Survey: Logistic regression
Number of strata = 105 Number of obs = 70,190
Number of PSUs = 214 Population size = 306,545,053
Subpop. no. obs = 21,020
Subpop. size = 141,524,461
Design df = 109
F( 6, 104) = 209.59
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_abn | Odds Ratio Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_IMJ2_1 | 1.016611 .0806805 0.21 0.836 .8686463 1.189779
_IMJ2_2 | 1.076644 .1762755 0.45 0.653 .7782926 1.489365
gndr | 2.312864 .091466 21.20 0.000 2.138504 2.501441
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.232537 .1356927 1.90 0.060 .9909169 1.533072
_IMJ2Xgndr_2 | 1.116694 .2305348 0.53 0.594 .7417129 1.68125
ridageyr | 1.047681 .0020636 23.65 0.000 1.043599 1.051779
_cons | .093324 .0075496 -29.32 0.000 .0794989 .1095534
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)
xi: svy,subpop(if include==1): logit BP_abn i.MJ2 gndr i.MJ2*gndr ridageyr i.AL_cat bmxbmi, or pweight: wtmec12yr
VCE: linearized
Single unit: missing
Strata 1: sdmvstra
SU 1: sdmvpsu
FPC 1: <zero>
i.MJ2 _IMJ2_0-2 (naturally coded; _IMJ2_0 omitted)
i.MJ2*gndr _IMJ2Xgndr_# (coded as above)
i.AL_cat _IAL_cat_0-2 (naturally coded; _IAL_cat_0 omitted)
(running logit on estimation sample)
note: _IMJ2_1 omitted because of collinearity
note: _IMJ2_2 omitted because of collinearity
note: gndr omitted because of collinearity
Survey: Logistic regression
Number of strata = 105 Number of obs = 67,991
Number of PSUs = 214 Population size = 292,944,878
Subpop. no. obs = 18,821
Subpop. size = 127,924,286
Design df = 109
F( 9, 101) = 169.03
Prob > F = 0.0000
------------------------------------------------------------------------------
| Linearized
BP_abn | Odds Ratio Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
_IMJ2_1 | .9648459 .0902933 -0.38 0.703 .8015037 1.161476
_IMJ2_2 | .9536203 .1635884 -0.28 0.782 .6787614 1.339781
gndr | 2.53853 .109732 21.55 0.000 2.3301 2.765603
_IMJ2_1 | 1 (omitted)
_IMJ2_2 | 1 (omitted)
gndr | 1 (omitted)
_IMJ2Xgndr_1 | 1.433386 .1731692 2.98 0.004 1.128168 1.821178
_IMJ2Xgndr_2 | 1.408374 .3229539 1.49 0.138 .8940011 2.218698
ridageyr | 1.04789 .0022721 21.57 0.000 1.043396 1.052403
_IAL_cat_1 | 1.162273 .061677 2.83 0.005 1.04624 1.291175
_IAL_cat_2 | 1.378395 .0932491 4.74 0.000 1.205433 1.576175
bmxbmi | 1.077433 .0039443 20.37 0.000 1.069644 1.085279
_cons | .0090766 .0012913 -33.05 0.000 .0068465 .0120333
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.