why it's not enough to scan values of (simple) Rhat blindly

73 views
Skip to first unread message

Marc Kery

unread,
Jul 3, 2025, 6:58:08 AMJul 3
to 'Google Groups' via hmecology: Hierarchical Modeling in Ecology
Dear all,

here's my cautionary note for the day ...

Sometimes the simple Rhat, as computed by many software including JAGS when running it with jagsUI, can be really misleading as a convergence criterion, as these plots show:


Solution: always check out at least some traceplots. Or (as Matt H. would argue), use more refined versions of Rhat such as split-chain Rhat ...

Best regards  --- Marc




Matthijs Hollanders

unread,
Jul 8, 2025, 6:40:06 AMJul 8
to Marc Kery, 'Google Groups' via hmecology: Hierarchical Modeling in Ecology
Hey Marc, thanks for sharing. How did these traceplots come about?

--
*** Three hierarchical modeling email lists ***
(1) unmarked: for questions specific to the R package unmarked
(2) SCR: for design and Bayesian or non-bayesian analysis of spatial capture-recapture
(3) HMecology (this list): for everything else, especially material covered in the books by Royle & Dorazio (2008), Kéry & Schaub (2012), Kéry & Royle (2016, 2021) and Schaub & Kéry (2022)
---
You received this message because you are subscribed to the Google Groups "hmecology: Hierarchical Modeling in Ecology" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hmecology+...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/hmecology/ZR0P278MB0869FE324D4A65FD2A574310EB43A%40ZR0P278MB0869.CHEP278.PROD.OUTLOOK.COM.

Marc Kery

unread,
Jul 8, 2025, 7:44:42 AMJul 8
to Matthijs Hollanders, 'Google Groups' via hmecology: Hierarchical Modeling in Ecology
Dear Matt,

to be honest, I can't remember :( I think I have never seen this sort of pathology, but it does illustrate that something simple forms of Rhat may get it spectacularly wrong.

Best regards  --- Marc




From: Matthijs Hollanders <matthijs....@gmail.com>
Sent: Tuesday, July 8, 2025 12:39
To: Marc Kery <marc...@vogelwarte.ch>
Cc: 'Google Groups' via hmecology: Hierarchical Modeling in Ecology <hmec...@googlegroups.com>
Subject: Re: why it's not enough to scan values of (simple) Rhat blindly
 
Reply all
Reply to author
Forward
0 new messages