Function for applying R inverse to u in REML EM algorithm

Rinv.u(G, AtG, GtAAtG, GtAu, u, tau, sigma)

G | Matrix of genotype markers (size n x L) |
---|---|

AtG | AtG from precomputation |

GtAAtG | GtAAtG from precomputation |

GtAu | GtAu from precomputation |

u | u=Aty from REML EM |

tau | Variance component from G main effect |

sigma | Variance component from model noise epsilon |

Vector resulting from left multiplication of Rinv with input vector u