[Rd] R binomial RNG stuck at 1 via Fortran call

Gilles GUILLOT gilles.guillot at inapg.inra.fr
Wed Nov 30 10:51:22 CET 2005


I have some trouble with the result of a fortran function calling the R
binomial RNG:

The C function rbinom is wrapped as in the file attached.

My main fortran program starts as
call rndstart()

and ends as
call rndend()  

I happen to call the binomial RNG within a loop as 

b = ggrbinom(1.d0,0.5d0)
write(*,*) 'b=',b 

In certain cases, after a few iterations in the loop,
b get stuck at 1

Any hint to explain that would help.


Gilles Guillot 
INRA - MIA Paris 

Currently working from Göteborg Stochastic Centre
Eklandagatan 86 - Rum1439 
Chalmers University of Technology
SE 412-96 Göteborg Sweden 
Phone +46 31 772 3514 / Fax +46 31772 3508


More information about the R-devel mailing list