model_i_zeta function

Computes the overall association theta and the row and column effects zeta