The first thing you are doing wrong is that you are not including a copy of cs for us to see ;).
Based on what you have written, I speculate that cs does not use the index correctly. if so then a simple, although inefficient, workaround is to rewrite cs: cs <- function(dataframe, index) { dataframe <- dataframe[index,] ... } Good luck, Andrew On Tue, Jan 17, 2006 at 10:02:49PM -0800, Ben Ridenhour wrote: > Hello, > I am new to using R and I am having problems get boot() to work properly. > Here is what I am trying to do: > > I have statistic called "cs". cs takes a data matrix (154 x 5) and > calculates 12 different scores for me. cs outputs the data as a vector (12 x > 1). cs doesn't really use weights, per se, however I have included this as > one of the 2 arguments cs can take. > > I try performing a bootstrap by issuing: > myout<-boot(data, cs,R=999) > I have tried other versions where I specify stype="w", etc... > > The problem I get is that the dataset does not seem to be resampled. I end > up with 999 replicates that have the exact same value of the output of cs. > > In the end I have something like > > Bootstrap Statistics : > original bias std. error > t1* 0.865122275 1.698641e-14 0 > t2* -0.005248414 -9.627715e-17 0 > t3* -0.052833740 -8.812395e-16 0 > t4* 0.807040121 1.287859e-14 0 > t5* 0.542082588 -9.103829e-15 0 > t6* -0.018617838 -7.285839e-17 0 > t7* 0.006409704 1.422473e-16 0 > t8* 0.529874453 8.104628e-15 0 > t9* 0.074804390 2.359224e-16 0 > t10* -0.007153634 1.301043e-16 0 > t11* -0.018241243 -2.359224e-16 0 > t12* 0.049409513 -1.200429e-15 0 > > Clearly the bootstrap is not working. What am I doing wrong? > > Thanks, > Ben > > > ------------------------------- > Benjamin Ridenhour > School of Biological Sciences > Washigton State University > P.O. Box 644236 > Pullman, WA 99164-4236 > Phone (509)335-7218 > -------------------------------- > "Nothing in biology makes sense except in the light of evolution." > -T. Dobzhansky > > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Andrew Robinson Department of Mathematics and Statistics Tel: +61-3-8344-9763 University of Melbourne, VIC 3010 Australia Fax: +61-3-8344-4599 Email: [EMAIL PROTECTED] http://www.ms.unimelb.edu.au ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html