plot.hov | R Documentation |
Oneway analysis of variance makes the assumption that the variances of
the groups are equal.
Brown and Forsyth, 1974 present the recommended test of this assumption.
The Brown and Forsyth test statistic is the F
statistic resulting
from an ordinary one-way analysis of variance on the
absolute deviations from the median. The hovPlot
function
graphs the components of the Brown and Forsyth test statistic.
hovPlot(x, data=NULL, method = "bf", ## x is a formula
transpose = TRUE, ...)
## users will normally use the formula above and will not call the
## method directly.
hovPlot.bf(x, group, ## x is the response variable
y.name = deparse(substitute(x)),
group.name = deparse(substitute(group)),
transpose = TRUE, ...)
## users will normally use the formula above and will not call the
## panel function directly.
panel.hov(..., transpose = TRUE)
x |
Formula appropriate for oneway anova in |
data |
data.frame |
method |
Character string defining method. At this time the only recognized method is "bf" for the Brown-Forsyth method. |
transpose |
Always |
group |
factor. |
y.name |
name of response variable, defaults to variable name in formula. |
group.name |
name of factor, defaults to variable name in formula. |
... |
additional arguments. |
"trellis"
object with three panels containing boxplots for each
group: The observed data "y"
, the data with the median
subtracted "y-med(y)"
, and the absolute deviations from the
median "abs(y-med(y))"
The Brown and Forsyth test statistic is
the F
statistic resulting from an ordinary one-way analysis of
variance on the data points in the third panel.
Richard M. Heiberger <rmh@temple.edu>
Heiberger, Richard M. and Holland, Burt (2015). Statistical Analysis and Data Display: An Intermediate Course with Examples in R. Second Edition. Springer-Verlag, New York. https://link.springer.com/book/10.1007/978-1-4939-2122-5
Brown, M.~B. and Forsyth, A.~B. (1974). Robust tests for equality of variances. Journal of the American Statistical Association, 69:364–367.
aov
, hov
data(turkey)
hov(wt.gain ~ diet, data=turkey)
hovPlot(wt.gain ~ diet, data=turkey)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.