Multiple-trait quantitative trait locus mapping with incomplete phenotypic data
© Guo and Nelson; licensee BioMed Central Ltd. 2008
Received: 28 March 2008
Accepted: 05 December 2008
Published: 05 December 2008
Conventional multiple-trait quantitative trait locus (QTL) mapping methods must discard cases (individuals) with incomplete phenotypic data, thereby sacrificing other phenotypic and genotypic information contained in the discarded cases. Under standard assumptions about the missing-data mechanism, it is possible to exploit these cases.
We present an expectation-maximization (EM) algorithm, derived for recombinant inbred and F2 genetic models but extensible to any mating design, that supports conventional hypothesis tests for QTL main effect, pleiotropy, and QTL-by-environment interaction in multiple-trait analyses with missing phenotypic data. We evaluate its performance by simulations and illustrate with a real-data example.
The EM method affords improved QTL detection power and precision of QTL location and effect estimation in comparison with case deletion or imputation methods. It may be incorporated into any least-squares or likelihood-maximization QTL-mapping approach.
Statistical methods for identifying and mapping genes controlling complex traits, commonly known as quantitative trait loci or QTL, have been developed to a high degree. The primary focus has been on methods for single traits ([1–8] and many others). It was proposed [9, 10] that multi-trait QTL mapping methods that consider simultaneously several correlated phenotypic traits, or a single trait measured in several environments, offer increased detection power and precision of location and effect estimation over single-trait QTL mapping. This is because trait-by-trait QTL-searching neglects information contained in the data about the common influence of a QTL on more than one trait, in more than one environment, or  at more than one developmental stage. Multi-trait (MT) QTL mapping allows a formal test of pleiotropy of a QTL for multiple traits or QTL-by-environment interaction for a single trait measured across multiple environments. The enhancement by MT of QTL-detection power is greatest when the QTL induces covariation between the tested traits in the direction opposite to that from "background" sources [12, 13]. These advantages have been exploited in animal [14–16] and plant  studies.
With the promise of increased power from a multivariate approach comes an interesting problem: what to do when some of the multivariate data are missing.
Two main statistical approaches have been elaborated for multi-trait QTL analysis: regression [10, 18–21] and maximum likelihood or ML . Regression QTL-mapping methods, though easier to implement and faster to compute, give biased parameter estimates with sparse markers  or when QTLs interact or are closely linked , while ML methods are free of these defects . It has also been proposed to transform multiple traits into canonical variates so that conventional univariate interval QTL mapping can be applied [18, 24, 25], but interpretation of the results is difficult.
Though QTL-mapping data are often incomplete, information-recovery methods are at present applied only to genotypic data. For incompletely informative marker-genotype data, posterior distributions are readily estimated from flanking markers in the same individual . For unknown QTL genotypes at tested positions in map intervals, ML methods estimate posterior distributions simultaneously with the parameters of a phenotypic mixture distribution , while regression methods  replace missing QTL genotypes with their expectations given flanking markers. Variations based on sampling include multiple imputation (MI) as described by  and  and Bayesian approaches (e.g. [5–7, 28, 29]).
In contrast to genotypic data, missing phenotypic data for any trait results in discarding all cases (individuals) lacking even one value, sacrificing all other phenotypic and genotypic information available for these cases. The problem was recognized by , but they provided no solution, nor does conventional QTL-mapping software offer an alternative to this "casewise"  deletion. Is there one?
Completion of incomplete multivariate data may be done by imputation (single or multiple), by EM algorithm, or by Bayesian approaches. Single imputation typically replaces missing data with three kinds of values: a value drawn from a model-based distribution, a mean of other observations of the same variable, or a conditional mean calculated by least-squares regression on predictors. MI [31, 32] fills in missing data by imputing multiple (e.g. 3–5) times to produce several complete datasets, with parameter estimates calculated as the average over the results from these datasets. The defect of imputation methods, in analyses such as QTL mapping where we want ML estimates of statistics, is that bias is introduced by maximization of the likelihood over both original and imputed data. In contrast, the EM algorithm as described by  focuses not on replacing a missing value with its expectation, but on using the information available in the original dataset. In the framework of EM, missing data imputed are in effect integrated out of the complete-data log likelihood by iterative refinement of their expectation. An EM method described by  addressed the problem of missing genotype or phenotype data in single-marker QTL analyses by the use of flanking-marker genotypes. While free of the dependence on recombination estimates to which interval-mapping methods are subject, the method accommodated only single traits.  provided an EM algorithm for incomplete multivariate data and extended it to accommodate multiple regression with missing responses. A Bayesian approach developed by  for joint mapping with multiple traits in outbred populations employed an identity-by-descent (IBD)-based variance-components model and reversible-jump Markov-chain Monte Carlo (MCMC) estimation, but did not consider missing phenotypic and genotypic data. It would be possible to derive an MCMC algorithm to sample missing phenotypic entries from their posterior distribution, though Bayesian approaches are computationally intensive and often criticized for lacking a test statistic.
Here we describe an adaptation of the EM method of  to the case of multi-trait QTL mapping with incomplete phenotypic data. For simplicity we have limited our scope of mating design to biparental crosses between inbreds. We show that the tests for QTL main effects may be constructed as in , and we describe the properties and behavior of the test statistics and QTL effect and position estimates based on simulation studies and a real example.
All the methods gave similar QTL-detection specificity of 0.98 to 1.00, except with sample size 100 and missing proportion 0.40, where CaD gave specificity as low as 0.93.
Accuracy and precision of QTL effect estimation
Real data analyses
QTL statistics for analysis of two traits in a 325-line doubled-haploid rice population with 10% simulated missing data, estimated by several approaches.
While any substitution of unavailable data by their expectations based on trends in available data will increase the precision of parameter estimates over those obtained by CaD, the EM-based multi-trait QTL mapping method we propose here is superior to MS and CMS for several reasons. MS underestimates phenotypic variation and QTL effect due to fill-in of missing data with a single value, resulting in decreased power compared with our method especially when amounts of missing data are relatively large. The same trend can be observed for CMS, which, as a precursor of the EM algorithm, is closely related to a single EM iteration . Although CMS gave better estimates of QTL effect than MS, it still underestimates variance .
While we did not include MI [31, 37] in the simulation study, we doubt its potential utility for multi-trait QTL mapping with missing trait data. We investigated MI by filling in missing trait data with values sampled from their conditional distributions under the null and alternative hypotheses given the observed trait values. Resulting logarithm-of-odds (LOD) profiles were sawtoothed (not shown here) due to random sampling, and a different profile could be obtained with each analysis even with many imputations (e.g. 100 compared with 3–5 in regular MI) performed at each QTL test position. For these reasons, apart from the high computational cost, we did not pursue this method further.
For MS, CMS, and even MI, the effects on QTL mapping of introducing imputed data need further study. Although simulation results showed specificities close to those of our method, complete-data analysis, and CaD, the bias imposed on the LOD test statistic by introduction of these "artificial" data remains unknown. Interestingly, in the real-data example we chose for illustration, CMS, besides apparently biasing the location of the QTL, gave higher background LOD scores than CoD in regions away from the QTL, while the EM and other approaches did not (Fig. 5). In fact, imputation of missing data is also performed in the E step of our EM algorithm. But this kind of imputation only furnishes a pivot to facilitate parameter estimation and is actually not involved in the likelihood calculation. Thus, theoretically, the EM-based method does not bias QTL detection and parameter estimation as may imputation methods.
The information gain of our method over CaD, MS, and CMS depends on the amount of missing trait data. The reason is readily explained by the following example for CaD. Consider a sample of 200 individuals with missing proportion 0.1 for each of two traits independently. The average number of individuals available for CaD is 162 and that for EM 198, and the difference is 36. This difference expands to 96 with a missing proportion of 0.4. In other words, power is lost more slowly with data loss when the information-recovering EM method is applied.
Cofactor markers too may lack genotype data. In our model, these are replaced by expectations given flanking markers, computed by the method of . Prediction of missing trait data employs these imputed cofactor genotypes, with a resulting potential for error in parameter estimation. While such error could be minimized by extension of our method to include missing cofactor data, we suspect that the improvement would hardly justify the computation, and remark that the issue is shared by all composite interval-mapping methods.
While the EM method should give more power than the information-discarding alternatives regardless of the chosen QTL acceptance threshold, we do not know how to find the optimum threshold for multiple-trait mapping when our algorithm is applied. Indeed, the question has not to our knowledge been satisfactorily answered to date even for the case of complete data. Thresholds in our simulation experiment were based on analyses of simulated populations lacking QTLs but with complete trait data conforming to the same variance structure as the QTL population – an option not available to an analyst in practice. A working method might be to adopt a threshold lying between the overly conservative one calculated from permutation (of individual records for all traits) applied to the complete data remaining after CaD and the insufficiently conservative one that would be obtained from permutation applied under our reconstruction algorithm. For nonpathological data sets our simulation results suggest that this range will be relatively narrow. The quick and approximate method of  might serve as well as any for establishing thresholds.
Some extensions of the EM method are promising. First, we have derived the EM calculation of the hypothesis test for QTL main effect. By following the procedure of , one may derive specific EM implementations for other hypothesis tests including for QTL-by-environment interaction, pleiotropy, and pleiotropy vs. close linkage. Second, the EM method may be extended to multiple interval mapping  with multiple traits and incomplete phenotypic data. Third, mixed-model QTL mapping as recommended by  can now be applied to incomplete trait data as an alternative method for multi-trait QTL mapping. When multiple traits are actually different expressions of a single trait in different environments (locations or years), a mixed model allows treating environmental effect as a random and QTL effect as a fixed factor [39, 40]. One of the advantages of the mixed model is in accommodating both balanced and unbalanced data structure.
In principle, our EM approach may be used to handle multi-trait QTL mapping with any proportion of missing phenotypic data. But if the data may not be MAR, as is especially likely when the missing proportion is large, a more prudent course of analysis is to find out why not, and work out an appropriate statistical method for QTL mapping.
The method we have presented requires more computing time than the conventional EM or ECM interval-mapping algorithm. There are two reasons for this. First, to obtain parameter estimates, the EM algorithm must be applied under both null and alternative hypotheses, because the trait data are missing in both cases. In contrast, conventional methods require EM iteration only under the alternative hypothesis. Second, our EM algorithm is used to complete both QTL genotype and phenotype in the case of ML-based QTL mapping, while the conventional method must complete only QTL genotype. The computing load increases with the proportion of missing data, but the extreme amounts of missing data we have simulated are unusual in real experiments.
Missing-data mechanism is ignored
Several kinds of "missingness" have been defined . Under MAR, "missing at random", the probability of missing phenotypic data within any genotype class is unrelated to the phenotypic value. Either for MAR or the stronger assumption, MCAR or "missing completely at random" (missingness also independent of genotype), estimation methods need not model a missing-data mechanism. Either assumption seems reasonable in conventional mapping practice and is accommodated by the method described here.
Multivariate regression with incomplete data
Consider the linear model
Yn × m= Xn × pBp × m+ En × m, (1)
where y i is the response and x i the predictor vector of the ith individual, YT = [y1, y2, ..., y n ], and XT = [x1, x2, ..., x n ]. Parameters B and V are estimated by maximization of (3).
where is the observed part of y i , is the part of the predictor vector associated with , b obs contains the regression coefficients associated with , and is the submatrix of V representing the variance-covariance matrix of the traits for which the ith individual has complete data. Since missing entries vary among individuals, the log likelihood is a logarithm sum of multivariate normal probabilities of varying dimensions.
Given these partitions, the MLEs of parameters in model (1) are obtained as follows.
ALGORITHM 1: Starting with random initial values , iterate the following two steps until convergence.
Multi-trait QTL mapping with incomplete phenotypic data by regression
We now describe our multi-trait QTL mapping method with incomplete data. Though the method given is based on a recombinant inbred line (RIL) population, it is readily extended to other mating designs, as we show for the F2. According to the statistical model for multiple-trait analysis [9, 10, 19] based on complete phenotypic data, the model for incomplete phenotypic data is written as
Yn × m= zn × 1a(1 × m)+ xn × (p+1)b(p+1) × m+ En × m (13)
where Y = [y1, y2,..., y n ]' is a matrix of phenotypic data for n lines and m traits, and y1, y2,...,y n are 1 × m vectors composed of observations and missing trait data; z is a matrix of QTL genotypes represented as 2 for QQ and 0 for qq; a is a matrix of additive effects of a putative QTL at a tested position; x is a matrix of genotypes of p cofactor markers with the first column ones; b is a matrix of cofactor marker effects; and E is a matrix of residual errors e ij (i = 1, 2, ..., n; j = 1, 2,..., m), which are assumed to be correlated between traits and to follow a multivariate normal distribution with means zero and covariance matrix as in (2). Equation (13) is readily seen to be a variant of (1).
In this model, QTL genotype is replaced with its conditional expectation given flanking-marker genotypes [1, 26]. Least-squares estimates of the parameters can then be obtained by multiple regression based on ALGORITHM 1. If considering a F2 population, we may instead use the model
Yn × m= zn × 1a1 × m+ wn × 1d1 × m+ xn × (p+1)b(p+1) × m+ En × m, (14)
where w is a matrix of QTL genotypes represented as 1 for Qq and 0 for qq and QQ, and d is a matrix of dominance effects of a putative QTL at a tested position. Since w is unobservable, it is also replaced by its probability conditional on flanking markers.
Multitrait QTL mapping with incomplete phenotypic data by ECM
Instead of replacing a missing QTL genotype with its expectation given flanking markers, ECM (expectation/conditional maximization) treats QTL genotype as missing data included in model (14) and estimates parameters at a QTL position by repeatedly updating the posterior probability of QTL genotype given both flanking-marker genotypes and phenotypes. Since we now have two types of missing data in model (14), QTL genotype and phenotype, we may extend the ECM method of  for multi-trait QTL mapping as follows:
iterate the following two steps until convergence (< 10-7 change in log likelihood between two iterations).
where , and are the predicted phenotypic means of QTL genotypes QQ, Qq, and qq given cofactor markers, and V is updated by (25).
In the ECM algorithm, incautious selection of initial parameter values may lead to convergence on local maxima. We used additive and dominance effects of 0, and the estimates of μ and V under the null hypothesis H0: a = 0, d = 0.
Hypothesis tests for QTL effects with missing phenotypic data
Hypothesis tests for QTL main effects, pleiotropy effects and close linkage vs. pleiotropy are constructed according to  and can be tested by ALGORITHM 1 if regression is chosen or ALGORITHM 2 if the ECM method is used. As test statistic the likelihood ratio (LR) or its transformation to a logarithm-of-odds (LOD) are commonly used. For example, to test main QTL effects in a two-trait example, the hypotheses can be formulated as H0: a1 = 0, a2 = 0 and H1: at least one a ≠ 0. For the regression method, parameters under H0 or H1 are estimated by ALGORITHM 1 (Equations 8–12) depending on whether or not QTL effects are included in model (13). If the ECM method is used, first these quantities are estimated under H0 by ALGORITHM 1 without inclusion of QTL effect and then those of the full model under H1 are obtained by ALGORITHM 2 (Equations 16–25 for RIL or 16–19, 26–32, and 25 for F2). Then the LR is obtained as LR = -2(ℓ reduced - ℓ full ), where ℓ reduced is the log likelihood of the reduced model under H0, and ℓ full is that of the full model under H1 . Both are calculated from (4) and a LOD score is calculated as LR/(2 ln 10).
QTL effects and variances for two traits used for simulation of multi-trait QTL mapping.
2. 9 (4.6%)
Total genetic variance
Trait values of each line were calculated as the sum of QTL effects plus a random vector of environmental effects with means zero and variance given in Table 2. Then a specified proportion (0.05, 0.10, 0.20, or 0.40) of values for each trait independently was set to missing. Lines lacking data for both traits were dropped. Analyses were performed on 500 replicates.
In the QTL analyses, the calculation interval (step size) used was 1 cM. Cofactor markers for each trait were selected by forward stepwise regression at a significance level of 0.01 and combined for multi-trait analysis. Cofactors lying within 10 cM of a QTL testing position were dropped from the model.
Though the LR at each test position asymptotically follows a chi-square distribution with degrees of freedom determined by the corresponding hypothesis test, an acceptance threshold applying over all test positions must be found. Genomewide LOD thresholds of 3.71, 3.54 and 3.43 for n = 100, 200, and 300 at significance level 0.05 were calculated from 5000 simulations with no missing data under the null hypothesis of no QTL . When sample size or heritability is relatively small, the effect of a QTL may extend to adjacent intervals due to sampling error. Rather than including heritability as an experimental parameter for investigation, we chose the cautious expedient of declaring a QTL if a LOD peak higher than threshold was found within the interval containing the simulated QTL and the intervals on either side of the QTL interval. Power of QTL detection was calculated as the number of correctly declared ("true positive") QTLs divided by the number of actual QTLs simulated, while specificity was calculated as the number of true positive QTLs divided by the total number declared.
Real data analysis
We applied the EM method to a population of 325 doubled-haploid lines (unpublished data) tested for rice sheath-blight disease in field and greenhouse studies and genotyped with 114 codominant markers. The traits we chose for analysis were sheath-blight score and heading date measured in Stuttgart, Arkansas in 2006 and correlated at r = -0.57. The preliminary SIM and CIM analyses indicated that a QTL on chromosome 12 influenced both traits, and suggested a multi-trait analysis.
Since there were no missing trait data, we generated a new sample with 10% randomly missing phenotype scores. Genomewide LOD thresholds 3.73, 3.77, 3.90, 3.80, and 3.83 for EM, CoD, CaD, MS, and CMS at significance level 0.05 were calculated from 1000 permutations, based on shuffling the phenotypic records for both traits at once in order to preserve their correlation structure.
We thank the reviewers for close attention leading to improvement of the article. Support for this work has been provided by National Science Foundation NSF-DBI grant 0109879 and U. S. Department of Agriculture Cooperative State Research, Education and Extension Service – National Research Initiative – Applied Plant Genomics Program grant 2004-35317-14867 entitled "RiceCAP: A coordinated research, education, and extension project for the application of genomic discoveries to improve rice in the United States". This is contribution KAES 08-378-J from the Kansas Agriculture Experiment Station.
- Haley CS, Knott SA: A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity. 1992, 69: 315-324.View ArticlePubMedGoogle Scholar
- Jansen RC: Interval mapping of multiple quantitative trait loci. Genetics. 1993, 135: 205-211.PubMed CentralPubMedGoogle Scholar
- Kao C-H, Zeng Z-B, Teasdale RD: Multiple interval mapping for quantitative trait loci. Genetics. 1999, 152: 1203-1216.PubMed CentralPubMedGoogle Scholar
- Lander ES, Botstein D: Mapping Mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989, 121: 185-199.PubMed CentralPubMedGoogle Scholar
- Satagopan JM, Yandell BS, Newton MA, Osborn TC: A Bayesian approach to detect quantitative trait loci using Markov chain Monte Carlo. Genetics. 1996, 144: 805-816.PubMed CentralPubMedGoogle Scholar
- Wang H, Zhang Y-M, Li X, Masinde GL, Mohan S, Baylink DJ, Xu S: Bayesian shrinkage estimation of quantitative trait loci parameters. Genetics. 2005, 170: 465-480. 10.1534/genetics.104.039354.PubMed CentralView ArticlePubMedGoogle Scholar
- Yi N, Xu S: Bayesian mapping of quantitative trait loci under complicated mating designs. Genetics. 2001, 157: 1759-1771.PubMed CentralPubMedGoogle Scholar
- Zeng Z-B: Precision mapping of quantitative trait loci. Genetics. 1994, 136: 1457-1468.PubMed CentralPubMedGoogle Scholar
- Jiang C, Zeng Z-B: Multiple trait analysis of genetic mapping for quantitative trait loci. Genetics. 1995, 140: 1111-1127.PubMed CentralPubMedGoogle Scholar
- Korol AB, Ronin YI, Kirzhner VM: Interval mapping of quantitative trait loci employing correlated trait complexes. Genetics. 1995, 140: 1137-1147.PubMed CentralPubMedGoogle Scholar
- Ma CX, Casella G, Wu RL: Functional mapping of quantitative trait loci underlying the character process: A theoretical framework. Genetics. 2002, 161: 1751-1762.PubMed CentralPubMedGoogle Scholar
- Allison DB, Thiel B, St Jean P, Elston RC, Infante MC, Schork NJ: Multiple phenotype modeling in gene-mapping studies of quantitative traits: Power advantages. Am J Hum Genet. 1998, 63: 1190-1201. 10.1086/302038.PubMed CentralView ArticlePubMedGoogle Scholar
- Evans DM: The power of multivariate quantitative-trait loci linkage analysis is influenced by the correlation between variables. Am J Hum Genet. 2002, 70 (6): 1599-1602. 10.1086/340850.PubMed CentralView ArticlePubMedGoogle Scholar
- Kučerová J, Lund MS, Sorensen P, Sahana G, Guldbrandtsen B, Nielsen VH, Thomsen B, Bendixen C: Multitrait quantitative trait loci mapping for milk production traits in Danish Holstein cattle. J Dairy Sci. 2006, 89 (6): 2245-2256.View ArticlePubMedGoogle Scholar
- Neuschl C, Brockmann GA, Knott SA: Multiple-trait QTL mapping for body and organ weights in a cross between NMRI8 and DBA/2 mice. Genet Res. 2007, 89 (1): 47-59. 10.1017/S001667230700852X.View ArticlePubMedGoogle Scholar
- Thomasen J: Quantitative trait loci affecting calving traits in Danish Holstein cattle. J Dairy Sci. 2008, 91 (5): 2098-2105. 10.3168/jds.2007-0602.View ArticlePubMedGoogle Scholar
- Malosetti M: A multi-trait multi-environment QTL mixed model with an application to drought and nitrogen stress trials in maize (Zea mays L.). Euphytica. 2008, 161 (1–2): 241-257. 10.1007/s10681-007-9594-0.View ArticleGoogle Scholar
- Calinski T, Kaczmarek Z, Krajewski P, Frova C, Sari-Gorla M: A multivariate approach to the problem of QTL localization. Heredity. 2000, 84 (3): 303-310. 10.1046/j.1365-2540.2000.00675.x.View ArticlePubMedGoogle Scholar
- Hackett CA, Meyer RC, Thomas WTB: Multi-trait QTL mapping in barley using multivariate regression. Genet Res. 2001, 77: 95-106. 10.1017/S0016672300004869.View ArticlePubMedGoogle Scholar
- Knott SA, Haley CS: Multitrait least squares for quantitative trait loci detection. Genetics. 2000, 156: 899-911.PubMed CentralPubMedGoogle Scholar
- Korol AB, Ronin YI, Nevo E, Hayes PM: Multi-interval mapping of correlated trait complexes: simulation analysis and evidence from barley. Heredity. 1998, 80: 273-284. 10.1046/j.1365-2540.1998.00253.x.View ArticleGoogle Scholar
- Xu S: A comment on the simple regression method for interval mapping. Genetics. 1995, 141: 1657-1659.PubMed CentralPubMedGoogle Scholar
- Kao C-H: On the differences between maximum likelihood and regression interval mapping in the analyis of quantitative trait loci. Genetics. 2000, 156: 855-865.PubMed CentralPubMedGoogle Scholar
- Mangin B, Thoquet P, Grimsley N: Pleiotropic QTL analysis. Biometrics. 1998, 54: 88-99. 10.2307/2533998.View ArticleGoogle Scholar
- Weller JI, Wiggans GR, VanRaden PM, Ron M: Application of a canonical transformation to detection of quantitative trait loci with the aid of genetic markers in a multi-trait experiment. Theor Appl Genet. 1996, 92: 998-1002. 10.1007/BF00224040.View ArticlePubMedGoogle Scholar
- Jiang C, Zeng Z-B: Mapping quantitative trait loci with dominant and missing markers in various crosses from two inbred lines. Genetica. 1997, 101: 47-58. 10.1023/A:1018394410659.View ArticlePubMedGoogle Scholar
- Sen S, Churchill GA: A statistical framework for quantitative trait mapping. Genetics. 2001, 159: 371-387.PubMed CentralPubMedGoogle Scholar
- Sillanpaa MJ, Arjas E: Bayesian mapping of multiple quantitative trait loci from incomplete outbred offspring data. Genetics. 1999, 151 (4): 1605-1619.PubMed CentralPubMedGoogle Scholar
- Sillanpää MJ, Arjas E: Bayesian mapping of multiple quantitative trait loci from incomplete inbred line cross data. Genetics. 1998, 148: 1373-1388.PubMed CentralPubMedGoogle Scholar
- Allison PD: Missing Data. 2002, Thousand Oaks, Calif.: Sage PublicationsGoogle Scholar
- Rubin DB: Multiple Imputation for Nonresponse in Surveys. 1987, New York: WileyView ArticleGoogle Scholar
- Rubin DB: Multiple imputation after 18+ years. Journal of the American Statistical Association. 1996, 91: 473-489. 10.2307/2291635.View ArticleGoogle Scholar
- Dempster AP, Laird NM, Rubin DB: Maximum likelihood from incomplete data via the EM algorithm. J Royal Stat Soc B. 1977, 39 (1): 1-38.Google Scholar
- Niu T, Ding AA, Kreutz R, Lindpaintner K: An expectation-maximization-likelihood-ratio test for handling missing data: application in experimental crosses. Genetics. 2005, 169 (2): 1021-1031. 10.1534/genetics.103.019752.PubMed CentralView ArticlePubMedGoogle Scholar
- Little RJA, Rubin DB: Statistical Analysis with Missing Data. 2001, Hoboken, New Jersey: John Wiley & SonsGoogle Scholar
- Liu JF, Liu YJ, Liu XG, Deng HW: Bayesian mapping of quantitative trait loci for multiple complex traits with the use of variance components. Am J Hum Genet. 2007, 81 (2): 304-320. 10.1086/519495.PubMed CentralView ArticlePubMedGoogle Scholar
- Rubin DB: Inference and missing data. Biometrika. 1976, 63: 581-592. 10.1093/biomet/63.3.581.View ArticleGoogle Scholar
- Piepho HP: A quick method for computing approximate thresholds for quantitative trait loci detection. Genetics. 2001, 157: 425-432.PubMed CentralPubMedGoogle Scholar
- Piepho H-P: A mixed-model approach to mapping quantitative trait loci in barley on the basis of multiple environment data. Genetics. 2000, 156: 2043-2050.PubMed CentralPubMedGoogle Scholar
- Wang DL, Zhu J, Li ZK, Paterson AH: Mapping QTLs with epistatic effects and QTL × environment interactions by mixed linear model approaches. Theor Appl Genet. 1999, 99: 1255-1264. 10.1007/s001220051331.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.