[R] Bootstrapping in R

Christoph Puschmann c.puschmann at student.unsw.edu.au
Sat Oct 1 12:36:58 CEST 2016


Dear Rui,

You can insert a “formula” argument in the code. For example, if you boot a regression, you can insert the formula in the command. Though I just realised that it is not necessary to do.

All the best,

Christoph


> On 1 Oct 2016, at 19:49, ruipbarradas at sapo.pt wrote:
> 
> Sorry, but what formula? formula is not a ?boot argument.
> To the OP: Michael is probably right, if you reset the seed each time, you'll get equal values, otherwise you should get different results due to randomization.
> 
> Rui Barradas
> 
> 
> Quoting Christoph Puschmann <c.puschmann at student.unsw.edu.au>:
> 
>> Dear Bryan,
>> 
>> Did you try to include formula in the boot command? like:
>> 
>> results <- boot(data, statistic, R, formula)
>> 
>> All the best,
>> 
>> Christoph
>> 
>>> On 1 Oct 2016, at 19:24, Michael Dewey <lists at dewey.myzen.co.uk> wrote:
>>> 
>>> Dear Bryan
>>> 
>>> You are not resetting the seed each time by any chance?
>>> 
>>> Michael
>>> 
>>> On 01/10/2016 02:44, Bryan Mac wrote:
>>>> Hi,
>>>> 
>>>> I have read the help page and it was helpful but, I am having concerns because each time I run this code I get the same value.
>>>> I expected that each time I run the code, I will get different values due to random sampling.
>>>> 
>>>> How do I get this randomization? The values shouldn’t be the same each time the code is run, correct?
>>>> 
>>>> 
>>>> result <- boot(n_data, statistic = DataSummary, R = 100).
>>>> 
>>>> 
>>>> Best,
>>>> 
>>>> Bryan Mac
>>>> bryanmac.24 at gmail.com
>>>> 
>>>> 
>>>> 
>>>>> On Sep 29, 2016, at 12:16 PM, ruipbarradas at sapo.pt wrote:
>>>>> 
>>>>> Hello,
>>>>> 
>>>>> Read the help page ?boot::boot.
>>>>> For instance, try the following.
>>>>> 
>>>>> 
>>>>> library(boot)
>>>>> 
>>>>> x <- rnorm(100)
>>>>> stat <- function(x, f) mean(x[f])
>>>>> boot(x, stat, R = 100)
>>>>> 
>>>>> Hope this helps,
>>>>> 
>>>>> Rui Barradas
>>>>> 
>>>>> 
>>>>> 
>>>>> Citando bryan.mac24 <bryan.mac24 at gmail.com>:
>>>>> 
>>>>>> Hi all,
>>>>>> I am wondering how to conduct bootstrapping in R. I need bootstrap 100 times. The trick I need to figure out is how to do get a random sample of 100 out of the total number of case.
>>>>>> Best,
>>>>>> BM
>>>>>> 	[[alternative HTML version deleted]]
>>>>>> 
>>>>>> ______________________________________________
>>>>>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>>>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>>>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>>>>> and provide commented, minimal, self-contained, reproducible code.
>>>>> 
>>>>> 
>>>>> 
>>>> 
>>>> 
>>>> 	[[alternative HTML version deleted]]
>>>> 
>>>> ______________________________________________
>>>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>>> and provide commented, minimal, self-contained, reproducible code.
>>>> 
>>> 
>>> --
>>> Michael
>>> http://www.dewey.myzen.co.uk/home.html
>>> 
>>> ______________________________________________
>>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>> https://stat.ethz.ch/mailman/listinfo/r-help
>>> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
> 
> 



More information about the R-help mailing list