The function calculates the consistency statistics k and corresponding p-values for each group (lab) according to Practice E 691 ASTM.
mandelkTest(x, ...)# S3 method for default
mandelkTest(x, g, ...)
# S3 method for formula
mandelkTest(formula, data, subset, na.action, ...)
a numeric vector of data values, or a list of numeric data vectors.
further arguments to be passed to or from methods.
a vector or factor object giving the group for the
corresponding elements of "x"
.
Ignored with a warning if "x"
is a list.
a formula of the form response ~ group
where
response
gives the data values and group
a vector or
factor of the corresponding groups.
an optional matrix or data frame (or similar: see
model.frame
) containing the variables in the
formula formula
. By default the variables are taken from
environment(formula)
.
an optional vector specifying a subset of observations to be used.
a function which indicates what should happen when
the data contain NA
s. Defaults to getOption("na.action")
.
A list with class "mandel"
containing the following components:
a character string indicating what type of test was performed.
a character string giving the name(s) of the data.
the p-value for the test.
the estimated quantiles of Mandel's statistic.
a character string describing the alternative hypothesis.
a character vector describing the levels of the groups.
the number of replicates for each group.
Practice E 691 (2005) Standard Practice for Conducting an Interlaboratory Study to Determine the Precision of a Test Method, ASTM International.
# NOT RUN {
data(Pentosan)
mandelkTest(value ~ lab, data=Pentosan, subset=(material == "A"))
# }
Run the code above in your browser using DataLab