Genetic evaluation by best linear unbiased prediction (BLUP) requires modeling genetic means, variances, and covariances. This paper presents theory to model means, variances, and covariances in a multibreed population, given marker and breed information, in the presence of gametic disequilibrium between the marker locus (ML) and linked quantitative trait locus (MQTL). Theory and algorithms are...