Displaying 2 results from an estimated 2 matches for "ggrgam".
2005 Nov 30
1
RNG stuck via Fortran call
...chain in Fortran interfaced with R.
Each loop of my Fortran calls various functions of the R RNG through
the wrapper given below.
In a run of 100 iterations of the Markov kernel,
after 20 iterations, the RNG seems like frozen.
For example, the first call to the RNG in my loop is:
rpostlamb = ggrgam(dble(m),1.d0)
write(*,*) 'rpostlamb=',rpostlamb
after the 20th iteration, it will return the same value
rpostlamb= 1.24634557
for all the following interations.
Is there any suggestion of explanation for this strange fact ?
Thanks in advance
Gilles
R Version 2.2.0 compiled unde...
2005 Nov 30
2
R binomial RNG stuck at 1 via Fortran call
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