bkreg = wk.BKRegression(XDATA, YDATA, method='wilson')