calcLSD: Approximate pairwise least significant interval from an...

In kohleth/kmisc: Kohleth Miscellaneous

Description

It approximate a pariwise LSD from the sed matrix from an asreml predict call. The goal is to ease plotting.

Usage

calcLSD(asrpred, weight = NULL)

plotLSD(asrpred, weight = NULL, order = c("no", "increasing", "decreasing"),
  ...)

Arguments

asrpred	Result of a predict.asreml(..., sed=TRUE)$prediction call.
weight	A weight matrix with the same dimension as the sed matrix. It feeds into the lm call which weighs the sed matrix.
order	One of "no"(default), "increasing", or "decreasing". Whether the resulting plot should be ordered by the predicted mean.

Examples

data(oats)
oats.asr <-
  asreml(yield ~ Variety * Nitrogen,
         random = ~ Blocks / Wplots,
         data = oats)
oats.pred <-
  predict(oats.asr, classify = "Variety", sed = TRUE)$prediction
plotLSD(oats.pred)
plotLSD(oats.pred, order = "inc")
plotLSD(oats.pred, order = "dec")
calcLSD(oats.pred)

kohleth/kmisc documentation built on May 20, 2019, 12:53 p.m.