thresh_load <- measEq.syntax(configural.model =
mod.cat,
data = dat,
ordered = paste0("R09", LETTERS[1:4]),
parameterization = "delta",
ID.cat = "Wu.Estabrook.2016",
group = "IDCNTRY",
group.equal = c("thresholds", "loadings", "means"))
model.thresh_load <- as.character(thresh_load)
fit.thresh_load <- cfa(model.thresh_load, data = dat, group = "IDCNTRY", ordered = paste0("R09", LETTERS[1:4]))