grr me demo.mo

import%20marimo%0A%0A__generated_with%20%3D%20%220.23.6%22%0Aapp%20%3D%20marimo.App(width%3D%22medium%22)%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%23%23%20Importing%20tools%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Aasync%20def%20_()%3A%0A%0A%20%20%20%20import%20micropip%0A%0A%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20import%20pyarrow%20as%20pa%0A%20%20%20%20except%20ModuleNotFoundError%3A%0A%20%20%20%20%20%20%20%20await%20micropip.install(%22pyarrow%22)%0A%20%20%20%20import%20pyarrow%20as%20pa%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20import%20polars%20as%20pl%0A%20%20%20%20import%20hvplot%0A%20%20%20%20import%20hvplot.polars%0A%20%20%20%20import%20hvplot.pandas%0A%20%20%20%20import%20statsmodels.formula.api%20as%20smf%0A%20%20%20%20import%20pandas%20as%20pd%0A%0A%20%20%20%20return%20mo%2C%20np%2C%20pd%2C%20pl%2C%20smf%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%23%23%20A%20function%20to%20generate%20synthetic%20data%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(np%2C%20pd)%3A%0A%20%20%20%20def%20gen_synthetic_grr_data(%0A%20%20%20%20%20%20%20%20n_parts%3D10%2C%0A%20%20%20%20%20%20%20%20n_appraisals%3D3%2C%0A%20%20%20%20%20%20%20%20n_replicates%3D12%2C%0A%20%20%20%20%20%20%20%20var_part%3D1.0%2C%0A%20%20%20%20%20%20%20%20var_tester%3D0.2%2C%0A%20%20%20%20%20%20%20%20var_interaction%3D0.03%2C%0A%20%20%20%20%20%20%20%20var_repeat%3D0.3%2C%0A%20%20%20%20%20%20%20%20grand_mean%3D10.0%2C%0A%20%20%20%20%20%20%20%20random_state%3D42%2C%0A%20%20%20%20%20%20%20%20csv_fname%3D'data.csv'%0A%20%20%20%20)%3A%0A%20%20%20%20%20%20%20%20rng%20%3D%20np.random.default_rng(random_state)%0A%0A%20%20%20%20%20%20%20%20parts%20%20%20%3D%20%5Bf%22P%7Bi%2B1%7D%22%20%20for%20i%20in%20range(n_parts)%5D%0A%20%20%20%20%20%20%20%20testers%20%3D%20%5Bf%22T%7Bi%2B1%7D%22%20for%20i%20in%20range(n_appraisals)%5D%0A%0A%20%20%20%20%20%20%20%20part_effects%20%20%20%20%3D%20%7Bp%3A%20rng.normal(0%2C%20np.sqrt(var_part))%20%20%20%20for%20p%20in%20parts%7D%0A%20%20%20%20%20%20%20%20tester_effects%20%20%3D%20%7Bt%3A%20rng.normal(0%2C%20np.sqrt(var_tester))%20%20for%20t%20in%20testers%7D%0A%20%20%20%20%20%20%20%20interact_effects%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20(t%2C%20p)%3A%20rng.normal(0%2C%20np.sqrt(var_interaction))%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20t%20in%20testers%20for%20p%20in%20parts%0A%20%20%20%20%20%20%20%20%7D%0A%0A%20%20%20%20%20%20%20%20rows%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20t%20in%20testers%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20p%20in%20parts%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20for%20_%20in%20range(n_replicates)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20y%20%3D%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%20grand_mean%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%2B%20part_effects%5Bp%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%2B%20tester_effects%5Bt%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%2B%20interact_effects%5B(t%2C%20p)%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%2B%20rng.normal(0%2C%20np.sqrt(var_repeat))%20%20%20%23%20measurement%20noise%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20rows.append(%7B%22appraisal%22%3A%20t%2C%20%22part%22%3A%20p%2C%20%22y%22%3A%20y%7D)%0A%0A%20%20%20%20%20%20%20%20pd.DataFrame(rows).to_csv(csv_fname)%0A%0A%0A%20%20%20%20return%20(gen_synthetic_grr_data%2C)%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%23%23%20Read%20the%20data%0A%20%20%20%20Read%20the%20data%20and%20define%20other%20important%20quantities!%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(gen_synthetic_grr_data%2C%20pd)%3A%0A%20%20%20%20gen_synthetic_grr_data()%0A%20%20%20%20df%20%3D%20pd.read_csv('data.csv')%0A%20%20%20%20value_col%20%3D%20%22y%22%0A%20%20%20%20appraisal_col%20%3D%20%22appraisal%22%0A%20%20%20%20part_col%20%3D%20%22part%22%0A%20%20%20%20tolerance%20%3D%205.0%0A%20%20%20%20return%20appraisal_col%2C%20df%2C%20part_col%2C%20tolerance%2C%20value_col%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%23%23%20Modify%20the%20data%0A%20%20%20%20Add%20new%20columns%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(appraisal_col%2C%20df%2C%20part_col%2C%20value_col)%3A%0A%20%20%20%20d%20%3D%20df%5B%5Bvalue_col%2C%20appraisal_col%2C%20part_col%5D%5D.dropna().copy()%0A%20%20%20%20d%20%3D%20d.rename(columns%3D%7B%0A%20%20%20%20%20%20%20%20value_col%3A%20%22y%22%2C%20appraisal_col%3A%20%22appraisal%22%2C%20part_col%3A%20%22part%22%0A%20%20%20%20%7D)%0A%20%20%20%20d%5B'appraisal'%5D%20%3D%20d%5B'appraisal'%5D.astype('category')%0A%20%20%20%20d%5B'part'%5D%20%3D%20d%5B'part'%5D.astype('category')%0A%20%20%20%20d%5B'appraisal_part'%5D%20%3D%20(d%5B'appraisal'%5D.astype('str')%20%2B%20%22%3A%22%20%2B%20d%5B'part'%5D.astype('str')).astype('category')%0A%20%20%20%20d%5B'_grp'%5D%20%3D%201%0A%20%20%20%20return%20(d%2C)%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%20Visualize%20the%20data!%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(d%2C%20mo)%3A%0A%20%20%20%20mo.vstack(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20d.hvplot.violin(y%3D'y'%2C%20by%3D'part')%2C%0A%20%20%20%20%20%20%20%20d.hvplot.violin(y%3D'y'%2C%20by%3D'appraisal')%2C%0A%20%20%20%20%5D%20)%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%23%23%20The%20model%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(np%2C%20pd%2C%20smf)%3A%0A%20%20%20%20def%20get_model_smf(d)%3A%0A%20%20%20%20%20%20%20%20model%20%3D%20smf.mixedlm(%0A%20%20%20%20%20%20%20%20%20%20%20%20%22y%20~%201%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20data%3Dd%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20groups%3Dd%5B%22_grp%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20re_formula%3D%220%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20vc_formula%3D%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22part%22%3A%20%220%20%2B%20C(part)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22appraisal%22%3A%20%220%20%2B%20C(appraisal)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22appraisal_part%22%3A%20%220%20%2B%20C(appraisal_part)%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%2C%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20return%20model%0A%0A%20%20%20%20def%20grr_table_from_fit(res%2C%20model%2C%20tol%3DNone)%3A%0A%20%20%20%20%20%20%20%20vcomp%20%3D%20pd.Series(res.vcomp%2C%20index%3Dlist(model.exog_vc.names))%0A%0A%20%20%20%20%20%20%20%20var_repeat%20%20%20%20%20%20%3D%20float(res.scale)%0A%20%20%20%20%20%20%20%20var_part%20%20%20%20%20%20%20%20%3D%20float(vcomp.get(%22part%22%2C%20np.nan))%0A%20%20%20%20%20%20%20%20var_appr%20%20%20%20%20%20%3D%20float(vcomp.get(%22appraisal%22%2C%20np.nan))%0A%20%20%20%20%20%20%20%20var_interaction%20%3D%20float(vcomp.get(%22appraisal_part%22%2C%20np.nan))%0A%0A%20%20%20%20%20%20%20%20var_reprod%20%3D%20var_appr%20%2B%20var_interaction%0A%20%20%20%20%20%20%20%20var_grr%20%20%20%20%3D%20var_repeat%20%2B%20var_reprod%0A%20%20%20%20%20%20%20%20var_total%20%20%3D%20var_grr%20%2B%20var_part%0A%0A%20%20%20%20%20%20%20%20out%20%3D%20pd.DataFrame(%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%22Variance%22%3A%20%7B%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Repeatability%22%3A%20%20%20%20var_repeat%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Appraisal%22%3A%20%20%20%20%20%20%20%20%20%20%20var_appr%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Appraisal%20x%20Part%22%3A%20%20%20%20var_interaction%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Reproducibility%22%3A%20%20var_reprod%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Part-to-Part%22%3A%20%20%20%20%20var_part%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Total%20Gage%20R%26R%22%3A%20%20%20var_grr%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Total%20Variation%22%3A%20%20var_total%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%7D%0A%20%20%20%20%20%20%20%20%7D)%0A%0A%20%20%20%20%20%20%20%20out%5B%22StdDev%22%5D%20%20%20%20%20%20%20%3D%20np.sqrt(out%5B%22Variance%22%5D)%0A%20%20%20%20%20%20%20%20out%5B%22StudyVar_6SD%22%5D%20%3D%206%20*%20out%5B%22StdDev%22%5D%0A%20%20%20%20%20%20%20%20out%5B%22%25Contribution%22%5D%20%3D%20100%20*%20out%5B%22Variance%22%5D%20%2F%20var_total%0A%20%20%20%20%20%20%20%20out%5B%22%25StudyVar%22%5D%20%20%20%20%3D%20100%20*%20out%5B%22StdDev%22%5D%20%2F%20np.sqrt(var_total)%0A%0A%20%20%20%20%20%20%20%20if%20tol%20is%20not%20None%20and%20tol%20%3E%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20out%5B%22%25Tolerance%22%5D%20%3D%20100%20*%20out%5B%22StudyVar_6SD%22%5D%20%2F%20tol%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20out%5B%22%25Tolerance%22%5D%20%3D%20np.nan%0A%0A%20%20%20%20%20%20%20%20return%20out%0A%0A%20%20%20%20return%20get_model_smf%2C%20grr_table_from_fit%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%20Run%20the%20actual%20model%20and%20create%20the%20GRR%20table%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_(d%2C%20get_model_smf%2C%20grr_table_from_fit%2C%20tolerance)%3A%0A%20%20%20%20model_0%20%3D%20get_model_smf(d)%0A%20%20%20%20res_0%20%3D%20model_0.fit(reml%3DTrue%2Cmethod%3D'lbfgs')%0A%20%20%20%20grr_table_0%20%3D%20grr_table_from_fit(res_0%2C%20model_0%2C%20tol%3Dtolerance)%0A%20%20%20%20return%20grr_table_0%2C%20res_0%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%20Bootstrapping%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20d%2C%0A%20%20%20%20get_model_smf%2C%0A%20%20%20%20grr_table_0%2C%0A%20%20%20%20grr_table_from_fit%2C%0A%20%20%20%20mo%2C%0A%20%20%20%20np%2C%0A%20%20%20%20pd%2C%0A%20%20%20%20pl%2C%0A%20%20%20%20res_0%2C%0A%20%20%20%20tolerance%2C%0A)%3A%0A%20%20%20%20n_bootstrap%20%3D%20100%0A%20%20%20%20random_state%20%3D%2042%0A%20%20%20%20boot_records%20%3D%20%5B%5D%0A%20%20%20%20rng%20%3D%20np.random.default_rng(seed%3Drandom_state)%0A%20%20%20%20cell_groups%20%3D%20list(d.groupby(%5B'appraisal'%2C%20'part'%5D%2C%20observed%3DTrue))%0A%20%20%20%20n_fail%3D0%0A%0A%20%20%20%20for%20boot_id%20in%20mo.status.progress_bar(range(n_bootstrap))%3A%0A%20%20%20%20%20%20%20%20resampled_cells%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20cell_key%2C%20cell_df%20in%20cell_groups%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20n_rows%20%3D%20len(cell_df)%0A%20%20%20%20%20%20%20%20%20%20%20%20row_ids%20%3D%20rng.integers(low%3D0%2Chigh%3Dn_rows%2Csize%3Dn_rows)%0A%20%20%20%20%20%20%20%20%20%20%20%20resampled_cell%20%3D%20cell_df.iloc%5Brow_ids%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20resampled_cells.append(resampled_cell)%0A%0A%20%20%20%20%20%20%20%20d_boot%20%3D%20pd.concat(resampled_cells%2C%20ignore_index%3DTrue)%0A%20%20%20%20%20%20%20%20d_boot%5B%22appraisal%22%5D%20%20%20%20%20%20%3D%20d_boot%5B%22appraisal%22%5D.astype(%22category%22)%0A%20%20%20%20%20%20%20%20d_boot%5B%22part%22%5D%20%20%20%20%20%20%20%20%3D%20d_boot%5B%22part%22%5D.astype(%22category%22)%0A%20%20%20%20%20%20%20%20d_boot%5B%22appraisal_part%22%5D%20%3D%20(%0A%20%20%20%20%20%20%20%20%20%20%20%20d_boot%5B%22appraisal%22%5D.astype(str)%20%2B%20%22%3A%22%20%2B%20d_boot%5B%22part%22%5D.astype(str)%0A%20%20%20%20%20%20%20%20).astype(%22category%22)%0A%20%20%20%20%20%20%20%20d_boot%5B%22_grp%22%5D%20%3D%201%0A%20%20%20%20%20%20%20%20try%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20model%20%3D%20get_model_smf(d%3Dd_boot)%0A%20%20%20%20%20%20%20%20%20%20%20%20res%20%3D%20model.fit(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20reml%3DTrue%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20method%3D'lbfgs'%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20start_params%3Dres_0.params%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20%20%20%20%20grr_table%20%3D%20grr_table_from_fit(model%3Dmodel%2Cres%3Dres%2Ctol%3Dtolerance)%0A%20%20%20%20%20%20%20%20%20%20%20%20boot_records.append(grr_table)%0A%20%20%20%20%20%20%20%20except%20Exception%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20n_fail%20%2B%3D%201%0A%0A%0A%20%20%20%20computed_quantities%20%3D%20list(grr_table_0)%0A%20%20%20%20bootstrapped_quantities_dict%20%3D%20%7B%7D%0A%20%20%20%20for%20quantity%20in%20computed_quantities%3A%0A%20%20%20%20%20%20%20%20quantity_series_list%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20tbl%20in%20boot_records%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20quantity_series_list.append(tbl%5Bquantity%5D)%0A%0A%20%20%20%20%20%20%20%20quantity_matrix%20%3D%20pd.concat(quantity_series_list%2C%20axis%3D1).T.reset_index(drop%3DTrue)%0A%20%20%20%20%20%20%20%20quantity_matrix%20%3D%20quantity_matrix.reset_index(names%3D%22bootstrap_id%22)%0A%20%20%20%20%20%20%20%20bootstrapped_quantities_dict%5Bquantity%5D%20%3D%20quantity_matrix%0A%0A%20%20%20%20bootstrapped_quantities_df%20%3D%20pd.concat(%0A%20%20%20%20%20%20%20%20%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20quantity_df.melt(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20id_vars%3D%22bootstrap_id%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20var_name%3D%22component%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20value_name%3D%22value%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20).assign(quantity%3Dquantity)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20quantity%2C%20quantity_df%20in%20bootstrapped_quantities_dict.items()%0A%20%20%20%20%20%20%20%20%5D%2C%0A%20%20%20%20%20%20%20%20ignore_index%3DTrue%2C%0A%20%20%20%20)%5B%5B%22bootstrap_id%22%2C%20%22quantity%22%2C%20%22component%22%2C%20%22value%22%5D%5D%0A%0A%20%20%20%20all_bootstrap%20%3D%20pl.from_pandas(bootstrapped_quantities_df)%0A%20%20%20%20return%20(all_bootstrap%2C)%0A%0A%0A%40app.cell%0Adef%20_(all_bootstrap%2C%20grr_table_0%2C%20mo%2C%20pl)%3A%0A%20%20%20%20import%20holoviews%20as%20hv%0A%20%20%20%20from%20functools%20import%20reduce%0A%0A%20%20%20%20%23%201.%20Initialize%20the%20HoloViews%20extension%20for%20Bokeh%0A%20%20%20%20hv.extension('bokeh')%0A%0A%20%20%20%20quantities%20%3D%20all_bootstrap%5B%22quantity%22%5D.unique().to_list()%0A%20%20%20%20components%20%3D%20all_bootstrap%5B%22component%22%5D.unique().to_list()%0A%0A%20%20%20%20component_tabs%20%3D%20%7B%7D%0A%20%20%20%20component_tabs%5B'Table'%5D%20%3D%20grr_table_0%0A%0A%20%20%20%20for%20c%20in%20components%3A%0A%20%20%20%20%20%20%20%20quantity_plots%20%3D%20%5B%5D%0A%20%20%20%20%20%20%20%20if%20c%20%3D%3D%20%22Total%20Variation%22%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%0A%20%20%20%20%20%20%20%20for%20q%20in%20quantities%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20sub_df%20%3D%20all_bootstrap.filter(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20(pl.col(%22component%22)%20%3D%3D%20c)%20%26%20(pl.col(%22quantity%22)%20%3D%3D%20q)%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20p%20%3D%20sub_df.hvplot.hist(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20y%3D%22value%22%2C%20bins%3D20%2C%20title%3Df%22%7Bc%7D%20-%20%7Bq%7D%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20width%3D300%2C%20height%3D300%2C%20normed%3DTrue%2C%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D'blue'%0A%20%20%20%20%20%20%20%20%20%20%20%20)%20*%20sub_df.hvplot.kde(y%3D'value')%0A%0A%20%20%20%20%20%20%20%20%20%20%20%20quantity_plots.append(p)%0A%0A%20%20%20%20%20%20%20%20plots_layout%20%3D%20hv.Layout(quantity_plots).cols(3)%20%0A%20%20%20%20%20%20%20%20plots_layout.opts(shared_axes%3DFalse)%0A%20%20%20%20%20%20%20%20filename%20%3D%20f%22data%2F%7Bc%7D_plot.html%22%0A%20%20%20%20%20%20%20%20%23hv.save(plots_layout%2C%20filename%2C%20fmt%3D'html'%2C%20resources%3D'inline')%0A%0A%20%20%20%20%20%20%20%20%23%20Add%20to%20marimo%20dictionary%0A%20%20%20%20%20%20%20%20component_tabs%5Bc%5D%20%3D%20plots_layout%0A%0A%20%20%20%20%23%20Create%20the%20final%20tabbed%20interface%0A%20%20%20%20mo.tabs(component_tabs)%0A%20%20%20%20return%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A