My model consists of a higher-order construct which has two first-order constructs. I am getting good fit in CFA, however on running path analysis, I am not getting such a good fit. Modification indices show that if I correlate error terms between my first-order constructs, I will achieve fit statistics. Can I do that?