Recently detected bugs fixed
model function rrdi(Exp) is treated as shorthand so that rrd2(Exp).grm(Geno) is replaced by rr2(Exp).grm(Geno) + diag(Exp).grm(Geno). This is a sparser formulation of a factor analytic model that the equivalent xfa2(Exp).grm(Geno).