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
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/multcompBoxplot.R
\name{multcompBoxplot}
\alias{multcompBoxplot}
\title{boxplot with multcomp graphics}
\usage{
multcompBoxplot(
formula,
data,
horizontal = TRUE,
compFn = "TukeyHSD",
sortFn = "mean",
decreasing = TRUE,
plotList = list(boxplot = list(fig = c(0, 0.75, 0, 1)), multcompTs = list(fig = c(0.7,
0.85, 0, 1)), multcompLetters = list(fig = c(0.87, 0.97, 0.03, 0.98), fontsize = 20,
fontface = "bold"))
)
}
\arguments{
\item{formula}{a two sided formula like "y~z", where both "y" and "z" are
columns of the data.frame "data", "y" is numeric, and "z" is a factor. This
will be passed as the first argument for both 'boxplot' and 'compFn', and so
must work in both contexts. NOTE: Any more complicated formula may produce
errors or unanticipated results.}
\item{data}{A data.frame for evaluating 'formula'.}
\item{horizontal}{TRUE for horizontal boxplots and vertical multcompTs and /
or multcompLetters; FALSE for the opposite.}
\item{compFn}{a function whose output will serve as the the only non-default
input to either 'multcompTs' or 'multcompLetters'. The default "TukeyHSD"
actually translates to 'TukeyHSD(aov(formula, data))[[1]][, "p adj"]'.}
\item{sortFn}{If sortFn is a function or a character string naming a
function, it is used to summarize the subset of y corresponding to each
level of z into a single number, which will then be used to sort the levels
of z according to the argument 'decreasing'. This step is skipped if
sortFn is NULL or NA or if it is neither a function nor a character string
that might name a function. If sortFn is a character string but a function
by that name is not found in the search path, multcompBoxplot stops with,
'Error in do.call(sortFn, list(x = x)) : could not find function ...'.}
\item{decreasing}{If the levels of z are to be sorted using the output of
'sortFn', this is uses as the 'decreasing' in 'order' to sort the levels of
z for plotting.}
\item{plotList}{A list with names in c("boxplot", "multcompTs",
"multcompLetters"). Replicates are allowed. If present, they produce,
e.g., multiple "multcompTs" side by side. This can be used to compare the
visual effects of different arguments to "plot.multcompTs". Each component
of 'plotList' is itself a list of arguments to pass to either "boxplot",
"plot.multcompTs" or "plot.multcompLetters". Placement can be controlled
via 'fig' arguments passed (indirectly) of the form 'c(x1, x2, y1, y2)'.
If(horizontal==TRUE), fig gives the coordinates of the figure region in the
display region of the plot device, as described on the \code{\link{par}}
help page; if(horizontal==FALSE), fig[c(3,4,1,2)] gives c(x1, x2, y1, y2)
for placement of that portion of the plot.}
}
\value{
This function invisibly returns a list with one component for each
component of plotList containing the output of the appropriate
"plot.multcomp" call plus the output of "compFn".
}
\description{
Create boxplots with multcompTs and / or multcompLetters
}
\details{
For formula = y~z, if 'sortFn' is a function or the name of a function,
'multcompBoxplot' starts by applying sortFn to the subsets of y
corresponding to each level of z, and then sorting those summaries in
increasing or decreasing order, per 'decreasing'. If 'sortFn' is NULL or
NA, this sort step is skipped.
'multcompBoxplot' then creates 'boxplot' as specified in 'plotList'. Next,
'compFn' is called to generate comparisons to feed to the functions
(\code{\link{multcompTs}} and / or \code{\link{multcompLetters}}, whose
output is then passed to (\code{\link{plot.multcomp}}) for plotting.
Components of the relevant sublists of 'plotList' are made available to
\code{\link{par}} or (for \code{\link{plot.multcompLetters}}) to
\code{\link[grid]{gpar}}.
}
\examples{
# Example from help("TukeyHSD")
multcompBoxplot(breaks~tension, data=warpbreaks)
# 'sortFn' can be either a function or a function name
# default order is 'decreasing=TRUE'
multcompBoxplot(breaks~tension, data=warpbreaks,
sortFn=median, decreasing=FALSE)
##################
library(multcomp)
data(recovery)
# Horizontal boxplots with both
# multcomp Ts and Letters on the right
# Using recovery{multcomp} data set
multcompBoxplot(minutes~blanket, recovery)
# Plotting boxes rather than letters and Ts
multcompBoxplot(minutes~blanket, recovery,
plotList=list(
boxplot=list(fig=c(0, 0.75, 0, 1), las=1,
cex.axis=1.5),
multcompTs=list(fig=c(0.7, 0.85, 0, 1),
type='boxes'),
multcompLetters=list(
fig=c(0.87, 0.97, 0.03, 0.98),
type='boxes') ) )
####################
# Vertical boxplots with both
# multcomp Ts and Letters on the top
multcompBoxplot(minutes~blanket, recovery,
horizontal=FALSE)
# Horizontal boxplots with 2 different
# displays of the "Ts" on the left
multcompBoxplot(minutes~blanket, recovery,
plotList=list(
boxplot=list(fig=c(0.3, 1, 0, 1)),
multcompTs=list(fig=c(0, 0.15, 0, 1),
orientation="reverse"),
multcompTs=list(fig=c(0.15, 0.3, 0, 1),
type="boxes", orientation="reverse",
mar=c(5,2, 4, 0)+.1) ) )
library(MASS)
anorx <-
multcompBoxplot(Postwt~Treat, data=anorexia)
\dontrun{
# Confirm than sortFn=NULL or NA
# leaves the order unchanged
library(multcomp)
data(cholesterol)
cholesterol$trt3 <- with(cholesterol, factor(
as.character(trt), levels=levels(trt)[c(5:4,1:3)]))
multcompBoxplot(response ~ trt3, cholesterol,
sortFn=NULL)
multcompBoxplot(response ~ trt3, cholesterol,
sortFn=NA)
}
}
\seealso{
\code{\link{boxplot}} \code{\link{multcompTs}}
\code{\link{multcompLetters}} \code{\link{plot.multcomp}}
\code{\link{TukeyHSD}} \code{\link{par}} \code{\link[grid]{gpar}}
}
\author{
Spencer Graves
}
\keyword{aplot}
|
