I wish to fit an error variance which varies smoothly with a time variable e.g.
a quadratic
var(t) = b0 + b1 * t + b2 * t * t
Is it possible to fit a function like this? I know I can divide my time variable
up into discrete intervals and get an error variance for each 'environment' but
a quadratic would fit fewer parameters.
A second question. When specifying an error which depends on the level of some
environmental factor (with unequal class frequencies) is it necessary to
specify, for the R structure, the numbers in each class and sort the data on
this factor or is there a simple way of just specifying 'let the error vary with
the levels of this factor' ? Thanks.
___________________________________________________________________________
THE FOLLOWING WARNINGS ATTACH TO THIS COMMUNICATION
This material is intended for the named recipient only and has been created
solely for the purpose and function of Livestock Improvement Corporation Limited
and may contain information that is subject to the Privacy Act / Confidentiality
/ Privilege / Copyright. Any person other than the named recipient is prohibited
from retaining / disclosing / copying / distributing / using all or part of this
material.
--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml