Fusarium head blight is one of the most important wheat diseases causing grain yield and quality losses as well as mycotoxin contamination all over the world. Since Fusarium cannot be reliably controlled with fungicides, breeding has become a favorable tool to decrease the infection severity. In most cases, selection for Fusarium resistance is done by artificial infection in the field. However, there is a risk in preferring late heading genotypes, because heading of wheat is negatively correlated to head blight severity. Because an indirect selection for late maturity is not intended, we considered a statistical approach to avoid this problem. In this paper, we propose a mixed model to analyze extensive Fusarium head blight rating in resistance breeding experiments of wheat. The objective of the analysis was to select for Fusarium resistance, while at the same time ensuring that late heading genotypes, which show less head blight over the shorter vegetation period, are not preferred. Thus, selection was to be done such that genetic variability for heading date was retained. Therefore, the statistical model contained a covariate to adjust for differences in the heading date. The use of covariate adjustment is an easily handled alternative to a bivariate analysis. Covariate adjustment will in practice often work almost equally well as bivariate analysis. Any statistical software with powerful mixed model analysis tools can be used for this type of analysis. We propose an ad hoc method to obtain heritability estimates and a form of LSD (least significance difference) as a measure of accuracy on the basis of the proposed model and under special consideration of the experimental design. The ad hoc LSD was used as a rough measure to judge rankings of genotypic means (BLUPs). Friedman's super smoother was used to compare smoothed rank estimates for adjusted and unadjusted genotypes against increasing smoothed heading dates. Traits were transformed to meet the model assumptions, especially homogeneity of errors and normality, and back-transformation of means and standard errors was conducted by using the delta method.