Approximating prediction error variances and reliabilities in multiple-trait genomic prediction model using Monte Carlo sampling

Genomic prediction models, such as genomic best linear unbiased prediction (GBLUP), use genomic data to improve the accuracy of estimated genetic values. As the number of genotypes and traits increases, the exact calculation of prediction error variances (PEVs) and reliabilities becomes computationally infeasible due to the need to invert the coefficient matrix of the mixed model equations, whose dimension increases directly with the number of individuals and traits. The objective of this study was to evaluate the applicability of the Monte Carlo (MC) sampling method to approximate PEVs and reliabilities in a multiple-trait GBLUP framework relevant to hybrid breeding. The MC method avoids direct matrix inversion by repeatedly sampling genetic values from their assumed distributions to approximate PEVs. We applied the MC method using four previously published formulas to approximate PEVs and reliabilities. All formulas produced consistent estimates of PEVs and reliabilities, with convergence rates depending on the formula, the level of reliability, and the MC sample size.