I am not sure how to specify thresholds when > 1, I was hoping for some help. I have run the ACE models with one threshold (binary response data) successfully and have tried to incorporate the correct code for thresholds >1 into my script from several different scripts but have had no luck.
It seems that I am not putting the right information into dimnames? My data contains one variable with 6 levels, for each twin. I'm also not sure how the specified matrices combine and any tips on how that occurs would be greatly appreciated.
I'll attach my script for the saturated model (which I'll then adapt for the model with parameters specified).
p.s. the script may be a little disorganised - I have been playing around with it
|ord5ThrModel.R ||1.88 KB|