[R] Bootstrapping in R

Michael Dewey lists at dewey.myzen.co.uk
Sat Oct 1 11:24:54 CEST 2016


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



More information about the R-help mailing list