Adjusted model

Model 2; from model_building3:

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.

Table 3: e/evated/Stage 1/ Stage2

Model 1: Crude

Women

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.

Men

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.

Overall

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.

Model 2: adjusting for confounders

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.

Women

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.

Men

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.

Overall

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.

Model 3: fully adjusted (including mediators)

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     =       68,614
Number of PSUs     =       214                Population size   =  296,797,030
                                              Subpop. no. obs   =       19,444
                                              Subpop. size      =  131,776,438
                                              Design df         =          109
                                              F(  21,     89)   =        88.29
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
     _IMJ2_1 |   1.205777   .1003403     2.25   0.027      1.02244    1.421988
     _IMJ2_2 |   1.303003   .1732425     1.99   0.049     1.001157    1.695854
        gndr |    2.40363    .123965    17.00   0.000     2.170076    2.662321
    ridageyr |   1.026969   .0030686     8.91   0.000     1.020905    1.033069
  _IAL_cat_1 |   1.081904   .0766249     1.11   0.269     .9402134    1.244948
  _IAL_cat_2 |   1.310025   .1081971     3.27   0.001     1.112214    1.543018
      bmxbmi |   1.064217   .0049431    13.40   0.000     1.054465    1.074059
       _cons |    .011041    .001912   -26.02   0.000     .0078334    .0155621
-------------+----------------------------------------------------------------
Stage_1_HTN  |
     _IMJ2_1 |   1.171018   .1020712     1.81   0.073     .9852266    1.391846
     _IMJ2_2 |   1.074954   .1348259     0.58   0.566     .8383583    1.378321
        gndr |   2.592831   .1429113    17.29   0.000     2.324509    2.892127
    ridageyr |   1.048616   .0027539    18.08   0.000     1.043172    1.054088
  _IAL_cat_1 |    1.12873   .0774353     1.77   0.080     .9852323    1.293127
  _IAL_cat_2 |   1.300216   .0995582     3.43   0.001     1.117139    1.513297
      bmxbmi |   1.077195   .0042692    18.76   0.000     1.068766    1.085689
       _cons |   .0039603    .000705   -31.07   0.000     .0027829     .005636
-------------+----------------------------------------------------------------
Stage_2_HTN  |
     _IMJ2_1 |   1.244438   .1362472     2.00   0.048      1.00169    1.546013
     _IMJ2_2 |   1.252059   .2226305     1.26   0.209      .880181    1.781055
        gndr |   3.096627   .2266328    15.44   0.000     2.678506    3.580018
    ridageyr |   1.088977   .0032248    28.78   0.000     1.082605    1.095388
  _IAL_cat_1 |   1.293558   .0992258     3.36   0.001     1.111115    1.505957
  _IAL_cat_2 |   1.635608   .1629647     4.94   0.000     1.342509    1.992697
      bmxbmi |   1.104262   .0056514    19.38   0.000     1.093118     1.11552
       _cons |   .0001461   .0000337   -38.31   0.000     .0000925    .0002307
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.

Women

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, 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,498
Number of PSUs     =       214                Population size   =  302,374,280
                                              Subpop. no. obs   =        9,987
                                              Subpop. size      = 66,159,736.7
                                              Design df         =          109
                                              F(  18,     92)   =        50.52
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
     _IMJ2_1 |   1.057479   .1663734     0.36   0.723     .7741917    1.444425
     _IMJ2_2 |    .940194   .2318249    -0.25   0.803     .5767389    1.532695
    ridageyr |   1.052513   .0044564    12.09   0.000     1.043718    1.061383
  _IAL_cat_1 |   1.019527   .0997199     0.20   0.844     .8398625    1.237626
  _IAL_cat_2 |   1.257752   .1824334     1.58   0.117     .9435061    1.676662
      bmxbmi |   1.070391   .0056648    12.85   0.000     1.059222    1.081677
       _cons |    .003589   .0008772   -23.03   0.000      .002211    .0058256
-------------+----------------------------------------------------------------
Stage_1_HTN  |
     _IMJ2_1 |   1.200665   .1657115     1.33   0.188     .9133205    1.578413
     _IMJ2_2 |   1.121593    .255824     0.50   0.616     .7136825    1.762647
    ridageyr |   1.070265   .0049779    14.60   0.000     1.060445    1.080177
  _IAL_cat_1 |   1.093994   .1108211     0.89   0.377     .8949945    1.337239
  _IAL_cat_2 |   1.151269    .173583     0.93   0.352     .8538783    1.552237
      bmxbmi |   1.073228   .0051134    14.83   0.000     1.063142    1.083411
       _cons |   .0019576   .0005675   -21.51   0.000     .0011021    .0034773
-------------+----------------------------------------------------------------
Stage_2_HTN  |
     _IMJ2_1 |    1.05934    .178823     0.34   0.733     .7581138    1.480255
     _IMJ2_2 |   1.288758   .3670241     0.89   0.375     .7328868    2.266238
    ridageyr |   1.120796   .0051528    24.81   0.000      1.11063    1.131055
  _IAL_cat_1 |   1.157486   .1347632     1.26   0.212     .9189666    1.457913
  _IAL_cat_2 |   1.721043   .2775155     3.37   0.001     1.250247    2.369123
      bmxbmi |   1.090813   .0075302    12.59   0.000      1.07599    1.105841
       _cons |   .0000635    .000022   -27.84   0.000     .0000319    .0001263
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.

Men

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, 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,306
Number of PSUs     =       214                Population size   =  300,967,803
                                              Subpop. no. obs   =        9,457
                                              Subpop. size      =   65,616,701
                                              Design df         =          109
                                              F(  18,     92)   =        33.84
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
     _IMJ2_1 |   1.264938   .1165922     2.55   0.012     1.053735    1.518473
     _IMJ2_2 |   1.379079   .2322852     1.91   0.059     .9876595    1.925621
    ridageyr |   1.004976   .0036403     1.37   0.173     .9977867    1.012217
  _IAL_cat_1 |   1.169623   .1080269     1.70   0.093     .9739708    1.404578
  _IAL_cat_2 |   1.355378   .1352987     3.05   0.003     1.112081    1.651902
      bmxbmi |   1.064844    .007682     8.71   0.000     1.049727    1.080179
       _cons |   .0569171   .0146403   -11.14   0.000     .0341851     .094765
-------------+----------------------------------------------------------------
Stage_1_HTN  |
     _IMJ2_1 |   1.177737   .1347072     1.43   0.155     .9388505    1.477406
     _IMJ2_2 |   1.041642   .1827003     0.23   0.817      .735774    1.474661
    ridageyr |   1.028865   .0030745     9.52   0.000     1.022789    1.034976
  _IAL_cat_1 |   1.166109   .0964313     1.86   0.066     .9898256    1.373787
  _IAL_cat_2 |   1.362106     .13169     3.20   0.002     1.124583    1.649795
      bmxbmi |   1.089893   .0066267    14.16   0.000     1.076838    1.103106
       _cons |   .0149519     .00329   -19.10   0.000      .009667    .0231259
-------------+----------------------------------------------------------------
Stage_2_HTN  |
     _IMJ2_1 |   1.359583   .1903247     2.19   0.030     1.030174    1.794325
     _IMJ2_2 |   1.229514   .2835473     0.90   0.372     .7784446    1.941955
    ridageyr |   1.065788   .0043725    15.53   0.000     1.057157     1.07449
  _IAL_cat_1 |   1.429549   .1600167     3.19   0.002     1.145118     1.78463
  _IAL_cat_2 |   1.624723   .1901672     4.15   0.000     1.288342    2.048932
      bmxbmi |   1.127731   .0087038    15.58   0.000     1.110612    1.145114
       _cons |   .0005553   .0001779   -23.40   0.000     .0002943    .0010477
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.

Overall

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1): mlogit BP_cat i.MJ2 ridageyr gndr 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     =       68,614
Number of PSUs     =       214                Population size   =  296,797,030
                                              Subpop. no. obs   =       19,444
                                              Subpop. size      =  131,776,438
                                              Design df         =          109
                                              F(  21,     89)   =        88.29
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
     _IMJ2_1 |   1.205777   .1003403     2.25   0.027      1.02244    1.421988
     _IMJ2_2 |   1.303003   .1732425     1.99   0.049     1.001157    1.695854
    ridageyr |   1.026969   .0030686     8.91   0.000     1.020905    1.033069
        gndr |    2.40363    .123965    17.00   0.000     2.170076    2.662321
  _IAL_cat_1 |   1.081904   .0766249     1.11   0.269     .9402134    1.244948
  _IAL_cat_2 |   1.310025   .1081971     3.27   0.001     1.112214    1.543018
      bmxbmi |   1.064217   .0049431    13.40   0.000     1.054465    1.074059
       _cons |    .011041    .001912   -26.02   0.000     .0078334    .0155621
-------------+----------------------------------------------------------------
Stage_1_HTN  |
     _IMJ2_1 |   1.171018   .1020712     1.81   0.073     .9852266    1.391846
     _IMJ2_2 |   1.074954   .1348259     0.58   0.566     .8383583    1.378321
    ridageyr |   1.048616   .0027539    18.08   0.000     1.043172    1.054088
        gndr |   2.592831   .1429113    17.29   0.000     2.324509    2.892127
  _IAL_cat_1 |    1.12873   .0774353     1.77   0.080     .9852323    1.293127
  _IAL_cat_2 |   1.300216   .0995582     3.43   0.001     1.117139    1.513297
      bmxbmi |   1.077195   .0042692    18.76   0.000     1.068766    1.085689
       _cons |   .0039603    .000705   -31.07   0.000     .0027829     .005636
-------------+----------------------------------------------------------------
Stage_2_HTN  |
     _IMJ2_1 |   1.244438   .1362472     2.00   0.048      1.00169    1.546013
     _IMJ2_2 |   1.252059   .2226305     1.26   0.209      .880181    1.781055
    ridageyr |   1.088977   .0032248    28.78   0.000     1.082605    1.095388
        gndr |   3.096627   .2266328    15.44   0.000     2.678506    3.580018
  _IAL_cat_1 |   1.293558   .0992258     3.36   0.001     1.111115    1.505957
  _IAL_cat_2 |   1.635608   .1629647     4.94   0.000     1.342509    1.992697
      bmxbmi |   1.104262   .0056514    19.38   0.000     1.093118     1.11552
       _cons |   .0001461   .0000337   -38.31   0.000     .0000925    .0002307
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.

Table 3: N Tabulation

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta gndr MJ2 if include==1
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

           |      Cannabis Use Collapsed
    Gender | Never/Pas  Light-Mod  Frequent  |     Total
-----------+---------------------------------+----------
    Female |     9,532        826        321 |    10,679 
      Male |     8,475      1,271        595 |    10,341 
-----------+---------------------------------+----------
     Total |    18,007      2,097        916 |    21,020 

Table 3 Abnormal BP

Model 1: Crude

Women

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): 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   =       10,679
                                              Subpop. size      = 70,330,509.5
                                              Design df         =          109
                                              F(   2,    108)   =         4.13
                                              Prob > F          =       0.0187

------------------------------------------------------------------------------
             |             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
       _cons |    .630271   .0197932   -14.70   0.000     .5922374    .6707472
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Men

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==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   =       10,341
                                              Subpop. size      = 71,193,951.1
                                              Design df         =          109
                                              F(   2,    108)   =         0.09
                                              Prob > F          =       0.9181

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |    .996259    .075491    -0.05   0.961     .8573316    1.157699
     _IMJ2_2 |   .9545312   .1077371    -0.41   0.681     .7631979    1.193832
       _cons |   1.376755   .0486795     9.04   0.000     1.283577    1.476697
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Overall

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.

Model 2: adjusting for confounders

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.

Women

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 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   =       10,679
                                              Subpop. size      = 70,330,509.5
                                              Design df         =          109
                                              F(   3,    107)   =       210.39
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.140678   .0985242     1.52   0.130     .9612064     1.35366
     _IMJ2_2 |   1.186826   .2077724     0.98   0.330     .8388764    1.679098
    ridageyr |   1.071178   .0029664    24.83   0.000     1.065315    1.077074
       _cons |   .0367337   .0044348   -27.37   0.000     .0289165    .0466641
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Men

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 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   =       10,341
                                              Subpop. size      = 71,193,951.1
                                              Design df         =          109
                                              F(   3,    107)   =        43.27
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.138015   .0880995     1.67   0.098     .9761408    1.326733
     _IMJ2_2 |   1.092278   .1285448     0.75   0.455     .8650377    1.379214
    ridageyr |    1.02792   .0025521    11.09   0.000     1.022874     1.03299
       _cons |   .4581415   .0501147    -7.14   0.000     .3688452    .5690563
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Overall

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.

Model 3: fully adjusted (including mediators)

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     =       68,614
Number of PSUs     =       214                Population size   =  296,797,030
                                              Subpop. no. obs   =       19,444
                                              Subpop. size      =  131,776,438
                                              Design df         =          109
                                              F(   7,    103)   =       225.00
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.202959   .0733552     3.03   0.003     1.066013    1.357497
     _IMJ2_2 |   1.200407   .1243141     1.76   0.081     .9776615    1.473902
        gndr |   2.605168   .1086162    22.97   0.000     2.398549    2.829587
    ridageyr |   1.047848   .0022484    21.78   0.000     1.043401    1.052314
  _IAL_cat_1 |   1.138895   .0584702     2.53   0.013      1.02871    1.260883
  _IAL_cat_2 |    1.36091   .0891629     4.70   0.000     1.195185    1.549615
      bmxbmi |   1.077406   .0037993    21.14   0.000     1.069902    1.084963
       _cons |   .0091305   .0012481   -34.36   0.000     .0069636    .0119717
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Women

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 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     =       69,498
Number of PSUs     =       214                Population size   =  302,374,280
                                              Subpop. no. obs   =        9,987
                                              Subpop. size      = 66,159,736.7
                                              Design df         =          109
                                              F(   6,    104)   =       114.68
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.126048   .1130626     1.18   0.240     .9228504    1.373987
     _IMJ2_2 |   1.083061   .1872574     0.46   0.645     .7688314     1.52572
    ridageyr |   1.072018    .003441    21.67   0.000     1.065219    1.078859
  _IAL_cat_1 |   1.075934   .0780864     1.01   0.315     .9317854    1.242383
  _IAL_cat_2 |   1.282675   .1337174     2.39   0.019     1.043239    1.577066
      bmxbmi |   1.075201   .0042838    18.20   0.000     1.066744    1.083724
       _cons |   .0039128   .0007979   -27.18   0.000     .0026119    .0058617
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Men

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 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     =       69,306
Number of PSUs     =       214                Population size   =  300,967,803
                                              Subpop. no. obs   =        9,457
                                              Subpop. size      =   65,616,701
                                              Design df         =          109
                                              F(   6,    104)   =        64.55
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.250669    .102945     2.72   0.008     1.062409    1.472289
     _IMJ2_2 |   1.218662   .1688295     1.43   0.156     .9260544    1.603724
    ridageyr |   1.027107    .002808     9.78   0.000     1.021557    1.032688
  _IAL_cat_1 |   1.212714   .0768931     3.04   0.003     1.069501    1.375103
  _IAL_cat_2 |   1.405375   .1175936     4.07   0.000     1.190608    1.658881
      bmxbmi |   1.087982   .0057816    15.87   0.000     1.076583    1.099501
       _cons |   .0377906   .0071692   -17.27   0.000     .0259473    .0550397
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Overall

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1): logit BP_abn i.MJ2 ridageyr gndr 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     =       68,614
Number of PSUs     =       214                Population size   =  296,797,030
                                              Subpop. no. obs   =       19,444
                                              Subpop. size      =  131,776,438
                                              Design df         =          109
                                              F(   7,    103)   =       225.00
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.202959   .0733552     3.03   0.003     1.066013    1.357497
     _IMJ2_2 |   1.200407   .1243141     1.76   0.081     .9776615    1.473902
    ridageyr |   1.047848   .0022484    21.78   0.000     1.043401    1.052314
        gndr |   2.605168   .1086162    22.97   0.000     2.398549    2.829587
  _IAL_cat_1 |   1.138895   .0584702     2.53   0.013      1.02871    1.260883
  _IAL_cat_2 |    1.36091   .0891629     4.70   0.000     1.195185    1.549615
      bmxbmi |   1.077406   .0037993    21.14   0.000     1.069902    1.084963
       _cons |   .0091305   .0012481   -34.36   0.000     .0069636    .0119717
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Table 2

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta BP_cat MJ_cat if include==1,col 

svy, subpop(if include==1): ta BP_cat MJ_cat,col percent
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

            |      Cannabis Use Category
BP Category | Never Use   Past Use  Current U |     Total
------------+---------------------------------+----------
     Normal |     5,149      4,193      1,539 |    10,881 
            |     53.07      50.49      51.08 |     51.76 
------------+---------------------------------+----------
   Elevated |     1,478      1,311        555 |     3,344 
            |     15.23      15.79      18.42 |     15.91 
------------+---------------------------------+----------
Stage 1 HTN |     1,938      1,760        577 |     4,275 
            |     19.97      21.19      19.15 |     20.34 
------------+---------------------------------+----------
Stage 2 HTN |     1,138      1,040        342 |     2,520 
            |     11.73      12.52      11.35 |     11.99 
------------+---------------------------------+----------
      Total |     9,703      8,304      3,013 |    21,020 
            |    100.00     100.00     100.00 |    100.00 

(running tabulate on estimation sample)

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

--------------------------------------------------
BP        |         Cannabis Use Category         
Category  | Never Us  Past Use   Current     Total
----------+---------------------------------------
   Normal |    52.66     51.49     51.37     51.94
 Elevated |    15.63     16.04      18.8     16.27
  Stage 1 |    20.84     21.16     19.68     20.82
  Stage 2 |    10.87     11.31     10.15     10.97
          | 
    Total |      100       100       100       100
--------------------------------------------------
  Key:  column percentage

  Pearson:
    Uncorrected   chi2(6)         =   70.3954
    Design-based  F(5.56, 605.89) =    1.9542     P = 0.0758

Men

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta BP_cat MJ_cat if include==1 & gndr==1,col 

svy, subpop(if include==1 & gndr==1): ta BP_cat MJ_cat,col percent
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

            |      Cannabis Use Category
BP Category | Never Use   Past Use  Current U |     Total
------------+---------------------------------+----------
     Normal |     1,770      1,739        798 |     4,307 
            |     42.85      40.03      42.77 |     41.65 
------------+---------------------------------+----------
   Elevated |       753        809        410 |     1,972 
            |     18.23      18.62      21.97 |     19.07 
------------+---------------------------------+----------
Stage 1 HTN |     1,009      1,128        412 |     2,549 
            |     24.43      25.97      22.08 |     24.65 
------------+---------------------------------+----------
Stage 2 HTN |       599        668        246 |     1,513 
            |     14.50      15.38      13.18 |     14.63 
------------+---------------------------------+----------
      Total |     4,131      4,344      1,866 |    10,341 
            |    100.00     100.00     100.00 |    100.00 

(running tabulate on estimation sample)

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

--------------------------------------------------
BP        |         Cannabis Use Category         
Category  | Never Us  Past Use   Current     Total
----------+---------------------------------------
   Normal |    43.49     41.02     42.51     42.15
 Elevated |    18.13     18.98     22.83     19.37
  Stage 1 |     24.7     26.02     22.55     24.94
  Stage 2 |    13.68     13.98      12.1     13.54
          | 
    Total |      100       100       100       100
--------------------------------------------------
  Key:  column percentage

  Pearson:
    Uncorrected   chi2(6)         =  189.5565
    Design-based  F(5.54, 604.35) =    2.4934     P = 0.0251

Women

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta BP_cat MJ_cat if include==1 & gndr==0,col 

svy, subpop(if include==1 & gndr==0): ta BP_cat MJ_cat,col percent
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

            |      Cannabis Use Category
BP Category | Never Use   Past Use  Current U |     Total
------------+---------------------------------+----------
     Normal |     3,379      2,454        741 |     6,574 
            |     60.64      61.97      64.60 |     61.56 
------------+---------------------------------+----------
   Elevated |       725        502        145 |     1,372 
            |     13.01      12.68      12.64 |     12.85 
------------+---------------------------------+----------
Stage 1 HTN |       929        632        165 |     1,726 
            |     16.67      15.96      14.39 |     16.16 
------------+---------------------------------+----------
Stage 2 HTN |       539        372         96 |     1,007 
            |      9.67       9.39       8.37 |      9.43 
------------+---------------------------------+----------
      Total |     5,572      3,960      1,147 |    10,679 
            |    100.00     100.00     100.00 |    100.00 

(running tabulate on estimation sample)

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

--------------------------------------------------
BP        |         Cannabis Use Category         
Category  | Never Us  Past Use   Current     Total
----------+---------------------------------------
   Normal |    59.89     62.82     66.03     61.85
 Elevated |    13.66     12.85     12.12     13.13
  Stage 1 |    17.79     15.91     14.93     16.65
  Stage 2 |    8.654     8.417     6.927     8.361
          | 
    Total |      100       100       100       100
--------------------------------------------------
  Key:  column percentage

  Pearson:
    Uncorrected   chi2(6)         =  131.1846
    Design-based  F(5.63, 613.53) =    2.1713     P = 0.0481

Alternate Table 2 based on collapsed MJ variable

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta BP_cat MJ2 if include==1,col 

svy, subpop(if include==1): ta BP_cat MJ2,col percent
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

            |      Cannabis Use Collapsed
BP Category | Never/Pas  Light-Mod  Frequent  |     Total
------------+---------------------------------+----------
     Normal |     9,342      1,067        472 |    10,881 
            |     51.88      50.88      51.53 |     51.76 
------------+---------------------------------+----------
   Elevated |     2,789        369        186 |     3,344 
            |     15.49      17.60      20.31 |     15.91 
------------+---------------------------------+----------
Stage 1 HTN |     3,698        414        163 |     4,275 
            |     20.54      19.74      17.79 |     20.34 
------------+---------------------------------+----------
Stage 2 HTN |     2,178        247         95 |     2,520 
            |     12.10      11.78      10.37 |     11.99 
------------+---------------------------------+----------
      Total |    18,007      2,097        916 |    21,020 
            |    100.00     100.00     100.00 |    100.00 

(running tabulate on estimation sample)

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

--------------------------------------------------
BP        |         Cannabis Use Collapsed        
Category  | Never/Pa  Light-Mo  Frequent     Total
----------+---------------------------------------
   Normal |    52.04     51.82     50.38     51.94
 Elevated |    15.85     17.81     20.99     16.27
  Stage 1 |    21.01      20.1     18.74     20.82
  Stage 2 |     11.1     10.27     9.887     10.97
          | 
    Total |      100       100       100       100
--------------------------------------------------
  Key:  column percentage

  Pearson:
    Uncorrected   chi2(6)         =   78.0180
    Design-based  F(5.64, 614.81) =    2.3247     P = 0.0348

Men

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta BP_cat MJ2 if include==1 & gndr==1,col 

svy, subpop(if include==1 & gndr==1): ta BP_cat MJ2,col percent
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

            |      Cannabis Use Collapsed
BP Category | Never/Pas  Light-Mod  Frequent  |     Total
------------+---------------------------------+----------
     Normal |     3,509        535        263 |     4,307 
            |     41.40      42.09      44.20 |     41.65 
------------+---------------------------------+----------
   Elevated |     1,562        264        146 |     1,972 
            |     18.43      20.77      24.54 |     19.07 
------------+---------------------------------+----------
Stage 1 HTN |     2,137        294        118 |     2,549 
            |     25.22      23.13      19.83 |     24.65 
------------+---------------------------------+----------
Stage 2 HTN |     1,267        178         68 |     1,513 
            |     14.95      14.00      11.43 |     14.63 
------------+---------------------------------+----------
      Total |     8,475      1,271        595 |    10,341 
            |    100.00     100.00     100.00 |    100.00 

(running tabulate on estimation sample)

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

--------------------------------------------------
BP        |         Cannabis Use Collapsed        
Category  | Never/Pa  Light-Mo  Frequent     Total
----------+---------------------------------------
   Normal |    42.07     42.17     43.21     42.15
 Elevated |    18.62     21.81     24.93     19.37
  Stage 1 |    25.46     23.45     20.71     24.94
  Stage 2 |    13.85     12.57     11.15     13.54
          | 
    Total |      100       100       100       100
--------------------------------------------------
  Key:  column percentage

  Pearson:
    Uncorrected   chi2(6)         =  180.2785
    Design-based  F(5.49, 598.30) =    2.4796     P = 0.0263

Women

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta BP_cat MJ2 if include==1 & gndr==0,col 

svy, subpop(if include==1 & gndr==0): ta BP_cat MJ2,col percent
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

            |      Cannabis Use Collapsed
BP Category | Never/Pas  Light-Mod  Frequent  |     Total
------------+---------------------------------+----------
     Normal |     5,833        532        209 |     6,574 
            |     61.19      64.41      65.11 |     61.56 
------------+---------------------------------+----------
   Elevated |     1,227        105         40 |     1,372 
            |     12.87      12.71      12.46 |     12.85 
------------+---------------------------------+----------
Stage 1 HTN |     1,561        120         45 |     1,726 
            |     16.38      14.53      14.02 |     16.16 
------------+---------------------------------+----------
Stage 2 HTN |       911         69         27 |     1,007 
            |      9.56       8.35       8.41 |      9.43 
------------+---------------------------------+----------
      Total |     9,532        826        321 |    10,679 
            |    100.00     100.00     100.00 |    100.00 

(running tabulate on estimation sample)

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

--------------------------------------------------
BP        |         Cannabis Use Collapsed        
Category  | Never/Pa  Light-Mo  Frequent     Total
----------+---------------------------------------
   Normal |    61.34     66.66     64.39     61.85
 Elevated |    13.26     11.67     13.28     13.13
  Stage 1 |    16.86     14.94      14.9     16.65
  Stage 2 |    8.537     6.735     7.424     8.361
          | 
    Total |      100       100       100       100
--------------------------------------------------
  Key:  column percentage

  Pearson:
    Uncorrected   chi2(6)         =   72.8007
    Design-based  F(5.53, 602.54) =    1.1777     P = 0.3176

Table 4

N Tabulation

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

ta gndr MJ2 if include==1,col 

svy, subpop(if include==1): ta gndr MJ2,col percent
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>


+-------------------+
| Key               |
|-------------------|
|     frequency     |
| column percentage |
+-------------------+

           |      Cannabis Use Collapsed
    Gender | Never/Pas  Light-Mod  Frequent  |     Total
-----------+---------------------------------+----------
    Female |     9,532        826        321 |    10,679 
           |     52.93      39.39      35.04 |     50.80 
-----------+---------------------------------+----------
      Male |     8,475      1,271        595 |    10,341 
           |     47.07      60.61      64.96 |     49.20 
-----------+---------------------------------+----------
     Total |    18,007      2,097        916 |    21,020 
           |    100.00     100.00     100.00 |    100.00 

(running tabulate on estimation sample)

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

--------------------------------------------------
          |         Cannabis Use Collapsed        
   Gender | Never/Pa  Light-Mo  Frequent     Total
----------+---------------------------------------
   Female |    51.72     39.41     33.86     49.69
     Male |    48.28     60.59     66.14     50.31
          | 
    Total |      100       100       100       100
--------------------------------------------------
  Key:  column percentage

  Pearson:
    Uncorrected   chi2(2)         =  708.2258
    Design-based  F(1.87, 204.04) =   80.5623     P = 0.0000

M1

Overall

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.

Male

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==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   =       10,341
                                              Subpop. size      = 71,193,951.1
                                              Design df         =          109
                                              F(   2,    108)   =         0.09
                                              Prob > F          =       0.9181

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |    .996259    .075491    -0.05   0.961     .8573316    1.157699
     _IMJ2_2 |   .9545312   .1077371    -0.41   0.681     .7631979    1.193832
       _cons |   1.376755   .0486795     9.04   0.000     1.283577    1.476697
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Female

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): 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   =       10,679
                                              Subpop. size      = 70,330,509.5
                                              Design df         =          109
                                              F(   2,    108)   =         4.13
                                              Prob > F          =       0.0187

------------------------------------------------------------------------------
             |             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
       _cons |    .630271   .0197932   -14.70   0.000     .5922374    .6707472
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

M2

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.

Male

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 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   =       10,341
                                              Subpop. size      = 71,193,951.1
                                              Design df         =          109
                                              F(   3,    107)   =        43.27
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.138015   .0880995     1.67   0.098     .9761408    1.326733
     _IMJ2_2 |   1.092278   .1285448     0.75   0.455     .8650377    1.379214
    ridageyr |    1.02792   .0025521    11.09   0.000     1.022874     1.03299
       _cons |   .4581415   .0501147    -7.14   0.000     .3688452    .5690563
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Female

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 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   =       10,679
                                              Subpop. size      = 70,330,509.5
                                              Design df         =          109
                                              F(   3,    107)   =       210.39
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.140678   .0985242     1.52   0.130     .9612064     1.35366
     _IMJ2_2 |   1.186826   .2077724     0.98   0.330     .8388764    1.679098
    ridageyr |   1.071178   .0029664    24.83   0.000     1.065315    1.077074
       _cons |   .0367337   .0044348   -27.37   0.000     .0289165    .0466641
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

M3

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     =       68,614
Number of PSUs     =       214                Population size   =  296,797,030
                                              Subpop. no. obs   =       19,444
                                              Subpop. size      =  131,776,438
                                              Design df         =          109
                                              F(   7,    103)   =       225.00
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.202959   .0733552     3.03   0.003     1.066013    1.357497
     _IMJ2_2 |   1.200407   .1243141     1.76   0.081     .9776615    1.473902
        gndr |   2.605168   .1086162    22.97   0.000     2.398549    2.829587
    ridageyr |   1.047848   .0022484    21.78   0.000     1.043401    1.052314
  _IAL_cat_1 |   1.138895   .0584702     2.53   0.013      1.02871    1.260883
  _IAL_cat_2 |    1.36091   .0891629     4.70   0.000     1.195185    1.549615
      bmxbmi |   1.077406   .0037993    21.14   0.000     1.069902    1.084963
       _cons |   .0091305   .0012481   -34.36   0.000     .0069636    .0119717
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Male

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 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     =       69,306
Number of PSUs     =       214                Population size   =  300,967,803
                                              Subpop. no. obs   =        9,457
                                              Subpop. size      =   65,616,701
                                              Design df         =          109
                                              F(   6,    104)   =        64.55
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.250669    .102945     2.72   0.008     1.062409    1.472289
     _IMJ2_2 |   1.218662   .1688295     1.43   0.156     .9260544    1.603724
    ridageyr |   1.027107    .002808     9.78   0.000     1.021557    1.032688
  _IAL_cat_1 |   1.212714   .0768931     3.04   0.003     1.069501    1.375103
  _IAL_cat_2 |   1.405375   .1175936     4.07   0.000     1.190608    1.658881
      bmxbmi |   1.087982   .0057816    15.87   0.000     1.076583    1.099501
       _cons |   .0377906   .0071692   -17.27   0.000     .0259473    .0550397
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Female

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 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     =       69,498
Number of PSUs     =       214                Population size   =  302,374,280
                                              Subpop. no. obs   =        9,987
                                              Subpop. size      = 66,159,736.7
                                              Design df         =          109
                                              F(   6,    104)   =       114.68
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.126048   .1130626     1.18   0.240     .9228504    1.373987
     _IMJ2_2 |   1.083061   .1872574     0.46   0.645     .7688314     1.52572
    ridageyr |   1.072018    .003441    21.67   0.000     1.065219    1.078859
  _IAL_cat_1 |   1.075934   .0780864     1.01   0.315     .9317854    1.242383
  _IAL_cat_2 |   1.282675   .1337174     2.39   0.019     1.043239    1.577066
      bmxbmi |   1.075201   .0042838    18.20   0.000     1.066744    1.083724
       _cons |   .0039128   .0007979   -27.18   0.000     .0026119    .0058617
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

T5

Adjusted OR of abnormal blood pressure according to cannabis use: Effect of Adjustment for possible Mediators

Overall

Base Model

Model 2

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.

BMI

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 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)
(running logit on estimation sample)

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       70,080
Number of PSUs     =       214                Population size   =  305,959,182
                                              Subpop. no. obs   =       20,910
                                              Subpop. size      =  140,938,589
                                              Design df         =          109
                                              F(   5,    105)   =       297.74
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.286639    .078038     4.16   0.000     1.140905    1.450988
     _IMJ2_2 |   1.284134   .1277675     2.51   0.013     1.054308    1.564059
        gndr |   2.515776   .0967975    23.98   0.000     2.331059     2.71513
    ridageyr |    1.04516   .0021152    21.83   0.000     1.040976     1.04936
      bmxbmi |   1.077584    .003698    21.77   0.000      1.07028    1.084938
       _cons |   .0110468   .0015006   -33.17   0.000     .0084394    .0144596
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

HEI 2015

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 hei2015, 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     =       68,641
Number of PSUs     =       214                Population size   =  297,338,173
                                              Subpop. no. obs   =       19,471
                                              Subpop. size      =  132,317,580
                                              Design df         =          109
                                              F(   5,    105)   =       229.38
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.186144     .07235     2.80   0.006     1.051077    1.338567
     _IMJ2_2 |   1.157676   .1254377     1.35   0.179      .933945    1.435003
        gndr |   2.272486   .0847978    22.00   0.000     2.110484    2.446924
    ridageyr |   1.050607   .0022368    23.19   0.000     1.046183    1.055049
     hei2015 |    .989016   .0014472    -7.55   0.000     .9861519    .9918884
       _cons |   .1484859     .01377   -20.57   0.000     .1235555    .1784466
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Alcohol Use

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 AL_cat, 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     =       68,710
Number of PSUs     =       214                Population size   =  297,316,373
                                              Subpop. no. obs   =       19,540
                                              Subpop. size      =  132,295,780
                                              Design df         =          109
                                              F(   5,    105)   =       241.11
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |    1.07985   .0632877     1.31   0.193     .9614271     1.21286
     _IMJ2_2 |   1.077968   .1064205     0.76   0.449     .8863986    1.310939
        gndr |   2.461356   .0998585    22.20   0.000     2.271187    2.667447
    ridageyr |   1.050699   .0022082    23.53   0.000     1.046332    1.055085
      AL_cat |   1.178286   .0386243     5.00   0.000     1.104167    1.257379
       _cons |   .0735804   .0066442   -28.90   0.000      .061523    .0880009
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Tobacco use

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 SMK_cat, 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,139
Number of PSUs     =       214                Population size   =  306,306,164
                                              Subpop. no. obs   =       20,969
                                              Subpop. size      =  141,285,572
                                              Design df         =          109
                                              F(   5,    105)   =       254.28
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.127442   .0664951     2.03   0.044     1.003062    1.267245
     _IMJ2_2 |   1.110426   .1061853     1.10   0.276     .9187111    1.342147
        gndr |   2.361468   .0887486    22.86   0.000     2.191963    2.544082
    ridageyr |   1.047309   .0020738    23.34   0.000     1.043206    1.051427
     SMK_cat |    1.03314   .0195549     1.72   0.088     .9951009    1.072634
       _cons |   .0916514    .007393   -29.63   0.000       .07811    .1075404
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Male

Base Model

Model 2

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==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)

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   =       10,341
                                              Subpop. size      = 71,193,951.1
                                              Design df         =          109
                                              F(   3,    107)   =        43.27
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.138015   .0880995     1.67   0.098     .9761408    1.326733
     _IMJ2_2 |   1.092278   .1285448     0.75   0.455     .8650377    1.379214
        gndr |          1  (omitted)
    ridageyr |    1.02792   .0025521    11.09   0.000     1.022874     1.03299
       _cons |   .4581415   .0501147    -7.14   0.000     .3688452    .5690563
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

BMI

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 gndr ridageyr 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)
(running logit on estimation sample)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       70,145
Number of PSUs     =       214                Population size   =  306,263,641
                                              Subpop. no. obs   =       10,296
                                              Subpop. size      = 70,912,538.8
                                              Design df         =          109
                                              F(   4,    106)   =        82.74
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.325917   .1061475     3.52   0.001     1.131377    1.553907
     _IMJ2_2 |   1.277555   .1669605     1.87   0.064     .9860281    1.655274
        gndr |          1  (omitted)
    ridageyr |   1.023479   .0026718     8.89   0.000     1.018197    1.028788
      bmxbmi |   1.088338   .0060264    15.29   0.000      1.07646    1.100348
       _cons |   .0471976   .0090845   -15.86   0.000     .0322288    .0691189
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

HEI 2015

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 gndr ridageyr hei2015, 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)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       69,287
Number of PSUs     =       214                Population size   =  301,034,783
                                              Subpop. no. obs   =        9,438
                                              Subpop. size      = 65,683,681.8
                                              Design df         =          109
                                              F(   4,    106)   =        30.58
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.163569   .0945698     1.86   0.065     .9904523    1.366944
     _IMJ2_2 |   1.093711   .1408378     0.70   0.488     .8473492    1.411702
        gndr |          1  (omitted)
    ridageyr |   1.029383   .0027225    10.95   0.000     1.024001    1.034793
     hei2015 |   .9912399    .002232    -3.91   0.000      .986826    .9956736
       _cons |   .6674628    .093011    -2.90   0.004     .5063845    .8797791
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Alcohol Use

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 gndr ridageyr AL_cat, 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)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       69,343
Number of PSUs     =       214                Population size   =  301,203,657
                                              Subpop. no. obs   =        9,494
                                              Subpop. size      = 65,852,555.7
                                              Design df         =          109
                                              F(   4,    106)   =        40.17
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.071611   .0868686     0.85   0.395     .9125588    1.258384
     _IMJ2_2 |    1.04393   .1315131     0.34   0.734     .8132688    1.340012
        gndr |          1  (omitted)
    ridageyr |   1.031874   .0026697    12.13   0.000     1.026596    1.037178
      AL_cat |   1.218727   .0489625     4.92   0.000     1.125448    1.319737
       _cons |   .3537031   .0419791    -8.76   0.000     .2795634    .4475047
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Tobacco use

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==1): logit BP_abn i.MJ2 gndr ridageyr SMK_cat, 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)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       70,160
Number of PSUs     =       214                Population size   =  306,422,650
                                              Subpop. no. obs   =       10,311
                                              Subpop. size      = 71,071,548.1
                                              Design df         =          109
                                              F(   4,    106)   =        32.53
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.129273   .0915413     1.50   0.137     .9616653    1.326092
     _IMJ2_2 |   1.075755    .129635     0.61   0.546     .8472023    1.365966
        gndr |          1  (omitted)
    ridageyr |   1.027758   .0025881    10.87   0.000     1.022641      1.0329
     SMK_cat |   1.011251   .0261313     0.43   0.666     .9607632    1.064391
       _cons |    .457494    .049951    -7.16   0.000     .3684721    .5680234
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Female

Base Model

Model 2

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): 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)

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   =       10,679
                                              Subpop. size      = 70,330,509.5
                                              Design df         =          109
                                              F(   3,    107)   =       210.39
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.140678   .0985242     1.52   0.130     .9612064     1.35366
     _IMJ2_2 |   1.186826   .2077724     0.98   0.330     .8388764    1.679098
        gndr |          1  (omitted)
    ridageyr |   1.071178   .0029664    24.83   0.000     1.065315    1.077074
       _cons |   .0367337   .0044348   -27.37   0.000     .0289165    .0466641
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

BMI

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 gndr ridageyr 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)
(running logit on estimation sample)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       70,125
Number of PSUs     =       214                Population size   =  306,240,594
                                              Subpop. no. obs   =       10,614
                                              Subpop. size      = 70,026,050.6
                                              Design df         =          109
                                              F(   4,    106)   =       185.85
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.215196   .1148793     2.06   0.042     1.007567    1.465611
     _IMJ2_2 |   1.226002   .2101288     1.19   0.237     .8728968    1.721946
        gndr |          1  (omitted)
    ridageyr |   1.071233   .0032017    23.02   0.000     1.064906    1.077598
      bmxbmi |   1.075628    .004072    19.26   0.000     1.067588    1.083729
       _cons |   .0042092   .0008249   -27.92   0.000     .0028544    .0062069
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

HEI 2015

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 gndr ridageyr hei2015, 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)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       69,544
Number of PSUs     =       214                Population size   =  302,848,442
                                              Subpop. no. obs   =       10,033
                                              Subpop. size      = 66,633,898.5
                                              Design df         =          109
                                              F(   4,    106)   =       150.71
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.175673   .1093945     1.74   0.085     .9776751     1.41377
     _IMJ2_2 |   1.200887   .2273775     0.97   0.336     .8251354    1.747749
        gndr |          1  (omitted)
    ridageyr |   1.075446   .0031842    24.57   0.000     1.069153    1.081775
     hei2015 |   .9864173   .0018019    -7.49   0.000     .9828525    .9899951
       _cons |   .0646525   .0093425   -18.95   0.000     .0485516    .0860929
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Alcohol Use

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 gndr ridageyr AL_cat, 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)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       69,557
Number of PSUs     =       214                Population size   =  302,657,768
                                              Subpop. no. obs   =       10,046
                                              Subpop. size      = 66,443,224.6
                                              Design df         =          109
                                              F(   4,    106)   =       137.09
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.065488   .0981535     0.69   0.493     .8876771    1.278917
     _IMJ2_2 |   1.044505   .1848386     0.25   0.806     .7355094    1.483312
        gndr |          1  (omitted)
    ridageyr |   1.072271   .0032528    23.00   0.000     1.065843    1.078737
      AL_cat |   1.117404   .0551134     2.25   0.026      1.01334    1.232154
       _cons |   .0328599    .004608   -24.36   0.000     .0248862    .0433883
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.

Tobacco use

use "data\NHANES0518_new.dta", clear
svyset sdmvpsu [pw=wtmec12yr], strata(sdmvstra)

xi: svy,subpop(if include==1 & gndr==0): logit BP_abn i.MJ2 gndr ridageyr SMK_cat, 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)

note: gndr omitted because of collinearity

Survey: Logistic regression

Number of strata   =       105                Number of obs     =       70,169
Number of PSUs     =       214                Population size   =  306,428,567
                                              Subpop. no. obs   =       10,658
                                              Subpop. size      = 70,214,023.5
                                              Design df         =          109
                                              F(   4,    106)   =       158.67
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_abn | Odds Ratio   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
     _IMJ2_1 |   1.086909   .0954157     0.95   0.345     .9133354    1.293468
     _IMJ2_2 |   1.095308   .1909757     0.52   0.603     .7752755    1.547448
        gndr |          1  (omitted)
    ridageyr |   1.070537    .002987    24.43   0.000     1.064633    1.076473
     SMK_cat |    1.06419   .0250602     2.64   0.009     1.015663    1.115036
       _cons |   .0361999   .0043938   -27.34   0.000     .0284598    .0460451
------------------------------------------------------------------------------
Note: _cons estimates baseline odds.