Estimation of individual animal SNP-BLUP reliability using full Monte Carlo sampling

Calculation of individual animal reliability of estimated genomic breeding value by SNP-BLUP requires inversion of the mixed model equations (MME). When the SNP-BLUP model includes a residual polygenic (RPG) effect, the size of the MME will be at least the number of genotyped animals (n) plus the number of SNP markers (m). Inversion of the MME in SNP-BLUP involves computations proportional to the cube of the MME size; that is, (n + m)(3), which can present a considerable computational burden. We introduce a full Monte Carlo (MC) sampling-based method for approximating reliability in the SNP-BLUP model and compare its performance to the genomic BLUP (GBLUP) model. The performance of the full MC approach was evaluated using 2 data sets, including 19,757 and 222,619 genotyped animals selected from populations with 231,186 and 13.35 million pedigree animals, respectively. Genotypes were available in the data sets for 11,729 and 50,240 SNP markers. An advantage of the full MC approximation method was its low computational demand. A drawback was its tendency to overestimate reliability for animals with low reliability, especially when the weight of the RPG effect was high. The overestimation can be lessened by increasing the number of MC samples.