h2(x) = phi(x)' C phi(x) / V_P(x), with C the covariance of the reaction-norm group and
phi(x) the Legendre basis evaluated at x. The phenotypic variance V_P(x) is evaluated at x
too: every other random regression on the same basis (the permanent environment of a
test-day model) enters as phi(x)' P phi(x), each scalar term with its variance, and the
residual. A covariance between two scalar terms enters with coefficient 1, as in h2(),
and the residual of a fit with weights is read at the harmonic mean weight of the
records, as in h2() (weights w and 10 w give the same curve).
Every random regression of the fit is evaluated on the Legendre basis with these limits.
Arguments
- fit
result of
model()with an rn() term- group
name of the reaction-norm group; by default the group of the direct genetic term, as in
h2()- limits
the same ones passed to
legendre()when preparing the data- points
where to evaluate the curve