| linear.approx {boot} | R Documentation |
This function takes a bootstrap object and for each bootstrap replicate it calculates the linear approximation to the statistic of interest for that bootstrap sample.
linear.approx(boot.out, L=NULL, index=1, type=NULL, t0=NULL,
t=NULL, ...)
boot.out |
An object of class "boot" representing a nonparametric bootstrap. It will
usually be created by the function boot.
|
L |
A vector containing the empirical influence values for the statistic of
interest. If it is not supplied then L is calculated through a call
to empinf.
|
index |
The index of the variable of interest within the output of
boot.out$statistic.
|
type |
This gives the type of empirical influence values to be calculated. It is
not used if L is supplied. The possible types of empirical influence
values are described in the helpfile for empinf.
|
t0 |
The observed value of the statistic of interest. The input value is used only
if one of t or L is also supplied. The default value is
boot.out$t0[index]. If t0 is supplied but neither t nor L are supplied
then t0 is set to boot.out$t0[index] and a warning is generated.
|
t |
A vector of bootstrap replicates of the statistic of interest. If t0 is
missing then t is not used, otherwise it is used to calculate the empirical
influence values (if they are not supplied in L).
|
... |
Any extra arguments required by boot.out$statistic. These are needed if
L is not supplied as they are used by empinf to calculate empirical
influence values.
|
The linear approximation to a bootstrap replicate with frequency vector f
is given by t0 + sum(L * f)/n in the one sample with an easy extension
to the stratified case. The frequencies are found by calling boot.array.
A vector of length boot.out$R with the linear approximations to the
statistic of interest for each of the bootstrap samples.
Davison, A.C. and Hinkley, D.V. (1997) Bootstrap Methods and Their Application. Cambridge University Press.
# Using the city data let us look at the linear approximation to the
# ratio statistic and its logarithm. We compare these with the
# corresponding plots for the bigcity data
ratio <- function(d, w)
sum(d$x * w)/sum(d$u * w)
city.boot <- boot(city, ratio, R=499, stype="w")
bigcity.boot <- boot(bigcity, ratio, R=499, stype="w")
par(pty="s")
par(mfrow=c(2,2))
# The first plot is for the city data ratio statistic.
city.lin1 <- linear.approx(city.boot)
lim <- range(c(city.boot$t,city.lin1))
plot(city.boot$t, city.lin1, xlim=lim,ylim=lim,
main="Ratio; n=10", xlab="t*", ylab="tL*")
abline(0,1)
# Now for the log of the ratio statistic for the city data.
city.lin2 <- linear.approx(city.boot,t0=log(city.boot$t0),
t=log(city.boot$t))
lim <- range(c(log(city.boot$t),city.lin2))
plot(log(city.boot$t), city.lin2, xlim=lim, ylim=lim,
main="Log(Ratio); n=10", xlab="t*", ylab="tL*")
abline(0,1)
# The ratio statistic for the bigcity data.
bigcity.lin1 <- linear.approx(bigcity.boot)
lim <- range(c(bigcity.boot$t,bigcity.lin1))
plot(bigcity.lin1,bigcity.boot$t, xlim=lim,ylim=lim,
main="Ratio; n=49", xlab="t*", ylab="tL*")
abline(0,1)
# Finally the log of the ratio statistic for the bigcity data.
bigcity.lin2 <- linear.approx(bigcity.boot,t0=log(bigcity.boot$t0),
t=log(bigcity.boot$t))
lim <- range(c(log(bigcity.boot$t),bigcity.lin2))
plot(bigcity.lin2,log(bigcity.boot$t), xlim=lim,ylim=lim,
main="Log(Ratio); n=49", xlab="t*", ylab="tL*")
abline(0,1)
par(mfrow=c(1,1))