From gilles.guillot at inapg.inra.fr  Wed Nov 30 10:51:22 2005
From: gilles.guillot at inapg.inra.fr (Gilles GUILLOT)
Date: Wed, 30 Nov 2005 10:51:22 +0100
Subject: [Rd] R binomial RNG stuck at 1 via Fortran call
Message-ID: <200511301051.22383.gilles.guillot@inapg.inra.fr>

Hi, 

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 


-- 
???????????????????????????????????????????????????????????????????
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

www.inapg.inra.fr/ens_rech/mathinfo/personnel/guillot/welcome.html
???????????????????????????????????????????????????????????????????



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