indirect effect with multiple groups

[KwonHyun Kim]

  First of all, thank you for the wonderful package!

  Anyhow, I found that when I do multiple group analysis with FIML missing with indirect effect,

  indirect effects are only reported once rather than for each groups...

  It seems like there's a little bug...

[Terrence Jorgensen]

That doesn't seem likely.  Note that you would need to label each group's relevant direct effects in your model syntax, then define the indirect separately for each group in your model syntax.  If that is what you are doing and you don't see all of the indirect effects you defined in your syntax when you print results, please post your R script (and perhaps some data) that reproduces this.

Terrence D. Jorgensen

Postdoctoral Researcher, Methods and Statistics

Research Institute for Child Development and Education, the University of Amsterdam

UvA web page: http://www.uva.nl/profile/t.d.jorgensen

[KwonHyun Kim]

  Can I label each group's relevant direct effects? 

  In the manual, https://dornsife.usc.edu/assets/sites/210/docs/GC3/lavaan_tutorial.pdf

  I couldn't find how to do it...

  Here's my code and data

df <- read.csv( file="sample.csv")

theModel <- "

bidToEst ~ b2*n.bidder + est + b4*DM

CO ~ b1*bidToEst + est + DM

indirect1 := b1*b4

"

theFitMissing <- sem(theModel, data=df, group="PT1", meanstructure = T, 

                     missing="fiml", estimator="MLR")

summary(theFitMissing)

The result is like this

lavaan (0.5-20) converged normally after 135 iterations

  Number of observations per group         

  Building                                          21

  Civil                                             29

  Number of missing patterns per group     

  Building                                           1

  Civil                                              2

  Estimator                                         ML      Robust

  Minimum Function Test Statistic                3.746       2.807

  Degrees of freedom                                 2           2

  P-value (Chi-square)                           0.154       0.246

  Scaling correction factor                                  1.335

    for the Yuan-Bentler correction

Chi-square for each group:

  Building                                       3.641       2.728

  Civil                                          0.104       0.078

Parameter Estimates:

  Information                                 Observed

  Standard Errors                   Robust.huber.white

Group 1 [Building]:

Regressions:

                   Estimate  Std.Err  Z-value  P(>|z|)

  bidToEst ~                                          

    n.bidder  (b2)    0.019    0.009    2.169    0.030

    est               0.678    1.653    0.410    0.682

    DM        (b4)  -13.616    2.214   -6.149    0.000

  CO ~                                                

    bidToEst  (b1)   -0.453    0.457   -0.990    0.322

    est               0.907    3.026    0.300    0.764

    DM                5.353    6.249    0.857    0.392

Intercepts:

                   Estimate  Std.Err  Z-value  P(>|z|)

    bidToEst        105.451    4.100   25.721    0.000

    CO               40.795   45.749    0.892    0.373

Variances:

                   Estimate  Std.Err  Z-value  P(>|z|)

    bidToEst         13.506    4.381    3.083    0.002

    CO               97.516   41.811    2.332    0.020

Group 2 [Civil]:

Regressions:

                   Estimate  Std.Err  Z-value  P(>|z|)

  bidToEst ~                                          

    n.bidder         -0.000    0.025   -0.000    1.000

    est              -2.422    1.991   -1.216    0.224

    DM               -2.728    3.891   -0.701    0.483

  CO ~                                                

    bidToEst         -0.828    0.779   -1.063    0.288

    est             -14.404    6.902   -2.087    0.037

    DM               -7.946    9.829   -0.808    0.419

Intercepts:

                   Estimate  Std.Err  Z-value  P(>|z|)

    bidToEst         88.487    7.893   11.211    0.000

    CO              104.284   72.177    1.445    0.149

Variances:

                   Estimate  Std.Err  Z-value  P(>|z|)

    bidToEst         25.771   14.807    1.740    0.082

    CO              272.193   52.292    5.205    0.000

Defined Parameters:

                   Estimate  Std.Err  Z-value  P(>|z|)

    indirect1         6.162    6.084    1.013    0.311

Since I can't upload my file, here's the data.

"DM","bidToEst","PT1","CO","n.bidder","est","group"

"DB",83.8953850027579,"Building",11.8538013768767,26.1489973189028,0.610091689742867,"Building DB"

"DBB",85.3986892001144,"Building",14.7828761370761,300.74812145789,0.358698995687694,"Building DBB"

"DBB",81.4429960392939,"Building",3.69912499787114,168.988241973462,0.227192814980091,"Building DBB"

"DB",81.4691296306825,"Civil",3.6008547008547,39.3874249090285,1.12220523334596,"Civil DB"

"DB",96.4267422236228,"Building",5.27016469657253,1.27557004722181,1.55021448978163,"Building DB"

"DB",83.6838769092819,"Civil",0.375725167477131,21.3783094701085,1.27069633291945,"Civil DB"

"DBB",77.6528829233968,"Civil",-3.49505205799899,30.4859385024828,0.409371551856459,"Civil DBB"

"DB",90.4453237947309,"Building",9.66206337628657,26.3950884328478,0.374388311489585,"Building DB"

"DBB",82.0169705117793,"Building",45.0667485596846,113.430471167304,0.475531091690759,"Building DBB"

"DBB",81.4922750509868,"Building",12.7132212441329,65.4506811040261,0.513118122887595,"Building DBB"

"DB",88.0044545928245,"Civil",21.6655768901235,43.930824875474,1.14929769499773,"Civil DB"

"DB",91.6506353899002,"Civil",-20.3344021836164,30.8578076046398,1.10326436108454,"Civil DB"

"DBB",83.2791298477013,"Civil",31.3219190878873,201.509602136067,0.522390979653338,"Civil DBB"

"DB",62.7506072492357,"Civil",54.1566868440535,NA,0.837059273523532,"Civil DB"

"DBB",79.1138025285149,"Building",6.92847920239302,50.9888597277118,1.17781177702804,"Building DBB"

"DBB",82.8909740338334,"Civil",4.88811939040394,12.3742038982423,0.440354387088696,"Civil DBB"

"DBB",75.4416115055647,"Civil",3.72232261027225,NA,1.69519785901629,"Civil DBB"

"DBB",72.7231225242743,"Building",23.8761673493541,23.7328245523964,1.13818953837208,"Building DBB"

"DBB",85.452581383254,"Civil",12.9427664267905,43.3511078808967,0.312992616826135,"Civil DBB"

"DBB",84.5166549792916,"Civil",2.48604072586593,109.605515107659,0.175815874614995,"Civil DBB"

"DBB",83.1935676988753,"Civil",3.76612630420366,45.5852729783663,0.396483636313663,"Civil DBB"

"DBB",79.4104657312587,"Civil",8.25138689696803,95.4548741276217,0.428862066782892,"Civil DBB"

"DBB",83.8666209608296,"Building",6.15781824745712,40.6061102498214,0.241494781574357,"Building DBB"

"DB",88.5333740887958,"Civil",48.2564102564103,NA,0.488208640707848,"Civil DB"

"DBB",78.1612407806485,"Civil",7.37194913493186,NA,0.263184964326562,"Civil DBB"

"DBB",81.2931118428235,"Civil",11.0075039999968,NA,0.143895060097667,"Civil DBB"

"DBB",83.0443707634046,"Building",24.1340748537767,244.102768507199,0.498400231651987,"Building DBB"

"DB",82.4747191053291,"Civil",0.114130316069985,50.8118147943856,1.87886021948522,"Civil DB"

"DB",93.4240434496653,"Building",6.69664323470111,24.5523012002782,0.524033184662273,"Building DB"

"DBB",81.8010456658452,"Civil",-4.05022307591037,19.6870500610748,0.380006170929469,"Civil DBB"

"DBB",87.5424033629156,"Civil",24.9712000600781,26.6711814441543,0.469490913674489,"Civil DBB"

"DBB",87.644666119422,"Civil",40.2013424737863,53.1727755063689,0.434995292624726,"Civil DBB"

"DBB",79.7131959743441,"Building",27.9440308063223,144.833618756418,0.386578045864907,"Building DBB"

"DBB",84.6392509343134,"Civil",8.35116973246018,NA,0.324890904413535,"Civil DBB"

"DBB",84.7762072759827,"Civil",28.4848644882796,28.6037893152554,0.444905758863617,"Civil DBB"

"DBB",80.8927374493548,"Civil",24.0588489118928,60.1373387930844,1.0197094350605,"Civil DBB"

"DBB",83.5744013343398,"Building",5.39992101028155,31.5805351705382,0.601814572270448,"Building DBB"

"DB",93.2345282781678,"Building",9.48750914867786,40.312865270657,1.67099665997445,"Building DB"

"DB",98.1786561595878,"Building",-11.3930542857143,-21.6563695583504,0.757422268819872,"Building DB"

"DB",92.0588370267413,"Building",6.68201221542356,14.9416750349154,0.36348330527428,"Building DB"

"DBB",80.4147086158586,"Civil",45.60993523261,-4.51043924242669,0.324382785535711,"Civil DBB"

"DBB",82.5290448293158,"Civil",24.1176920768758,34.2781431011458,-0.0125577377351569,"Civil DBB"

"DBB",77.0037589192996,"Civil",8.74646642385534,71.6072597683236,0.665455326960194,"Civil DBB"

"DB",84.7825448116935,"Civil",0,-3.22962594449218,1.29910370425686,"Civil DB"

"DB",91.1758969729986,"Building",-1.74236513899324,24.4369446590635,0.592718858854658,"Building DB"

"DBB",73.202675368757,"Building",6.30737296805729,34.1619121158615,0.832680723680093,"Building DBB"

"DB",85.1542985510818,"Civil",0.496335776149234,-1.80671419818071,0.482281951312482,"Civil DB"

"DBB",79.3471478635322,"Civil",-6.71220548901901,121.724942073048,0.541510035066437,"Civil DBB"

"DB",96.8377542322004,"Building",8.12321439480449,53.8208794998502,1.80015946499683,"Building DB"

"DBB",80.6953154147242,"Building",10.1128568281245,0.581357942534066,0.58231535735652,"Building DBB"

2016년 8월 9일 화요일 오후 10시 29분 17초 UTC+9, Terrence Jorgensen 님의 말:

[Terrence Jorgensen]

Can I label each group's relevant direct effects?  In the manual, https://dornsife.usc.edu/assets/sites/210/docs/GC3/lavaan_tutorial.pdf

  I couldn't find how to do it...

Look at the bottom of p. 18

[KwonHyun Kim]

Oh I remember reading this part but I thought it was irrelevant...

Thank you million times!