Two-Way factorial Mixed Model usinf Mixed Model Procedure of Henderson
I´ve been trying to duplicate the analysis of the classic machine x operator example given in the paper by McLean, Sanders and Stroup 1991, A unified approach to Mixed Linear Models , The American Statistician, 45 (1): 54-63). Unfortunately, I could not get the solution vector B of fixed effects (machines) and U of random effects (person), and neither the estiamtes of variance components for the Operator, machine x operator and residual variances by REML. Apparently, the current module of Mixed Linear Model of JASP, does not suppport the analysis for such an example of a 2-way factorial model with one randdom effect, using the Mixed Model Procedure (Henderson1963, 1973).
Is there any way that we can program a package for such type of statistical problems?.