def runFit(th,dth,i): """ runs a defult SLR fit using a lasso with bs to estimate penalty. """ y = dth[:,i] X = np.sin(th.T-th[:,i]).T enm = SLR2_2.select(X,y) return(enm)
def runFitnEst(th,dth,i): """ runs a defult SLR fit using a lasso, also estimates std err of coef and error of model (ie full solution). """ y = dth[:,i] X = np.sin(th.T-th[:,i]).T solution,enm = SLR2_2.estModel(X,y) return(solution,enm)