Attachment | Size |
---|---|
joint_ord_cont_model.R [6] | 13.22 KB |
Hi everyone,
I am trying to run a joint ordinal-continuous model.
In my model I have a continuous variable (varA) and an ordinal variable with 1 threshold (varB).
There are 5 zygosity groups (MZM, DZM, MZF, DZF, DOS) and a covariate age which
I modeled separately for men and women.
The model is actually running perfectly fine (no errors) and the estimates for the thresholds for varB
look good and are what I would have expected.
However, the estimates for the means are all way too low when I compare them with the raw data.
Yet the pattern of the means between the groups is similar to the raw data (higher means in males than females).
Does anyone have an idea as to why my means are lower than they should be?
Attached is my script.
Thank you very much in advance!
Jorien