Hi All,
I need to calculate variance explained for my endogenous latent variables after fitting the pooled data to a structural model in Stage 2 (using tssem2). I have specified a model (all observed variables) with 3 dependent variables (two of which are mediators, with one key "final" dependent variable).
(I guess I'll drop in the A matrix so you know what it looks like:)
KSSE Prestige EnjHlpg BossExp ExtrMotiv Recip TrustB ATS ITS KS KSSE "0" "0" "0" "0" "0" "0" "0" "0" ".3*KSSE_ITS" ".3*KSSE_KS" Prestige "0" "0" "0" "0" "0" "0" "0" ".3*Prest_ATS" "0" "0" EnjHlpg "0" "0" "0" "0" "0" "0" "0" ".3*EnjHlp_ATS" "0" "0" BossExp "0" "0" "0" "0" "0" "0" "0" "0" ".3*BossExp_ITS" "0" ExtrMotiv "0" "0" "0" "0" "0" "0" "0" ".3*ExtrMotiv_ATS" "0" "0" Recip "0" "0" "0" "0" "0" "0" "0" ".3*Recip_ATS" "0" "0" TrustB "0" "0" "0" "0" "0" "0" "0" ".3*Trst_ATS" "0" "0" ATS "0" "0" "0" "0" "0" "0" "0" "0" ".3*ATS_ITS" "0" ITS "0" "0" "0" "0" "0" "0" "0" "0" "0" ".3*ITS_KS" KS "0" "0" "0" "0" "0" "0" "0" "0" "0" "0"
The S matrix looks (in part) like this:
$labels KSSE Prestige EnjHlpg BossExp ExtrMotiv Recip TrustB ATS ITS KS KSSE "e1" "cov" "cov" "cov" "cov" "cov" "cov" NA NA NA Prestige "cov" "e2" "cov" "cov" "cov" "cov" "cov" NA NA NA EnjHlpg "cov" "cov" "e3" "cov" "cov" "cov" "cov" NA NA NA BossExp "cov" "cov" "cov" "e4" "cov" "cov" "cov" NA NA NA ExtrMotiv "cov" "cov" "cov" "cov" "e5" "cov" "cov" NA NA NA Recip "cov" "cov" "cov" "cov" "cov" "e6" "cov" NA NA NA TrustB "cov" "cov" "cov" "cov" "cov" "cov" "e7" NA NA NA ATS NA NA NA NA NA NA NA "e8" NA NA ITS NA NA NA NA NA NA NA NA "e9" NA KS NA NA NA NA NA NA NA NA NA "e10"
In the output of the stage 2 analysis, I get parameter estimates for all the paths as expected, and I get values for the error terms e1-e10:
Estimate Std.Error lbound ubound z value Pr(>|z|) KSSE_KS 0.142375 NA 0.038505 0.240741 NA NA KSSE_ITS 0.400807 NA 0.309218 0.495398 NA NA Prest_ATS 0.422256 NA 0.355886 0.488867 NA NA EnjHlp_ATS 0.522933 NA 0.443224 0.602747 NA NA BossExp_ITS 0.437108 NA 0.376581 0.497560 NA NA ExtrMotiv_ATS 0.296702 NA 0.219903 0.373522 NA NA Recip_ATS 0.505664 NA 0.448948 0.562794 NA NA Trst_ATS 0.609971 NA 0.548550 0.671707 NA NA ATS_ITS 0.742128 NA 0.686109 0.799690 NA NA ITS_KS 0.592750 NA 0.540885 0.644920 NA NA e8 0.449247 NA 0.360477 0.529136 NA NA e4 0.808936 NA 0.752439 0.858183 NA NA e3 0.726541 NA 0.636709 0.803560 NA NA e5 0.911968 NA 0.860482 0.951646 NA NA e9 0.648647 NA 0.584078 0.707454 NA NA e10 1.000000 NA NA NA NA NA e1 0.751432 NA 0.701149 0.796394 NA NA e2 0.821700 NA 0.760925 0.873350 NA NA e6 0.744304 NA 0.683333 0.798446 NA NA e7 0.627935 NA 0.548822 0.699112 NA NA
I need r-square values for all three dependent variables (ATS, ITS, and KS). These are not output by default, but I think they can be calculated from some of the output. Except that the key dependent variable (for which I'd really like an r-squared value) is KS (associated with e10), which is being constrained automatically by the metaSEM package.
I can't imagine I'm the first to want to calculate r-square values using metaSEM, but I didn't see any reference to their calculation in the documentation or here in the forums. Any help you can provide will be most useful.
Thanks!