setvarianceterms.call {asremlPlus} | R Documentation |
allows the setting of bounds and initial values for terms in the
random
and residual
arguments of an asreml
call, with
the resulting call being evaluated.
Description
Takes an unevaluated call
and evaluates the call
after setting the bounds and initial values for the terms
specified in terms
. The elements of terms
are matched
with those generated by asreml and used, for example, in the
varcomp
component of a summary.asreml
object.
These names generally include descriptive suffices. To match an
element of terms
that includes such a suffix, set
ignore.suffices
to FALSE
so that a literal match
between the element and the assigned names is sought.
Usage
## S3 method for class 'call'
setvarianceterms(call, terms, ignore.suffices = TRUE,
bounds = "P", initial.values = NA, ...)
Arguments
call |
an unevaluated |
terms |
A character vector specifying the terms
that are to have bounds and/or initial values specified.
The names must match those in the |
ignore.suffices |
A logical vector specifying whether the suffices of the
|
bounds |
A The codes used by ASReml are:
|
initial.values |
A character vector specifying the initial values for
the terms specified in |
... |
additional arguments to be added to the call, or arguments in the call with changed values. |
Value
An asreml
object.
Author(s)
Chris Brien
References
Butler, D. G., Cullis, B. R., Gilmour, A. R., Gogel, B. J. and Thompson, R. (2023). ASReml-R Reference Manual Version 4.2. VSN International Ltd, https://asreml.kb.vsni.co.uk/.
See Also
newfit.asreml
, update.asreml
Examples
## Not run:
m1.call <- call("asreml",
fixed = Height ~ (Block + Irrig)*csDay.num,
random= ~ spl(csDay.num)/(Irrig+Block)
+ dev(csDay.num)
+ str(~Block:Plot/csDay.num, ~us(2):id(20))
+ Block:Plot:spl(csDay.num),
data=quote(dat)) ##use quote to stop evaluation of dat here
terms <- c("Block:Plot+Block:Plot:csDay.num!us(2).2:1", "R!variance")
m1.asreml <- setvarianceterms(m1.call, terms, bounds=c("U","P"),
initial=c(NA,3), ignore.suffices=c(FALSE,TRUE))
summary(m1.asreml)
## End(Not run)