Figure 2.63 shows the report obtained for a fit to the Investment Castings.jmp sample data using the REML method. Run the script Model - REML, and then fit the model. Note that Casting is a random effect and is nested within Temperature.
Figure 2.63 Fit Least Squares Report for REML Method
For each term in the model, this report gives an empirical estimate of its best linear unbiased predictor (BLUP) and a test for whether the corresponding coefficient is zero.
The empirical estimate of the best linear unbiased predictor (BLUP) for each random effect. See Best Linear Unbiased Predictors.
The t ratio for testing that the effect is zero. The t ratio is obtained by dividing the BLUP by its standard error.
The p-value for the test.
The term best linear unbiased predictor (BLUP) refers to an estimator of a random effect. Specifically, it is an estimator that, among all unbiased estimators, minimizes mean square prediction error. The Random Effect Predictions report gives estimates of the BLUPs, or empirical BLUPs. These are empirical because the BLUPs depend on the values of the variance components, which are unknown. The estimated values of the variance components are substituted into the formulas for the BLUPs, resulting in the estimates shown in the report.
The p-value for the test that the covariance parameter is equal to zero. This column appears only when you have selected Unbounded Variance Components in the Fit Model launch window.
The Kackar-Harville correction. See Kackar-Harville Correction. This column appears only if you right-click in the report and select Columns > Norm KHC.
The estimates of σ2 and the variance components in G are obtained by maximizing a residual log-likelihood function that depends on only these parameters. An iterative procedure attempts to maximize the residual log-likelihood function, or equivalently, to minimize twice the negative of the residual log-likelihood (–2LogLike). The Iterations report provides details about this procedure.
The convergence criterion is based on the gradient, with a default tolerance of 10-8. You can change the criterion in the Fit Model launch window by selecting the option Convergence Settings > Convergence Limit and specifying the desired tolerance.
The p-value for the effect test.
When you use the REML method, six additional options appear in the Save Columns menu. These option names start with the adjective Conditional. This prefix indicates that the calculations for these columns use the predicted values for the terms associated with the random effects, rather than their expected values of zero.

Help created on 7/12/2018