In `loo_shortcut(X, y)` function, acorrding to the ITSL, `mse = ((y[row] - ypred[row]) / (1 - hi))**2` should be right