# to see contents of create_plots function, look at the top-most algorithm's code
baseline, params = baseline_fitter.pspline_asls(y, lam=(1e3, 1e0), p=0.005)
create_plots(y, baseline)