Almost every regression project I've seen in the last decade starts the same way. Load the data, df. corr() , plot the heatmap, scan for red squares, drop one variable from every pair above 0.
Source: [Dev.to](https://dev.to/mohit_modi_e86a932fb11e61/correlation-is-pairwise-multicollinearity-isnt-2970)