Itzik posted on Thursday, March 05, 2015 - 5:43 am
Dear Drs. I have conducted a bi-factor geomin efa, with 3 factors. I wanted to ask you how can I use the output to calculate hierarchical omega coefficient(for the general factor). It is my understanding that the formula is something like this: (squared sum of unstandardized loadings on G)/(squared sum of unstandardized loadings on G + squared sum of unstandardized loadings on S1 + squared sum of unstandardized loadings on S2 + error variance for this item). Is that correct? can I use the 'estimated residual variances' seen in the output as the error for the item?
The rotation should be orthogonal, otherwise you have to include covariance terms among the factors in computing the total variance in the denominator.
Itzik posted on Saturday, March 07, 2015 - 1:10 am
I see. Where can I find the covariance terms among factors in the output? Do I have to also calculate their squared sum?
Itzik posted on Saturday, March 07, 2015 - 2:58 am
Let me be more specific - in the oblique rotation, the correlation between the two specific factors is -0.169. The program does not provide correlation with the general factor. Can I use this correlation to caclulate the covariance and add it to the denominator?
And one last question, if its okay - when adding the residual variances to the denominator - should I calculate the squared sum of variances as given, or should I calculate the squared sum of the square root of residual variances?
Ok. I think another approach would help me more. I can simply calculate the denominator by summing the correlation/covariance matrix, right? So the only question is whether I should use the correlation or covariance matrix? In any case - how can I ask Mplus to produce the matrices used for the EFA (geomin bi-factor efa with WSLMV estimator).
I'm sorry for my troubles understanding, but I need a method to calculate the denominator which is the total variance in standardized units (as the actual V(x) is in unstandardized units, being close to the sum of the covariance matrix). So how can I get an output of the correlation matrix which is used by the WSLMV estimator? Alternatively, how can I use the formula you wrote above, to calculate the variance based on standardized loadings?
I am using bi-factor CFA models for a recent analysis and need to compute omega, omegaH and omegaS. In this thread a user posted about using the unstandardized factor loadings to compute omega coefficients, but other examples I have seen all use standardized loadings in the computation (e.g., Reise, Bonifay, & Haviland, 2012; Gignac, 2014; Watkins, 2013).
Is it appropriate to use either loading, or should standardized loadings be used?
Using the unstandardized is correct but perhaps you get the same result using standardized.
Daniel Lee posted on Tuesday, June 06, 2017 - 12:56 pm
Hi Dr. Muthen,
In a Bifactor CFA model (of an anxiety scale), if the OmegaHS was .87 for the general factor and ranged from .01 -.06 (across 4 subtypes), ECV was close to .80, and the percent of uncontaminated correlation was close to 100, does that suggest that the scale is unidimensional and that I should interpret it as such?
If this scale is truly unidimensional, why might the model fit for the one-factor model be so poor (e.g., RMSEA > .15)? Also items loaded well on subtype factors (factor loading > .3) even after loading on the general factor, which further confuses me as to whether this scale is uni- or multidimensional.
As always, thank you for sharing your knowledge and wisdom!