Advertisement · 728 × 90

Posts by Rob Davies

Preview
Models as Prediction Machines: How to Convert Confusing Coefficients Into Clear Quantities - Julia M. Rohrer, Vincent Arel-Bundock, 2026 Psychological researchers usually make sense of regression models by interpreting coefficient estimates directly. This works well enough for simple linear model...

Good news everyone 🥳 Our (w @vincentab.bsky.social) primer on models as prediction machines (with the marginaleffects package) is finally officially published!>

journals.sagepub.com/doi/10.1177/...

13 hours ago 264 104 8 6
Close up with a beautiful nudibranch on a rocky reef. The lower half of the image is a palette of swatches made from the nudi's vibrant colours.

Close up with a beautiful nudibranch on a rocky reef. The lower half of the image is a palette of swatches made from the nudi's vibrant colours.

Once upon a time I just admired nudis and sea slugs.

But it was not enough. Now apparently I'm creating an R package to celebrate their colour palettes? 😅

First up, my Sydney fave, Hypselodoris bennetti.

#rstats #nudibranch #dataviz 🦑🐙🧪 #marinelife #invertebrates

3 weeks ago 1342 324 39 60

Nice clear presentation of commonplace meta-analysis failure mode: Pooling coefficients that mean different things, because original models had different adjustment sets. A coefficient gets its meaning from the whole model, not just from the predictor variable it multiplies.

1 week ago 82 17 1 3
Top down view of a Bullina bubble snail exploring algal turf on a shallow rocky reef. Two tiny black dots for eye-spots behind its headshield. Neon blue edge, a pale blue body, cherry-red stripes on a tiny white shell. What's not to love. Lower half of the image is swatches of the Bullina's colours.

Top down view of a Bullina bubble snail exploring algal turf on a shallow rocky reef. Two tiny black dots for eye-spots behind its headshield. Neon blue edge, a pale blue body, cherry-red stripes on a tiny white shell. What's not to love. Lower half of the image is swatches of the Bullina's colours.

Last of the palettes included in the R package, but certainly not least, the beautiful Bullina bubble-baby. Tiny, friend-shaped, perfect.

#rstats #dataviz 🦑🐙🧪 #marinelife #invertebrates

20 hours ago 67 14 0 1

Brief fun survey from Jessica, Andrew & myself:

If you are a faculty member, research scientist, postdoc, or senior Ph.D. student in any area of science, please take five minutes and fill it out. We’ll share the results widely along with some reflections.

docs.google.com/forms/d/e/1F...

1 week ago 52 52 6 5

Super interesting reading the definitions for the candidate epistemic virtues for science inquiry

(have narrowed my picks but still ranking them...& not sure when this survey closes...)

2 days ago 2 5 1 0

holy shit

2 days ago 462 66 9 0
Advertisement
Post image Post image

Bird nests made from anti-bird spikes?! 🤯

Hi, I'm a nest researcher 👋 and new here on BlueSky, sharing the craziest #bird nests I've ever found. 👀 Today, I’m sharing my discovery of rebellious birds that build nests out of anti-bird spikes. And honestly, it's like telling a joke...

A thread. 🧵

1 year ago 2373 780 62 210
Magpies’ nest made from anti-bird spikes and a strip of the spikes (bottom right). Photograph: Auke-Florian Hiemstra / Naturalis Biodiversity Center
It is an impressively large nest

Magpies’ nest made from anti-bird spikes and a strip of the spikes (bottom right). Photograph: Auke-Florian Hiemstra / Naturalis Biodiversity Center It is an impressively large nest

Magpies are absolute superstars for this artform:

Magpies’ nest made from anti-bird spikes and a strip of the spikes (bottom right). Photograph: Auke-Florian Hiemstra / Naturalis Biodiversity Center

3 days ago 70 11 3 2
Preview
The decision of a Prime Minister Only one person decided to appoint Lord Mandelson

'When the current Prime Minister of the United Kingdom makes a bad decision it always seems that others must take the blame.' Excellent summary of the increasing political weakness of Starmer's position over Mandelson @davidallengreen.bsky.social emptycity.substack.com/p/the-decisi...

3 days ago 59 22 9 1
Preview
Copyediting 101 [a homily]

"filth"

Leo FTW

5 days ago 144 20 4 0
Wildlife trade drives animal-to-human pathogen transmission over 40 years

Wildlife trade drives animal-to-human pathogen transmission over 40 years

New in @science.org ‼️ In the most comprehensive study to date, we show that wildlife trade is driving animal-to-human zoonotic spillover at a planetary scale, with +1 spillover per host every 10 years. Live animal markets and illegal trade pose even greater risks. 🔓 www.science.org/doi/10.1126/...

1 week ago 685 353 10 23
Video

TRUMP, throwing insults at Pope Leo:
“weak on crime”
“terrible for foreign policy”
“catering to the Radical Left”
“owes his position to me”
“sit down and mind your own business”

POPE LEO, unbothered:
“We are called to love.”

6 days ago 3358 965 95 83

I'll send the link when I'm at my desk 👍. It's really not my field but people in the disability sphere have got me thinking (better late than never I suppose). The thing that really hit home is the adoption of copilot at my university in the name of "productivity". I mean, define productivity!

6 days ago 3 1 2 0
Advertisement
An interactive OJS playground demonstrating a linear congruential generator (LCG) using the formula X_n = (aX_{n-1} + c) mod m. Controls on the left set modulus (m=8), multiplier (a=5), increment (c=3), seed (X_0=1), and numbers to generate (12). A table on the right shows the resulting sequence of X values, intermediate calculations, mod m results, and normalized values X_n/m, with the final "random" numbers highlighted in yellow.

An interactive OJS playground demonstrating a linear congruential generator (LCG) using the formula X_n = (aX_{n-1} + c) mod m. Controls on the left set modulus (m=8), multiplier (a=5), increment (c=3), seed (X_0=1), and numbers to generate (12). A table on the right shows the resulting sequence of X values, intermediate calculations, mod m results, and normalized values X_n/m, with the final "random" numbers highlighted in yellow.

Excerpt from the blog post with R code that tests all seeds from 1 to 10,000 to find which ones produce 10 heads in a row when simulating coin flips. The possible_seeds data frame is filtered to show 10 seeds (614, 1667, 3212, 4166, 4580, 5527, 5824, 7365, 7468, 8975) that meet this criterion. The post notes that seed 614 actually produces 13 heads in a row, confirmed with a withr::with_seed(614, ...) call below.

Excerpt from the blog post with R code that tests all seeds from 1 to 10,000 to find which ones produce 10 heads in a row when simulating coin flips. The possible_seeds data frame is filtered to show 10 seeds (614, 1667, 3212, 4166, 4580, 5527, 5824, 7365, 7468, 8975) that meet this criterion. The post notes that seed 614 actually produces 13 heads in a row, confirmed with a withr::with_seed(614, ...) call below.

R console output demonstrating that set.seed(1234) produces reproducible results. The first block calls runif(5) and returns five values: 0.1137, 0.6223, 0.6093, 0.6234, 0.8609. The second block uses the same seed but splits the draw into runif(2) then runif(3), returning the same five values in the same order, showing that the sequence is preserved regardless of how many numbers are drawn at a time.

R console output demonstrating that set.seed(1234) produces reproducible results. The first block calls runif(5) and returns five values: 0.1137, 0.6223, 0.6093, 0.6234, 0.8609. The second block uses the same seed but splits the draw into runif(2) then runif(3), returning the same five values in the same order, showing that the sequence is preserved regardless of how many numbers are drawn at a time.

Table of contents for the post:

Introduction
Seeds and reproducible randomness
My (somewhat incorrect) mental model of how seeds work
Making “random” numbers with an equation
    Live interactive playground
    Cycles and fancier algorithms
Why does it matter if “random” numbers aren’t actually random?
    You’re limiting yourself to narrow, known universes
    You can seed hack and get any values you want
    Real world bad things can happen because of pseudorandom numbers
Can computers even create true randomness?
    Moving a mouse around
    Lava lamps
    Atmospheric noise
How I use true randomness in my own work
“…as an ook cometh of a litel spyr…”

Table of contents for the post: Introduction Seeds and reproducible randomness My (somewhat incorrect) mental model of how seeds work Making “random” numbers with an equation Live interactive playground Cycles and fancier algorithms Why does it matter if “random” numbers aren’t actually random? You’re limiting yourself to narrow, known universes You can seed hack and get any values you want Real world bad things can happen because of pseudorandom numbers Can computers even create true randomness? Moving a mouse around Lava lamps Atmospheric noise How I use true randomness in my own work “…as an ook cometh of a litel spyr…”

I've been using random seeds for years but I have no idea how they work. Seeds somehow(?) make the same random numbers?

So I figured it out! New post includes an interactive PRNG generator, lava lamps, lottery fraud, @random.org, Chaucer, and Minecraft #rstats

www.andrewheiss.com/blog/2026/04...

1 week ago 100 23 6 3

Looks fantastically useful

1 week ago 1 0 0 0
Join Lexical Norms to Your Word List Your word list should ideally be nested in Tidy format (one word per row, within one column of a dataframe). Your word vector should NOT be a factor but a chr. Set up your word list like this. You can split/unlist a language sample to get it in this format also.

How to join zillions of lexical norms to each word in your language sample the easy way: a quick tutorial and demo reilly-lab.github.io/Jamie_JoinLe...

1 week ago 17 9 1 0

If your Bayesian solution ends up with the same uncertainty under both scenarios you are maybe not taking your generative mechanism seriously enough. I (by which I mean Lauren) wrote about this once in a paper that didn’t get as much love as the worse one I wrote before arxiv.org/abs/1905.10341

1 week ago 7 1 1 1
cover of the book "Bayesian Workflow" by Gelman, Vehtari, et al. Coming out later this year, in the summer probably.

cover of the book "Bayesian Workflow" by Gelman, Vehtari, et al. Coming out later this year, in the summer probably.

I would have preferred to have the "draw the rest of the owl" meme on the cover, but this will do. Seems like it is on schedule, and we'll leave some typos so you know we didn't write it with AI.

2 weeks ago 376 57 12 8

More details about the Bayesian Workflow book and case studies now available on the book web site avehtari.github.io/Bayesian-Wor... (but you still need to wait a bit for the book)

2 weeks ago 98 28 2 0

- Replicability is not estimated consistently across studies. Different estimates have different properties and interpretations. Estimates vary greatly depending on the choice of replicability metric.
- Replications do not track truth but I've written a lot about this elsewhere so I'll stop here.

10 months ago 44 1 1 0
Preview
Scientific discovery in a model-centric framework: Reproducibility, innovation, and epistemic diversity Consistent confirmations obtained independently of each other lend credibility to a scientific result. We refer to results satisfying this consistency as reproducible and assume that reproducibility i...

I could recommend these four papers of ours for starters:
1. journals.plos.org/plosone/arti...
2. royalsocietypublishing.org/doi/abs/10.1...
3. royalsocietypublishing.org/doi/abs/10.1...
4. philsci-archive.pitt.edu/24738/

10 months ago 24 8 1 1
Advertisement

OK for fun: top funniest tweet ever? I am torn between "moon's haunted" and "Denise I was at your wedding."

2 weeks ago 6463 773 1593 4415
Video

I built a Quarto extension that adds accessibility features to Reveal.js presentations: keyboard nav, screen reader support, dyslexia-friendly fonts, high contrast, slide change visual/audio cue and more.

Feedback more than welcome.

github.com/mcanouil/qua...

#Quarto #RevealJS #Accessibility

3 weeks ago 30 6 1 1

Part 2 of my shrinkage estimator series is out! Part 1 covered the univariate case, but now we dive into multivariate shrinkage 🤓

We cover Spearman's classic correlation disattenuation formula, multivariate James-Stein estimators, and hierarchical methods too

haines-lab.com/post/how-to-...

3 weeks ago 42 15 2 3

This is basically how my personal development meetings go. I say to my students: there’s room in the world for everyone. There should be

3 weeks ago 2 0 0 0
Photograph of a car with a squirrel inside of it, perched on the steering wheel. The car is the color of champaign at a beige convention, and the squirrel is the color of squirrels. The squirrel is holding a package of crackers in it's mouth. They are the kind like you get at a restaurant, where you get two crackers wrapped in plastic.

The driver's side window is slightly cracked. This is how the squirrel got in, and how it got out. It threw the crackers out first, and then climbed out after them. Everything in this operation suggested that this was not the squirrel's first rodeo.

Photograph of a car with a squirrel inside of it, perched on the steering wheel. The car is the color of champaign at a beige convention, and the squirrel is the color of squirrels. The squirrel is holding a package of crackers in it's mouth. They are the kind like you get at a restaurant, where you get two crackers wrapped in plastic. The driver's side window is slightly cracked. This is how the squirrel got in, and how it got out. It threw the crackers out first, and then climbed out after them. Everything in this operation suggested that this was not the squirrel's first rodeo.

A closeup of the squirrel sitting on the steering wheel. The squirrel deserves a name, so we'll call her Anjeloma, and she's what you might call a winner. She is still squirrel colored. The crackers are white, and labeled "Zest." As if Anjeloma needed more zest. Squirrel, please.

You can't see much of the car, but you can see smudges of grunge at the edges of the windshield, where the wipers have cast aside the debris of previous rains and pollen-falls.

A closeup of the squirrel sitting on the steering wheel. The squirrel deserves a name, so we'll call her Anjeloma, and she's what you might call a winner. She is still squirrel colored. The crackers are white, and labeled "Zest." As if Anjeloma needed more zest. Squirrel, please. You can't see much of the car, but you can see smudges of grunge at the edges of the windshield, where the wipers have cast aside the debris of previous rains and pollen-falls.

The world is stupid, but I just watched a squirrel break into a car in the parking lot below me, steal a package of crackers, and escape to a nearby tree. So at least somebody is winning.

3 weeks ago 9068 1593 149 126

Great to see this! Learn more about the different ways to easily reuse our data, with our Chart Data API and enhanced data downloads: ourworldindata.org/easier-to-re...

3 weeks ago 45 7 0 2
Video

Code Window 1.0.0 for Quarto is here.

Add window-style decorations to your code blocks: macOS traffic lights, Windows buttons, or a plain title bar.

Hotfixes: Filename and code-annotations for Typst/PDF in Quarto 1.9.

github.com/mcanouil/quarto-code-window

#Quarto #Typst #TechnicalWriting

3 weeks ago 15 4 0 0

For those deterred by the (free) registration link, here's the full piece:

3 weeks ago 236 105 7 2
Advertisement