comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » monte carlo analysis

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: monte carlo analysis [message #18343 is a reply to message #18279]

Tue, 21 December 1999 00:00

Andy Loughe
Messages: 174
Registered: November 1995

Senior Member

William Thompson wrote:

<snip>

> There was some discussion a
> while back about the properties of the random number generators in IDL. You
> might check out www.dejanews.com, for example

One concern is generating the initial seed.
Here is one simple method using the system clock...

seed = long((systime(1) - long(systime(1))) * 1.e8)

; Then for a Normal distribution:
num_rans = 100000
ran_nums = randomn(seed, num_rans)
stats = moment(ran_nums)
plot, histogram(ran_nums, binsize=.1), thick=3
print, '1)MEAN, 2)VARIANCE, 3)SKEWNESS, 4)KURTOSIS :', stats

Then there is no need to set the seed again, unless you
delvar, seed

--
| Andrew F. Loughe | email: afl@cdc.noaa.gov
| NOAA-CIRES CDC | web : www.cdc.noaa.gov/~afl
| 325 Broadway Mail Code R/E/CD1 | voice: (303)497-6211
| Boulder, CO 80303-3328 | fax : (303)497-7013
------------------------------------------------------------ -----
"We must believe in free will, we have no choice"-Isaac B. Singer

Report message to a moderator

[Message index]

		monte carlo analysis By: charliesolomon on Fri, 17 December 1999 00:00
		Re: monte carlo analysis By: htonishi on Wed, 22 December 1999 00:00
		Re: monte carlo analysis By: Andy Loughe on Tue, 21 December 1999 00:00
		Re: monte carlo analysis By: thompson on Sat, 18 December 1999 00:00

Previous Topic:	Marketplace for Technical Expertise
Next Topic:	Re: idl 5.3 runtime on idl 5.2 system work?

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Sat Nov 29 14:57:38 PST 2025

Total time taken to generate the page: 1.43859 seconds