Key message We evaluated four existing Monte Carlo-based methods for approximating prediction error variances and reliabilities in a multiple-trait genomic prediction model. All four methods yielded consistent approximations of the exact values. Abstract 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 geno...