#A | A matrix which stores ZZ' for all variance components |
#phi_start | Staring values for the variance components |
#M_phi | Matrix with VC estimates at each iteration |
#phi | VC estimate from the latest iteration |
#DL | Gradient of the restricted log-likelihood |
#FS | Fisher's Information matrix |
#V | Variance matrix of y |
#P | The projection matrix |
#llh | Restricted log-likelihood |
#beta_hat | Estimates of fixed effects |
#conv_test | Binary variable equal to 1 if the algorithm converges within n_maxiter iterations |