[R] inconsistencies between HLM and clmm results for the signs for the coefficients

Thierry Onkelinx th|erry@onke||nx @end|ng |rom |nbo@be
Fri Oct 21 10:02:20 CEST 2022


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 using 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 using 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 using 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]]



More information about the R-help mailing list