For each pair of methods in
data, a regression of the differences on
the averages between methods is made and a linear relationship between
methods with prediction standard deviations is derived.
1 2 3 4 5 6 7 8
A character string, or a list of two functions, each
other's inverse. The measurements are transformed by this before analysis.
Possibilities are: "exp", "log", "logit", "pctlogit" (transforms percentages
by the logit), "sqrt", "sq" (square), "cll" (complementary log-minus-log),
"ll" (log-minus-log). For further details see
The tolerance used to check whether the supplied
transformation and its inverse combine to the identity. Only used if
Should the results be printed?
If methods really are a random selection of raters,
neither intercept nor slope different from 0 are sensible, so if this is
If this is TRUE, a slope of the differences in the verages is estimated, otherwise the relationship is assumed constant.
If the input object contains replicate measurements these are taken as separate items in the order they appear in the dataset.
DA.reg returns a
MethComp object, i.e. a list
with three components,
Conv is a three-dimensional array, with dimensions
From (both with levels equal to the methods in
data) and an
unnamed dimension with levels
"beta=1", referring to the linear relationship of
"sd(t-f)", referring to the regression of the differences on the
referring to the regression of the absoulte residuals on the averages, and
LoA-hi, the limits of agreement.
Converting from method l to method k using
with prediction standard deviation
σ, just requires the entries
[k,l,c("alpha","beta","sd.pred")], if we assume the s.d. is constant.
The next entry is the p-values for the hypothesis β=1, intercept
and slope of the SD of the differences as a linear function of the average
and finally p-value of the hypothesis that standard errors are constant over
the range. The latter three are derived by regressing the absolute values of
the residuals on the averages, and can be used to produce LoA where the s.d.
increases (or decreases) by the mean, using the function
VarComp element of the list is
NULL, and only present for
compatibility with the print method for
data element is the input dataframe. The measurements in
are left un-transformed, even if data are transformed (i.e. if the
Transform attribute of the object is non-null).
DA2y returns a 2 by 3 matrix with rownames
c("int","slope","sd"), calculated under the
assumption that the differences were formed as
D <- y1 - y2.
y2DA returns a 3-component vector with names
c("DA-int","DA-slope","DA-sd"), referring to differences
D=y1-y2 as a linear function of
B. Carstensen: Comparing methods of measurement: Extending the LoA by regression. Stat Med, 29:401-410, 2010.
1 2 3 4 5 6 7 8
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.