s-news
[Top] [All Lists]

bootstrap2

To: s-news@lists.biostat.wustl.edu
Subject: bootstrap2
From: Karin <twentestat@yahoo.com>
Date: Mon, 20 Dec 2004 02:59:59 -0800 (PST)
Comment: DomainKeys? See http://antispam.yahoo.com/domainkeys
Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=s1024; d=yahoo.com; b=WClr7GRhGiluE+uvoTKWpXecHXeHeo9IrS6YB1AiUB508StiowV/k4iHmcgXBVDdVrOqTF6oiKMEpvd7Fn4/4KRuOsbS2b7rDzDCWmwPV5DI8dcs/sDrkanhi28nMfA/dbn4Jm8JxcbStmsV3DpBuL+tW5OXi80grbC14ZrOsew= ;
Hello,
 
I want to use the function bootstrap2 of the library resample from T. Hesterberg.
When I am trying the two examples of the help sections of bootstrap2 I got an error-message (see below):
Am I doing something wrong?
 
Thanks in advance,
 
Karin
--------------------------------
 
the first example:

set.seed(0)

x <- matrix(rnorm(15*3), 15)

treatment <- rep(c(T,F), length=15)

bootstrap2(x,statistic=colMeans,treatment=treatment,seed = 1)

the first example:

and the error-message:

 

Problem in bootstats: argument bootstrap.objects= not matched: bootstats(replicates = boot[[1]]$replicates - boot[[2]]$
replicates, observed
  = boot[[1]]$observed - boot[[2]]$observed, n = n, call = this.call,
 seed.start = boot[[1]]$seed.start, seed.end = boot[[2]]$seed.end,
 dim.obs = boot[[1]]$dim.obs, treatment = if(save.treatment) treatment,
 group = if(save.group) group, subject = if(save.subject) subject,
 parent.frame = parent.frame, label = boot[[1]]$label, defaultLabel =
 resampMakeLabel("bootstrap", this.call$data, this.call$statistic,
 treatmentNames = u), bootstrap.objects = boot)
Use traceback() to see the call stack
Problem occured while combining results for two samples.
 Saving results in .bootstrap2.partial.results


Do you Yahoo!?
All your favorites on one personal page ? Try My Yahoo!
<Prev in Thread] Current Thread [Next in Thread>
  • bootstrap2, Karin <=