# Independent Variables
X_1 = day_weather_df[["TMAX", "TMIN", "PRCP", "SNOW", "SNWD"]].copy()
# Target Variable
y = day_weather_df[["COUNT(*)"]].copy()

# Ordinary least squared linear regression.
sklearnmodel = LinearRegression(fit_intercept = True)
# Fit linear model
sklearnmodel.fit(X_1, y)
# Returns coefficient of determination
r_squared = sklearnmodel.score(X_1, y)