def METHOD_NAME(self, X, y):
    """Solve linear regression analytically."""
    Xb = np.hstack((np.ones((X.shape[0], 1)), X))
    w = np.zeros(X.shape[1])
    z = np.linalg.inv(np.dot(Xb.T, Xb))
    params = np.dot(z, np.dot(Xb.T, y))
    b, w = np.array([params[0]]), params[1:].reshape(X.shape[1], 1)
    return b, w