[R] boot vs. bootstrap summary

Andrew Beckerman a.p.beckerman at stir.ac.uk
Thu Nov 21 10:24:02 CET 2002


Dear listers -

Thanks to all for the insight and help on defining the statistic term in 
the boot() function from the boot library in R - it "obviously" differs 
from the bootstrap function in SPLUS - hence my question about translating 
the syntax from the SPLUS example to R - but what i did not realize is that 
the boot library exisits in SPLUS as well.  Many thanks to those who took 
the time to explain the syntax to me... the list, as usual, was most 
helpful (even for the incesant reminders to try ?boot (how could I have 
known about the statisitic portion of boot without looking at that?).

Prof Ripley suggested reading the help more closely and examples in 
Venables and Ripley's MASS book which I did not realize covers the boot 
library as well... my oversight - its page 175 in the 3rd edition for 
linear models.

John Fox suggested his R and S-PLUS Companion which has an on-line appendix 
on bootstrapping regression models that might be of help. You can find it 
(and other appendices) at 
<http://www.socsci.mcmaster.ca/jfox/Books/Companion/appendix.html>.

Angelo Cany provided this advice
...statistic must be an R function (not an expression) and must
take 2 arguments, the original data and a vector of indices defining
the bootstrap replicate.  Your second error is that the number of
bootstrap replicates does not have a default (deliberately so) and
you must specify the number required.

Your example could be done as follows assuming that you want row
resampling.   Let data.model be the original lm fit.
bootcoef <- function(data, i, model) {
     coef(update(model, data=data[i,]))
}
boot(data.frame, statistic=bootcoef,R=999, model=data.model)

And Roger Peng, noting as all others did my confusion (inability to 
uderstand the ?boot text - hence my post) with the statistic term in 
boot(), said:

I wrote a very simple package for bootstrapping in a few classical
situations, including linear models.  You can resample rows and residuals.
It's available at http://department.stat.ucla.edu/~rpeng/uclaboot.  The
interface is (in my opinion, of course!) a little easier to use than the
boot library.

Hope these are of help to others.

cheers
andrew
----------------------------------------------------------
Dr. Andrew Beckerman
Institute of Biological Science
University of Stirling
Stirling FK9 4LA
+44 (0)1786 then wk-467808 fx-464994
(APB is not responsible for anything below this)

-- 
The University of Stirling is a university established in Scotland by
charter at Stirling, FK9 4LA.  Privileged/Confidential Information may
be contained in this message.  If you are not the addressee indicated
in this message (or responsible for delivery of the message to such
person), you may not disclose, copy or deliver this message to anyone
and any action taken or omitted to be taken in reliance on it, is
prohibited and may be unlawful.  In such case, you should destroy this
message and kindly notify the sender by reply email.  Please advise
immediately if you or your employer do not consent to Internet email
for messages of this kind.  Opinions, conclusions and other
information in this message that do not relate to the official
business of the University of Stirling shall be understood as neither
given nor endorsed by it.

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list