comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » "bootstrap" statistics
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
"bootstrap" statistics [message #30849] Mon, 20 May 2002 12:41 Go to previous message
wmconnolley is currently offline  wmconnolley
Messages: 106
Registered: November 2000
Senior Member
Hello group. I want to do what I think of as "bootstrap" statistics, viz
given a timeseries I take a random subsample (with, say, half the number
of elements), compute some statistic (say, then mean); then take another
random subsample; then again lots of times (say 1000 or 10000) and end
up with a distribution of the statistic concerned.

So: to do this I need a means to generate n/2 random (non-repeating)
indices from 0...n-1. At the moment I do this by:
once I have m indices I generate one more at
random; see if its in the list of m; if not, good; if it is, generate another
one. This is hideously inefficient and slow: there *must* be a better way.

I have found:

http://www.astro.washington.edu/deutsch-bin/getpro/library14 .html?PERMUTE

which is a somewhat better way, but still slow. Is there a no-loops version?

-W.

--
William M Connolley | wmc@bas.ac.uk | http://www.nerc-bas.ac.uk/icd/wmc/
Climate Modeller, British Antarctic Survey | Disclaimer: I speak for myself
I'm a .signature virus! copy me into your .signature file & help me spread!
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Re: Job offer for IDL Programmers near Frankfurt/Germany
Next Topic: peculiar things with Z-device

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Wed Oct 08 13:55:13 PDT 2025

Total time taken to generate the page: 0.00591 seconds