Announcement

Collapse
No announcement yet.
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • gsem mlogit: set base category

    Hi, I've estimated a latent profile analysis in Stata 17.0 with three continuous dependent variables and four latent classes. I've also used gsem to run a multinomial logit predicting the classes, with two predictors (control1 and 2). The code structure is:

    Code:
    gsem (dv1 dv2 dv3 <-,)(C <- control1 control2), lclass(C 4) vce(cluster pid)
    I would like to change the base category for the multinomial logit, but the usual factor variable notation doesn't work. I'm able to run the following, but there's no different in output, ie. base category is still class 1:

    Code:
    gsem (dv1 dv2 dv3 <-,)(b4.C <- control1 control2), lclass(C 4) vce(cluster pid)
    Wonder if anyone can help? Thank you!

  • #2
    It’s actually in the manual, although the manual is long.

    https://www.statalist.org/forums/for...94#post1554494
    Be aware that it can be very hard to answer a question without sample data. You can use the dataex command for this. Type help dataex at the command line.

    When presenting code or results, please use the code delimiters format them. Use the # button on the formatting toolbar, between the " (double quote) and <> buttons.

    Comment


    • #3
      Thank you so much Weiwen! This works for me. I will try the dataex next time I post.

      Comment


      • #4
        Hi Weiwen, using the recommended code causes the class membership to change, and also the category order to be shifted; any idea why? Thank you!

        Data:
        Code:
        * Example generated by -dataex-. For more info, type help dataex
        clear
        input float(meetingx meetingsd meetinggd) int page byte(pmale pedu) float pmarr byte(pwork metro econ health pkids) float youngest long pid
                3         0          0 73 0 1 0 0 1 50 50 4 46  11
              2.5  .8660254   .6666667 75 0 1 0 0 1 90 80 4 48  11
             3.25 1.0897248  2.3333333 77 0 1 0 0 1 80 60 4 50  11
                3         0          0 79 0 1 0 0 1 80 60 4 52  11
             2.75  .4330127  .33333325 81 0 1 0 0 1 90 50 4 54  11
                5 1.8708287          0 83 0 1 0 0 1 80 80 4 56  11
             5.25  .4330127  -.3333335 85 0 1 0 0 1 80 40 4 58  11
                3 1.4142135          3 70 0 0 1 0 1  0 30 3 41  41
         3.333333  1.885618          4 72 0 0 1 0 1 50 40 3 43  41
        2.3333333  .4714045          1 74 0 0 1 0 1 80 80 3 45  41
        1.6666666   .942809          2 76 0 0 1 0 1 60 70 3 47  41
                2         0          0 78 0 0 1 0 1 30 50 3 49  41
         2.666667  .4714045         -1 80 0 0 1 0 1 30 50 3 51  41
         6.666667   .942809         -1 82 0 0 1 0 1 70 70 3 53  41
                3 1.4142135          3 71 1 1 1 0 1 50 80 3 41  42
         3.333333  1.885618          4 73 1 1 1 0 1 50 80 3 43  42
        2.3333333  .4714045          1 75 1 1 1 0 1 60 50 3 45  42
        1.6666666   .942809          2 77 1 1 1 0 1 50 60 3 47  42
                2         0          0 79 1 1 1 0 1 30 80 3 49  42
         2.666667  .4714045         -1 81 1 1 1 0 1 30 70 3 51  42
        4.6666665 2.3570225       -2.5 83 1 1 1 0 1 60 60 3 53  42
              6.5       3.5          7 67 1 0 1 0 1 60 70 3 33  61
         6.666667  2.867442        5.5 69 1 0 1 0 1 30 60 3 35  61
         7.666667 3.2998316          7 71 1 0 1 0 1 70 50 3 37  61
                6 2.9439204        4.5 73 1 0 1 0 1 70 60 3 39  61
         6.666667  2.867442        5.5 75 1 0 1 0 1 30 30 3 41  61
                6   2.44949        4.5 77 1 0 1 0 1 30 50 3 43  61
        4.3333335   .942809         -1 79 1 0 1 0 1 70 60 3 45  61
              6.5       3.5          7 62 0 0 1 1 1 30 70 3 33  62
         6.666667  2.867442        5.5 64 0 0 1 1 1 40 50 3 35  62
         7.666667 3.2998316          7 66 0 0 1 0 1 70 60 3 37  62
                6 2.9439204        4.5 68 0 0 1 0 1 50 60 3 39  62
         7.333333  3.091206        6.5 70 0 0 1 0 1 20 50 3 41  62
                6 2.1602468        4.5 72 0 0 1 0 1 50 50 3 43  62
         5.333333  .4714045         .5 74 0 0 1 0 1 60 60 3 45  62
             4.75  3.269174          7 87 0 0 0 0 1 80 70 4 42  91
              4.5  3.201562   7.333333 89 0 0 0 0 1 60 60 4 44  91
                5  3.082207         -4 91 0 0 0 0 1 60 80 4 46  91
                7 2.1602468          . 93 0 0 0 0 1 80 80 4 48  91
             4.75  3.031089          7 95 0 0 0 0 1 50 50 4 50  91
             3.75  1.299038          3 69 0 0 0 1 1 70 70 4 42 101
              4.5  3.201562   7.333333 71 0 0 0 1 1 70 60 4 44 101
             4.75  3.031089          7 73 0 0 0 1 1 80 60 4 46 101
                3         0          . 75 0 0 0 1 1 70 70 4 48 101
                3         0          0 77 0 0 0 1 1 70 80 4 50 101
                3         0          . 79 0 0 0 1 1 50 50 4 52 101
              6.5       1.5  3.3333335 81 0 0 0 1 1 70 70 4 54 101
         6.333333  .4714045         .5 64 1 0 1 0 1 70 70 3 32 111
                4 1.4142135       -1.5 66 1 0 1 1 1 60 70 3 34 111
         5.333333 3.2998316        3.5 68 1 0 1 1 1 50 80 3 36 111
         5.333333 3.2998316        3.5 70 1 0 1 1 1 80 80 3 38 111
         6.666667   .942809          1 72 1 0 1 1 1 60 60 3 40 111
                7         0          0 61 0 0 1 0 1 60 50 3 32 112
                4 1.4142135       -1.5 63 0 0 1 0 1 60 20 3 34 112
         5.333333 3.2998316        3.5 65 0 0 1 0 1 50 40 3 36 112
         6.333333  2.867442         .5 67 0 0 1 0 1 70 50 3 38 112
         7.333333   .942809          2 69 0 0 1 0 1 70 80 3 40 112
              5.8 3.4292855 -1.1666665 78 0 0 0 0 1 60 70 5 43 141
              6.4 3.1368775  -.1666665 80 0 0 0 0 1 70 90 5 45 141
              6.6   2.87054       -1.5 82 0 0 0 0 1 50 50 5 47 141
             8.25  2.046338        1.5 71 1 0 1 1 1 70 70 4 34 201
              6.5       3.5          0 73 1 0 1 0 1 30 60 4 36 201
        end
        Do file:
        Code:
        gsem (meetingx meetingsd meetinggd <-,)(C <- page pmale pedu pmarr pwork metro econ health pkids youngest), lclass(C 4, base(1)) vce(cluster pid) from(B, skip)
        estat lcprob
        gsem (meetingx meetingsd meetinggd <-,)(C <- page pmale pedu pmarr pwork metro econ health pkids youngest), lclass(C 4, base(2)) vce(cluster pid) from(B, skip)
        estat lcprob
        Output (I have deleted the iteration records to stay within character limit):
        Code:
        . gsem (meetingx meetingsd meetinggd <-,)(C <- page pmale pedu pmarr pwork metro econ health pkids y
        > oungest), lclass(C 4, base(1)) vce(cluster pid) from(B, skip)
        
        Fitting class model:
        
        Iteration 0:   (class) log likelihood = -43781.058  
        Iteration 8:   (class) log likelihood = -6692.8873  
        
        Fitting outcome model:
        
        Iteration 0:   (outcome) log likelihood = -281792.23  
        Iteration 7:   (outcome) log likelihood = -153716.21  
        
        Refining starting values:
        
        Iteration 0:   (EM) log likelihood = -159639.55
        Iteration 20:  (EM) log likelihood = -156672.26
        note: EM algorithm reached maximum iterations.
        
        Fitting full model:
        
        Iteration 0:   log pseudolikelihood = -146317.02  (not concave)
        Iteration 8:   log pseudolikelihood = -142326.08  
        
        Generalized structural equation model                   Number of obs = 26,678
        Log pseudolikelihood = -142326.08
        
         ( 1)  [/]var(e.meetingx)#1bn.C - [/]var(e.meetingx)#4.C = 0
         ( 2)  [/]var(e.meetingx)#2.C - [/]var(e.meetingx)#4.C = 0
         ( 3)  [/]var(e.meetingx)#3.C - [/]var(e.meetingx)#4.C = 0
         ( 4)  [/]var(e.meetingsd)#1bn.C - [/]var(e.meetingsd)#4.C = 0
         ( 5)  [/]var(e.meetingsd)#2.C - [/]var(e.meetingsd)#4.C = 0
         ( 6)  [/]var(e.meetingsd)#3.C - [/]var(e.meetingsd)#4.C = 0
         ( 7)  [/]var(e.meetinggd)#1bn.C - [/]var(e.meetinggd)#4.C = 0
         ( 8)  [/]var(e.meetinggd)#2.C - [/]var(e.meetinggd)#4.C = 0
         ( 9)  [/]var(e.meetinggd)#3.C - [/]var(e.meetinggd)#4.C = 0
        
                                        (Std. err. adjusted for 5,849 clusters in pid)
        ------------------------------------------------------------------------------
                     |               Robust
                     | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -------------+----------------------------------------------------------------
        1.C          |  (base outcome)
        -------------+----------------------------------------------------------------
        2.C          |
                page |  -.0032349   .0096245    -0.34   0.737    -.0220986    .0156288
               pmale |  -.0461375   .0937949    -0.49   0.623    -.2299722    .1376972
                pedu |   .2901108   .0898664     3.23   0.001     .1139759    .4662457
               pmarr |   .3027186   .0975886     3.10   0.002     .1114485    .4939887
               pwork |   .0025648   .0822975     0.03   0.975    -.1587353    .1638649
               metro |   .4637131   .0787419     5.89   0.000     .3093818    .6180444
                econ |   .0075708   .0018843     4.02   0.000     .0038776     .011264
              health |  -.0004807   .0019266    -0.25   0.803    -.0042567    .0032953
               pkids |  -.5721241   .0501857   -11.40   0.000    -.6704863    -.473762
            youngest |  -.0526812   .0111648    -4.72   0.000    -.0745639   -.0307986
               _cons |   2.774327   .4343398     6.39   0.000     1.923037    3.625617
        -------------+----------------------------------------------------------------
        3.C          |
                page |   .0030391   .0095167     0.32   0.749    -.0156133    .0216914
               pmale |  -.0810071   .0970496    -0.83   0.404    -.2712208    .1092066
                pedu |   .3187712   .0929806     3.43   0.001     .1365325    .5010099
               pmarr |   .0811879   .0855257     0.95   0.342    -.0864394    .2488152
               pwork |  -.1229037   .0837606    -1.47   0.142    -.2870715     .041264
               metro |   .1158267   .0791939     1.46   0.144    -.0393904    .2710438
                econ |  -.0058754   .0017177    -3.42   0.001    -.0092419   -.0025088
              health |   .0031695   .0017584     1.80   0.071     -.000277     .006616
               pkids |  -.1619267   .0391078    -4.14   0.000    -.2385765   -.0852769
            youngest |  -.0261519   .0106378    -2.46   0.014    -.0470015   -.0053022
               _cons |   .7027791   .4037749     1.74   0.082    -.0886052    1.494163
        -------------+----------------------------------------------------------------
        4.C          |
                page |  -.0275709   .0078581    -3.51   0.000    -.0429726   -.0121693
               pmale |    .204673   .0762452     2.68   0.007     .0552352    .3541108
                pedu |   .1564602   .0716772     2.18   0.029     .0159755    .2969449
               pmarr |   .7047544   .0720235     9.79   0.000      .563591    .8459178
               pwork |   .3199581   .0658001     4.86   0.000     .1909923    .4489238
               metro |  -.4476219   .0635297    -7.05   0.000    -.5721377    -.323106
                econ |  -.0019328   .0013445    -1.44   0.151     -.004568    .0007024
              health |  -.0019929   .0012467    -1.60   0.110    -.0044364    .0004506
               pkids |   .0203683   .0292182     0.70   0.486    -.0368983    .0776349
            youngest |   .0348195   .0085603     4.07   0.000     .0180415    .0515974
               _cons |   .5232246   .3289114     1.59   0.112    -.1214298    1.167879
        ------------------------------------------------------------------------------
        
        Class:    1        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   5.807185   .0242754   239.22   0.000     5.759606    5.854763
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |   2.504744   .0136736   183.18   0.000     2.477944    2.531544
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |   3.376923    .050161    67.32   0.000     3.278609    3.475236
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|    1.45987     .02127                      1.418772     1.50216
         var(e.meetingsd)|   .3822204   .0047181                      .3730841    .3915805
         var(e.meetinggd)|   2.259825   .0408712                      2.181122    2.341368
        ----------------------------------------------------------------------------------
        
        Class:    2        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   7.767166   .0554482   140.08   0.000      7.65849    7.875843
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |   .8260735   .0228727    36.12   0.000     .7812439    .8709031
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |   .1984991   .0415538     4.78   0.000      .117055    .2799431
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|    1.45987     .02127                      1.418772     1.50216
         var(e.meetingsd)|   .3822204   .0047181                      .3730841    .3915805
         var(e.meetinggd)|   2.259825   .0408712                      2.181122    2.341368
        ----------------------------------------------------------------------------------
        
        Class:    3        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   5.958739   .0355316   167.70   0.000     5.889099     6.02838
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |   2.440381   .0213646   114.23   0.000     2.398508    2.482255
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |  -3.366734   .0803316   -41.91   0.000    -3.524181   -3.209287
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|    1.45987     .02127                      1.418772     1.50216
         var(e.meetingsd)|   .3822204   .0047181                      .3730841    .3915805
         var(e.meetinggd)|   2.259825   .0408712                      2.181122    2.341368
        ----------------------------------------------------------------------------------
        
        Class:    4        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   4.224788   .0277109   152.46   0.000     4.170475      4.2791
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |   .6458396   .0129199    49.99   0.000     .6205171    .6711621
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |   .0361678   .0136705     2.65   0.008     .0093741    .0629615
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|    1.45987     .02127                      1.418772     1.50216
         var(e.meetingsd)|   .3822204   .0047181                      .3730841    .3915805
         var(e.meetinggd)|   2.259825   .0408712                      2.181122    2.341368
        ----------------------------------------------------------------------------------
        
        . estat lcprob
        
        Latent class marginal probabilities                     Number of obs = 26,678
        
        --------------------------------------------------------------
                     |            Delta-method
                     |     Margin   std. err.     [95% conf. interval]
        -------------+------------------------------------------------
                   C |
                  1  |   .2917769   .0057049      .2807228     .303083
                  2  |    .153538   .0060884      .1419808    .1658541
                  3  |   .1382465   .0044964      .1296672    .1472974
                  4  |   .4164386   .0072529      .4022963    .4307197
        --------------------------------------------------------------
        
        . gsem (meetingx meetingsd meetinggd <-,)(C <- page pmale pedu pmarr pwork metro econ health pkids y
        > oungest), lclass(C 4, base(2)) vce(cluster pid) from(B, skip)
        
        Fitting class model:
        
        Iteration 0:   (class) log likelihood = -41534.777  
        Iteration 8:   (class) log likelihood = -6692.8873  
        
        Fitting outcome model:
        
        Iteration 0:   (outcome) log likelihood = -281792.23   
        Iteration 7:   (outcome) log likelihood = -153716.21  
        
        Refining starting values:
        
        Iteration 0:   (EM) log likelihood = -162932.36
        Iteration 20:  (EM) log likelihood = -150678.08
        note: EM algorithm reached maximum iterations.
        
        Fitting full model:
        
        Iteration 0:   log pseudolikelihood = -143473.13  (not concave)
        Iteration 5:   log pseudolikelihood = -142950.18  
        
        Generalized structural equation model                   Number of obs = 26,678
        Log pseudolikelihood = -142950.18
        
         ( 1)  [/]var(e.meetingx)#1bn.C - [/]var(e.meetingx)#4.C = 0
         ( 2)  [/]var(e.meetingx)#2.C - [/]var(e.meetingx)#4.C = 0
         ( 3)  [/]var(e.meetingx)#3.C - [/]var(e.meetingx)#4.C = 0
         ( 4)  [/]var(e.meetingsd)#1bn.C - [/]var(e.meetingsd)#4.C = 0
         ( 5)  [/]var(e.meetingsd)#2.C - [/]var(e.meetingsd)#4.C = 0
         ( 6)  [/]var(e.meetingsd)#3.C - [/]var(e.meetingsd)#4.C = 0
         ( 7)  [/]var(e.meetinggd)#1bn.C - [/]var(e.meetinggd)#4.C = 0
         ( 8)  [/]var(e.meetinggd)#2.C - [/]var(e.meetinggd)#4.C = 0
         ( 9)  [/]var(e.meetinggd)#3.C - [/]var(e.meetinggd)#4.C = 0
        
                                        (Std. err. adjusted for 5,849 clusters in pid)
        ------------------------------------------------------------------------------
                     |               Robust
                     | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -------------+----------------------------------------------------------------
        1.C          |
                page |   .0204243   .0074725     2.73   0.006     .0057785    .0350701
               pmale |  -.1040803   .0717653    -1.45   0.147    -.2447377     .036577
                pedu |  -.1637761   .0672238    -2.44   0.015    -.2955324   -.0320198
               pmarr |  -.5096934   .0698017    -7.30   0.000    -.6465023   -.3728845
               pwork |  -.3000703    .063975    -4.69   0.000    -.4254591   -.1746816
               metro |    .342702   .0589389     5.81   0.000     .2271839    .4582202
                econ |   .0006217   .0013159     0.47   0.637    -.0019574    .0032008
              health |   .0014688   .0012492     1.18   0.240    -.0009795    .0039171
               pkids |    .253729   .0315213     8.05   0.000     .1919483    .3155096
            youngest |   -.017477   .0079516    -2.20   0.028    -.0330619   -.0018922
               _cons |  -2.085031   .3304226    -6.31   0.000    -2.732648   -1.437415
        -------------+----------------------------------------------------------------
        2.C          |  (base outcome)
        -------------+----------------------------------------------------------------
        3.C          |
                page |   .0236474    .008472     2.79   0.005     .0070426    .0402521
               pmale |  -.1973168   .0848559    -2.33   0.020    -.3636313   -.0310022
                pedu |    .091123    .079884     1.14   0.254    -.0654467    .2476928
               pmarr |  -.5409875   .0810001    -6.68   0.000    -.6997449   -.3822302
               pwork |  -.3734373   .0711102    -5.25   0.000    -.5128107    -.234064
               metro |   .3647285   .0693064     5.26   0.000     .2288905    .5005665
                econ |  -.0071846    .001516    -4.74   0.000    -.0101559   -.0042133
              health |    .004535   .0015079     3.01   0.003     .0015796    .0074905
               pkids |  -.0877064   .0360988    -2.43   0.015    -.1584588    -.016954
            youngest |  -.0423358   .0094112    -4.50   0.000    -.0607814   -.0238902
               _cons |  -.4074428   .3539598    -1.15   0.250    -1.101191    .2863056
        -------------+----------------------------------------------------------------
        4.C          |
                page |   .0241856   .0106706     2.27   0.023     .0032717    .0450995
               pmale |  -.1805619   .1029905    -1.75   0.080    -.3824194    .0212957
                pedu |    -.23138   .0987172    -2.34   0.019    -.4248621   -.0378978
               pmarr |  -.8478983   .0924016    -9.18   0.000    -1.029002   -.6667945
               pwork |   -.148535   .0889952    -1.67   0.095    -.3229623    .0258924
               metro |   .1746956   .0872895     2.00   0.045     .0036115    .3457798
                econ |  -.0050352   .0020023    -2.51   0.012    -.0089597   -.0011108
              health |   .0023525   .0019866     1.18   0.236    -.0015411    .0062461
               pkids |  -.3197573   .0496235    -6.44   0.000    -.4170175   -.2224971
            youngest |  -.0228336   .0123386    -1.85   0.064    -.0470167    .0013495
               _cons |  -.4845144   .4533944    -1.07   0.285    -1.373151    .4041223
        ------------------------------------------------------------------------------
        
        Class:    1        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   5.809119   .0307645   188.83   0.000     5.748821    5.869416
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |    2.10436   .0191687   109.78   0.000      2.06679     2.14193
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |   2.173597   .0588696    36.92   0.000     2.058215     2.28898
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|   2.791409   .0411317                      2.711945    2.873201
         var(e.meetingsd)|   .3364498   .0041352                      .3284419     .344653
         var(e.meetinggd)|   1.465292   .0286423                      1.410216    1.522519
        ----------------------------------------------------------------------------------
        
        Class:    2        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   5.030891   .0325866   154.39   0.000     4.967023     5.09476
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |   .5618507   .0106153    52.93   0.000     .5410451    .5826563
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |  -.0755626   .0115269    -6.56   0.000     -.098155   -.0529702
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|   2.791409   .0411317                      2.711945    2.873201
         var(e.meetingsd)|   .3364498   .0041352                      .3284419     .344653
         var(e.meetinggd)|   1.465292   .0286423                      1.410216    1.522519
        ----------------------------------------------------------------------------------
        
        Class:    3        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   5.992643   .0343397   174.51   0.000     5.925338    6.059948
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |   2.376467   .0209966   113.18   0.000     2.335314     2.41762
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |  -3.472029   .0672514   -51.63   0.000     -3.60384   -3.340219
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|   2.791409   .0411317                      2.711945    2.873201
         var(e.meetingsd)|   .3364498   .0041352                      .3284419     .344653
         var(e.meetinggd)|   1.465292   .0286423                      1.410216    1.522519
        ----------------------------------------------------------------------------------
        
        Class:    4        
        
        Response: meetingx                                      Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetingsd                                     Number of obs = 26,678
        Family:   Gaussian 
        Link:     Identity 
        
        Response: meetinggd                                     Number of obs = 26,600
        Family:   Gaussian 
        Link:     Identity 
        
                                            (Std. err. adjusted for 5,849 clusters in pid)
        ----------------------------------------------------------------------------------
                         |               Robust
                         | Coefficient  std. err.      z    P>|z|     [95% conf. interval]
        -----------------+----------------------------------------------------------------
        meetingx         |
                   _cons |   5.990273   .0338807   176.80   0.000     5.923868    6.056678
        -----------------+----------------------------------------------------------------
        meetingsd        |
                   _cons |   2.931292   .0244568   119.86   0.000     2.883358    2.979226
        -----------------+----------------------------------------------------------------
        meetinggd        |
                   _cons |   5.316938   .0925858    57.43   0.000     5.135474    5.498403
        -----------------+----------------------------------------------------------------
          var(e.meetingx)|   2.791409   .0411317                      2.711945    2.873201
         var(e.meetingsd)|   .3364498   .0041352                      .3284419     .344653
         var(e.meetinggd)|   1.465292   .0286423                      1.410216    1.522519
        ----------------------------------------------------------------------------------
        
        . estat lcprob
        
        Latent class marginal probabilities                     Number of obs = 26,678
        
        --------------------------------------------------------------
                     |            Delta-method
                     |     Margin   std. err.     [95% conf. interval]
        -------------+------------------------------------------------
                   C |
                  1  |   .2709957     .00536      .2606191    .2816281
                  2  |     .49138   .0062657      .4791072    .5036633
                  3  |   .1422319   .0040713      .1344374    .1503997
                  4  |   .0953924   .0054951      .0851526    .1067199
        --------------------------------------------------------------

        Comment

        Working...
        X