2-way Bayesian repeated measures ANOVA

Discuss statistics related things
Post Reply
Philip2023
Posts: 3
Joined: Mon May 08, 2023 5:37 pm

2-way Bayesian repeated measures ANOVA

Post by Philip2023 »

Hi,

I noticed strange results with a Bayesian RM ANOVA with two factors (2 and 11 levels, respectively). The attachment shows the results, compared with those obtained with JASP, and the are major discrepancies. It looks like something is not working properly in the JSQ module in Jamovi. :thinking:

Cheers,

Phil
Attachments
jamovi-jasp.jpg
jamovi-jasp.jpg (384.24 KiB) Viewed 3469 times
User avatar
jonathon
Posts: 2625
Joined: Fri Jan 27, 2017 10:04 am

Re: 2-way Bayesian repeated measures ANOVA

Post by jonathon »

hi, can you attach a .omv file? (you may need to zip it up first).

with thanks
Philip2023
Posts: 3
Joined: Mon May 08, 2023 5:37 pm

Re: 2-way Bayesian repeated measures ANOVA

Post by Philip2023 »

Hi, I'm attaching the .omv file (zipped). The Bayesian ANOVA is at the bottom of the Results section. IRRPRO should have a BF10 around 2 or so, but it's 10^7 ...

Thanks much for looking into this.

-Phil
Attachments
CK-300-350frontal.zip
(71.78 KiB) Downloaded 136 times
User avatar
jonathon
Posts: 2625
Joined: Fri Jan 27, 2017 10:04 am

Re: 2-way Bayesian repeated measures ANOVA

Post by jonathon »

ah, wait, i remember this one. the jasp guys changed the RM ANOVA implementation. we use the RM ANOVA implementation from the BayesFactor R package. the two implementations are different, and i take the fact that BayesFactor isn't updating it's implementation as a sign that they think the original procedure is sound.

so i think this might be different stats folks disagreeing on which procedure is the best for a given design (which is fairly standard).

jonathon
User avatar
reason180
Posts: 276
Joined: Mon Jul 24, 2017 4:56 pm

Re: 2-way Bayesian repeated measures ANOVA

Post by reason180 »

Yes, it may be fairly standard that "different stats folks [disagree] on which procedure is the best for a given design," but while there is strong consensus about how to implement a classic t test, ANOVA, or regression (you tend to get exactly the same results regardless of which stats package you use), this is not the case for Bayesian analysis. In my opinion, the reason is that the application of Bayesian analysis to non-binomial data and vague-ish hypotheses is somewhat forced, leading to lots of complexity, and lots of freedom to shop around for an implementation that happens to provide a desirable results. For this reason, while I somewhat-trust Bayesian versions of *simple* classical analyses such as t test, my trust breaks down for more complicated analyses such as ANOVAs.
Last edited by reason180 on Thu May 11, 2023 4:54 pm, edited 1 time in total.
Philip2023
Posts: 3
Joined: Mon May 08, 2023 5:37 pm

Re: 2-way Bayesian repeated measures ANOVA

Post by Philip2023 »

Very helpful, and something that I hadn't heard much about, not being a Bayesian analysis expert. I had no idea that different methods could lead to such a wild range of results: from a BF=2 (consistent with the standard NHT ANOVA results) to a BF>10^7 :/

Is the moral that we just shouldn't rely on Bayesian analyses for RM ANOVAs (at least, until there is some consensus)?

Thanks!

-Phil
User avatar
jonathon
Posts: 2625
Joined: Fri Jan 27, 2017 10:04 am

Re: 2-way Bayesian repeated measures ANOVA

Post by jonathon »

that is quite a large discrepancy ... maybe as a sanity check, compare descriptive plots to make sure you're comparing apples with apples.

jonathon
Post Reply