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.
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.
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. > >