Debbie Hahs-Vaughn
2022-Oct-20 01:51 UTC
[R] inconsistencies between HLM and clmm results for the signs for the coefficients
I've generated a very simple multilevel ordinal model using HLM software and using the clmm function from the ordinal package in R. The threshold values and coefficients are comparable, however the signs for the coefficients (and only the coefficients; not the thresholds) are opposite (e.g., positive in HLM, negative in clmm). This then impacts the odds ratios for the predictors. HLM reports odds greater than one and clmm reports odds less than one. I anticipate there is a very easy reason for what's happening, and it's just escaping me. Can someone please advise on what I've overlooked? [[alternative HTML version deleted]]
Thierry Onkelinx
2022-Oct-21 08:02 UTC
[R] inconsistencies between HLM and clmm results for the signs for the coefficients
Dear Debbie, Have a look at the vignette of the ordinal package. Look for the equation that defines \eta. And note the minus sign associated with \beta. You'll need to find the equation used in HLM and compare it with the equation from ordinal. Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance thierry.onkelinx at inbo.be Havenlaan 88 bus 73, 1000 Brussel www.inbo.be /////////////////////////////////////////////////////////////////////////////////////////// To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher The plural of anecdote is not data. ~ Roger Brinner The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey /////////////////////////////////////////////////////////////////////////////////////////// <https://www.inbo.be> Op do 20 okt. 2022 om 09:22 schreef Debbie Hahs-Vaughn <debbie at ucf.edu>:> I've generated a very simple multilevel ordinal model using HLM software > and using the clmm function from the ordinal package in R. The threshold > values and coefficients are comparable, however the signs for the > coefficients (and only the coefficients; not the thresholds) are opposite > (e.g., positive in HLM, negative in clmm). This then impacts the odds > ratios for the predictors. HLM reports odds greater than one and clmm > reports odds less than one. I anticipate there is a very easy reason for > what's happening, and it's just escaping me. Can someone please advise on > what I've overlooked? > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]