comp.lang.idl-pvwave archive: archive

Home » Public Forums » archive » Re: Help needed!!

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

Re: Help needed!! [message #70506 is a reply to message #70504]

Fri, 16 April 2010 09:40

bala murugan
Messages: 23
Registered: April 2010

Junior Member

On Apr 16, 10:33 am, pp <pp.pente...@gmail.com> wrote:
> On Apr 16, 1:18 pm, bala murugan <bala2...@gmail.com> wrote:
>
>
>
>
>
>> My aim is to do the following,
>
>> To write an IDL routine that takes as input the mean of a Poisson
>> distribution (fLambda = a floating point number >= 0.0) and a number
>> of realizations (N), and generates N samples from the Poisson
>> distribution with mean fLambda. The output samples should be provided
>> as a vector (list) of integers (each >= 0). I will also want to plot
>> a histogram of the samples (IDL probably has a built-in histogramming
>> routine).
>
>> My subroutine/program will look something like this...
>
>> piSamples = PoissonDist(fLambda, N)
>> piSamples = vector of N integer samples returned by the routine
>> fLambda = mean of the Poisson distribution
>> N = number of samples to generate
>
>> This routine will need to loop n = 1..N Each time through the loop,
>> it will need to call RANDOMU to get a random number between 0.0 and
>> 1.0. It will then need to call IMSL_POISSONCDF (probably repeatedly)
>> to determine which integer in the Poisson distribution corresponds to
>> the random number gotten from RANDOMU. This integer is then placed in
>> the output vector, etc.
>
>> I am not using IMSL_POISSONCDF. Rather I am calculating the poisson
>> probability myself.
>
>> Can you please help me with this?
>
> Better now, but you still did not say what the problem is. One problem
> that is obvious is the line
>
> if (x EQ r) THEN a[i]=j
>
> since neither a or i exist at that point. And that condition is
> probably never going to be true.

@pp

Can you please go through the summary of what I want to accomplish and
suggest me some method to do it?

Report message to a moderator

[Message index]

		Re: Help needed!! By: R.G.Stockwell on Sat, 17 April 2010 10:07
		Re: Help needed!! By: R.G.Stockwell on Sat, 17 April 2010 10:04
		Re: Help needed!! By: jeanh on Sat, 17 April 2010 07:57
		Re: Help needed!! By: Craig Markwardt on Fri, 16 April 2010 20:59
		Re: Help needed!! By: bala murugan on Fri, 16 April 2010 12:46
		Re: Help needed!! By: bala murugan on Fri, 16 April 2010 12:19
		Re: Help needed!! By: Gray on Fri, 16 April 2010 12:13
		Re: Help needed!! By: bala murugan on Fri, 16 April 2010 12:11
		Re: Help needed!! By: R.G.Stockwell on Fri, 16 April 2010 11:59
		Re: Help needed!! By: bala murugan on Fri, 16 April 2010 11:06
		Re: Help needed!! By: Gray on Fri, 16 April 2010 10:53
		Re: Help needed!! By: R.G.Stockwell on Fri, 16 April 2010 10:48
		Re: Help needed!! By: Gray on Fri, 16 April 2010 10:47
		Re: Help needed!! By: bala murugan on Fri, 16 April 2010 09:40
		Re: Help needed!! By: bala murugan on Fri, 16 April 2010 09:39
		Re: Help needed!! By: penteado on Fri, 16 April 2010 09:33
		Re: Help needed!! By: Craig Markwardt on Fri, 16 April 2010 09:23
		Re: Help needed!! By: bala murugan on Fri, 16 April 2010 09:18
		Re: Help needed!! By: penteado on Fri, 16 April 2010 09:10
		Re: Help needed!! By: bala murugan on Mon, 19 April 2010 10:16

Previous Topic:	Re: GridData Conundrum
Next Topic:	Comparison operators and floating-point errors

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Wed Oct 08 16:17:49 PDT 2025

Total time taken to generate the page: 0.00416 seconds