Added nonparametric version of bootstrapped_splines

Work in progress, but a step in the right direction. TODO: Account for
decline in N at points in the trial, perhaps by using bootstrapped CI’s
instead of SE’s over smooth.spline() fits. Right now, visuals don’t
align perfectly with the bootstrapped p-values as far as divergences go.