ACCLAB / DABEST-python

Data Analysis with Bootstrapped ESTimation
https://acclab.github.io/DABEST-python/
Apache License 2.0
337 stars 46 forks source link

Possibility to do mixed model statistics ? #143

Closed dhuzard closed 1 year ago

dhuzard commented 1 year ago

Hi, thanks a lot for this tool, I love it. But I'm considering adding my GLMM statistics and would like to know if it is something that could be added easily ? (I have experimental design with strong 'cage-effect' among my animals, and would like to control for it as a 'random effect' in the GLMM, also in a repeated-measure design experiment). Thank you very much!

adamcc commented 1 year ago

Hi @dhuzard Thanks for your kind words. We will discuss this in our forthcoming developer meeting. Also, please feel free to fork the repo and do your own assessment of how hard this will be. Cheers

adamcc commented 1 year ago

Hi @dhuzard ,

We just discussed this, and decided that its not something that we could easily do at the moment.

Additional notes from @hwchoi912 "Thank you for your inquiry, and for your recommendation of incorporating linear mixed effects modeling into the DABEST repertoire. It is unfortunately difficult to reconcile the two worlds mainly because DABEST embodies calculation and visualization of effect sizes and their variability based on bootstrapping. The other challenge there is that random effects reported from mixed effects regression models are predicted values (best unbiased linear predictors, or BLUPS), not statistical estimates. Therefore, bootstrapping will not be applicable for the characterization of the “random effect” distributions, whether it be computed at the cage level or at the individual animal level."

Kind regards, Adam