Amy D Willis
@amydwillis.bsky.social
530 followers 120 following 62 posts
Biodiversity-loving, error bar-needing statistics nerd; Associate Professor @UWBiostat. Methods & software for #microbiome & #biodiversity data. She/her.
Posts Media Videos Starter Packs
Pinned
amydwillis.bsky.social
I am restarting the StatDivLab blog! Goal is rigorous but accessible statistics for microbiome science. Here's the new scope & welcome: 😽✉️📝

statdivlab.github.io/blog/index.h...

PS. I see you, friends. What a disaster. Please keep resisting, and encourage those around you to do the same.
StatDivLab Blog
This is a way to host a blog for the Statistical Diversity Lab.
statdivlab.github.io
amydwillis.bsky.social
Wisdom from @titus.idyll.org on the final day of #STAMPS2025: "There's free as in beer, and there's free as in kittens. #Bioinformatics software is free as in kittens. You have to love and care for them or they... well... yeah."
amydwillis.bsky.social
The question was "...*on the same sequencing run*?"
amydwillis.bsky.social
Given that a Poisson regression with log link targets the same parameter as your NB regression, I'd be curious to see the coverage of robust Wald CI's. `rigr` wraps this, so does `raoBust`, so should be easy to add.

(Sorry -- I'm at a workshop today or I'd do it myself)
amydwillis.bsky.social
raoBust doesn't invert score tests @nlaroy.bsky.social, but it does implement (model-misspecification) robust score tests, which are amazing for inference.

Feel free to open a feature request. We'll see what we can do.

github.com/statdivlab/r...
GitHub - statdivlab/raoBust: Generalized Linear Models with robust and non-robust Wald and Rao (score) tests
Generalized Linear Models with robust and non-robust Wald and Rao (score) tests - statdivlab/raoBust
github.com
amydwillis.bsky.social
The StatDivLab is getting excited for #STAMPS2025 and working on our lectures for 🤩 Stats Day 😻.

We ran this think-pair-share activity last year and, wow, the ensuing discussion was *very* engaged 😹😅

Reach out if you'd like to attend! Woods Hole, MA, July 14-24 2025. #microbiome #dataanalysis
amydwillis.bsky.social
Agreed that (eg) MAG assembly vs taxonomic estimation makes a huge difference to your answer. Also, please definitely read the Conflict of Interest statement for the shallow shotgun paper and note how much of its claims rely on bioinfomatic subsampling *and not actual shallow sequencing data*
Reposted by Amy D Willis
ilafores.bsky.social
Best workshop you can ever follow
amydwillis.bsky.social
We had a last minute cancellation for STAMPS -- a phenomenal course on microbiome data analysis. Woods Hole, MA, July 14-24 2025.

If you or *anyone* you know is interested to attend, please email Titus &/or me... and we'll do what we can!!!

Thx for sharing widely! 🤞❤️
titus.idyll.org
Hi folks, we had a last minute cancellation at MBL STAMPS, the microbiome data analysis course in Woods Hole that runs from July 14th to July 24th. We can provide housing space for a woman. If you or anyone you know is interested, pls drop me an e-mail at [email protected]. Thx!
amydwillis.bsky.social
That's wonderful to hear, Isabelle!! Thanks for sharing, too. Hope to see you there in 2026.
amydwillis.bsky.social
We had a last minute cancellation for STAMPS -- a phenomenal course on microbiome data analysis. Woods Hole, MA, July 14-24 2025.

If you or *anyone* you know is interested to attend, please email Titus &/or me... and we'll do what we can!!!

Thx for sharing widely! 🤞❤️
titus.idyll.org
Hi folks, we had a last minute cancellation at MBL STAMPS, the microbiome data analysis course in Woods Hole that runs from July 14th to July 24th. We can provide housing space for a woman. If you or anyone you know is interested, pls drop me an e-mail at [email protected]. Thx!
amydwillis.bsky.social
The StatDivLab is *fully* reliant on our NIGMS #MIRA R35 to bring you top-quality statistical methods for microbiome research. #NIH

They are quietly taking MIRAs away from this year's applicants. It's sneaky and scary.

Fight back! Call your reps and tell them to protect science!
peiferlabunc.bsky.social
The fantastically successful NIGMS MIRA R35 "expired" on May 17, after the most recent deadline--these removals are the most aggressive attack on American science in the current series of devastating actions 🧪
gserratomarks.bsky.social
NIH moved up the expiration date for some funding opportunities, effectively *canceling* the next application date without removing it from the standard table. I can't even imagine how many people are still prepping these grants.

👀 List of impacted opportunities:
grants.nih.gov/grants/guide...
Reposted by Amy D Willis
ethanwhite.weecology.org
It only takes 2 minutes to help defend our amazing program officers (and many other federal employees) and allow them to keep serving all Americans instead of responding to the whims of a single person
ljrissler.bsky.social
You have 2 days left to comment (by May 23 at 11:59 pm)!

Go to the link if you don’t want your NSF or NIH program officers who took an oath to the constitution (not a political party) to instead be subject to the whims of a president and fired if they don’t “conform”.
Regulations.gov
www.regulations.gov
amydwillis.bsky.social
How do I describe data with a lot of zeroes? Sparse.

How do I describe data with a lot of variance? High-variance.

How do I describe data where the totals convey complex information about an unknown quantity I care about? (abundance)

I don't. I just state my assumptions.
GitHub - statdivlab/radEmu
Contribute to statdivlab/radEmu development by creating an account on GitHub.
github.com
amydwillis.bsky.social
I hope you're all having a better week than me.

😽😽 7/6
amydwillis.bsky.social
I leave you with the StatDivLab mantra:

1. choose something meaningful to estimate
2. choose a sensible way to estimate it
3. choose tests that control Type 1 error

That's what we will keep doing, even if anonymous reviewers insist on buzzwords.

6/6
amydwillis.bsky.social
Saying "microbiome data is zero-inflated" leads people to seek out "zero-inflated models." Usually, these are bad methods with bad properties. Stay away. 5/6
amydwillis.bsky.social
Why does this matter? This sort of thinking leads biologists to trust estimators based on highly-parametrised parametric models that are (1) surely misspecified and (2) have terrible properties under misspecification. 4/6
amydwillis.bsky.social
Don't get me started on overdispersed, let alone compositional. Microbiome data is none of these, and I'm not new to this field. 3/6
amydwillis.bsky.social
In fact, if you look at blanks and other control data, you see a lot of incorrect detections. There's better evidence that microbiome data is NON-ZERO inflated than zero-inflated.

2/6
amydwillis.bsky.social
Please, #microbiome and #sequencing data are NOT zero-inflated. Let's stop repeating this nonsense. Zero-inflated compared to what?? Those zeroes carry important information about abundance and sequencing depth, and are not "inflated" in any sense. 1/6
amydwillis.bsky.social
How are DOGE & Trump planning to permanently decrease the NIH budget, and how do we fight back? -- by a courageous NIH whistleblower.

Thank you to this brave person for their clear analysis and constructive suggestions. #NIH #DOGE

donmoynihan.substack.com/p/the-nih-bu...
The NIH budget is on a fast track to disaster
An NIH insider explains what Republicans are likely to do next, and what we can do
donmoynihan.substack.com
Reposted by Amy D Willis
loimai.bsky.social
#EBAME 10 Microbial Ecogenomics Workshop! Brest🇫🇷, Oct 11-25, 2025.
A two-week workshop of lectures & tutorials to learn about omics data analysis for microbial ecology and evolution, all in the beautiful Brest Bay!
Apply here until June 1: maignienlab.gitlab.io/ebame
amydwillis.bsky.social
Please take 2 minutes to log that the ACA should be protected and expanded, and that trans healthcare is both the right thing to do and a great investment in Americans and America (or however you want to put it!) 💙💓🤍🤎🖤
threnody.bsky.social
and if you're like "but trump doesn't listen to courts" or "it'll eventually pass," i have this to say:

every single day a trans person can access the lifesaving, essential healthcare they need is worth fighting for.

it's worth court battles & protests & especially 2 minutes of your time.🧵
amydwillis.bsky.social
absolutely!!!!!!! 🦤❤️🦠
amydwillis.bsky.social
We've just released radEmu v2.0.0 🥳🦤😻

`remotes::install_github("statdivlab/radEmu")`

A huge thanks to users for sharing their requests and questions, and to the maintenance team (Sarah and @davidandacat.bsky.social ) for their time and commitment!

Release notes: github.com/statdivlab/r...
A screenshot of the release notes for radEmu v2. Text available at https://github.com/statdivlab/radEmu/releases