von_Eye_diagonal_linear_by_linear function

Fits the diagonal effects model, where each category has its own parameter delta[k], while also incorporating a linear-by-linear term.