Re: [Rd] Problem with rnorm ?

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Sat, 12 May 2007 12:27:33 +0100 (BST)

  1. You are linking an executable against -lR. You can't do that as here, for R is never initialized (and I very much doubt you will be able to give a reference that recommended you to do that). You need to link against standalone libRmath, or start up R properly if you really need R.
  2. The arguments to rnorm() are doubles, not integers: did you not get a warning? I would expect automatic conversion to occur, but it is better to get this right.

On Fri, 11 May 2007, Tong Wang wrote:

> Hi,
> I couldn't get the rnorm() work right in C, for example, the following code produce strange results.
>
> #include <stdio.h>
> #include <Rmath.h>
>
> int main(void){
> double x[3];
> for(int i=0;i<3;i++) x[i]=rnorm(0,1);
> printf("%lf,%lf,%lf",x[0],x[1],x[2]);
> return 0;
> }
>
> output : -8.773321,-8.773321,-8.773321
> compiling script: gcc -std=gnu99 -Wall -IC:/PROGRA~1/R/R-24~1.1/include test.c -LC:/PROGRA~1/R/R-24~1.1/bin -lR -o test
>
> Could someone tell me what's wrong with my code ? Thanks a lot.
>
> ______________________________________________
> R-devel_at_r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>

-- 
Brian D. Ripley,                  ripley_at_stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Received on Sat 12 May 2007 - 11:30:14 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 13 May 2007 - 04:33:04 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.