Introduction to GLMS

import%20marimo%0A%0A__generated_with%20%3D%20%220.18.2%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22%2C%20app_title%3D%22Introduction%20to%20GLMS%22)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20pandas%20as%20pd%0A%20%20%20%20import%20altair%20as%20alt%0A%20%20%20%20import%20numpy%20as%20np%0A%0A%20%20%20%20import%20statsmodels.api%20as%20sm%0A%20%20%20%20import%20statsmodels.formula.api%20as%20smf%0A%20%20%20%20from%20sklearn.model_selection%20import%20train_test_split%0A%20%20%20%20return%20alt%2C%20mo%2C%20np%2C%20pd%2C%20sm%2C%20train_test_split%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(text%3D%20%22%23%20Generalized%20Linear%20Models%20-%20Seoul%20Rented%20Bikes%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(text%3D%22%22%22%0A%20%20%20%20The%20objective%20of%20this%20notebook%20application%20is%20to%20introduce%20students%20to%20generalized%20linear%20models%20in%20python%0A%0A%20%20%20%20Data%20source%3A%20%3Ccode%3Ehttps%3A%2F%2Farchive.ics.uci.edu%2Fdataset%2F560%2Fseoul%2Bbike%2Bsharing%2Bdemand%3C%2Fcode%3E%0A%20%20%20%20%22%22%22).callout(kind%3D%22info%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20label1%20%3D%20mo.md(text%3D%22**1.%20Click%20the%20select%20to%20see%20the%20variables%20descriptions**%22)%0A%20%20%20%20return%20(label1%2C)%0A%0A%0A%40app.cell%0Adef%20_(column_description)%3A%0A%20%20%20%20column_description%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(%22%22%22%0A%20%20%20%20**2.%20Inspect%20the%20dataframe%20interacting%20with%20the%20ui%20elements**%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20datafetch(mo%2C%20pd)%3A%0A%20%20%20%20%23%23%20fetching%20data%0A%20%20%20%20%40mo.cache()%0A%20%20%20%20def%20fetch_bike_data()%3A%0A%20%20%20%20%20%20%20%20url%20%3D%20%22https%3A%2F%2Fraw.githubusercontent.com%2Fdatagus%2Fglm_exercise_marimo%2Fmain%2Fdata%2FSeoulBikeData.csv%22%0A%20%20%20%20%20%20%20%20bike_df%20%3D%20pd.read_csv(url%2C%20encoding%3D'latin-1')%0A%0A%20%20%20%20%20%20%20%20%23%20Remove%20units%20from%20all%20columns%0A%20%20%20%20%20%20%20%20bike_df.columns%20%3D%20bike_df.columns.str.replace(r'%5C(.*%3F%5C)'%2C%20''%2C%20regex%3DTrue).str.strip()%0A%20%20%20%20%20%20%20%20return%20bike_df%0A%0A%20%20%20%20bike_df%20%3D%20fetch_bike_data()%0A%20%20%20%20bike_df%0A%20%20%20%20return%20(bike_df%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20add_description(bike_df%2C%20mo)%3A%0A%20%20%20%20%23creating%20dictionary%20for%20column%20descriptions%0A%20%20%20%20bike_dict%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Date%22%3A%20%22Date%20the%20bike%20was%20rented%22%2C%0A%20%20%20%20%20%20%20%20%22Rented%20Bike%20Count%22%3A%20%22Count%20of%20bikes%20rented%20at%20each%20hour%22%2C%0A%20%20%20%20%20%20%20%20%22Hour%22%3A%20%22Hour%20of%20the%20day%22%2C%0A%20%20%20%20%20%20%20%20%22Temperature%22%3A%20%22Temperature%20in%20Celsius%22%2C%0A%20%20%20%20%20%20%20%20%22Humidity%22%3A%20%22Percentage%20humidity%22%2C%0A%20%20%20%20%20%20%20%20%22Wind%20speed%22%3A%20%22Wind%20speed%20in%20m%2Fs%22%2C%0A%20%20%20%20%20%20%20%20%22Visibility%22%3A%20%22Visibility%20in%20units%20of%2010%20meters%22%2C%0A%20%20%20%20%20%20%20%20%22Dew%20point%20temperature%22%3A%20%22Dew%20point%20temperature%20in%20Celsius%22%2C%0A%20%20%20%20%20%20%20%20%22Solar%20Radiation%22%3A%20%22Solar%20radiation%20in%20MJ%2Fm%C2%B2%22%2C%0A%20%20%20%20%20%20%20%20%22Rainfall%22%3A%20%22Rainfall%20in%20mm%22%2C%0A%20%20%20%20%20%20%20%20%22Snowfall%22%3A%20%22Snowfall%20in%20cm%22%2C%0A%20%20%20%20%20%20%20%20%22Seasons%22%3A%20%22Season%20(Winter%2C%20Spring%2C%20Summer%2C%20Autumn)%22%2C%0A%20%20%20%20%20%20%20%20%22Holiday%22%3A%20%22Holiday%20vs.%20no%20holiday%22%2C%0A%20%20%20%20%20%20%20%20%22Functioning%20Day%22%3A%20%22NoFunc%20(Non-functional%20hours)%20or%20Fun%20(Functional%20hours)%22%0A%20%20%20%20%7D%0A%0A%20%20%20%20%23%20Create%20dropdown%20with%20all%20columns%20from%20bike_df%0A%20%20%20%20column_selector%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dbike_df.columns.tolist()%2C%0A%20%20%20%20%20%20%20%20value%3Dbike_df.columns%5B0%5D%2C%0A%20%20%20%20%20%20%20%20full_width%3DTrue%0A%20%20%20%20)%0A%20%20%20%20return%20bike_dict%2C%20column_selector%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20selector_description(bike_dict%2C%20column_selector%2C%20label1%2C%20mo)%3A%0A%20%20%20%20%23%20creating%20ui%20to%20check%20the%20column%20descriptions%0A%20%20%20%20selected_column%20%3D%20column_selector.value%0A%0A%20%20%20%20des%20%3D%20mo.md(bike_dict%5Bselected_column%5D)%0A%0A%20%20%20%20column_description%20%3D%20mo.hstack(%5Blabel1%2C%20mo.vstack(%5Bcolumn_selector%2C%20des%5D%2C%20gap%3D0.05)%5D)%0A%20%20%20%20return%20(column_description%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(bike_df%2C%20pd)%3A%0A%20%20%20%20%23changing%20the%20date%20format%20and%20splitting%20them%0A%0A%20%20%20%20bike_df%5B%22Date%22%5D%20%3D%20pd.to_datetime(bike_df%5B%22Date%22%5D%2C%20format%3D%22%25d%2F%25m%2F%25Y%22)%0A%20%20%20%20bike_df%5B'year'%5D%20%3D%20bike_df%5B%22Date%22%5D.dt.year%0A%20%20%20%20bike_df%5B'month'%5D%20%3D%20bike_df%5B%22Date%22%5D.dt.month%0A%20%20%20%20bike_df%5B'month_name'%5D%20%3D%20bike_df%5B%22Date%22%5D.dt.month_name()%0A%20%20%20%20bike_df%5B'weekday'%5D%20%3D%20bike_df%5B%22Date%22%5D.dt.dayofweek%20%2B%201%0A%20%20%20%20bike_df%5B'weekday_name'%5D%20%3D%20bike_df%5B%22Date%22%5D.dt.day_name()%0A%20%20%20%20bike_df%5B'day'%5D%20%3D%20bike_df%5B%22Date%22%5D.dt.day%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(bike_df%2C%20np)%3A%0A%20%20%20%20%23converting%20circular%20data%0A%0A%20%20%20%20%23%20Hour%20(0%E2%80%9323)%0A%20%20%20%20bike_df%5B%22hour_sin%22%5D%20%3D%20np.sin(2%20*%20np.pi%20*%20bike_df%5B%22Hour%22%5D%20%2F%2024)%0A%20%20%20%20bike_df%5B%22hour_cos%22%5D%20%3D%20np.cos(2%20*%20np.pi%20*%20bike_df%5B%22Hour%22%5D%20%2F%2024)%0A%0A%20%20%20%20%23%20Day%20of%20week%20(0%3DMon%20...%206%3DSun)%0A%20%20%20%20bike_df%5B%22dow_sin%22%5D%20%3D%20np.sin(2%20*%20np.pi%20*%20bike_df%5B%22weekday%22%5D%20%2F%207)%0A%20%20%20%20bike_df%5B%22dow_cos%22%5D%20%3D%20np.cos(2%20*%20np.pi%20*%20bike_df%5B%22weekday%22%5D%20%2F%207)%0A%0A%20%20%20%20%23%20Month%20(1%E2%80%9312)%0A%20%20%20%20bike_df%5B%22month_sin%22%5D%20%3D%20np.sin(2%20*%20np.pi%20*%20bike_df%5B%22month%22%5D%20%2F%2012)%0A%20%20%20%20bike_df%5B%22month_cos%22%5D%20%3D%20np.cos(2%20*%20np.pi%20*%20bike_df%5B%22month%22%5D%20%2F%2012)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(bike_df%2C%20mo)%3A%0A%20%20%20%20%23creating%20year%20and%20month%20selectors%0A%20%20%20%20year_selector%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dbike_df%5B%22year%22%5D.unique().tolist()%2C%0A%20%20%20%20%20%20%20%20value%3D2017%2C%0A%20%20%20%20%20%20%20%20full_width%3DTrue%0A%20%20%20%20)%0A%0A%20%20%20%20month_selector%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dbike_df%5B%22month_name%22%5D.unique().tolist()%2C%0A%20%20%20%20%20%20%20%20value%3D%22December%22%0A%20%20%20%20)%0A%0A%20%20%20%20day_selector%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20options%3Dbike_df%5B%22day%22%5D.unique().tolist()%2C%0A%20%20%20%20%20%20%20%20value%3D%221%22%2C%0A%20%20%20%20%20%20%20%20full_width%3DTrue%0A%20%20%20%20)%0A%20%20%20%20return%20day_selector%2C%20month_selector%2C%20year_selector%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(day_selector%2C%20mo%2C%20month_selector%2C%20year_selector)%3A%0A%20%20%20%20%23stacking%20the%20selector%0A%0A%20%20%20%20selectors%20%3D%20mo.hstack(%5Byear_selector%2C%20month_selector%2Cday_selector%5D%2C%20justify%3D%22center%22)%0A%20%20%20%20return%20(selectors%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(text%3D%22%23%23%20**Exploration%20charts**%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20**3.%20Select%20the%20date%20you%20want%20to%20visualize**%0A%20%20%20%20On%20the%20left%20plot%20you%20see%20the%20count%20of%20bikes%20rented%20per%20hour%20in%20a%20single%20day.%20On%20the%20right%20you%20see%20the%20average%20count%20of%20bikes%20rented%20in%20a%20specific%20month%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(bike_df%2C%20day_selector%2C%20month_selector%2C%20year_selector)%3A%0A%20%20%20%20%23creating%20subselected%20dataframe%20based%20on%20the%20selectors%0A%20%20%20%20single_df%20%3D%20(%0A%20%20%20%20%20%20%20%20bike_df.loc%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20(bike_df%5B%22year%22%5D%3D%3Dyear_selector.value)%20%26%20%0A%20%20%20%20%20%20%20%20%20%20%20%20(bike_df%5B%22month_name%22%5D%3D%3Dmonth_selector.value)%20%26%0A%20%20%20%20%20%20%20%20%20%20%20%20(bike_df%5B%22day%22%5D%3D%3Dday_selector.value)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%0A%20%20%20%20)%0A%20%20%20%20return%20(single_df%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(alt%2C%20single_df)%3A%0A%20%20%20%20%23creating%20time%20series%20chart%20based%20on%20the%20selectors%0A%0A%20%20%20%20chart%20%3D%20alt.Chart(single_df).mark_line(point%3DTrue).encode(%0A%20%20%20%20%20%20%20%20x%3Dalt.X(%22Hour%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20axis%3Dalt.Axis(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20labelAngle%3D-45%2C%20%20%23%20Rotate%20labels%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20title%3D'Time'%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20))%2C%23%20%3AT%20tells%20Altair%20this%20is%20a%20temporal%20(time)%20field%0A%20%20%20%20%20%20%20%20y%3D%22Rented%20Bike%20Count%3AQ%22%2C%20%20%23%20%3AQ%20tells%20Altair%20this%20is%20a%20quantitative%20field%0A%20%20%20%20%20%20%20%20tooltip%3D%5B%22Temperature%22%2C%22Solar%20Radiation%22%2C%20%22Rainfall%22%2C%22Snowfall%22%5D%0A%20%20%20%20).properties(%0A%20%20%20%20%20%20%20%20width%3D350%2C%0A%20%20%20%20%20%20%20%20height%3D300%2C%0A%20%20%20%20%20%20%20%20title%3D%22Count%20of%20Rented%20Bikes%20for%20the%20specified%20Date%22%0A%20%20%20%20).configure_axis(%0A%20%20%20%20%20%20%20%20grid%3DFalse%0A%20%20%20%20)%0A%20%20%20%20return%20(chart%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(bike_df)%3A%0A%20%20%20%20%23grouping%20dataframe%20to%20get%20average%20counts%20per%20day%0A%0A%20%20%20%20grouped_df%20%3D%20(%0A%20%20%20%20%20%20%20%20bike_df.copy()%0A%20%20%20%20%20%20%20%20%20%20.groupby(%5B%22Date%22%2C%22year%22%2C%20%22Seasons%22%2C%22month_name%22%2C%20%22day%22%2C%22weekday%22%2C%20%22weekday_name%22%2C%22Holiday%22%5D)%0A%20%20%20%20%20%20%20%20%20%20.mean(numeric_only%3DTrue)%20%20%20%20%20%20%20%23%20mean%20for%20all%20numeric%20columns%0A%20%20%20%20%20%20%20%20%20%20.round(2)%0A%20%20%20%20%20%20%20%20%20%20.reset_index()%0A%20%20%20%20)%0A%0A%0A%20%20%20%20drop_columns%20%3D%20%5B%22Hour%22%2C%20%22month%22%5D%0A%20%20%20%20grouped_df%20%3D%20grouped_df.drop(drop_columns%2C%20axis%3D1)%0A%0A%20%20%20%20discrete_columns%20%3D%20%5B%22Rented%20Bike%20Count%22%2C%22Humidity%22%2C%20%22Visibility%22%5D%0A%20%20%20%20grouped_df%5Bdiscrete_columns%5D%20%3D%20grouped_df%5Bdiscrete_columns%5D.round(0).astype(int)%0A%20%20%20%20return%20(grouped_df%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(grouped_df%2C%20month_selector)%3A%0A%20%20%20%20%23subselecting%20grouped_df%20based%20on%20the%20month%20selector%0A%0A%20%20%20%20grouped_month%20%3D%20(%0A%20%20%20%20%20%20%20%20grouped_df.loc%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20(grouped_df%5B%22month_name%22%5D%3D%3Dmonth_selector.value)%0A%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20)%0A%20%20%20%20return%20(grouped_month%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(alt%2C%20grouped_month%2C%20month_selector)%3A%0A%20%20%20%20%23getting%20the%20timeseries%20average%20counts%20plot%20%0A%0A%20%20%20%20chart_month%20%3D%20alt.Chart(grouped_month).mark_line(point%3DTrue).encode(%0A%20%20%20%20%20%20%20%20x%3Dalt.X(%22Date%3AT%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20axis%3Dalt.Axis(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20labelAngle%3D-45%2C%20%20%23%20Rotate%20labels%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20title%3D'Time'%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20))%2C%23%20%3AT%20tells%20Altair%20this%20is%20a%20temporal%20(time)%20field%0A%20%20%20%20%20%20%20%20y%3D%22Rented%20Bike%20Count%3AQ%22%2C%20%20%23%20%3AQ%20tells%20Altair%20this%20is%20a%20quantitative%20field%0A%20%20%20%20%20%20%20%20color%3Dalt.value(%22purple%22)%2C%0A%20%20%20%20%20%20%20%20tooltip%3D%5B%22Rented%20Bike%20Count%22%2C%22Temperature%22%2C%22Solar%20Radiation%22%2C%20%22Rainfall%22%2C%22Snowfall%22%5D%0A%20%20%20%20).properties(%0A%20%20%20%20%20%20%20%20width%3D350%2C%0A%20%20%20%20%20%20%20%20height%3D300%2C%0A%20%20%20%20%20%20%20%20title%3Df%22Average%20count%20of%20rented%20bikes%20for%20%7Bmonth_selector.value%7D%22%0A%20%20%20%20).configure_axis(%0A%20%20%20%20%20%20%20%20grid%3DFalse%0A%20%20%20%20)%0A%20%20%20%20return%20(chart_month%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(chart%2C%20chart_month%2C%20mo%2C%20selectors)%3A%0A%20%20%20%20mo.hstack(%5B%0A%20%20%20%20%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5Bselectors%2C%20chart%5D%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20heights%3D%22equal%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5B%22%5Cn%22%2Cchart_month%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20align%3D%22center%22%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%5D%2Cwidths%3D%22equal%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20%23mo.md(text%3D%22%23%23%20Splitting%20the%20data%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(bike_df%2C%20pd%2C%20sm%2C%20train_test_split)%3A%0A%20%20%20%20df%20%3D%20bike_df.copy()%0A%0A%20%20%20%20df%20%3D%20df.rename(columns%3D%7B%22Rented%20Bike%20Count%22%3A%22bike_count%22%2C%22Wind%20speed%22%3A%22Wind_speed%22%2C%22Dew%20point%20temperature%22%3A%22Dew_point_temperature%22%2C%22Solar%20Radiation%22%3A%22Solar_Radiation%22%7D)%0A%0A%20%20%20%20y%20%3D%20df%5B%5B%22bike_count%22%5D%5D%0A%20%20%20%20X%20%3D%20df%5B%5B%22Date%22%2C%22Hour%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Temperature%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Humidity%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Wind_speed%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Visibility%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Dew_point_temperature%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Solar_Radiation%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Rainfall%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Snowfall%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22year%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22month%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22month_name%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22weekday%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22day%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22hour_sin%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22hour_cos%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22dow_sin%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22dow_cos%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%22month_sin%22%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%22month_cos%22%0A%20%20%20%20%20%20%20%20%20%20%20%5D%5D%0A%0A%20%20%20%20X%20%3D%20sm.add_constant(X)%0A%20%20%20%20X_train%2C%20X_test%2C%20y_train%2C%20y_test%20%3D%20train_test_split(X%2Cy%2C%20test_size%20%3D%200.25%2C%20random_state%3D42)%20%0A%0A%20%20%20%20train%20%3D%20pd.concat(%5BX_train%2C%20y_train%5D%2C%20axis%3D1)%0A%20%20%20%20test%20%3D%20pd.concat(%5BX_test%2C%20y_test%5D%2C%20axis%3D1)%0A%20%20%20%20return%20test%2C%20train%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20%23mo.md(text%3D%22%23%23%20Putting%20interactive%20predictors%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20%23%20Create%20checkboxes%20for%20variable%20selection%0A%20%20%20%20predictor_names%20%3D%20%5B'Hour'%2C%20'Temperature'%2C%20'Humidity'%2C%20'Wind_speed'%2C%20'Visibility'%2C%20'Dew_point_temperature'%2C%20'Solar_Radiation'%2C%20'Rainfall'%2C%20'Snowfall'%2C%20'year'%2C%20'month'%2C%20'weekday'%2C%20'day'%2C%20'hour_sin'%2C'hour_cos'%2C'dow_sin'%2C'dow_cos'%2C'month_sin'%2C%20'month_cos'%5D%0A%0A%20%20%20%20predictor_selector%20%3D%20mo.ui.array(%5B%0A%20%20%20%20%20%20%20%20mo.ui.checkbox(label%3Dvar%2C%20value%3DFalse)%20for%20var%20in%20predictor_names%0A%20%20%20%20%5D)%0A%20%20%20%20return%20predictor_names%2C%20predictor_selector%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20predictor_selector)%3A%0A%20%20%20%20predictors%20%3D%20mo.vstack(%5B%0A%20%20%20%20%20%20%20%20mo.md(%22**4.Select%20Variables%20for%20Model**%22)%2C%0A%20%20%20%20%20%20%20%20mo.hstack(%0A%20%20%20%20%20%20%20%20%20%20%20%20%5Bpredictor_selector%5Bi%5D%20for%20i%20in%20range(0%2C%20len(predictor_selector))%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20widths%3D%22equal%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20wrap%3DTrue%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%5D)%0A%20%20%20%20return%20(predictors%2C)%0A%0A%0A%40app.cell%0Adef%20_(predictor_names%2C%20predictor_selector)%3A%0A%20%20%20%20%23%20Now%20this%20will%20react%20to%20changes%0A%20%20%20%20def%20put_name(number)%3A%0A%20%20%20%20%20%20%20if%20predictor_selector%5Bnumber%5D.value%3A%0A%20%20%20%20%20%20%20%20%20%20%20a%20%3D%20predictor_names%5Bnumber%5D%20%2B%20%22%20%2B%20%22%0A%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20a%20%3D%20%22%22%0A%20%20%20%20%20%20%20return%20a%0A%20%20%20%20return%20(put_name%2C)%0A%0A%0A%40app.cell%0Adef%20_(predictor_selector%2C%20put_name)%3A%0A%20%20%20%20for%20i%20in%20range(0%2Clen(predictor_selector))%3A%0A%20%20%20%20%20%20%20if%20i%20%3D%3D%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20formula_pred%20%3D%20put_name(i)%0A%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20formula_pred%20%3D%20formula_pred%20%2B%20put_name(i)%0A%0A%20%20%20%20formula_pred%20%3D%20formula_pred%5B%3A-3%5D%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(predictor_names%2C%20predictor_selector)%3A%0A%20%20%20%20coli%20%3D%20%5Bpredictor_names%5Bj%5D%20for%20j%20in%20range(0%2Clen(predictor_selector))%20if%20predictor_selector%5Bj%5D.value%5D%0A%20%20%20%20return%20(coli%2C)%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20%23run_model_btn%20%3D%20mo.ui.run_button(label%3D%22Fit%20Model%22)%0A%20%20%20%20%23run_model_btn%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20%23run_model_btn.value%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(sm%2C%20train)%3A%0A%20%20%20%20def%20glm_model(coli)%3A%0A%0A%20%20%20%20%20%20%20%20y%20%3D%20train%5B%22bike_count%22%5D%0A%20%20%20%20%20%20%20%20X%20%3D%20train%5Bcoli%5D%0A%20%20%20%20%20%20%20%20X%20%3D%20sm.add_constant(X)%0A%20%20%20%20%20%20%20%20%23formula%20%3D%20(f%22%22%22bike_count%20~%20%7Bformula_pred%7D%22%22%22)%0A%20%20%20%20%20%20%20%20family%20%3D%20sm.families.Poisson()%0A%20%20%20%20%20%20%20%20%23model%20%3D%20smf.GLM(formula%2C%20train%2C%20family%3Dfamily)%0A%20%20%20%20%20%20%20%20model%20%3Dsm.GLM(y%2CX%2Cfamily%3Dfamily)%0A%20%20%20%20%20%20%20%20result%20%3D%20model.fit()%0A%20%20%20%20%20%20%20%20return%20result%0A%20%20%20%20return%20(glm_model%2C)%0A%0A%0A%40app.cell%0Adef%20_(coli%2C%20glm_model%2C%20mo%2C%20pd)%3A%0A%20%20%20%20result%20%3D%20glm_model(coli)%0A%0A%20%20%20%20%23%20Create%20tabs%20for%20different%20model%20outputs%0A%20%20%20%20model_tabs%20%3D%20mo.ui.tabs(%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Model%20Statistics%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20pd.DataFrame(%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'Metric'%3A%20%5B'Deviance'%2C%20'Null%20Deviance'%2C%20'Deviance%20Explained'%2C%20'Log-Likelihood'%2C%20'AIC'%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'Value'%3A%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20round(result.deviance%2C2)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20round(result.null_deviance%2C2)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%201%20-%20(result.deviance%20%2F%20result.null_deviance)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20round(result.llf%2C2)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20round(result.aic%2C2)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D)%2C%0A%20%20%20%20%20%20%20%20%22Coefficients%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20pd.DataFrame(%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'Coefficient'%3A%20result.params%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'Std%20Error'%3A%20result.bse%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'z'%3A%20result.tvalues%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'P%3E%7Cz%7C'%3A%20result.pvalues%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'%5B0.025'%3A%20result.conf_int()%5B0%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20'0.975%5D'%3A%20result.conf_int()%5B1%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D)%0A%20%20%20%20%7D)%0A%20%20%20%20return%20model_tabs%2C%20result%0A%0A%0A%40app.cell%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(text%3D%22%23%23%20**Add%20variables%20to%20the%20model**%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(mo%2C%20model_tabs%2C%20predictors%2C%20residuals_histogram%2C%20time_series_chart)%3A%0A%20%20%20%20final_ui%20%3D%20mo.hstack(%0A%20%20%20%20%20%20%20%20%5B%0A%0A%20%20%20%20%20%20%20%20mo.vstack(%5Bpredictors%2C%20model_tabs%5D)%2C%0A%20%20%20%20%20%20%20%20mo.vstack(%5Btime_series_chart%2Cresiduals_histogram%5D)%0A%20%20%20%20%5D%2C%20%20%20%20%20widths%3D%5B1%2C%201%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20final_ui%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20%23%23mo.md(text%3D%22%23%23%20Predicting%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(coli%2C%20result%2C%20sm%2C%20test)%3A%0A%20%20%20%20X_test_with_const%20%3D%20sm.add_constant(test%5Bcoli%5D)%0A%20%20%20%20y_predict%20%3D%20result.predict(X_test_with_const)%0A%20%20%20%20return%20(y_predict%2C)%0A%0A%0A%40app.cell%0Adef%20_(pd%2C%20test%2C%20y_predict)%3A%0A%20%20%20%20predictions_comparison%20%3D%20pd.DataFrame(%7B%0A%20%20%20%20%20%20%20%20'Date'%3A%20test%5B%22Date%22%5D%2C%0A%20%20%20%20%20%20%20%20'Hour'%3A%20test%5B%22Hour%22%5D%2C%0A%20%20%20%20%20%20%20%20'year'%3A%20test%5B%22year%22%5D%2C%0A%20%20%20%20%20%20%20%20'month'%3A%20test%5B%22month_name%22%5D%2C%0A%20%20%20%20%20%20%20%20'day'%3A%20test%5B%22day%22%5D%2C%0A%20%20%20%20%20%20%20%20'Actual'%3A%20test%5B%22bike_count%22%5D%2C%0A%20%20%20%20%20%20%20%20'Predicted'%3A%20y_predict%2C%0A%20%20%20%20%20%20%20%20'Residual'%3A%20test%5B%22bike_count%22%5D%20-%20y_predict%0A%20%20%20%20%7D)%0A%20%20%20%20return%20(predictions_comparison%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(alt%2C%20predictions_comparison)%3A%0A%20%20%20%20%23%20Reshape%20data%20for%20Altair%20(melt%20to%20long%20format)%0A%20%20%20%20plot_data%20%3D%20predictions_comparison.melt(%0A%20%20%20%20%20%20%20%20id_vars%3D%5B'Date'%5D%2C%20%0A%20%20%20%20%20%20%20%20value_vars%3D%5B'Actual'%2C%20'Predicted'%5D%2C%0A%20%20%20%20%20%20%20%20var_name%3D'Type'%2C%0A%20%20%20%20%20%20%20%20value_name%3D'Bike%20Count'%0A%20%20%20%20)%0A%0A%20%20%20%20%23%20Create%20the%20chart%20with%20monthly%20aggregation%0A%20%20%20%20time_series_chart%20%3D%20alt.Chart(plot_data).mark_point().encode(%0A%20%20%20%20%20%20%20%20x%3Dalt.X('yearmonth(Date)%3AT'%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20axis%3Dalt.Axis(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20labelAngle%3D-45%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20title%3D'Month'%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20))%2C%0A%20%20%20%20%20%20%20%20y%3Dalt.Y('mean(Bike%20Count)%3AQ'%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20title%3D'Average%20Bike%20Count'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20axis%3Dalt.Axis(format%3D'.0f'))%2C%0A%20%20%20%20%20%20%20%20color%3Dalt.Color('Type%3AN'%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20scale%3Dalt.Scale(domain%3D%5B'Actual'%2C%20'Predicted'%5D%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20range%3D%5B'%231f77b4'%2C%20'%23ff7f0e'%5D)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20legend%3Dalt.Legend(title%3D'Value%20Type'))%2C%0A%20%20%20%20%20%20%20%20tooltip%3D%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip('yearmonth(Date)%3AT'%2C%20title%3D'Month')%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip('Type%3AN'%2C%20title%3D'Type')%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip('mean(Bike%20Count)%3AQ'%2C%20title%3D'Average%20Count'%2C%20format%3D'.0f')%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20).properties(%0A%20%20%20%20%20%20%20%20width%3D350%2C%0A%20%20%20%20%20%20%20%20height%3D200%2C%0A%20%20%20%20%20%20%20%20title%3D'Monthly%20Average%3A%20Actual%20vs%20Predicted%20Bike%20Rentals'%0A%20%20%20%20).configure_axis(%0A%20%20%20%20%20%20%20%20grid%3DFalse%0A%20%20%20%20)%0A%20%20%20%20return%20(time_series_chart%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(alt%2C%20predictions_comparison)%3A%0A%20%20%20%20residuals_histogram%20%3D%20alt.Chart(predictions_comparison).mark_bar(color%3D'red').encode(%0A%20%20%20%20%20%20%20%20x%3Dalt.X('Residual%3AQ'%2C%20bin%3Dalt.Bin(maxbins%3D50)%2C%20title%3D'Residual')%2C%0A%20%20%20%20%20%20%20%20y%3Dalt.Y('count()'%2C%20title%3D'Frequency')%2C%0A%20%20%20%20%20%20%20%20tooltip%3D%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip('Residual%3AQ'%2C%20bin%3Dalt.Bin(maxbins%3D50)%2C%20title%3D'Residual%20Range')%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20alt.Tooltip('count()'%2C%20title%3D'Count')%0A%20%20%20%20%20%20%20%20%5D%0A%20%20%20%20).properties(%0A%20%20%20%20%20%20%20%20width%3D350%2C%0A%20%20%20%20%20%20%20%20height%3D200%2C%0A%20%20%20%20%20%20%20%20title%3D'Distribution%20of%20Residuals%20(Actual%20-%20Predicted)'%0A%20%20%20%20)%0A%20%20%20%20return%20(residuals_histogram%2C)%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A