Here's the formula that the video ends with;
import patsy as ps
y, X = ps.dmatrices("n_born ~ wday + yday - 1 + C(month)", df_clean)
print(X[:5])
patsy: categories
Calmcode - patsy: Categories
load plz