bootMtrPairs-class {metRology} | R Documentation |
Object returned by bootMtrPairs
and associated methods.
Description
The object class returned by boot.mtr.pairwise
and associated
print
, summary
, and plotting classes.
Usage
## S3 method for class 'bootMtrPairs'
print(x, ...)
## S3 method for class 'bootMtrPairs'
plot(x, ...)
## S3 method for class 'bootMtrPairs'
barplot(height, ylab=NULL, names.arg=height$labels,
crit.vals=TRUE, lty.crit=c(2,1), col.crit=2, lwd.crit=c(1,2), ylim=NULL, ... )
## S3 method for class 'bootMtrPairs'
summary(object, p.adjust="none", ...)
## S3 method for class 'summary.bootMtrPairs'
print(x, digits=3, ...,
signif.stars = getOption("show.signif.stars"),
signif.legend=signif.stars)
Arguments
x |
An R object. For |
height |
An object of class |
object |
An object of class |
p.adjust |
Multiple correction method for calculated p-values, passed to
|
ylab |
Label for vertical axis, passed to |
names.arg |
Labels for individual bars in bar plot, passed to |
crit.vals |
If |
lty.crit , col.crit , lwd.crit |
Vectors of line style parameters for plotted critical values, passed to
|
ylim |
Limits for plot y range, passed to |
digits |
integer; passed to |
signif.stars |
logical; if |
signif.legend |
logical; if |
... |
Parameters passed to other methods. |
Details
The default plot
method is an alias for the barplot
method.
For the plot methods, quantiles for each point are taken directly from the quantiles
calulated by boot.mtr.pairwise
and retained in the returned object.
For the summary
method, p-values are initially calculated as the observed
proportion of simulated values exceeding the MSD value calculated by msd
. The
summary method additionally returns p-values after adjustment
for multiple comparisons using the adjustment method specified.
The print
method for the summary.bootMtrPairs
object prints the summary as a data
frame with columns for the calculated values (re-titled for the particular statstic),
data-specific upper quantiles (one column for each probability supplied to bootMtrPairs
and the p-values after adjustment for multiple comparisons based on the
proportion of simulated values exceeding the observed MSD.
Where that proportion is zero, the summary replaces the
raw zero proportion with 1/B
, corrects that proportion using the requested
adjustment method, and reports the p-value as less than ("<") the resulting
adjusted value.
Value
The print
method returns the object, invisibly.
The plot
and barplot
methods return the values at the midpoint of each bar.
The summary
method returns an object of class "summary.bootMtrPairs"
which
is a list with members:
t0 |
|
labels |
character vector of labels, by default taken from |
probs |
vector of probabilities supplied and used for quantiles |
critical.values |
matrix of quantiles. Each row corresponds to a probability
in |
pvals |
p-values estimated as the observed proportion of
simulated values exceeding the calculated values |
p.adjust |
Character value containing the name of the p-value adjustment method used. |
p.adj |
p-values adjusted using the given p-value adjustment method
specified by |
B |
Number of bootstrap replicates used. |
method |
The sampling method used by the parametric bootstrap. |
stat |
The statistic subjected to bootstrapping. Either |
t |
If |
.
Author(s)
S. L. R. Ellison s.ellison@lgcgroup.com
See Also
Examples
## Not run:
data(Pb)
pwch.Pb<-pdchisq(Pb$value, Pb$u) # Uses individual standard uncertainties
set.seed(1023)
boot.Pb <- boot.mtr.pairwise(pwch.Pb)
summary(boot.Pb)
# The default summary gives individual observation p-values. To
# avoid over-interpretation for the study as a whole,
# apply a sensible p-value adjustment:
summary(boot.Pb, p.adjust="holm")
plot(boot.Pb, crit=TRUE)
## End(Not run)