how to obtain 99% credible intervals for defined parameters

[follhim]

Hello! Title says it all. I can't seem to get 99% credible intervals for my defined parameters. How to do so?

> summary(ih_suicide_mgm_race_diffs)
blavaan 0.5.9 ended normally after 6000 iterations

  Estimator                                      BAYES
  Optimization method                             MCMC
  Number of model parameters                       104

  Number of observations per group:                   
    White                                          931
    Latine                                         158
    Other_Mult                                     249
    Black                                          180
  Number of missing patterns per group:               
    White                                            5
    Latine                                           6
    Other_Mult                                       6
    Black                                            6

  Statistic                                 MargLogLik         PPP
  Value                                             NA       1.000

Parameter Estimates:



Group 1 [White]:

Latent Variables:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  SUID =~                                                                      
    t1suicide         1.000                                                    
    t2suicide         1.000                                                    
    t3suicide         1.000                                                    
  STIGMA =~                                                                    
    w1feltstigma      1.000                                                    
    w2feltstigma      1.000                                                    
    w3feltstigma      1.000                                                    
  T1SUID =~                                                                    
    t1suicide         1.000                                                    
  T2SUID =~                                                                    
    t2suicide         1.000                                                    
  T3SUID =~                                                                    
    t3suicide         1.000                                                    
  T1STIGMA =~                                                                  
    w1feltstigma      1.000                                                    
  T2STIGMA =~                                                                  
    w2feltstigma      1.000                                                    
  T3STIGMA =~                                                                  
    w3feltstigma      1.000                                                    

Regressions:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T2SUID ~                                                                     
    T1SUID  (a1_1)    0.127    0.196   -0.353    0.396    1.001    normal(0,10)
    T1STIGM (b1_1)    0.182    0.290   -0.386    0.763    1.000    normal(0,10)
  T2STIGMA ~                                                                   
    T1SUID  (c1_1)    0.063    0.053   -0.036    0.173    1.000    normal(0,10)
    T1STIGM (d1_1)    0.131    0.133   -0.149    0.372    1.000    normal(0,10)
  T3SUID ~                                                                     
    T2SUID  (a2_1)    0.462    0.034    0.396    0.527    1.000    normal(0,10)
    T2STIGM (b2_1)    0.080    0.232   -0.379    0.539    1.000    normal(0,10)
  T3STIGMA ~                                                                   
    T2SUID  (c2_1)    0.021    0.025   -0.028    0.069    1.000    normal(0,10)
    T2STIGM (d2_1)    0.236    0.115   -0.004    0.444    1.001    normal(0,10)

Covariances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T1SUID ~~                                                                    
    T1STIGMA          0.096    0.047    0.005    0.191    1.001     lkj_corr(1)
 .T2SUID ~~                                                                    
   .T2STIGMA          0.165    0.153   -0.135    0.465    1.000     lkj_corr(1)
 .T3SUID ~~                                                                    
   .T3STIGMA         -0.035    0.153   -0.334    0.264    1.000     lkj_corr(1)
  SUID ~~                                                                      
    STIGMA            0.033    0.048   -0.059    0.129    1.000     lkj_corr(1)

Intercepts:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
   .t1suicide         0.649    0.034    0.583    0.715    1.000    normal(0,32)
   .t2suicide         3.254    0.080    3.100    3.410    1.000    normal(0,32)
   .t3suicide         2.632    0.081    2.472    2.794    1.000    normal(0,32)
   .w1feltstigma      2.633    0.031    2.573    2.694    1.000    normal(0,32)
   .w2feltstigma      2.464    0.052    2.363    2.565    1.000    normal(0,32)
   .w3feltstigma      2.443    0.075    2.298    2.592    1.000    normal(0,32)
    SUID              0.000                                                    
    STIGMA            0.000                                                    
    T1SUID            0.000                                                    
   .T2SUID            0.000                                                    
   .T3SUID            0.000                                                    
    T1STIGMA          0.000                                                    
   .T2STIGMA          0.000                                                    
   .T3STIGMA          0.000                                                    

Variances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
    SUID              0.163    0.145    0.000    0.486    1.001 gamma(1,.5)[sd]
    STIGMA            0.588    0.051    0.484    0.686    1.000 gamma(1,.5)[sd]
    T1SUID            0.903    0.150    0.574    1.125    1.001 gamma(1,.5)[sd]
    T1STIGMA          0.309    0.044    0.228    0.402    1.000 gamma(1,.5)[sd]
   .T2SUID            5.547    0.295    4.967    6.127    1.000 gamma(1,.5)[sd]
   .T2STIGMA          0.263    0.045    0.173    0.349    1.000 gamma(1,.5)[sd]
   .T3SUID            4.739    0.225    4.314    5.192    1.000 gamma(1,.5)[sd]
   .T3STIGMA          0.384    0.035    0.317    0.456    1.000 gamma(1,.5)[sd]
   .t1suicide         0.000                                                    
   .t2suicide         0.000                                                    
   .t3suicide         0.000                                                    
   .w1feltstigma      0.000                                                    
   .w2feltstigma      0.000                                                    
   .w3feltstigma      0.000                                                    


Group 2 [Latine]:

Latent Variables:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  SUID =~                                                                      
    t1suicide         1.000                                                    
    t2suicide         1.000                                                    
    t3suicide         1.000                                                    
  STIGMA =~                                                                    
    w1feltstigma      1.000                                                    
    w2feltstigma      1.000                                                    
    w3feltstigma      1.000                                                    
  T1SUID =~                                                                    
    t1suicide         1.000                                                    
  T2SUID =~                                                                    
    t2suicide         1.000                                                    
  T3SUID =~                                                                    
    t3suicide         1.000                                                    
  T1STIGMA =~                                                                  
    w1feltstigma      1.000                                                    
  T2STIGMA =~                                                                  
    w2feltstigma      1.000                                                    
  T3STIGMA =~                                                                  
    w3feltstigma      1.000                                                    

Regressions:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T2SUID ~                                                                     
    T1SUID  (a1_2)   -5.294    3.292  -12.734    0.026    1.002    normal(0,10)
    T1STIGM (b1_2)   -0.390    1.539   -3.513    2.654    1.000    normal(0,10)
  T2STIGMA ~                                                                   
    T1SUID  (c1_2)   -0.243    0.774   -1.849    1.267    1.000    normal(0,10)
    T1STIGM (d1_2)   -0.259    0.488   -1.339    0.466    1.001    normal(0,10)
  T3SUID ~                                                                     
    T2SUID  (a2_2)    0.586    0.084    0.418    0.750    1.000    normal(0,10)
    T2STIGM (b2_2)   -0.201    0.462   -1.177    0.651    1.000    normal(0,10)
  T3STIGMA ~                                                                   
    T2SUID  (c2_2)   -0.106    0.102   -0.304    0.096    1.000    normal(0,10)
    T2STIGM (d2_2)    0.056    0.412   -0.916    0.707    1.001    normal(0,10)

Covariances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T1SUID ~~                                                                    
    T1STIGMA         -0.006    0.100   -0.220    0.187    1.000     lkj_corr(1)
 .T2SUID ~~                                                                    
   .T2STIGMA          0.047    0.321   -0.588    0.795    1.000     lkj_corr(1)
 .T3SUID ~~                                                                    
   .T3STIGMA         -0.171    0.361   -0.897    0.540    1.000     lkj_corr(1)
  SUID ~~                                                                      
    STIGMA            0.100    0.103   -0.091    0.317    1.000     lkj_corr(1)

Intercepts:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
   .t1suicide         0.633    0.087    0.464    0.803    1.000    normal(0,32)
   .t2suicide         2.082    0.199    1.694    2.473    1.000    normal(0,32)
   .t3suicide         1.568    0.204    1.169    1.969    1.000    normal(0,32)
   .w1feltstigma      2.733    0.075    2.585    2.882    1.000    normal(0,32)
   .w2feltstigma      2.424    0.221    1.988    2.858    1.000    normal(0,32)
   .w3feltstigma      2.779    0.344    2.109    3.467    1.000    normal(0,32)
    SUID              0.000                                                    
    STIGMA            0.000                                                    
    T1SUID            0.000                                                    
   .T2SUID            0.000                                                    
   .T3SUID            0.000                                                    
    T1STIGMA          0.000                                                    
   .T2STIGMA          0.000                                                    
   .T3STIGMA          0.000                                                    

Variances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
    SUID              0.956    0.237    0.234    1.290    1.004 gamma(1,.5)[sd]
    STIGMA            0.476    0.151    0.099    0.748    1.002 gamma(1,.5)[sd]
    T1SUID            0.232    0.223    0.035    0.960    1.005 gamma(1,.5)[sd]
    T1STIGMA          0.424    0.154    0.159    0.783    1.002 gamma(1,.5)[sd]
   .T2SUID            1.883    1.860    0.002    5.962    1.002 gamma(1,.5)[sd]
   .T2STIGMA          0.265    0.176    0.003    0.635    1.001 gamma(1,.5)[sd]
   .T3SUID            3.690    0.433    2.940    4.628    1.000 gamma(1,.5)[sd]
   .T3STIGMA          0.486    0.174    0.150    0.859    1.000 gamma(1,.5)[sd]
   .t1suicide         0.000                                                    
   .t2suicide         0.000                                                    
   .t3suicide         0.000                                                    
   .w1feltstigma      0.000                                                    
   .w2feltstigma      0.000                                                    
   .w3feltstigma      0.000                                                    


Group 3 [Other_Mult]:

Latent Variables:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  SUID =~                                                                      
    t1suicide         1.000                                                    
    t2suicide         1.000                                                    
    t3suicide         1.000                                                    
  STIGMA =~                                                                    
    w1feltstigma      1.000                                                    
    w2feltstigma      1.000                                                    
    w3feltstigma      1.000                                                    
  T1SUID =~                                                                    
    t1suicide         1.000                                                    
  T2SUID =~                                                                    
    t2suicide         1.000                                                    
  T3SUID =~                                                                    
    t3suicide         1.000                                                    
  T1STIGMA =~                                                                  
    w1feltstigma      1.000                                                    
  T2STIGMA =~                                                                  
    w2feltstigma      1.000                                                    
  T3STIGMA =~                                                                  
    w3feltstigma      1.000                                                    

Regressions:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T2SUID ~                                                                     
    T1SUID  (a1_3)    0.107    0.326   -0.450    0.508    1.002    normal(0,10)
    T1STIGM (b1_3)    0.008    0.512   -0.864    0.965    1.001    normal(0,10)
  T2STIGMA ~                                                                   
    T1SUID  (c1_3)   -0.014    0.120   -0.171    0.196    1.003    normal(0,10)
    T1STIGM (d1_3)    0.192    0.286   -0.463    0.588    1.003    normal(0,10)
  T3SUID ~                                                                     
    T2SUID  (a2_3)    0.404    0.083    0.243    0.572    1.000    normal(0,10)
    T2STIGM (b2_3)   -1.139    0.434   -2.130   -0.414    1.001    normal(0,10)
  T3STIGMA ~                                                                   
    T2SUID  (c2_3)    0.163    0.051    0.065    0.266    1.000    normal(0,10)
    T2STIGM (d2_3)    0.091    0.301   -0.608    0.597    1.002    normal(0,10)

Covariances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T1SUID ~~                                                                    
    T1STIGMA          0.137    0.087   -0.025    0.319    1.001     lkj_corr(1)
 .T2SUID ~~                                                                    
   .T2STIGMA          0.193    0.416   -0.627    1.009    1.000     lkj_corr(1)
 .T3SUID ~~                                                                    
   .T3STIGMA          0.188    0.332   -0.436    0.864    1.000     lkj_corr(1)
  SUID ~~                                                                      
    STIGMA            0.016    0.072   -0.141    0.165    1.001     lkj_corr(1)

Intercepts:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
   .t1suicide         0.819    0.070    0.683    0.958    1.000    normal(0,32)
   .t2suicide         2.401    0.172    2.064    2.742    1.000    normal(0,32)
   .t3suicide         1.372    0.146    1.084    1.662    1.000    normal(0,32)
   .w1feltstigma      2.568    0.064    2.445    2.691    1.000    normal(0,32)
   .w2feltstigma      2.447    0.170    2.113    2.781    1.000    normal(0,32)
   .w3feltstigma      1.871    0.294    1.286    2.440    1.000    normal(0,32)
    SUID              0.000                                                    
    STIGMA            0.000                                                    
    T1SUID            0.000                                                    
   .T2SUID            0.000                                                    
   .T3SUID            0.000                                                    
    T1STIGMA          0.000                                                    
   .T2STIGMA          0.000                                                    
   .T3STIGMA          0.000                                                    

Variances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
    SUID              0.085    0.116    0.000    0.408    1.001 gamma(1,.5)[sd]
    STIGMA            0.496    0.125    0.219    0.730    1.002 gamma(1,.5)[sd]
    T1SUID            1.118    0.154    0.766    1.387    1.001 gamma(1,.5)[sd]
    T1STIGMA          0.479    0.122    0.252    0.753    1.002 gamma(1,.5)[sd]
   .T2SUID            7.308    0.714    6.020    8.755    1.000 gamma(1,.5)[sd]
   .T2STIGMA          0.387    0.117    0.156    0.614    1.002 gamma(1,.5)[sd]
   .T3SUID            3.743    0.427    2.930    4.607    1.000 gamma(1,.5)[sd]
   .T3STIGMA          0.350    0.140    0.064    0.644    1.002 gamma(1,.5)[sd]
   .t1suicide         0.000                                                    
   .t2suicide         0.000                                                    
   .t3suicide         0.000                                                    
   .w1feltstigma      0.000                                                    
   .w2feltstigma      0.000                                                    
   .w3feltstigma      0.000                                                    


Group 4 [Black]:

Latent Variables:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  SUID =~                                                                      
    t1suicide         1.000                                                    
    t2suicide         1.000                                                    
    t3suicide         1.000                                                    
  STIGMA =~                                                                    
    w1feltstigma      1.000                                                    
    w2feltstigma      1.000                                                    
    w3feltstigma      1.000                                                    
  T1SUID =~                                                                    
    t1suicide         1.000                                                    
  T2SUID =~                                                                    
    t2suicide         1.000                                                    
  T3SUID =~                                                                    
    t3suicide         1.000                                                    
  T1STIGMA =~                                                                  
    w1feltstigma      1.000                                                    
  T2STIGMA =~                                                                  
    w2feltstigma      1.000                                                    
  T3STIGMA =~                                                                  
    w3feltstigma      1.000                                                    

Regressions:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T2SUID ~                                                                     
    T1SUID  (a1_4)   -0.272    1.224   -3.824    1.122    1.000    normal(0,10)
    T1STIGM (b1_4)   -2.498    2.853   -9.437    2.112    1.001    normal(0,10)
  T2STIGMA ~                                                                   
    T1SUID  (c1_4)    0.052    0.489   -0.665    1.200    1.002    normal(0,10)
    T1STIGM (d1_4)   -1.512    1.682   -5.710    0.295    1.001    normal(0,10)
  T3SUID ~                                                                     
    T2SUID  (a2_4)    0.513    0.065    0.383    0.637    1.000    normal(0,10)
    T2STIGM (b2_4)    0.216    0.250   -0.275    0.717    1.000    normal(0,10)
  T3STIGMA ~                                                                   
    T2SUID  (c2_4)   -0.046    0.076   -0.201    0.101    1.001    normal(0,10)
    T2STIGM (d2_4)   -0.040    0.255   -0.605    0.401    1.000    normal(0,10)

Covariances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
  T1SUID ~~                                                                    
    T1STIGMA          0.108    0.112   -0.105    0.329    1.000     lkj_corr(1)
 .T2SUID ~~                                                                    
   .T2STIGMA         -0.192    0.475   -1.142    0.808    1.000     lkj_corr(1)
 .T3SUID ~~                                                                    
   .T3STIGMA         -0.036    0.229   -0.487    0.416    1.000     lkj_corr(1)
  SUID ~~                                                                      
    STIGMA            0.149    0.115   -0.055    0.379    1.000     lkj_corr(1)

Intercepts:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
   .t1suicide         0.600    0.083    0.438    0.762    1.000    normal(0,32)
   .t2suicide         1.785    0.185    1.423    2.158    1.000    normal(0,32)
   .t3suicide         1.318    0.166    0.991    1.643    1.000    normal(0,32)
   .w1feltstigma      2.864    0.066    2.735    2.997    1.000    normal(0,32)
   .w2feltstigma      2.548    0.319    1.924    3.191    1.000    normal(0,32)
   .w3feltstigma      2.799    0.293    2.225    3.372    1.000    normal(0,32)
    SUID              0.000                                                    
    STIGMA            0.000                                                    
    T1SUID            0.000                                                    
   .T2SUID            0.000                                                    
   .T3SUID            0.000                                                    
    T1STIGMA          0.000                                                    
   .T2STIGMA          0.000                                                    
   .T3STIGMA          0.000                                                    

Variances:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
    SUID              0.270    0.250    0.001    0.932    1.000 gamma(1,.5)[sd]
    STIGMA            0.615    0.106    0.407    0.826    1.000 gamma(1,.5)[sd]
    T1SUID            0.978    0.271    0.298    1.407    1.000 gamma(1,.5)[sd]
    T1STIGMA          0.177    0.097    0.028    0.393    1.001 gamma(1,.5)[sd]
   .T2SUID            4.739    1.380    0.965    6.861    1.000 gamma(1,.5)[sd]
   .T2STIGMA          0.301    0.254    0.001    0.888    1.000 gamma(1,.5)[sd]
   .T3SUID            3.044    0.336    2.455    3.775    1.000 gamma(1,.5)[sd]
   .T3STIGMA          0.498    0.133    0.266    0.789    1.000 gamma(1,.5)[sd]
   .t1suicide         0.000                                                    
   .t2suicide         0.000                                                    
   .t3suicide         0.000                                                    
   .w1feltstigma      0.000                                                    
   .w2feltstigma      0.000                                                    
   .w3feltstigma      0.000                                                    

Defined Parameters:
                   Estimate  Post.SD pi.lower pi.upper     Rhat    Prior       
    T1STIGMA_T2SUI    0.572    1.572   -2.510    3.653                         
    T1STIGMA_T2SUI    0.174    0.589   -0.980    1.328                         
    T1STIGMA_T2SUI    2.680    2.870   -2.944    8.304                         
    T1STIGMA_T2SUI   -0.398    1.627   -3.587    2.791                         
    T1STIGMA_T2SUI    2.108    3.251   -4.263    8.480                         
    T1STIGMA_T2SUI    2.506    2.891   -3.161    8.173                         
    T1SUID_T2STIGM    0.306    0.776   -1.215    1.827                         
    T1SUID_T2STIGM    0.077    0.131   -0.179    0.333                         
    T1SUID_T2STIGM    0.011    0.492   -0.953    0.975                         
    T1SUID_T2STIGM   -0.229    0.785   -1.767    1.309                         
    T1SUID_T2STIGM   -0.295    0.915   -2.088    1.498                         
    T1SUID_T2STIGM   -0.066    0.504   -1.054    0.921                         
    T2STIGMA_T3SUI    0.281    0.517   -0.734    1.295                         
    T2STIGMA_T3SUI    1.219    0.494    0.252    2.186                         
    T2STIGMA_T3SUI   -0.136    0.341   -0.804    0.532                         
    T2STIGMA_T3SUI    0.938    0.637   -0.311    2.188                         
    T2STIGMA_T3SUI   -0.417    0.526   -1.448    0.615                         
    T2STIGMA_T3SUI   -1.355    0.502   -2.340   -0.370                         
    T2SUID_T3STIGM    0.127    0.105   -0.078    0.332                         
    T2SUID_T3STIGM   -0.142    0.057   -0.253   -0.031                         
    T2SUID_T3STIGM    0.068    0.080   -0.090    0.225                         
    T2SUID_T3STIGM   -0.270    0.113   -0.491   -0.048                         
    T2SUID_T3STIGM   -0.060    0.127   -0.309    0.190                         
    T2SUID_T3STIGM    0.210    0.092    0.029    0.390  

Here's what I've tried

> summary(ih_suicide_mgm_race_diffs,
+         prob = .99,
+         central.tendency = c("median"))
Error in .local(object, ...) : 
  unused arguments (prob = 0.99, central.tendency = "median")

> parameterEstimates(ih_suicide_mgm_race_diffs, 
+                    level = 0.99) %>%
+   filter(op == ":=")
                    lhs op         rhs block group                label    est
1  T1STIGMA_T2SUID_g1g2 := b1_g1-b1_g2     0     0 T1STIGMA_T2SUID_g1g2  0.572
2  T1STIGMA_T2SUID_g1g3 := b1_g1-b1_g3     0     0 T1STIGMA_T2SUID_g1g3  0.174
3  T1STIGMA_T2SUID_g1g4 := b1_g1-b1_g4     0     0 T1STIGMA_T2SUID_g1g4  2.680
4  T1STIGMA_T2SUID_g2g3 := b1_g2-b1_g3     0     0 T1STIGMA_T2SUID_g2g3 -0.398
5  T1STIGMA_T2SUID_g2g4 := b1_g2-b1_g4     0     0 T1STIGMA_T2SUID_g2g4  2.108
6  T1STIGMA_T2SUID_g3g4 := b1_g3-b1_g4     0     0 T1STIGMA_T2SUID_g3g4  2.506
7  T1SUID_T2STIGMA_g1g2 := c1_g1-c1_g2     0     0 T1SUID_T2STIGMA_g1g2  0.306
8  T1SUID_T2STIGMA_g1g3 := c1_g1-c1_g3     0     0 T1SUID_T2STIGMA_g1g3  0.077
9  T1SUID_T2STIGMA_g1g4 := c1_g1-c1_g4     0     0 T1SUID_T2STIGMA_g1g4  0.011
10 T1SUID_T2STIGMA_g2g3 := c1_g2-c1_g3     0     0 T1SUID_T2STIGMA_g2g3 -0.229
11 T1SUID_T2STIGMA_g2g4 := c1_g2-c1_g4     0     0 T1SUID_T2STIGMA_g2g4 -0.295
12 T1SUID_T2STIGMA_g3g4 := c1_g3-c1_g4     0     0 T1SUID_T2STIGMA_g3g4 -0.066
13 T2STIGMA_T3SUID_g1g2 := b2_g1-b2_g2     0     0 T2STIGMA_T3SUID_g1g2  0.281
14 T2STIGMA_T3SUID_g1g3 := b2_g1-b2_g3     0     0 T2STIGMA_T3SUID_g1g3  1.219
15 T2STIGMA_T3SUID_g1g4 := b2_g1-b2_g4     0     0 T2STIGMA_T3SUID_g1g4 -0.136
16 T2STIGMA_T3SUID_g2g3 := b2_g2-b2_g3     0     0 T2STIGMA_T3SUID_g2g3  0.938
17 T2STIGMA_T3SUID_g2g4 := b2_g2-b2_g4     0     0 T2STIGMA_T3SUID_g2g4 -0.417
18 T2STIGMA_T3SUID_g3g4 := b2_g3-b2_g4     0     0 T2STIGMA_T3SUID_g3g4 -1.355
19 T2SUID_T3STIGMA_g1g2 := c2_g1-c2_g2     0     0 T2SUID_T3STIGMA_g1g2  0.127
20 T2SUID_T3STIGMA_g1g3 := c2_g1-c2_g3     0     0 T2SUID_T3STIGMA_g1g3 -0.142
21 T2SUID_T3STIGMA_g1g4 := c2_g1-c2_g4     0     0 T2SUID_T3STIGMA_g1g4  0.068
22 T2SUID_T3STIGMA_g2g3 := c2_g2-c2_g3     0     0 T2SUID_T3STIGMA_g2g3 -0.270
23 T2SUID_T3STIGMA_g2g4 := c2_g2-c2_g4     0     0 T2SUID_T3STIGMA_g2g4 -0.060
24 T2SUID_T3STIGMA_g3g4 := c2_g3-c2_g4     0     0 T2SUID_T3STIGMA_g3g4  0.210

[ecmerkle]

Those intervals for the defined parameters are obtained in a different manner as compared to the intervals for other parameters, which partly explains why things are not working as you may expect.

I also believe there is a bug here that never got fixed. See this post to the google group for a workaround that involves obtaining the posterior samples and then applying the operations to obtain posterior distributions of defined parameters. And I will try to look at this in the near future.

[follhim]

Super helpful!

In case you or others are curious, here's a function I had AI create me that works almost exactly (I verified comparing my 95% calculations with the results from summary).

library(rstan)
library(dplyr)
library(coda)

#' Extract HPD Intervals for Defined Parameters from blavaan Model
#'
#' @param blavaan_model A fitted blavaan model object
#' @param prob The probability level for credible intervals (default = 0.99)
#' @return A data frame with parameter names, means, SDs, and HPD intervals
#' @export
#'
#' @details 
#' This function extracts posterior samples for defined parameters (those specified 
#' with :=) and computes Highest Posterior Density (HPD) intervals. 
#' 
#' The function:
#' 1. Extracts base parameter samples from Stan
#' 2. Computes defined parameter values from these samples
#' 3. Calculates summary statistics and HPD intervals
#' 
#' Uses posterior means (not medians) to match blavaan's summary() output.
#' 
#' @examples
#' # Get 99% HPD intervals
#' ci_99 <- get_defined_params_hpd(model, prob = 0.99)
#' 
#' # Get 95% HPD intervals
#' ci_95 <- get_defined_params_hpd(model, prob = 0.95)
#' 
get_defined_params_hpd <- function(blavaan_model, prob = 0.99) {
  
  # Validate inputs
  if (!grepl("stan", blavaan_model@Options$target)) {
    stop("This function currently only supports Stan backend models.")
  }
  
  if (prob <= 0 || prob >= 1) {
    stop("prob must be between 0 and 1")
  }
  
  # Extract Stan samples
  stan_samples <- extract(blavaan_model@external$mcmcout, permuted = TRUE)
  
  # Get parameter table
  pt_df <- data.frame(
    label = blavaan_model@ParTable$label,
    pxnames = blavaan_model@ParTable$pxnames,
    lhs = blavaan_model@ParTable$lhs,
    rhs = blavaan_model@ParTable$rhs,
    op = blavaan_model@ParTable$op,
    stringsAsFactors = FALSE
  )
  
  # Helper: Extract Stan parameter index
  get_stan_index <- function(pxname) {
    if(is.na(pxname)) return(NA)
    as.numeric(gsub(".*\\[(\\d+)\\].*", "\\1", pxname))
  }
  
  # Helper: Get samples for a labeled parameter
  get_param_samples <- function(label_name) {
    row <- pt_df[pt_df$label == label_name & !is.na(pt_df$label), ]
    if(nrow(row) == 0) return(NULL)
    
    idx <- get_stan_index(row$pxnames[1])
    if(is.na(idx)) return(NULL)
    
    return(stan_samples$bet_sign[, idx])
  }
  
  # Get all defined parameters
  defined_params <- pt_df %>%
    filter(op == ":=", label != "")
  
  if(nrow(defined_params) == 0) {
    stop("No defined parameters found in the model.")
  }
  
  # Calculate HPD intervals for each defined parameter
  results_list <- list()
  
  for(i in 1:nrow(defined_params)) {
    param_label <- defined_params$label[i]
    param_def <- defined_params$rhs[i]
    
    # Parse definition (format: param1 - param2)
    parts <- strsplit(gsub("\\s+", "", param_def), "-")[[1]]
    
    if(length(parts) == 2) {
      # Get samples
      samples1 <- get_param_samples(parts[1])
      samples2 <- get_param_samples(parts[2])
      
      if(!is.null(samples1) && !is.null(samples2)) {
        # Calculate difference
        diff_samples <- samples1 - samples2
        
        # Calculate HPD interval
        hpd <- HPDinterval(mcmc(diff_samples), prob = prob)
        
        # Store results
        results_list[[param_label]] <- data.frame(
          parameter = param_label,
          mean = mean(diff_samples),
          sd = sd(diff_samples),
          hpd_lower = hpd[1, "lower"],
          hpd_upper = hpd[1, "upper"],
          stringsAsFactors = FALSE
        )
      }
    }
  }
  
  # Combine results
  if(length(results_list) == 0) {
    stop("No valid defined parameters could be calculated.")
  }
  
  ci_results <- do.call(rbind, results_list)
  rownames(ci_results) <- NULL
  
  # Add formatted CI column
  ci_results <- ci_results %>%
    mutate(
      HPD = sprintf("[%.3f, %.3f]", hpd_lower, hpd_upper)
    ) %>%
    select(parameter, mean, sd, hpd_lower, hpd_upper, HPD)
  
  # Add metadata
  attr(ci_results, "prob") <- prob
  attr(ci_results, "n_samples") <- length(diff_samples)
  
  return(ci_results)
}

[ecmerkle]

Thanks for following up! And I will leave this issue open to remind me about it.