This function is called internally by plmm()
estimate_eta(n, s, U, y, incpt_flag)a numeric value with the estimated value of eta, the variance parameter
The number of observations
The non-zero eigenvalues of K, the realized relationship matrix
The eigenvectors of K associated with s
Continuous outcome vector
Logical: Does the model require fitting an intercept?