I am fitting a logistic model where the dependent variable is categorical and has say, 5 levels such as 5 countries. If I use dummy variables to identify countries 2-5, I obtain estimates for the threshold for country 1 plus changes in the threshold for each of the remaining 4 countries. Can I fit the same model in a way equivalent to say, the noconstant option in STATA, where instead I get estimates/SE of the 5 country-specific thresholds as part of the output?
I'm thinking more along the lines of the following example:
each subject is coded one if they are ill, 0 otherwise ('outcome'), with a dummy variable dummyUS=1 if from the US (0 otherwise), and dummyUK if from the UK (0 otherwise).
If I fit the logistic model outcome ON dummyUK
the threshold will be -the estimated log odds of being from the US, and the coefficient for dummyUK will the change in the estimated log odds of being ill between UK and the US.
Can I fit a model analogous to the STATA commmand:
logit outcome dummyUS dummyUK, noconstant
so that there will be a coefficient for dummyUS (the estimated log odds of being ill in the US), dummyUK (the estimated log odds of being ill in the UK), and with no additional 'threshold' estimated in the model?
*** FATAL ERROR THE SAMPLE COVARIANCE MATRIX FOR THE INDEPENDENT VARIABLES IN THE MODEL CANNOT BE INVERTED. THIS CAN OCCUR IF A VARIABLE HAS NO VARIATION OR IF TWO VARIABLES ARE PERFECTLY CORRELATED. CHECK YOUR DATA.
Has anyone had any further thoughts on other ways that could get round this in Mplus?
I encounter the same problem as mentioned in the post above (fatal error: covariance matrix cannot be inverted) when including several dummies. However, none of my dummies are perfectly correlated and I have no other independent variables. Do you know what else could be going wrong, or how I could make sure I can run the model without error?