similar to: pnorm, relative accuracy in the tails

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

Derar R-people I have troubles understanding what the function sort.list() is doing. On the homepage it says that it returns a permutation which rearranges a vector into ascending or descending order (like order() but on a vector instead of a sequence). > sort.list(c(0, 2, 10, 11, 4)) [1] 1 2 5 3 4 which does not make sense to me. In fact I am getting the same (non-sensical) result using

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-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

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"

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 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

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.

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

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, 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

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",

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

Dear R-people Does someone have a routine to ngenerate inverse-gaussian random numbers. I am thinking of something similar to rinvgauss, pinvgauss etc. in S-plus. best regards Helgi -- Helgi Tomasson FAX: 354-552-6806 University of Iceland PHONE:354-525-4571 Faculty of Economics and Business Administration

Dear R People, I getting some curious behaviour with the integrate function. Consider the function foo <- function(x) { if(x==0) 4 else 0 } I get the error > integrate(foo,0,1) Error in integrate(foo, 0, 1) : evaluation of function gave a result of wrong length But now consider foo <- function(x) { ifelse(x==0,4,0) } > integrate(foo,0,1) 0 with absolute error

Dear R Community, We have recently installed R 1.5.0 on several workstations. Ever since, if I try to extract certain numbers or lines from a file (which contains many NA's), I don't only get my desired line but all other lines which contain NA's. In a file with 15000 rows and 7500 NA's this is not very practical. In the earlier R-Version this worked when I typed in the

I'm a PhD student and I'm working with covariance function. I'm interested to know if exist some packages in R to calculate and plot the bidimensional Autocovariance Function. the input matrix is a matrix that describe a spatial location over a 2-D space and I want to use it in the same way I can use a time serie in the 1-D acf. Thanks, Nicola.

I'm a PhD student and I'm working with covariance function. I'm interested to know if exist some packages in R to calculate and plot the bidimensional Autocovariance Function. the input matrix is a matrix that describe a spatial location over a 2-D space and I want to use it in the same way I can use a time serie in the 1-D acf. Thanks, Nicola.

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,