Implement (in R) the EM algorithm for interval mapping with a
single-QTL model, under the assumption that the residual
variation is homoscedastic and normally distributed, for a
backcross. Consider just one pair of markers at a
recombination fraction, r; assuming complete data at the
markers and no crossover interference.