I am using the Environmental Motives Scale in a new population. My sample size is 263.
The results of my exploratory factor analysis showed 2 factors (Eigenvalue>1 and with loadings >0.3) - Biospheric Factor and Human Factor
Cronbach alpha was high for both factors (>0.8)
However, unexpectedly, confirmatory factor analysis showed that the model did not fit well:
RMSEA= 0.126, TLI =0.872 and SRMR = 0.063, AIC = 6786
After a long time on Youtube, I then checked the residual matrix and found that the standardized covariance residuals between two of the items in the Biospheric factor was 7.480. From what I understand if values are >3, it indicates that there may be additional factor/s that are accounting for correlation besides the named factor. I therefore tried covarying the error terms of those two items and rechecked the model fit using CFA.
Results of this model show much better model fit.
RMSEA = 0.083, TLI = 0.945, SRMR = 0.043, AIC = 6731 (not as much difference as I thought there would be)
The questions I am now left with (which google does not seem to have the answer to) are:
1. Is it acceptable to covary the error terms to improve model fit?
2. How does covarying error terms impact on the scoring of the individual scales? Can I still add up the items to measure biospheric vs human scales as I would have without the covarying terms?
I would be so grateful for any insight or assistance.
Thank you
Tabitha