similar to: Calling qnorm from Fortran

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() call dblepr("X was", 5, x, 1) call rndend() end but I don't understand

I would like to call C routines from Fortran 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() call dblepr("X was", 5, x, 1) call rndend() end but I don't understand the purpose

Hi all, Is there a way that I can import R functions into Fortran? Especially, I want to generate random numbers from some not-so-common distributions (e.g. inverted chi square) but did not find any routines written in Fortran that deal with distributions other than uniform and normal. Thanks. Anny [[alternative HTML version deleted]]

Dear R Developers, I should have had the Design library running in R by now but have kept putting off changing some calls to LINPACK routines to use those builtin to R. Specifically I call dqrsl1 and dqr. Would it be an easy task to put those in the next release of R? If not I'll finally bite the bullet and get back into reading LINPACK documentation (which I have but haven't examined

I found the following strange behavior using qnorm() and pnorm(): > x<-8.21;x-qnorm(pnorm(x)) [1] 0.0004638484 > x<-8.22;x-qnorm(pnorm(x)) [1] 0.01046385 > x<-8.23;x-qnorm(pnorm(x)) [1] 0.02046385 > x<-8.24;x-qnorm(pnorm(x)) [1] 0.03046385 > x<-8.25;x-qnorm(pnorm(x)) [1] 0.04046385 > x<-8.26;x-qnorm(pnorm(x)) [1] 0.05046385 > x<-8.27;x-qnorm(pnorm(x))

I am trying to replicate part of the function of the S-Plus terms.inner function. I don't need an entire model terms object but R expressions or character strings containing the innermost variable name given a vector of character strings or expressions. Here are some example inputs that each should result in the string or expression "x": x; x^2; g(x); h(g(x)); h(g(x^3));

Great news! Having spent a significant amount of time porting/mantaining some packages on S4/Splus 5.x and R, I'm interested in your experience as you port Hmisc and Design; I'd like to know of your progress in order to help improve available R/Splus portability tools both at the C and the S language code. These tools include the package SLanguage written by John Chambers, and

I am using read.spss in the foreign package to read an SPSS save file. For date variables I get huge values such as 11489990400. Does anyone know how to convert these values to R POSIXct date objects? Thanks in advance -Frank platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor

how R implement qnorm() I wonder anyone knows the mathematical process that R calculated the quantile? The reason I asked is soly by curiosity. I know the probability of a normal distribution is calculated through integrate the Gaussian function, which can be implemented easily (see code), while the calculation of quantile (or Zα) in R is a bit confusing as it requires inverse error function (X

The Hmisc and Design libraries have been updated respectively to versions 1.4-2 and 1.1-1. New versions for Linux/Unix/Windows may be obtained from http://hesweb1.med.virginia.edu/biostat/s/library/r . Web sites for the libraries are http://hesweb1.med.virginia.edu/biostat/s/Hmisc.html and http://hesweb1.med.virginia.edu/biostat/s/Design.html . Thanks to Xiao Gang Fan for porting the libraries

The Hmisc and Design libraries have been updated respectively to versions 1.4-2 and 1.1-1. New versions for Linux/Unix/Windows may be obtained from http://hesweb1.med.virginia.edu/biostat/s/library/r . Web sites for the libraries are http://hesweb1.med.virginia.edu/biostat/s/Hmisc.html and http://hesweb1.med.virginia.edu/biostat/s/Design.html . Thanks to Xiao Gang Fan for porting the libraries

Hello all, I have a question about basic statistics. Given a PDF value of 0.328161, how can I find out the value of -0.625 in R? It is like reversing the dnorm function but I do not know how to do it in R. > pdf.xb <- dnorm(-0.625) > pdf.xb [1] 0.328161 > qnorm(pdf.xb) [1] -0.444997 > pnorm(pdf.xb) [1] 0.628605 Many thanks, Edwin -- View this message in context:

Trenkler, Dietrich said: > > I found the following strange behavior using qnorm() and pnorm(): > > > x<-8.21;x-qnorm(pnorm(x)) > [1] 0.0004638484 > > x<-8.28;x-qnorm(pnorm(x)) > [1] 0.07046385 > > x<-8.29;x-qnorm(pnorm(x)) > [1] 0.08046385 > > x<-8.30;x-qnorm(pnorm(x)) > [1] -Inf > qnorm(1-.Machine$double.eps) [1] 8.12589

The Design library has been fully ported to R except for Cox proportional hazards regression modeling (using Therneau's survival package) which will be available in about two weeks. It will take much longer to make all the example code executable, is it currently contains many examples for which data are not provided. Thanks to Xiao Gang Fan <xiao.gang.fan1 at libertysurf.fr> who

The Design library has been fully ported to R except for Cox proportional hazards regression modeling (using Therneau's survival package) which will be available in about two weeks. It will take much longer to make all the example code executable, is it currently contains many examples for which data are not provided. Thanks to Xiao Gang Fan <xiao.gang.fan1 at libertysurf.fr> who

In S-Plus the ... argument for tapply is before simplify. In R (1.2.3) ... is after simplify, causing the user to have to specify simplify every time ... is used even though the default is usually OK. Unless there was a benefit to reordering the arguments, it might be a good idea to reconsider the current function definition. ... is probably used more than simplify (e.g., tapply(rainfall,

On Mon, 16 Jun 2003 09:56:07 +0200 Patrick Hausmann <p.hausmann at mr-regionalberatung.de> wrote: > Dear Frank, > > I can reproduce the problem from Philippe Glaziou (see attachment). > Maybe this can help. > Best, > Patrick > > Am Wollelager 11 > 27749 Delmenhorst > Tel. 04221 96373-0 > Fax 04221 96373-29 > > http://www.mr-regionalberatung.de >

In platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 7.1 year 2003 month 06 day 16 language R I get > sort(c(3,1,NA)) [1] 1 3 Shouldn't NAs be retained by default? Thanks -Frank ---

Hello, could You please help: I am looking for a way to formulate test accuracy measures such as test sensitivity, specificity, predictive values, and correct classification rate using p/d/qnorm. The tests' primary values follow a bimodal distribution, which is modelled by a mixture of two normal distributions: p * dnorm ((x - u1) / s1) / s1 + (1 - p) * dnorm ((x - u2) / s2) / s2)

The Hmisc and Design libraries have had major improvements. New versions for Linux/Unix/Windows may be obtained from http://hesweb1.med.virginia.edu/biostat/s/library/r Thanks to Xiao Gang Fan for porting the libraries to Windows once again. Change logs may be obtained from http://hesweb1.med.virginia.edu/biostat/s/library/common Several bugs have been fixed in Design, and many bugs have been