Re: [Rd] package Geneland / Rgui under windows

2007-10-15 Thread Gilles GUILLOT
I forgot to say that I was the package author. I suspected a bug of R as Geneland worked fine for two years and the problem popped up with the release of R 2.6.0 And I can't see any explanation. So, any clue would help. gilles > Gilles GUILLOT wrote: >> Hi, >> I experience

[Rd] package Geneland / Rgui under windows

2007-10-15 Thread Gilles GUILLOT
Hi, I experienced a problem with the package Geneland under R 2.6.0 with windows XP professional. The commands below should simulate a dataset, then make an MCMC simulation stored in tempdir(). It works with R 2.5.1 (both GUI and command line) It works with the command line of R 2.6.0 but not wi

[Rd] package Geneland / Rgui under windows (PR#9964)

2007-10-10 Thread guillot
Full_Name: Gilles Guillot Version: 2.6.0 OS: windows XP professional Submission from: (NULL) (129.240.88.50) This sequence of command does not work in the R gui 2.6.0: library(Geneland) set.seed(1) data <- simdata(nindiv=200, coord.lim=c(0,1,

[Rd] RNG stuck via Fortran call

2005-11-30 Thread Gilles GUILLOT
Having not much success with my previous question I try to reformulate it: I'm simulating a Markow 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 iterati

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

2005-11-30 Thread Gilles GUILLOT
wrapper now as part as the message, sorry about previous mail #include #include void F77_SUB(rndstart)(void) { GetRNGstate(); } void F77_SUB(rndend)(void) { PutRNGstate(); } double F77_SUB(ggrnorm)(double *mu, double *sigma) { return rnorm(*mu, *sigma); } double F77_SUB(ggrexp)(double *

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

2005-11-30 Thread Gilles GUILLOT
wrapper now attached -- ÅÅÅ 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

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

2005-11-30 Thread Gilles GUILLOT
(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 fro

[Rd] R build under mandriva 10.2

2005-11-08 Thread Gilles GUILLOT
After upgrading from mandrake 10.1 to mandriva 10.2 I can't build shared archive with R 2.2.0 [EMAIL PROTECTED] src]$ R CMD SHLIB main.f sub.f wrapper.c gcc -shared -L/usr/local/lib -o main.so main.o sub.o wrapper.o -lg2c -lm -lgcc_s /usr//bin/ld: cannot find -lg2c collect2: ld returned 1 exit

Re: [Rd] pb with dyn.load - fortran code now attached

2005-10-27 Thread Gilles GUILLOT
> > is.loaded("supsmu_") > > [1] FALSE > > > is.loaded("supsmu") > > [1] TRUE > > That is a Fortran entry point, and it complies with the description > quoted. > > What is.loaded() needs depends on how the symbol would be found > and so it is not much use. OK, thanks. But still > is.loaded("fst

[Rd] pb with dyn.load - fortran code now attached

2005-10-26 Thread Gilles GUILLOT
# QUESTION 1: [EMAIL PROTECTED] guillot]$ R R : Copyright 2005, The R Foundation for Statistical Computing Version 2.2.0 (2005-10-06 r35749) ISBN 3-900051-07-0 >system("R CMD SHLIB ~/tmp/test1.f") g77 -fPIC -g -O2 -c /home/guillot/tmp/test1.f -o /home/guillo

[Rd] pb with dyn.load

2005-10-26 Thread Gilles GUILLOT
# QUESTION 1: [EMAIL PROTECTED] guillot]$ R R : Copyright 2005, The R Foundation for Statistical Computing Version 2.2.0 (2005-10-06 r35749) ISBN 3-900051-07-0 >system("R CMD SHLIB ~/tmp/test1.f") g77 -fPIC -g -O2 -c /home/guillot/tmp/test1.f -o /home/guillo

[Rd] calling fortran functions CHOL and DPOTRF form Fortran

2005-07-05 Thread Gilles GUILLOT
? Thanks in advance, Gilles -- _________ Gilles GUILLOT INRA -Département Mathématiques et Informatique Appliquées Unité Mixte de Recherche INRA - INAPG - ENGREF Institut National Agronomique de Paris-Grignon 16, rue Claude Bernard 752

Re: [Rd] Calling C from Fortran

2005-06-15 Thread Gilles GUILLOT
at > Fortran <-> C conventions. > > Reid Huntsinger > > -Original Message- > From: Gilles GUILLOT [mailto:[EMAIL PROTECTED] > Sent: Wednesday, June 15, 2005 3:50 AM > To: Huntsinger, Reid > Subject: Re: [Rd] Calling C from Fortran > > > Thanks Reid! >

[Rd] Calling C function from Fortran

2005-06-15 Thread Gilles GUILLOT
Hi all, the example in the R doc and the hints from Shusong Jin , Ingmar Visser and Reid Huntsinger (thanks all three) refer to the case where the function does not have arguments. I'm still looking for a proper sequence of commands to call C functions with arguemnts from R. Imagine I want

[Rd] Calling C from Fortran

2005-06-14 Thread Gilles GUILLOT
I would like to call C routines from Fortran under linux as suggested in section 5.6 of the "Writing R extensions" documentation. I'm familiar with Fortran but not with C. I understand the example provided in Fortran: subroutine testit() double precision normrnd, x call rndstart() x = normrnd()