similar to: function sort.list()

Dear R people The function below should be decreasing, convex, and tend to zero when x tends to infinity. curve((1-pnorm(x))/dnorm(x),from=0, to=9) >From the plot we see that for x between 8.0 and 8.3 the function is fluctuating. As far as I understand, this is due to the function pnorm() not being sufficiently accurate in the tails. I am using pnorm() in a way that has probably not been

Dear R people We have two packages, where the first package (geoR) is required by the second (geoRglmm). Both packages have functions calling C-code via .C(). We would like to call C-functions included in the first package from within the C-code in the second package. Is this possible? An appropriate header file was included in the src directory of the second package but this alone did not

Dear R-devel During the last half a year I have several times encountered the following problem with optim() when using method= "L-BFGS-B". The function return a value which is clearly not the maximum (seen from printing the value each time the function is called). Some output is shown below. A few things I have observed (as I remember it): a. The problem seems to occur when the

Dear R-developers There is a discrepancy between the variables found by typing names(options()) and the variables described on the help page for the function ``options''. The following variables are in names(options()) but not described on the help page : "error.messages" "repositories" "show.coef.Pvalues" "help.try.all.packages"

Hi, Suppose I have created something like this in R: foo <- "myfoo" and I want to find out the number of character in foo (in other words, R should return 5 since "myfoo" has 5 charactors. How can I do it? I tried: length(foo) but it returned 1. Cheers, Kevin ------------------------------------------------------------------------------ Ko-Kang Kevin Wang

Dear R-list I am having troubles understanding how the function is.numeric() works. Any help appreciated. Some commands are given here (the output and my comments are given below) vec <- c(1.4, NA, NA, NA) sapply(vec,FUN=is.numeric) is.numeric(vec[2]) is.na(vec[2]) is.numeric(NA) vec <- c(TRUE, FALSE, NA, NA) sapply(vec,FUN=is.numeric) is.numeric(vec[2]) is.numeric(vec[3])

Dear friends, I want to sort a dataset according to one or two variables in the dataset, i thought sort could do it , but failed. e.g. x <- c(2, 9, 18, 3, 2) y<-c(2,5.6,5,9,8) z<-c(21,5,5,19,7) a<-cbind(x,y,z) a x y z [1,] 2 2.0 21 [2,] 9 5.6 5 [3,] 18 5.0 5 [4,] 3 9.0 19 [5,] 2 8.0 7 I want to sort dataset a according to ascending x and descending y, How to

Dear R-users, I wondered if anyone has extended GLMMGibbs to include correlated random effects, and if so, whether they would be willing to let me use their code? Jonathan Myles has no plans to extend glmm in this manner within the foreseeable future. With thanks, Patty -- -------------------------------------------------------------------------------- Assoc Prof Patty Solomon

Hi, Is it possible to run R in command line to evalute R expressions and return results to stdout, something like >R CMD -e "R.version$minor" then you got return >"8.1" Or do a simple calculation >R CMD -e "sin(1.2)" >0.932039 Thanks. -- Pingping Zheng Department of Mathematics and Statistics Fylde College Lancaster University Lancaster LA1 4YF

Dear all, I'm new on R and this mailing list. We work on spatial rainfall estimation with R and Grass. We have a problem with the krige function from the sgeostat package. We would like to limit the estimated area with a polygon limit. I use a 50 points polygon to describe my work area. The krige function work quiet well without limit. But if I use this option I have the following error

After "image" and "persp" would it be interesting to als have contour as a generic function? Cheers Paulo Justiniano Ribeiro Jr Dept Maths & Stats - Fylde College Lancaster University Lancaster LA1 4YF - U.K. e-mail: Paulo.Ribeiro@est.ufpr.br http://www.maths.lancs.ac.uk/~ribeiro -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Suppose I have a nx2 matrix of data, X, the following code generate density estimation for each column and plot them denlist <- apply(X, 2, density) par(mfrow=c(1,2)) lapply(denlist, plot) Does anyone know how to change the main title of each density plot to "var 1", "var 2" by passing optional argument "main"? I've tried lapply(denlist, plot,

Dear all, I wish use the "est.variogram" function under R. My question concerns the determination of number of lags. To make this calculation, i determined the biggest and the smallest distance between 2 pairs of points. I built the following tolerance : tol=mindistance/2 and i defined 1 lag as : 1lag=mindist+tol. Then i made a buckle the step of stop of which the maximal distance.

I have been trying the integrate function in R, a function which would be very useful for a current project of mine. But I am encountering errors integrating the one function I have tried. The function to be integrated is a product of a gamma demsity and a normal density: gamma.by.normal_function(y,x,shape,scale,stdev) return( dgamma(y,shape=shape,scale=scale)*

Hi, I have problems when calling fft_factor and fft_work directly in my C-code. Martin Maechler already told me that I should not do this -- anyway, the code works fine if the length of the vector is a multiple of 2,3,5. If the length of the vector has different form I get, in some cases, a segmentation fault when work and iwork are freed (I use C's malloc to allocate memory for work and

Hi R users: I want to know how can I easily sort in descending order a whole data frame by one or several variables. Thank you very much for your help. Kenneth Cabrera krcabrer at perseus.unalmed.edu.co krcabrer at epm.net.co Universidad Nacional de Colombia, Sede Medell?n -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello R-Users, I have I simple question that I could not solve: How to extract the elevation and the slope values from a linear model (lm) separately? coef(model.lm) gives both of them. Thanks, Antonio Olinto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Hi all, I face the following problem. I have two variables x,y. I have also calculated the surface z. When I am drawing the surface with function >levelplot (in library lattice) everything is OK. I have the legend etc ... Then I am trying to plot x,y on the same graph. So I am doing : >points(x,y) But, the points are not in the correct position. They are a little bit shifted to the

> Date: Thu, 9 Nov 2000 15:00:20 +0000 (GMT) > From: "Paulo J. Ribeiro Jr" <p.ribeiro at lancaster.ac.uk> > I'm trying to create a dll file for the windows distribution of a package. > > I've tried mingw32 and cygwin with the commands given below and it seems > that there is a problem related with C random() function. > (-mno-cygwin only included

I'm running a simulation that uses nls on each run. Sometimes, not unexpectedly, nls doesn't return parameter estimates (e.g., "singular gradient" error). Of course when this happens the simulation (in a loop) stops and no results are recorded. What I would like is for the simulation to simply record that nls couldn't obtain estimates on a particular run due to whatever