" x = np.c_[[1]*tmpdata.shape[0],tmpdata[:,0:p]]\n", " coef = np.matmul(np.matmul(np.linalg.inv(np.matmul(x.T,x)),x.T),tmpdata[:,p])\n", " epsilon = tmpdata[:,p] - np ...