def predict(tau,model,xT,yT): err = sp.zeros(tau.size) for j,t in enumerate(tau): yp = model.predict(xT,tau=t)[0] eq = sp.where(yp.ravel()==yT.ravel())[0] err[j] = eq.size*100.0/yT.size return err