I estimated measurement invariance using
measEq.syntax(
parameterization = "theta",
ordered = ...,
ID.cat = "lisrel",
group.equal = c("thresholds", "loadings", "intercepts")
)
By default, it fixes the factor variance in the first group (first year) to 1 for identification.
My question: is it possible to instead freely estimate the factor variances to compare them and identify the model by fixing the loading of the first item to 1 in each year?
In other words, can I change the identification strategy in measEq.syntax, or is it not supported?