








Study with the several resources on Docsity
Earn points by helping other students or get them with a premium plan
Prepare for your exams
Study with the several resources on Docsity
Earn points to download
Earn points by helping other students or get them with a premium plan
Community
Ask the community for help and clear up your study doubts
Discover the best universities in your country according to Docsity users
Free resources
Download our free guides on studying techniques, anxiety management strategies, and thesis advice from Docsity tutors
The application of Structural Equation Modeling (SEM) in neuroimaging research, focusing on creating connectivity models of brain regions, testing their goodness-of-fit, and finding a balance between model complexity, anatomical accuracy, and interpretability. The document also covers the estimation of path coefficients and their interpretation, as well as the comparison of different models.
What you will learn
Typology: Study Guides, Projects, Research
1 / 14
This page cannot be seen from the preview
Don't miss anything!
w
D c B A D yB zC C xA wB B vA A ψ ψ ψ ψ = + + = + + = + = For these equations, A, B, C and D are the known variables (measured covariances); v, w, x, y, and z are the unknown variables. For each region, a separate Ψ variable is included, and these represent the residual influences. Simply stated, this variable can be interpreted as the combined influences of areas outside the model and the influence of a brain region upon itself (McIntosh and Gonzalez-Lima, 1992). The path coefficients are normally computed using software packages such as AMOS, LISREL, and MX32. The starting values of the estimates are initially obtained using two-stage least squares, and they are iteratively modified using a maximum likelihood fit function (Joreskog and Sorbom, 1989). Minimizing the differences between observed and implied covariances is usually done with steepest-descent iterations. The structural equation modeling technique differs from other statistical approaches such as multiple regression or ANOVA where the regression coefficients are obtained from minimizing the sum squared differences between the predicted and observed dependent variables. In structural equation modeling, instead of considering individual observations (or variables) as with other usual statistical approaches, the covariance structure is emphasized. In the context of neural systems, the covariance measure corresponds to how much the neural activities of two or more brain regions are related. Applying structural equation modeling analysis to neuroimaging data has a particular advantage compared to applying it to economics, social sciences or psychology datasets, since the connections (or pathways) between the dependent variables (activity of brain areas) can be determined based on anatomical knowledge and the activity can also be measured directly. With applications in other fields, this is not always true: the models are sometimes hypothetical and cannot be measured directly. Goodness-of-fit Criteria Typically in SEM, statistical inference is used to measure: (1) the goodness of the overall fit of the model, i.e. how significantly different are the observed covariance structure and the covariance structure implied by the anatomical model, and (2) the difference between alternative models for modeling modulatory influence or experimental context by using the nested or stacked model approach. For the purpose of assessing the overall fit of the model, the χ^2 values relative to the degrees of freedom are most widely calculated. This is often referred to as the chi-square test and is an absolute test of model fit. If the p-value associated with the χ^2 value is below 0.05, the model is rejected in absolute fit sense. Because the χ^2 goodness-of-fit criterion is very sensitive to sample size and non-normality of the data, often other descriptive measures of fit are used in addition to the absolute χ^2 test. When the number of samples is greater than a few hundred, the χ^2 test has a high tendency to always show statistically significant results, ensuing in a rejected model. However other descriptive
Nested Model Comparison Since SEM is inherently linear, it cannot directly model non-linear changes in connection strength. However, to overcome this problem, two models can be constructed and these two models can be compared to test for non-linear changes. This is known as the nested (or stacked) model approach (McIntosh, 1998; Della-Maggiore et al., 2000). The first model defined in this approach is the restricted null model, in which the path coefficients are forced to be equal between all conditions and the second model is the corresponding alternate free model, in which the path coefficients are allowed to change between different conditions or subject groups. The χ^2 values are computed for both the null model and the alternative model with corresponding degrees of freedom. If the χ^2 value for the null model is statistically significantly larger than the alternative model, the null model is refuted and the alternative model is assumed to provide a better fit. In other words, different conditions within the free model are deemed to be significantly distinguishable in terms of their path connectivity, and one can infer that there is a statistically significant global difference in path coefficients between the conditions. The χ^2 diff is evaluated with the degrees of freedom equal to the difference in the degrees of freedom for the null and free model. If the χ^2 diff test for the null and free model is found to be statistically significant, one can also use pair-wise parameter comparisons (Arbuckle and Wothke, 1999) to determine which pairs of parameters are significantly different between the experimental conditions in the free model. For the pair-wise parameter comparison test, critical ratios for differences between two parameters in question are calculated by dividing the difference between the parameter estimates by an estimate of the standard error of the difference. Under appropriate assumptions and with a correct model, the critical ratios follow a standard normal distribution.
Results The structural model in Figure 1 was analyzed separately for the Hearing and the Deaf groups, and also separately for the left and the right hemispheres, resulting in four independent models – Hearing (left hemisphere), Hearing (right hemisphere), Deaf (left hemisphere), and Deaf (right hemisphere). In order to investigate if any of the path model’s connection strengths change between the CVCV Visual-only and the CVCV Audio-Visual conditions (i.e. to test for changes in the connection strengths between when the auditory speech information is absent vs. available), multi-group analyses were conducted with the nested models approach. The null (constrained) model’s parameters were restricted to be equal between the two conditions, whereas the free (unconstrained) model’s parameters were allowed to be different for the two separate conditions. Several indices for goodness-of-fit, as discussed in the Methods section, for the four nested models are listed in Table 11 along with their χ^2 statistics for model comparisons. The goodness-of-fit indices indicate that the anatomical model (Figure 1) adequately fits the experimental data for both subject groups and for both hemispheres especially when the models were unconstrained. This implies that our anatomical model is suitably represents a network of cortical regions that may underlie a audio-visual speech processing for both subject groups, while being sensitive enough to the changes in the availability of auditory speech. The χ^2 fit index for the Hearing (Right) model suggested that the absolute fit may not be acceptable (χ^2 (6) = 12.771, P = 0.047) as its p-value is near the borderline cut-off point of P > 0.05, but as stated in the Methods section, other descriptive fit statistics (RMR = 0.013, GFI = 0.998, AFGI = 0.986, RMSEA = 0.024) reflect a good overall fit, hence this model was not rejected in our analyses. The stability index (Fox, 1980; Bentler and Freeman, 1983) was also calculated for each model since our path model includes a nonrecursive subset of regions: AG, FG, STS, IFG, and PMC. As listed in Table 11, both the Hearing and the Deaf right hemisphere models’ estimates were found to be well below one and thus stable. However, the Deaf (Left) models’ stability indices were all greater than one (STI = 2.387, 2.387, 1.131). If the stability index value is greater than or equal to one for any of the nonrecursive subsets of a path model, the parameter estimates are known to yield an unstable system, producing results that are particularly difficult to interpret. Therefore, we decided not to present the parameter estimates from the Deaf (Left) model. All nested models except for the Hearing (Left) model (χ^2 diff = 23.995, df = 15, P = 0.065) showed statistically significant differences across unconstrained and constrained models. Since the Hearing (Left) model did not satisfy the conventional level of significance p < 0.05, its path coefficients should be interpreted with some caution.
X2 P RMR GFI AGFI RMSEA VO AV χ 2 diff df=
Tables 12 and 13 list the estimated path coefficients that minimize the difference between observed and model coefficients for the Hearing (Left and Right) and Deaf (Right) models along with their corresponding standard errors, critical ratios, and p-values for both the CVCV Visual-Only and CVCV Audio-Visual conditions. Here, the estimated path coefficients represent the strength of connections or the strength of the influence conveyed through that pathway. The last two columns of the table list critical ratios for pair-wise parameter differences between the two experimental tasks and their levels of significance. The estimated pathway connection strengths are also summarized graphically in Figures 10, 11 and 12 for the Hearing (Right), Hearing (Left) and the Deaf (Right) models respectively. In these figures, thicker arrows are used to represent pathways with statistically significant pair-wise differences in connection strengths across the two experimental tasks.
FG STS AG IFG FG STS AG IFG
path coefficient Standard Error Critical Ratio P- value Estimated path coefficient Standard Error Critical Ratio P- value Critical Ratio for Difference P- value Deaf (Right) V2 Î AG .392 .038 10.332 *** .294 .043 6.835 *** -2.310 * V2 Î FG .363 .032 11.227 *** .316 .036 8.911 *** -1. V2 Î STS -.101 .082 -1.223 .090 .091 .987 2.425 * FG Î AG .377 .052 7.239 *** .284 .039 7.289 *** -2.101 * AG Î FG .377 .052 7.239 *** .284 .039 7.289 *** -2.101 * AG Î STS .594 .120 4.951 *** .414 .182 2.274 * -1. STS Î AG -.249 .126 -1.977 * -.067 .156 -.430 2.308 * FG Î STS .417 .103 4.060 *** .164 .154 1.064 -1. STS Î FG -.154 .094 -1.629 .013 .129 .100 1. STS Î PMC .375 .144 2.607 ** .208 .206 1.008 -. PMC Î STS -.093 .200 -.467 .104 .298 .350. STS Î IFG .462 .042 11.089 *** .579 .069 8.374 *** 1. IFG Î STS -.203 .073 -2.775 ** -.328 .130 -2.515 * -. IFG Î PMC .302 .021 14.085 *** .248 .031 7.985 *** -1. PMC Î IFG .302 .021 14.085 *** .248 .031 7.985 *** -1. FG Î PMC .026 .058 .439 .046 .048 .960. AG Î PMC -.023 .076 -.298 .168 .090 1.881 2.124 *
V FG STS AG IFG PMC -. . . . . . . . . . . . -. . . . -. . . . -. . . . . . -. -. -. -.
V2 FG IFG PMC AG STS V2 FG IFG PMC AG STS