addto.test.summary | R Documentation |
test.summary
data.frame
.A row that summarizes the result of a proposed change to a model is added to a
test.summary
data.frame
. Only the values of those arguments for which there
are columns in test.summary
will be included in the row.
addto.test.summary(test.summary, terms, DF = 1, denDF = NA,
p = NA, AIC = NA, BIC = NA,
action = "Boundary")
test.summary |
A |
terms |
A |
DF |
A |
denDF |
A |
p |
A |
AIC |
A |
BIC |
A |
action |
A |
A data.frame
.
Chris Brien
asremlPlus-package
, asrtests.object
, print.test.summary
## Not run:
data(Wheat.dat)
## Fit an autocorrelation model
ar1.asr <- asreml(yield ~ Rep + WithinColPairs + Variety,
random = ~ Row + Column + units,
residual = ~ ar1(Row):ar1(Column),
data=Wheat.dat)
ar1.asrt <- as.asrtests(ar1.asr, NULL, NULL,
label = "Autocorrelation model")
ar1.asrt <- rmboundary.asrtests(ar1.asrt)
## Fit a tensor spline
Wheat.dat <- within(Wheat.dat,
{
cRow <- dae::as.numfac(Row)
cRow <- cRow - mean(unique(cRow))
cColumn <- dae::as.numfac(Column)
cColumn <- cColumn - mean(unique(cColumn))
})
ts.asr <- asreml(yield ~ Rep + cRow + cColumn + WithinColPairs +
Variety,
random = ~ spl(cRow) + spl(cColumn) +
dev(cRow) + dev(cColumn) +
spl(cRow):cColumn + cRow:spl(cColumn) +
spl(cRow):spl(cColumn),
residual = ~ Row:Column,
data=Wheat.dat)
ts.asrt <- as.asrtests(ts.asr, NULL, NULL,
label = "Tensor spline model")
ts.asrt <- rmboundary.asrtests(ts.asrt)
ar1.ic <- infoCriteria(ar1.asrt$asreml.obj)
ts.ic <- infoCriteria(ts.asrt$asreml.obj)
if (ar1.ic$AIC < ts.ic$AIC)
{
ic.diff <- ar1.ic - ts.ic
new.asrt <- ar1.asrt
new.asrt$test.summary <- addto.test.summary(ar1.asrt$test.summary,
terms = "Compare ar1 to ts",
DF = ic.diff$varDF,
AIC = ic.diff$AIC, BIC = ic.diff$BIC,
action = "Chose ar1")
} else
{
ic.diff <- ts.ic - ar1.ic
new.asrt <- ts.asrt
new.asrt$test.summary <- addto.test.summary(ts.asrt$test.summary,
terms = "Compare ar1 to ts",
DF = ic.diff$varDF,
AIC = ic.diff$AIC, BIC = ic.diff$BIC,
action = "Chose ts")
}
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.