# where the following magic happens (the red line)
mod_fit = lm_mod.fit(X, Y)