I was fitting an bivariate ordinal model, after the use of mxFactor(), I made the variables into ordered factors, like this
maxthresh <- 2
mzData$smoking1=mxFactor(mzData$smoking1, levels= c(0:maxthresh))
. After I ran the saturated (or ACE) model, the error appeared. It said that
In model 'Saturated_Fit' 'the thresholds in column 'smoking1' are not in sorted order.
Can anyone tell me what the reason is? Thank you in advance.