Ben Bolker's Avatar

Ben Bolker

@bbolker

Ecology, evolution, epidemiology, statistics (mixed models). McMaster University, Hamilton, Ontario https://math.mcmaster.ca/bolker

2,250
Followers
596
Following
817
Posts
07.10.2023
Joined
Posts Following

Latest posts by Ben Bolker @bbolker

Anyone up for this?

10.03.2026 21:00 πŸ‘ 12 πŸ” 2 πŸ’¬ 1 πŸ“Œ 0
Preview
Google Maps Find local businesses, view maps and get driving directions in Google Maps.

I remember such a sign on the road from Brisbane to Townsville (don't remember exactly where it was) πŸ‡¦πŸ‡Ί www.google.com/maps/dir/Bri...

09.03.2026 18:07 πŸ‘ 5 πŸ” 0 πŸ’¬ 0 πŸ“Œ 0
The Sakura Cherry Blossoms in High Park are firmly closed in brown buds. It is still too early to make a definitive prediction on when peak bloom will occur.

Note: Sakura Cherry Blossoms in High Park typically reach peak bloom anytime between late April to early May.

The Sakura Cherry Blossoms in High Park are firmly closed in brown buds. It is still too early to make a definitive prediction on when peak bloom will occur. Note: Sakura Cherry Blossoms in High Park typically reach peak bloom anytime between late April to early May.

Meanwhile in Toronto πŸ‡¨πŸ‡¦ highparknaturecentre.com/cherry-bloss...

09.03.2026 16:51 πŸ‘ 2 πŸ” 0 πŸ’¬ 1 πŸ“Œ 0

Shoulders of giants: Wahba β†’ Wood β†’ {Ross, Pedersen, Miller, Simpson} β†’ the rest of us ...

09.03.2026 16:02 πŸ‘ 3 πŸ” 0 πŸ’¬ 0 πŸ“Œ 0
Post image

While working full-time & raising a child, Grace Wahba earned advanced degrees, completing her PhD at Stanford and becoming the first female faculty member in statistics at the University of Wisconsin-Madison. #womenshistorymonth #statwomen magazine.amstat.org/blog/2026/03... #statssky

09.03.2026 13:19 πŸ‘ 22 πŸ” 11 πŸ’¬ 1 πŸ“Œ 1
For works first published prior to 1978, the term will vary depending on several factors. To determine the length of copyright protection for a particular work, consult chapter 3 of the Copyright Act (title 17 of the United States Code). More information on the term of copyright can be found in Circular 15a, Duration of Copyright, and Circular 1, Copyright Basics.

For works first published prior to 1978, the term will vary depending on several factors. To determine the length of copyright protection for a particular work, consult chapter 3 of the Copyright Act (title 17 of the United States Code). More information on the term of copyright can be found in Circular 15a, Duration of Copyright, and Circular 1, Copyright Basics.

(b) Copyrights in Their Renewal Term at the Time of the Effective Date of the Sonny Bono Copyright Term Extension Act7β€”Any copyright still in its renewal term at the time that the Sonny Bono Copyright Term Extension Act becomes effective shall have a copyright term of 95 years from the date copyright was originally secured.

(b) Copyrights in Their Renewal Term at the Time of the Effective Date of the Sonny Bono Copyright Term Extension Act7β€”Any copyright still in its renewal term at the time that the Sonny Bono Copyright Term Extension Act becomes effective shall have a copyright term of 95 years from the date copyright was originally secured.

Oops, it's more complicated (sorry). life+70 is for publications after 1978. Otherwise ... maybe have to ask someone who has this stuff at their fingertips ...

08.03.2026 19:59 πŸ‘ 2 πŸ” 0 πŸ’¬ 1 πŸ“Œ 0

Coming soon! (end of year of death of author [1957] + 70 years) = 2027-12-31 ?

08.03.2026 19:39 πŸ‘ 39 πŸ” 1 πŸ’¬ 2 πŸ“Œ 0
Preview
Ten Thousand

Congrats, it's useful (if obscure) obXKCD: xkcd.com/1053/ (OK, maybe & vs && isn't as exciting as diet Coke and mentos ...)

07.03.2026 23:56 πŸ‘ 2 πŸ” 0 πŸ’¬ 0 πŸ“Œ 0

Can you ask Claude to redo this on a scale where we don't get negative paylines? (e.g. something like a logistic fit with logit(percentile) as the preditor?)

07.03.2026 03:37 πŸ‘ 2 πŸ” 0 πŸ’¬ 0 πŸ“Œ 0

πŸŽ‰πŸŽ‰πŸŽ‰ Long-awaited: lme4 2.0-1 now handles some simple structured covariance matrices: diagonal ['true' - not just semantic splitting of `(x+y|f)` into `(x|f)` + `(y|f)`], compound symmetric, AR1 (no GPs/factor-analytic/reduced-rank/etc. yet ...)

05.03.2026 21:55 πŸ‘ 12 πŸ” 3 πŸ’¬ 0 πŸ“Œ 0
BMB peeves

it's on my "peeves" list too: bbolker.github.io/bbmisc/peeve...

05.03.2026 21:51 πŸ‘ 5 πŸ” 1 πŸ’¬ 1 πŸ“Œ 0
couverture du livre introduction Γ  la statistique bayΓ©sienne avec R

couverture du livre introduction Γ  la statistique bayΓ©sienne avec R

πŸ“˜ Nouveau livre !

Introduction Γ  la statistique bayΓ©sienne avec le logiciel R

Découvrez pas à pas la statistique bayésienne et sa mise en œuvre dans R : théorème de Bayes, MCMC, priors, régression, GLM/GLMM et comparaison/validation de modèles.

Aux Γ©ditions Quae πŸ™

πŸ“· CrΓ©dit photo : Yann Raulet

05.03.2026 15:32 πŸ‘ 52 πŸ” 17 πŸ’¬ 3 πŸ“Œ 3

What's your Erd"os-Bacon-Sabbath number? (Is it finite?)

05.03.2026 12:33 πŸ‘ 2 πŸ” 0 πŸ’¬ 1 πŸ“Œ 0
MCMC Interactive Gallery

Maybe not "statistical concept", but this is fun for MCMC algorithms chi-feng.github.io/mcmc-demo/ap... (references at chi-feng.github.io/mcmc-demo/)

05.03.2026 00:04 πŸ‘ 11 πŸ” 3 πŸ’¬ 1 πŸ“Œ 0
Preview
MixedModels.jl A Julia package for fitting (statistical) mixed-effects models

A quick plug for juliapackages.com/p/mixedmodels (by Doug Bates and phillipalday.com -- similar capabilities to lme4 but **much** faster

04.03.2026 23:32 πŸ‘ 14 πŸ” 5 πŸ’¬ 1 πŸ“Œ 0
Preview
What happens if the explanatory and response variables are sorted independently before regression? Suppose we have data set $(X_i,Y_i)$ with $n$ points. We want to perform a linear regression, but first we sort the $X_i$ values and the $Y_i$ values independently of each other, forming data set $...

For those who haven't yet seen this classic: stats.stackexchange.com/q/185507/2126 (someone whose manager was insisting that they do this ...)

04.03.2026 13:46 πŸ‘ 33 πŸ” 5 πŸ’¬ 1 πŸ“Œ 4
The Bayesian results imply much higher risk of early collapse than maximum likelihood methods. This difference is due to large probabilities of early collapse for certain parameter values that are plausible in light of the data. Because of simplifying assumptions, these results are not directly applicable to assessment. Nevertheless they imply that maximum likelihood and similar methods based upon point parameter estimates will grossly underestimate the risk of early collapse.

The Bayesian results imply much higher risk of early collapse than maximum likelihood methods. This difference is due to large probabilities of early collapse for certain parameter values that are plausible in light of the data. Because of simplifying assumptions, these results are not directly applicable to assessment. Nevertheless they imply that maximum likelihood and similar methods based upon point parameter estimates will grossly underestimate the risk of early collapse.

This is decision-theoretic rather than stats only, but: Ludwig, Donald. β€œUncertainty and the Assessment of Extinction Probabilities.” Ecological Applications 6, no. 4 (1996): 1067–76. doi.org/10.2307/2269....

03.03.2026 15:52 πŸ‘ 1 πŸ” 1 πŸ’¬ 0 πŸ“Œ 0

My thoughts on this aren't fully cooked, but: this presupposes that non-null results are more interesting than null results (which is mostly true given the way most scientists set up hypotheses but doesn't have to be). Is there a bias-interestingness tradeoff dial we could adjust?

02.03.2026 20:28 πŸ‘ 1 πŸ” 0 πŸ’¬ 1 πŸ“Œ 0

complete separation?

02.03.2026 14:33 πŸ‘ 1 πŸ” 0 πŸ’¬ 1 πŸ“Œ 0

It's good you asked, since I got to go down a rabbit hole/didn't know about him before.

02.03.2026 00:05 πŸ‘ 2 πŸ” 0 πŸ’¬ 0 πŸ“Œ 0
Preview
Ernest Everett Just - Wikipedia

I think your AI overlords are bullshitting/hallucinating. See the other part of thread on en.wikipedia.org/wiki/Ernest_... (who sounds like a true badass BTW ...) [off-topic: is there something I can read to update my mental model of BSky threading, which I don't get at all?]

02.03.2026 00:04 πŸ‘ 2 πŸ” 1 πŸ’¬ 2 πŸ“Œ 0

Sorry, missed alt-text on the second screenshot: "As Just remarked in the symposium this morning, he is interested more in the back than in the bristles on the back and more in eyes than in eye color"

01.03.2026 20:49 πŸ‘ 2 πŸ” 0 πŸ’¬ 0 πŸ“Œ 0
"Embryologist E. E. Just complained that genetics and selection could explain why populations of flies had more or fewer bristles on their backs, but it couldn't explain how a fly constructed its back in the first place (Harrison 1937: 372; Gilbert et al 1996: 361)

"Embryologist E. E. Just complained that genetics and selection could explain why populations of flies had more or fewer bristles on their backs, but it couldn't explain how a fly constructed its back in the first place (Harrison 1937: 372; Gilbert et al 1996: 361)

Post image

Quote is attributed from Amundson R (2005) The Changing role of the embryo in evolutionary thought: Roots of Evo-Devo.
Cambridge University Press, Cambridge (Google books screenshot β†’ Harrison, Ross G. 1937. β€œEmbryology and Its Relations.” Science 85 (2207): 369–74 (JSTOR screenshot)

01.03.2026 20:48 πŸ‘ 4 πŸ” 1 πŸ’¬ 2 πŸ“Œ 0

Think this is what you want? philarchive.org/archive/LIAT... quotes: "I am interested more in the fly’s back than the bristles on its back, and more in its eye than its eye color (E.E. Just)" . Paper has more details, plus see: en.wikipedia.org/wiki/Ernest_...

01.03.2026 20:40 πŸ‘ 3 πŸ” 0 πŸ’¬ 1 πŸ“Œ 0

I get the idea ... I thought the people I tagged might have an idea about the source (both fly people, although I think more on the evo of morphology side than strictly evo-devo ... [my characterization, not theirs])

01.03.2026 19:37 πŸ‘ 2 πŸ” 0 πŸ’¬ 4 πŸ“Œ 0

@idworkin.bsky.social ? @thelonglab.bsky.social ?

01.03.2026 19:14 πŸ‘ 1 πŸ” 0 πŸ’¬ 1 πŸ“Œ 0
Thinking inside the box Dirk Eddelbuettel, R, C++, Rcpp

"newbies" check (pkgs w/o prior releases on CRAN) have a particularly fussy human-administered set of checks (e.g. do all functions have explicitly documented return values?) Also, things like spell-check false positives (can be fixed via dirk.eddelbuettel.com/blog/2017/08... )

01.03.2026 17:54 πŸ‘ 1 πŸ” 1 πŸ’¬ 1 πŸ“Œ 0
Canada’s Fentanyl Czar - Privy Council Office - Canada.ca

don't forget the πŸ‡¨πŸ‡¦ fentanyl czar! www.canada.ca/en/privy-cou... (not really Canada's idea ...)

01.03.2026 00:01 πŸ‘ 6 πŸ” 0 πŸ’¬ 2 πŸ“Œ 0

... in the bottom of a locked file cabinet in a disused lavatory with a sign on the door saying "Beware of Leopard" ...

28.02.2026 00:24 πŸ‘ 4 πŸ” 0 πŸ’¬ 0 πŸ“Œ 0
results of `apropos("R.?[Vv]ersion", ignore.case=FALSE)`: c("getRversion", "R.version", "R.Version", "R.version.string")

results of `apropos("R.?[Vv]ersion", ignore.case=FALSE)`: c("getRversion", "R.version", "R.Version", "R.version.string")

argh.

28.02.2026 00:18 πŸ‘ 4 πŸ” 1 πŸ’¬ 1 πŸ“Œ 0