Internal core function used in SEL.HAP() to fit mixed models and
test the random effect of a haplotype design matrix z.
RANDOM(z, YFIX, KIN, PAR)If PAR is NULL, returns the estimated parameters (vector).
Otherwise returns a list with likelihood ratio test statistics and
parameter estimates.
Haplotype design matrix (n x r). When PAR is NULL,
z can be NULL and the function only estimates variance
components in the mixed model.
A data.frame or matrix with phenotype in the first column and fixed-effect covariates in the remaining columns.
A list of kinship matrices.
Optional vector of (log) variance components. If NULL, they are estimated by the internal mixed model.