boot and coefficients
Your error is based on a misunderstanding of what the 'statistic' should look like for the boot() function. The statistic needs to take the data as an argument and a vector of indices -- it is the indicies which get resampled. 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. -roger _______________________________ UCLA Department of Statistics rpeng at stat.ucla.edu http://www.stat.ucla.edu/~rpeng
On Wed, 20 Nov 2002, Andrew Beckerman wrote:
R1.61, WinNT Dear listers - I am trying to implement the SPLUS example for boostrapping coefficients from a linear model - boot(data.frame, coef(eval(model$call))) I get errors saying that my statistic is not defined and that the replicates are not defined. Changing this to - boot(recTTM, statisitic=coef(eval(mod$call)),R=999) still gives the error about the statistic. Clearly I do not understand what the statisitic is and how to make the coefficients a function as in the examples in the help... can someone guide me on syntax? 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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 _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._