patsy: categories

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])