predictiondiffs.asreml {asreml3Plus} | R Documentation |
Uses predictions and standard errors of pairwise differences from an
alldiffs
object to form, for those components not already present,
(i) a table of all pairwise differences of the predictions in an alldiffs
object, (ii) the p-values of each pairwise difference, and (iii) the minimum, mean and
maximum LSD values. Predictions that are aliased (or nonestimable) are removed
from the predictions
component of the alldiffs
object and
standarard errors of differences involving them are removed from the sed
component.
Each p-value is computed as the probability of a t-statistic as large as or larger than
the absolute value of the observed difference divided by its standard error. The
p-values are stored in the p.differences
component. The degrees of freedom of
the t-distribution is the degrees of freedom stored in the tdf
attribute of the
alldiffs
object. This t-distibrution is also used in calculating the
LSD statistics stored in the alldiffs
object.
predictiondiffs.asreml(classify, alldiffs.obj, x.num = NULL, x.fac = NULL, levels.length = NA, pairwise = TRUE, alpha = 0.05, inestimable.rm = TRUE)
classify |
a character string giving the variables that define the margins
of the multiway table to be predicted. Multiway tables are
specified by forming an interaction type term from the
classifying variables, that is, separating the variable names
with the |
alldiffs.obj |
An |
x.num |
A character string giving the name of the numeric covariate that
corresponds to |
x.fac |
A character string giving the name of the factor that corresponds to
|
levels.length |
The maximum number of characters from the the levels of factors to use in the row and column labels of the tables of pairwise differences and their p-values and standard errors. |
pairwise |
A logical indicating whether all pairwise differences of the
|
alpha |
The significance level for an LSD to compare a pair of predictions. |
inestimable.rm |
A |
An alldiffs
object that is a list with components
predictions
containing the predictions and
their standard errors, differences
containing all pairwise
differences between the predictions, p.differences
containing
p-values for all pairwise differences between the predictions,
sed
containing the standard errors of all pairwise differences
between the predictions, and,an LSD
containing the mean, minimum
and maximum LSDs.
asreml3Plus-package
, alldiffs
,
print.alldiffs
, predictionplot.asreml
,
predictparallel.asreml
, pred.present.asreml
## Not run: Var.pred <- predict(current.asr, classify="Variety", sed=TRUE)$predictions wald.tab <- current.asrt$wald.tab den.df <- wald.tab[match("Variety", rownames(wald.tab)), "denDF"] Var.diffs <- alldiffs(predictions = Var.pred$pvals, sed = Var.pred$sed, tdf = den.df) Var.diffs <- predictiondiffs.asreml(classify = "Variety", alldiffs.obj=Var.diffs) print.alldiffs(Var.diffs, which="differences") ## End(Not run)