similar to: na.omit

Dear all, There has been some time since I asked about interpolation in higher (>2) dimensions, and I must admit I failed to write a function to do this myself the last time, but eventually ended up doing it in MATLAB. I tried to translate the MATLAB code, but MATLAB code is so much more opaque than R (S) code, so I failed that too, mainly because I could only get one MATLAB session, I would

Dear all, I'm writing a lgrind definition for R. Lgrind is a program to produce pretty program listings for inclusion in LaTeX documents. Does anybody here have any experience with lgrindefs? My definition now looks like this: R:\ :pb=^\d?\p\d<-\dfunction\(\a\):\ :bb={:be=}:cb=\d#:ce=$sb=":se=\e":lb=':le=\e':id=.:\ :zb=@:ze=@:tb=%%:te=%%:mb=%\$:me=\$%:vb=%\|:ve=\|%:\

Dear all, I'm writing a lgrind definition for R. Lgrind is a program to produce pretty program listings for inclusion in LaTeX documents. Does anybody here have any experience with lgrindefs? My definition now looks like this: R:\ :pb=^\d?\p\d<-\dfunction\(\a\):\ :bb={:be=}:cb=\d#:ce=$sb=":se=\e":lb=':le=\e':id=.:\ :zb=@:ze=@:tb=%%:te=%%:mb=%\$:me=\$%:vb=%\|:ve=\|%:\

Dear all, I am interested by the various messages sent to the r-help mailing list these last few days about documentation on R. Yet, there are a lot of issues actually open on this topic, and it is not clear to me what is the "urgent need for documentation" on R. R comes with the 'Introduction to R' by Venables et al. I recommend the colleagues and students in my lab who are

Hi Rers, W/o trying to make a long story short (I hope you read that correctly), I have been working on a nice little R function ... Somewhere along the line I got an error that refers to a max function in my function (I think this is a vector vs. scalar issue): Error in max(..., na.rm = na.rm) : invalid "mode" of argument If that's not bad enough, even when I q() R and say

Dear all, I have been trying to understand what runif() is telling me. I am generating lots of numbers (billions and billions (wow, I''ve dreamed about saying that for many years... :-) ), for a distribution that has the following quantile function: 1 / (2 * sqrt(1 - p)) (that is, the distribution has a lower cutoff) As you can imagine, this has rather heavy upper tail. I was

Dear all, I think I must be going crazy.... If I do > mean(rpois(1000000, 14)) I get: [1] 13.50420 and again: > mean(rpois(1000000, 14)) [1] 13.49896 > mean(rpois(1000000, 14)) [1] 13.50161 > mean(rpois(1000000, 15)) [1] 14.49250 > mean(rpois(1000000, 15)) [1] 14.49897 > mean(rpois(1000000, 14.5)) [1] 13.99689 > mean(rpois(1000000, 14.5)) [1] 13.9963 This was on a platform

Full_Name: Kjetil Kjernsmo Version: 1.0.1 OS: osf1 Submission from: (NULL) (129.240.28.227) Ouch, another obscure bug... I just had a bus error crash, with a core dump, a small one this time. Usual <URL:http://www.astro.uio.no/~kjetikj/tmp/core>. I had done a few simulations, and done a chisq.test() on them. These simulations also write about 190MB's to file, but I have no idea whether

Dear all, Me and the Garbage Collector doesn't understand each other... I'm on R-1.0.1 (I have requested an upgrade to the right local people), an osf1 system with an alphaev6 CPU. I have fair amounts of memory to use, but it seems to be used up pretty fast nevertheless. I'm now trying to see what the garbage collector does, I have inserted a gc(v=T) call at a spot in my code. I

Hi I'm looking to a function in R equivalent to S-PLUS' "interp". Is there any implementation in R? If not, is there any similar interpolation function such that at the data points the interpolated values concides with the data? Thanks in advance P.J. Paulo Justiniano Ribeiro Jr Dept Maths & Stats - Fylde College Lancaster University Lancaster LA1 4YF - U.K.

Dear R members, Is there a way I can include more than one variable in a histogram? Many thanks in advance, Liqing, Liqing Zhang Department of Ecology & Evolutionary Biology 321 Steinhaus Hall U.C. Irvine Irvine, CA 92697-2525 Phone:(949)-824-7703 Fax: (949)-824-2181 Email: lzhang at uci.edu ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

How would you proceed with the following problem: I have an data frame containing a table relating raw scores on an iq test to 'decile-like' scores. data<-data.frame(age =c( 8, 9, 10, 11, 12, 13, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70), pers90=c(33, 39, 43, 47, 50, 54, 54, 52, 50, 49, 47, 46, 43, 42, 39, 37), perc75=c(27, 31, 37,

Hi I'm writing my MSC and I've used R to do the data analysis. Now I want to make a reference to the program but I don't know how ! I thought of something like: Anonymous. 2001. R - Data Analysis Environment. Version 1.2.3. www.r-project.org Can someone help me with this ? Thanks EJ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing

Dear all, I'm a bit confused about the output from functions that generates random values, e.g. rpois. I'm using 0.65.1 on Digital UNIX alphaev6. If I say > rpois(10,5) output is not unexpectedly: [1] 4 6 5 7 6 5 2 2 5 3 but I figured that if I go: > rpois(10,1:10) I would get a 10x10 matrix as output, with random values for each of the vector elements in the row vectors, but

Sorry, I pasted the wrong file in earlier... this is the correct one: pValCalculator(b, n=20, m=20) { ind <- 1:min(c(n,m)) prob <- (1-pnorm(b,sd=std*sqrt(ind))) prob1 <- sum((n-ind+1)*(m-ind+1)*prob) prob1 } inputData <- scan("/users/lvssso/projects/LAMA/output/pValLamaScore.tmp", list(block1 = "",block2 = "",width1 = 0,width2 = 0,alignment = 0,score

I am trying to install R on a unix workstation (Sun-SPARC) but I am not sure which of the binaries to use. The only binaries available are for linux, nextstep, osf, and windows. Which one of them will work for a SUN sparc station?? Thanks. PS. I tried to compile the source code but it looks like I have a probelm with fortran. So, if I could install one of the binaries it would be a lot

How do you obtain two plots on the same figure? for example plot(rnorm(100) plot(rnorm(100),type="l") -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To:

Dear all, Either I'm being totally blind now (as you know, the eyes are the first thing that is blinded :-)), or something is missing in "An Introduction to R". The section heading of section 5.7.1 is "Multiplication, inversion and solving linear equations", and since I was looking for the canonical way of inverting a matrix, this was the place to look. However, I found

Hi all, how can I remove objects from the workspace which starts with a certain pattern, e.g. lm (lm1, lm2, lm3 etc). Wildcards won?t work. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Dear all, I have a few files with binary data written by a C program a friend wrote. I allready have program to read these files, Thomas Lumley and Prof Brian D Ripley was kind enough to respond to a question to this list earlier with some code that works perfectly, so this is really no problem, it is more out curiousity. Most of my files consists of 40000 single precision floats, and I figured