Really Fast Bayesian Linear Regression
A little performance vignette. I measured the draw and update of a Bayesian linear regression against how fast one core can actually go, and found which scipy routines get there and which don’t.
6 items found
A little performance vignette. I measured the draw and update of a Bayesian linear regression against how fast one core can actually go, and found which scipy routines get there and which don’t.
My old churn model measured how surprising a user’s silence was against their own history. That’s a real question with a real answer, but it isn’t churn risk: it flagged regulars who took a week off and ignored light users who’d quietly stopped. It passed every backtest for two years; a confused user in a focus group is what caught it. The replacement is the table banks use for late accounts – bin users by how long they’ve been gone and how much they normally use the product, count who came back – with the same empirical Bayes shrinkage as before, pointed at the right quantity this time.
In my previous post, I covered running portable services—version-controlled config, atomic updates, zero-downtime restarts, all without a container runtime. This is only half of the container story, though. You also need to build the service images.
Update, August 2026: I no longer stand behind this approach. The model below measures how surprising a user’s silence is against their own history, and in production that turned out to flag the users least likely to leave. I’ve written up what went wrong and what replaced it in Churn Risk Is a Conditional Rate. The original post is left as written.
We were studying how microglia affect neuronal networks using a standard imaging experiment: 3 mice, 3 coverslips per condition, about 20 neurons measured per coverslip. Our question: Does LPS activation significantly increase PNA signal?
While hierarchical experimental designs are near-ubiquitous in neuroscience and biomedical research, researchers often do not take the structure of their datasets into account while performing statistical hypothesis tests. We present Hierarch, a Python package for analyzing nested experimental designs. Using a combination of permutation resampling and bootstrap aggregation, Hierarch can be used to perform hypothesis tests that maintain nominal Type I error rates and generate confidence intervals that maintain the nominal coverage probability without making distributional assumptions about the dataset of interest.