[Numpy-discussion] problem with poisson generators

[Flavio C. <fcc...@gm...>]

Hi,

I am having problems with the poisson random number generators of both=20
Numarray and Numeric.
I can't replicate it when calling the function from the python cosonle, but=
=20
it has consistently malfunctioned when used within one of my scripts.

What happens is that it frequently return a value greater than zero when=20
called with zero mean: poisson(0.0)

Unfortunately My program is too big to send attached but I have confirmed=
=20
the malfunction by printing both the mean and the result whenever it spits=
=20
out a wrong result.

This is very weird indeed, I have run poisson millions of times by itsel on=
=20
the python console, without any problems...

I hope it is some stupid mistake, but when I replace the poisson function=
=20
call within my program by the R equivalent command (rpois) via the rpy=20
wrapper, everything works just fine...=20

any Ideas?



--=20
Fl=E1vio Code=E7o Coelho
registered Linux user # 386432
---------------------------
Great spirits have always found violent opposition from mediocrities. The
latter cannot understand it when a man does not thoughtlessly submit to
hereditary prejudices but honestly and courageously uses his intelligence.
Albert Einstein