Re: [Rd] Identifying graphics files produced by R

2009-02-14 Thread Romain Francois
Hi, I don't know the answer but here is how I would try to get it from the source: $ cd R-devel/src $ grep -R "R Software" * Binary file library/grDevices/src/grDevices.so matches library/grDevices/src/devPS.c:fprintf(fp, "Creator: R Software\n"); Binary file library/grDevices/src/dev

Re: [Rd] [PATCH] typo in R-lang

2009-02-14 Thread Martin Maechler
Thank you, Peter; I just have committed the fix for the typo. (The inline patch was perfectly sufficient) Martin Maechler, ETH Zurich > "PC" == Peter Cowan > on Fri, 13 Feb 2009 20:02:28 -0800 writes: PC> Here is a patch for a small typo in the description of do.call() PC>

Re: [Rd] Generate random numbers in Fortran

2009-02-14 Thread Fabio Mathias
As I am wanting to generate a beta, then I created a function in C to generate a beta, but the problem appears when I go to compile My function in C is #include #include #include void F77_SUB(myrbeta)(double* px) {     GetRNGstate();     *px = rbeta(1.00,3.00);     PutRNGstate(); }

Re: [Rd] Generate random numbers in Fortran

2009-02-14 Thread Prof Brian Ripley
There are lots of invalid characters in your mail, but Fortran code starts in column 7 and that is what the compiler is telling you you have not done. (Some dialects will allow tabs instead, but I see no sing of those either.) If you are that unfamiliar with Fortran, why not just use C? On S

Re: [Rd] Generate random numbers in Fortran

2009-02-14 Thread Kjell Konis
Sorry, my example was just to give you a rough idea of how to do this. Here is an example that works. myrandom.c: #include #include void F77_SUB(fseedi)(void) { GetRNGstate(); } void F77_SUB(fseedo)(void) { PutRNGstate(); } void F77_SUB(myrbeta)(double* px, double* pa, double* pb)

Re: [Rd] proposed simulate.glm method

2009-02-14 Thread Nicholas Lewin-Koh
Hi, For extended glms such as gams, gnm or other distributions such as negative binomial, would there need to be a separate simulate method? Or, could the current framework, rather than stopping with an error look for the appropriate model matrix, coefficients, distribution function and family obje

Re: [Rd] proposed simulate.glm method

2009-02-14 Thread Martin Maechler
> "NicLK" == Nicholas Lewin-Koh > on Sat, 14 Feb 2009 08:34:45 -0800 writes: NicLK> Hi, For extended glms such as gams, gnm or other NicLK> distributions such as negative binomial, would there NicLK> need to be a separate simulate method? Not necessarily, as I said, th

Re: [Rd] proposed simulate.glm method

2009-02-14 Thread Nicholas Lewin-Koh
Hi, Well, my question wasn't that clear :-), but yes you mostly answered it. I guess the one case I would be concerned is in Heather's code, where the distribution to simulate from is chosen, that seemed to be hard coded. So if I built a family object, say for a model that assumes errors from a zi

Re: [Rd] Generate random numbers in Fortran

2009-02-14 Thread Ben Bolker
Wacek Kusnierczyk idi.ntnu.no> writes: > > you can always try to get hold of the extensive nag fortran libraries: > > http://www.nag.co.uk/numeric/fl/FLdescription.asp > > comsider also 'numerical recipes' by press et al., of which there are > fortran, c, and c++ editions (i think there was a

Re: [Rd] Generate random numbers in Fortran

2009-02-14 Thread Ben Bolker
Ben Bolker ufl.edu> writes: > > vQ > > If one can get the R routines to work, I think they have some > advantages over NAG and Numerical Recipes routines: > > * source code is freely redistributable > * I can't really claim expertise, but I believe there > has been some controversy (see t

Re: [Rd] proposed simulate.glm method

2009-02-14 Thread Prof Brian Ripley
On Sat, 14 Feb 2009, Nicholas Lewin-Koh wrote: Well, my question wasn't that clear :-), but yes you mostly answered it. I guess the one case I would be concerned is in Heather's code, where the distribution to simulate from is chosen, that seemed to be hard coded. Rather, all known glm famil

[Rd] SET_VECTOR_ELT and STRSXPs

2009-02-14 Thread Kasper Daniel Hansen
It seems (based on the NEWS file and on output from R CMD check) that we may no longer use SET_VECTOR_ELT on STRSXPs. So I guess that section 5.14 of R-extensions needs to be updated, the current phrasing is "By default a certain amount of misuse is allowed where the internal representation