newrcov.asrtests {asreml3Plus} | R Documentation |
asreml
.Fits a new rcov formula using asreml
. Any boundary terms resulting
after fitting the new rcov
model are removed using
rmboundary.asrtests
. A row is added to the test.summary
data.frame
using the supplied label
and the action
Set
.
newrcov.asrtests(terms=NULL, asrtests.obj, label = "R model", update = TRUE, allow.unconverged = TRUE, trace = FALSE, denDF="default", set.terms = NULL, ignore.suffices = TRUE, constraints = "P", initial.values = NA, ...)
terms |
A model for the |
asrtests.obj |
An |
label |
A character string to use as the label in |
update |
If |
allow.unconverged |
A |
trace |
If TRUE then partial iteration details are displayed when ASReml-R functions are invoked; if FALSE then no output is displayed. |
denDF |
Specifies the method to use in computing approximate denominator
degrees of freedom when |
set.terms |
A character vector specifying the terms that are to have constraints and/or initial values set prior to fitting. |
ignore.suffices |
A logical vector specifying whether the suffices of the
|
constraints |
A character vector specifying the constraints to be applied
to the terms specified in |
initial.values |
A character vector specifying the initial values for
the terms specified in |
... |
further arguments passed to |
An asrtests
object, which is a list containing:
asreml.obj
: an asreml
object containing the fit
after the term
has been omitted from the model;
wald.tab
: a 4-column data.frame
containing a
pseudo-anova table for the fixed terms produced by wald.asreml
;
test.summary
: a data.frame
with columns term
,
DF
, denDF
, p
and action
. A row is added to
it for each term
that is dropped, added or tested or a note that several terms have been
added or removed. A row contains the name of the term, the
DF, the p-value and the action taken. Possible codes are:
Dropped
, Retained
, Swapped
, Unswapped
,
Significant
, Nonsignificant
, Absent
, Added
,
Removed
and Boundary
. If the changed model did not
converge, Unconverged
will be added to the code.
Note that the logical asreml.obj$converge
also
reflects whether there is convergence.
If the term
is not in the model, then the supplied asreml
object will be returned. Also, reml.test
will have the likelihood
ratio and the p-value set to NA
and the degrees of freedom to zero.
Similarly, the row of test.summary
for the term
will have
its name, a p-value set to NA
, and action set to Absent.
asreml3Plus-package
, asrtests
, testrcov.asrtests
,
choose.model.asrtests
,
reml.lrt.asreml
, rmboundary.asrtests
,
newfit.asreml
, testswapran.asrtests
,
addrm.terms.asrtests
,
sig.devn.reparam.asrtests
## Not run: data(Wheat.dat) current.asr <- asreml(yield ~ Rep + WithinColPairs + Variety, random = ~ Row + Column + units, rcov = ~ ar1(Row):ar1(Column), data=Wheat.dat) current.asrt <- asrtests(current.asr, NULL, NULL) current.asrt <- rmboundary.asrtests(current.asrt) # Fit Row autocorrelation current.asrt <- newrcov.asrtests("~ Row:ar1(Column)", current.asrt, label="Row autocorrelation") print(current.asrt) ## End(Not run)