LeveneTest | R Documentation |
Computes Levene's test for homogeneity of variance across groups.
LeveneTest(x, ...)
## S3 method for class 'formula'
LeveneTest(formula, data, subset, na.action, ...)
## Default S3 method:
LeveneTest(x, g, center = median, ...)
x |
response variable for the default method, or a |
... |
arguments to be passed down, e.g., |
formula |
a formula of the form |
data |
an optional matrix or data frame (or similar: see
|
subset |
an optional vector specifying a subset of observations to be used. |
na.action |
a function which indicates what should happen
when the data contain NAs. Defaults to |
g |
factor defining groups. |
center |
The name of a function to compute the center of each group;
|
Let X_ij
be the jth observation of X for the ith group.
Let Z_ij = |X_ij - X_i|
, where X_i
is the mean of X in the ith group.
Levene’s test statistic is
W_0 = \frac{ \sum_i n_i (\bar{Z}_i - \bar{Z})^2 / (g - 1) }{ \sum_i
\sum_j (Z_{ij} - \bar{Z}_i)^2 / \sum_i (n_i - 1) }
where n_i
is the number of observations in group i and g is the number of
groups.
returns an object meant to be printed showing the results of the test.
This function is rewritten using common R standards based on car::leveneTest() using the same calculation logic.
andri.signorell andri@signorell.net; original version
written by John Fox jfox@mcmaster.ca based on a generic version
contributed by Derek Ogle
adapted from a response posted by Brian Ripley
to the r-help email list.
Fox, J. (2008) Applied Regression Analysis and Generalized Linear Models, Second Edition. Sage.
Fox, J. and Weisberg, S. (2011) An R Companion to Applied Regression, Second Edition, Sage.
fligner.test
for a rank-based (nonparametric)
k
-sample test for homogeneity of variances; mood.test
for another rank-based two-sample test for a difference in scale parameters;
var.test
and bartlett.test
for parametric tests
for the homogeneity in variance.
ansari_test
in package coin for exact
and approximate conditional p-values for the Ansari-Bradley test, as
well as different methods for handling ties.
## example from ansari.test:
## Hollander & Wolfe (1973, p. 86f):
## Serum iron determination using Hyland control sera
serum <- ToLong(data.frame(
ramsay=c(111, 107, 100, 99, 102, 106, 109, 108, 104, 99,
101, 96, 97, 102, 107, 113, 116, 113, 110, 98),
jung.parekh=c(107, 108, 106, 98, 105, 103, 110, 105, 104,
100, 96, 108, 103, 104, 114, 114, 113, 108, 106, 99)
))
LeveneTest(x ~ grp, data=serum)
LeveneTest(x ~ grp, data=serum, center=mean)
LeveneTest(x ~ grp, data=serum, center=mean, trim=0.1)
LeveneTest( c(rnorm(10), rnorm(10, 0, 2)), factor(rep(c("A","B"), each=10)) )
LeveneTest(Ozone ~ Month, data = airquality)
LeveneTest(count ~ spray, data = InsectSprays)
# Compare this to fligner.test() and bartlett.test()
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.