Description Usage Arguments Value Author(s) See Also Examples
This function does a radial plot based on the vector of letters resulted from pairwise comparisons.
1 2 3 4 5 6 7 8 |
cld |
Character vector with strings of letters that indicates which pair of treatment cells are not different. |
labels |
Vector of text to be annotated next each point. |
col |
Vector of colors to be used in the segments that joint points. |
means |
Numeric vector with the estimated means of treatment cells. It is used to place points at distances proportional to the differences on means. |
perim |
Logical value (default is |
legend |
Logical value (default is |
None is returned, only the plot is done.
Walmes Zeviani, walmes@ufpr.r.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 | set.seed(4321)
td <- data.frame(trt = rep(sample(1:20), each = 5))
td$y <- rnorm(nrow(td), mean = 0.15 * sort(td$trt), sd = 1)
plot(y ~ trt, data = td)
# Fit the model.
td$trt <- factor(td$trt)
m0 <- lm(y ~ trt, data = td)
anova(m0)
summary(m0)
library(multcomp)
library(doBy)
X <- LE_matrix(m0, effect = "trt")
rownames(X) <- levels(td$trt)
ci <- apmc(X, m0, focus = "trt", test = "fdr")
ci$cld <- with(ci, ordered_cld(cld, fit))
ci <- ci[order(ci$fit, decreasing = TRUE), ]
library(latticeExtra)
segplot(reorder(trt, fit) ~ lwr + upr,
centers = fit,
data = ci,
draw = FALSE,
cld = ci$cld) +
layer(panel.text(x = centers,
y = z,
labels = sprintf("%0.2f %s",
centers,
cld),
pos = 3))
radial_cld(cld = ci$cld)
radial_cld(cld = ci$cld, means = ci$fit, perim = TRUE)
radial_cld(cld = ci$cld, col = 1:3)
radial_cld(cld = ci$cld, col = 1:3)
radial_cld(cld = ci$cld, labels = sprintf("%0.2f %s", ci$fit, ci$cld))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.