For a few years now I’ve been chipping away at Regression Modeling Strategies by Frank Harrell. Between the exposition, the… Read more
Would you rather
Someone texted me the following image and asked for my response. So what would I rather do? Flip 3 coins… Read more
Fixing the p-value note in a HTML stargazer table
If you insert a stargazer table into R Markdown that outputs an HTML file, you may get a p-value note… Read more
Most useful tests for an ANCOVA model
In his book Regression Modeling Strategies, 2nd Ed, Frank Harrell provides a list of what he calls the “most useful… Read more
Using natural splines in linear modeling
Take a look at this scatterplot: It's clear there is a relationship between x and y, but the relationship is… Read more
Recreating a Geometric CDF plot from Casella and Berger
For reasons I no longer remember I decided a while back to reproduce Figure 1.5.2 from Casella and Berger (page… Read more
Revisiting some old R-code
One of my first blog posts on this site simulated rolling a die 6 times and observing if side i… Read more
Profile likelihood ratio confidence intervals
When you fit a generalized linear model (GLM) in R and call confint on the model object, you get confidence… Read more
Earthquake data and Benford’s Law
Much has been written about Benford's Law, that weird phenonmenon where if you have a naturally occuring set of numerical… Read more
Recreating a graph from a NCHS data brief
I was reading this sobering article about the rise of suicide rates in the United States and decided to read… Read more