Adjusted model

Model 2; from model_building2:

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

xi: svy,subpop(if include==1): mlogit BP_cat i.MJ gndr  ridageyr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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(  18,     92)   =        93.07
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |  -.0227526   .0547765    -0.42   0.679    -.1313178    .0858126
      _IMJ_2 |   .0797134   .0935508     0.85   0.396    -.1057013     .265128
      _IMJ_3 |   .4072004    .181229     2.25   0.027     .0480105    .7663903
      _IMJ_4 |   .2865253   .1303029     2.20   0.030     .0282692    .5447813
        gndr |   .7982729   .0478785    16.67   0.000     .7033792    .8931666
    ridageyr |   .0265302   .0027076     9.80   0.000     .0211638    .0318966
       _cons |  -2.599187   .1109633   -23.42   0.000    -2.819113   -2.379262
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |  -.0499135   .0560524    -0.89   0.375    -.1610075    .0611805
      _IMJ_2 |   .0921278   .0995805     0.93   0.357    -.1052375    .2894932
      _IMJ_3 |  -.0401177   .1842881    -0.22   0.828    -.4053708    .3251354
      _IMJ_4 |  -.0355165    .119337    -0.30   0.767    -.2720385    .2010055
        gndr |   .8633234   .0488686    17.67   0.000     .7664675    .9601793
    ridageyr |   .0470443   .0023972    19.62   0.000     .0422931    .0517955
       _cons |  -3.200029   .1069927   -29.91   0.000    -3.412085   -2.987973
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |  -.0548756   .0730604    -0.75   0.454    -.1996788    .0899277
      _IMJ_2 |   .1690277   .1130569     1.50   0.138    -.0550474    .3931028
      _IMJ_3 |   .1638104   .2305652     0.71   0.479    -.2931623    .6207831
      _IMJ_4 |   .0852999   .1726447     0.49   0.622    -.2568764    .4274762
        gndr |   .9855432    .067661    14.57   0.000     .8514412    1.119645
    ridageyr |   .0818487   .0025662    31.90   0.000     .0767627    .0869348
       _cons |  -5.475991   .1225016   -44.70   0.000    -5.718785   -5.233197
------------------------------------------------------------------------------

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.MJ, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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(  12,     98)   =         1.74
                                              Prob > F          =       0.0693

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   .8973761   .0661503    -1.47   0.145     .7753959    1.038545
      _IMJ_2 |   .6995585   .1080518    -2.31   0.023      .515079    .9501107
      _IMJ_3 |   1.155806   .3437167     0.49   0.627     .6410777    2.083816
      _IMJ_4 |    .904521   .2284574    -0.40   0.692     .5482956    1.492185
       _cons |   .2280199   .0109087   -30.90   0.000     .2073927    .2506987
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |   .8525828   .0682291    -1.99   0.049      .727534     .999125
      _IMJ_2 |   .7739365   .1111255    -1.78   0.077     .5822556     1.02872
      _IMJ_3 |   .6433332   .2198738    -1.29   0.200     .3267777    1.266542
      _IMJ_4 |   .7788413    .169243    -1.15   0.253     .5062971    1.198098
       _cons |   .2971067   .0156086   -23.10   0.000      .267727    .3297104
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   .9273116   .0817424    -0.86   0.394     .7786639    1.104336
      _IMJ_2 |   .6920029   .1041564    -2.45   0.016     .5135131    .9325333
      _IMJ_3 |   .7410859   .2980164    -0.75   0.458     .3339858    1.644406
      _IMJ_4 |   .7979358   .2167948    -0.83   0.408     .4656985    1.367197
       _cons |   .1444932   .0100806   -27.73   0.000     .1258334    .1659199
------------------------------------------------------------------------------
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.MJ, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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(  12,     98)   =         1.47
                                              Prob > F          =       0.1491

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   1.109859    .089886     1.29   0.201     .9452706    1.303105
      _IMJ_2 |   1.182803   .1327996     1.50   0.138     .9468269    1.477591
      _IMJ_3 |   1.457371   .3093195     1.77   0.079     .9569272    2.219533
      _IMJ_4 |   1.383918   .2196616     2.05   0.043     1.010386    1.895543
       _cons |   .4168928   .0228766   -15.94   0.000     .3739307    .4647909
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |   1.116623   .0899385     1.37   0.174     .9518685    1.309895
      _IMJ_2 |   1.016525   .1244252     0.13   0.894     .7975525    1.295618
      _IMJ_3 |   .8393784   .1818391    -0.81   0.421     .5463704    1.289521
      _IMJ_4 |   .8436662   .1414778    -1.01   0.313     .6050993    1.176291
       _cons |   .5680377   .0380708    -8.44   0.000     .4973795    .6487338
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   1.084126   .0962073     0.91   0.365     .9092737    1.292602
      _IMJ_2 |   .9773222   .1430684    -0.16   0.876     .7311947    1.306299
      _IMJ_3 |   .8386504   .2162303    -0.68   0.496     .5030959    1.398013
      _IMJ_4 |   .8203393   .1609687    -1.01   0.315     .5560244      1.2103
       _cons |   .3144748   .0228819   -15.90   0.000     .2722419    .3632592
------------------------------------------------------------------------------
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.MJ gndr ridageyr, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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(  18,     92)   =        93.07
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   .9775043   .0535442    -0.42   0.679      .876939    1.089602
      _IMJ_2 |   1.082977   .1013133     0.85   0.396     .8996934    1.303598
      _IMJ_3 |   1.502605   .2723156     2.25   0.027     1.049182    2.151984
      _IMJ_4 |   1.331792   .1735363     2.20   0.030     1.028673    1.724231
        gndr |   2.221701   .1063717    16.67   0.000     2.020569    2.442853
    ridageyr |   1.026885   .0027804     9.80   0.000     1.021389    1.032411
       _cons |    .074334   .0082483   -23.42   0.000     .0596588    .0926189
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |   .9513117   .0533233    -0.89   0.375     .8512857    1.063091
      _IMJ_2 |   1.096505   .1091905     0.93   0.357     .9001107     1.33575
      _IMJ_3 |   .9606764   .1770412    -0.22   0.828     .6667296    1.384218
      _IMJ_4 |   .9651068   .1151729    -0.30   0.767     .7618249    1.222631
        gndr |   2.371028   .1158687    17.67   0.000      2.15215    2.612165
    ridageyr |   1.048168   .0025127    19.62   0.000       1.0432     1.05316
       _cons |    .040761   .0043611   -29.91   0.000     .0329724    .0503895
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   .9466029   .0691592    -0.75   0.454     .8189937    1.094095
      _IMJ_2 |   1.184153   .1338767     1.50   0.138     .9464403    1.481571
      _IMJ_3 |   1.177991   .2716037     0.71   0.479     .7459011    1.860384
      _IMJ_4 |   1.089044   .1880177     0.49   0.622     .7734638    1.533383
        gndr |   2.679267    .181282    14.57   0.000     2.343021    3.063767
    ridageyr |   1.085292    .002785    31.90   0.000     1.079786    1.090826
       _cons |   .0041861   .0005128   -44.70   0.000     .0032837    .0053364
------------------------------------------------------------------------------
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.MJ ridageyr, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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(  15,     95)   =        65.56
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   .8921941   .0682114    -1.49   0.139     .7667459    1.038167
      _IMJ_2 |   .9350331   .1402465    -0.45   0.655     .6945783    1.258731
      _IMJ_3 |   1.434354   .4353485     1.19   0.237      .785965    2.617638
      _IMJ_4 |   1.142521   .2896618     0.53   0.600     .6912541    1.888387
    ridageyr |   1.053495   .0040032    13.71   0.000     1.045591    1.061459
       _cons |   .0279615   .0045801   -21.84   0.000       .02021     .038686
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |   .8479887   .0682946    -2.05   0.043     .7228814    .9947479
      _IMJ_2 |    1.11356   .1704099     0.70   0.484     .8222249    1.508123
      _IMJ_3 |   .8375161   .2688826    -0.55   0.582     .4432501    1.582477
      _IMJ_4 |   1.040011   .2298715     0.18   0.859     .6710989    1.611719
    ridageyr |   1.068497   .0042926    16.49   0.000     1.060023    1.077039
       _cons |   .0194879     .00372   -20.63   0.000     .0133493    .0284492
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   .9298373   .0905392    -0.75   0.457     .7666447    1.127768
      _IMJ_2 |   1.200483   .2093219     1.05   0.297     .8497089    1.696062
      _IMJ_3 |   1.081766   .4577725     0.19   0.853      .467612    2.502542
      _IMJ_4 |   1.216433   .3694505     0.65   0.520     .6662884    2.220825
    ridageyr |   1.115521   .0043328    28.15   0.000     1.106966    1.124141
       _cons |   .0012441   .0002491   -33.41   0.000     .0008366    .0018502
------------------------------------------------------------------------------
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.MJ ridageyr, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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(  15,     95)   =        25.69
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   1.104955   .0901474     1.22   0.224     .9399834    1.298881
      _IMJ_2 |   1.201143   .1349806     1.63   0.106     .9613147    1.500804
      _IMJ_3 |   1.490514   .3222008     1.85   0.068     .9711094    2.287726
      _IMJ_4 |    1.40665   .2272686     2.11   0.037     1.021212    1.937565
    ridageyr |   1.004048   .0033653     1.21   0.231     .9974003     1.01074
       _cons |   .3578418   .0496566    -7.41   0.000     .2717981    .4711246
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |   1.079369   .0883091     0.93   0.353     .9177975    1.269385
      _IMJ_2 |   1.132311   .1363701     1.03   0.304     .8918674    1.437578
      _IMJ_3 |   .9869854   .2206866    -0.06   0.953     .6336482    1.537352
      _IMJ_4 |   .9477825   .1577972    -0.32   0.748     .6813972    1.318309
    ridageyr |   1.029567   .0028833    10.40   0.000     1.023868    1.035297
       _cons |   .1804434   .0222421   -13.89   0.000     .1413321     .230378
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   1.007366   .0933502     0.08   0.937     .8383449    1.210463
      _IMJ_2 |   1.203217   .1700247     1.31   0.193     .9093097    1.592122
      _IMJ_3 |   1.156008   .3061477     0.55   0.585     .6839195    1.953965
      _IMJ_4 |   1.031802   .2246254     0.14   0.886     .6702046    1.588492
    ridageyr |    1.06312   .0037942    17.15   0.000     1.055627    1.070667
       _cons |     .02512   .0042299   -21.88   0.000     .0179919    .0350721
------------------------------------------------------------------------------
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.MJ gndr ridageyr i.AL_cat bmxbmi hei2015 SMK_cat, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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     =       66,571
Number of PSUs     =       214                Population size   =  284,502,347
                                              Subpop. no. obs   =       17,401
                                              Subpop. size      =  119,481,755
                                              Design df         =          109
                                              F(  33,     77)   =        67.02
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   .8994145   .0528021    -1.81   0.074     .8006215    1.010398
      _IMJ_2 |   1.045329   .1130524     0.41   0.683     .8436492    1.295221
      _IMJ_3 |   1.473469   .3060052     1.87   0.065     .9762966    2.223824
      _IMJ_4 |    1.22027   .1887768     1.29   0.201     .8980403    1.658119
        gndr |   2.409292   .1363836    15.53   0.000     2.153596    2.695346
    ridageyr |   1.028449   .0032319     8.93   0.000     1.022063    1.034874
  _IAL_cat_1 |   1.069678   .0856147     0.84   0.402     .9127668    1.253563
  _IAL_cat_2 |   1.252199   .1176576     2.39   0.018     1.039431    1.508519
      bmxbmi |   1.061583   .0052237    12.14   0.000      1.05128    1.071986
     hei2015 |   .9983394   .0023017    -0.72   0.473     .9937879    1.002912
     SMK_cat |   1.064639   .0322389     2.07   0.041     1.002622    1.130492
       _cons |   .0125239   .0027041   -20.29   0.000     .0081638    .0192127
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |   .8986767   .0636075    -1.51   0.134     .7810518    1.034016
      _IMJ_2 |   1.186498   .1300411     1.56   0.122     .9548326    1.474371
      _IMJ_3 |   1.024563   .2126323     0.12   0.907     .6790498    1.545879
      _IMJ_4 |   1.064499   .1648142     0.40   0.687     .7832055     1.44682
        gndr |   2.564286   .1410949    17.11   0.000     2.299349    2.859749
    ridageyr |   1.050281   .0032108    16.05   0.000     1.043937    1.056664
  _IAL_cat_1 |   1.115051   .0824399     1.47   0.144     .9630655    1.291023
  _IAL_cat_2 |   1.204248   .1049949     2.13   0.035     1.013139    1.431406
      bmxbmi |   1.076159   .0046918    16.84   0.000       1.0669    1.085498
     hei2015 |   .9945533   .0021367    -2.54   0.012     .9903275    .9987971
     SMK_cat |   1.022316   .0291191     0.77   0.440     .9662019    1.081689
       _cons |   .0052699   .0010959   -25.22   0.000     .0034898     .007958
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   .8820164   .0806483    -1.37   0.173      .735821    1.057258
      _IMJ_2 |    1.21905   .1776198     1.36   0.177     .9132845    1.627184
      _IMJ_3 |   .9440326   .2545762    -0.21   0.831     .5531813    1.611041
      _IMJ_4 |   1.044808   .2233101     0.21   0.838     .6840126    1.595912
        gndr |   3.081131   .2382477    14.55   0.000     2.643336    3.591435
    ridageyr |   1.090773   .0035634    26.60   0.000     1.083734    1.097859
  _IAL_cat_1 |   1.335062   .1097845     3.51   0.001     1.134279    1.571387
  _IAL_cat_2 |   1.666441   .1780087     4.78   0.000     1.348479    2.059377
      bmxbmi |   1.101779   .0062009    17.22   0.000     1.089558    1.114138
     hei2015 |   .9923819   .0026301    -2.89   0.005     .9871828    .9976084
     SMK_cat |   1.046439   .0379738     1.25   0.214     .9738192    1.124475
       _cons |   .0002159    .000063   -28.92   0.000      .000121     .000385
------------------------------------------------------------------------------
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.MJ gndr ridageyr i.AL_cat bmxbmi hei2015 SMK_cat, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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,570
Number of PSUs     =       214                Population size   =  296,912,292
                                              Subpop. no. obs   =        9,059
                                              Subpop. size      = 60,697,748.6
                                              Design df         =          109
                                              F(  30,     80)   =        38.00
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   .7949089   .0710841    -2.57   0.012     .6658016    .9490517
      _IMJ_2 |   .8260803   .1475029    -1.07   0.287     .5798659    1.176839
      _IMJ_3 |   1.097982    .348598     0.29   0.769     .5852125    2.060046
      _IMJ_4 |   .8437713   .2189417    -0.65   0.514     .5045166    1.411153
        gndr |          1  (omitted)
    ridageyr |   1.053451   .0048875    11.22   0.000     1.043809    1.063182
  _IAL_cat_1 |   1.037044   .1102129     0.34   0.733     .8400775    1.280192
  _IAL_cat_2 |   1.159048   .1765683     0.97   0.335     .8569872    1.567575
      bmxbmi |   1.067237   .0060512    11.48   0.000     1.055311    1.079298
     hei2015 |   .9987202   .0027401    -0.47   0.642     .9933042    1.004166
     SMK_cat |   1.128644   .0506613     2.70   0.008     1.032572    1.233655
       _cons |   .0041971     .00122   -18.83   0.000     .0023591    .0074672
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |   .7620305   .0755567    -2.74   0.007     .6260753    .9275088
      _IMJ_2 |   1.090053   .1814433     0.52   0.605     .7837393    1.516087
      _IMJ_3 |   .6329396   .2484922    -1.17   0.247     .2906891    1.378148
      _IMJ_4 |    .958582   .2533151    -0.16   0.873     .5677607    1.618428
        gndr |          1  (omitted)
    ridageyr |   1.072397   .0052039    14.40   0.000     1.062133    1.082761
  _IAL_cat_1 |   1.154032   .1197725     1.38   0.170     .9394704    1.417596
  _IAL_cat_2 |   1.136343   .1884004     0.77   0.442     .8180872    1.578407
      bmxbmi |   1.069907   .0057626    12.55   0.000     1.058546     1.08139
     hei2015 |   .9933766   .0029697    -2.22   0.028     .9875082    .9992798
     SMK_cat |   1.059519   .0464413     1.32   0.190     .9713591    1.155681
       _cons |   .0029992   .0010433   -16.70   0.000     .0015052    .0059762
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   .8640059   .1050208    -1.20   0.232     .6790327    1.099367
      _IMJ_2 |   1.095617    .232128     0.43   0.667     .7199299    1.667351
      _IMJ_3 |   .3178483   .0893315    -4.08   0.000     .1820972       .5548
      _IMJ_4 |    .912519   .3123359    -0.27   0.790     .4630447    1.798295
        gndr |          1  (omitted)
    ridageyr |   1.119427   .0055198    22.88   0.000      1.10854    1.130421
  _IAL_cat_1 |   1.190862   .1569827     1.33   0.188     .9170523    1.546425
  _IAL_cat_2 |   1.680285   .2992687     2.91   0.004      1.18053    2.391602
      bmxbmi |   1.086122   .0084733    10.59   0.000     1.069457    1.103046
     hei2015 |   .9898551   .0042158    -2.39   0.018     .9815348     .998246
     SMK_cat |   1.111808   .0676411     1.74   0.084      .985513    1.254288
       _cons |   .0001289   .0000603   -19.15   0.000      .000051    .0003257
------------------------------------------------------------------------------
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.MJ gndr ridageyr i.AL_cat bmxbmi hei2015 SMK_cat, rrr
      pweight: wtmec12yr
          VCE: linearized
  Single unit: missing
     Strata 1: sdmvstra
         SU 1: sdmvpsu
        FPC 1: <zero>

i.MJ              _IMJ_0-4            (naturally coded; _IMJ_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,191
Number of PSUs     =       214                Population size   =  294,135,108
                                              Subpop. no. obs   =        8,342
                                              Subpop. size      = 58,784,006.4
                                              Design df         =          109
                                              F(  30,     80)   =        19.69
                                              Prob > F          =       0.0000

------------------------------------------------------------------------------
             |             Linearized
      BP_cat |        RRR   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
Normal       |  (base outcome)
-------------+----------------------------------------------------------------
Elevated     |
      _IMJ_1 |   1.029215   .0954903     0.31   0.757     .8563378    1.236992
      _IMJ_2 |    1.23101   .1544341     1.66   0.100     .9600123    1.578506
      _IMJ_3 |    1.70292   .4283648     2.12   0.037     1.034363    2.803596
      _IMJ_4 |   1.420692    .278281     1.79   0.076     .9636008    2.094607
        gndr |          1  (omitted)
    ridageyr |   1.006308   .0037867     1.67   0.098     .9988309    1.013841
  _IAL_cat_1 |   1.110589   .1138247     1.02   0.308      .906429    1.360732
  _IAL_cat_2 |   1.301707   .1554497     2.21   0.029     1.027358    1.649319
      bmxbmi |   1.062824   .0085671     7.56   0.000     1.045979     1.07994
     hei2015 |   .9985276   .0033812    -0.44   0.664     .9918486    1.005252
     SMK_cat |   1.026822   .0426495     0.64   0.525     .9456777    1.114929
       _cons |    .060901   .0216853    -7.86   0.000     .0300698    .1233441
-------------+----------------------------------------------------------------
Stage_1_HTN  |
      _IMJ_1 |    1.05677   .1035199     0.56   0.574     .8702856    1.283215
      _IMJ_2 |   1.347343   .1805813     2.22   0.028      1.03303    1.757291
      _IMJ_3 |   1.290155   .3169866     1.04   0.302     .7927875    2.099554
      _IMJ_4 |   1.177875    .246408     0.78   0.436     .7780929    1.783065
        gndr |          1  (omitted)
    ridageyr |   1.029711   .0036452     8.27   0.000     1.022511    1.036961
  _IAL_cat_1 |   1.077236   .0951664     0.84   0.402     .9042089    1.283373
  _IAL_cat_2 |   1.214662   .1227036     1.93   0.057     .9942667    1.483912
      bmxbmi |   1.092266    .007082    13.61   0.000     1.078319    1.106392
     hei2015 |   .9957018    .003095    -1.39   0.169     .9895865    1.001855
     SMK_cat |   1.002658   .0383125     0.07   0.945     .9295285    1.081542
       _cons |   .0164899   .0041127   -16.46   0.000     .0100586    .0270331
-------------+----------------------------------------------------------------
Stage_2_HTN  |
      _IMJ_1 |   .9403952   .1050574    -0.55   0.583     .7536155    1.173467
      _IMJ_2 |   1.396226   .2578861     1.81   0.074     .9682153    2.013444
      _IMJ_3 |   1.324358   .4048688     0.92   0.360     .7225432    2.427433
      _IMJ_4 |   1.140741   .3121489     0.48   0.631     .6632115    1.962105
        gndr |          1  (omitted)
    ridageyr |   1.068822   .0049959    14.24   0.000     1.058966     1.07877
  _IAL_cat_1 |   1.432149   .1589488     3.24   0.002     1.149359    1.784516
  _IAL_cat_2 |   1.649894   .2017082     4.10   0.000     1.294864    2.102269
      bmxbmi |   1.128127   .0093439    14.56   0.000     1.109758    1.146799
     hei2015 |   .9942569   .0036443    -1.57   0.119     .9870601    1.001506
     SMK_cat |   1.009376   .0433372     0.22   0.828     .9270364     1.09903
       _cons |   .0006471    .000226   -21.03   0.000     .0003239    .0012928
------------------------------------------------------------------------------
Note: _cons estimates baseline relative risk for each outcome.

N Tabulation

For Table 3

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

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

           |                      Cannabis Use
    Gender | Never Use   Past Use  1-10 time  11-20 tim  21-30 tim |     Total
-----------+-------------------------------------------------------+----------
    Female |     5,572      3,960        698        128        321 |    10,679 
      Male |     4,131      4,344        993        278        595 |    10,341 
-----------+-------------------------------------------------------+----------
     Total |     9,703      8,304      1,691        406        916 |    21,020 

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