similar to: Package GSS for interpolation in more than 2D?

Full_Name: Kjetil Kjernsmo Version: 0.65.1 OS: Digital UNIX 4.0 Submission from: (NULL) (129.240.28.172) I'm experiencing problems with rpois. Occasionally, it draws a very high number. Yeah, I know, this is statistics, things like that does happen, but this really strange because a poisson distribution with a parameter of 3 shouldn't see the number 1932 very often, but the same,

Full_Name: Kjetil Kjernsmo Version: Version 1.0.0 OS: osf1 Submission from: (NULL) (129.240.28.227) Hello again! This is a follow-up to my message on R-help about a problem with match.arg() I have a little more on the topic, but not much really. First, it was not entirely true what I wrote in r-help that I call match.arg() in my ramp(), because ramp() just calls qamp() with runif() as

Full_Name: Kjetil Kjernsmo Version: 0.65.1 OS: Digital UNIX 4.0 Submission from: (NULL) (129.240.28.172) I have the following R code: ncloudsbin <- function(binno, ntotalclouds, numberofbins = 100, linewidth = numberofbins / 6) return(ntotalclouds * dnorm(binno, sd = linewidth)) photonsincidenttodetectorbin <- function(binno, ntotalclouds, avgphotonfromcloud, ampmode = 1, ampmean = 2,

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

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,

Dear all, I have done further modifications on the histogram functions that I reported earlier this year, and I hope this can be of use and perhaps included in the distribution. I have been using this stuff a couple of months myself, and while it is nothing sophisticated, it has it's applications. :-) I did a few small modifications today to make it a bit more compact. I have modified the

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=\|%:\

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

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

Dear Sir / Madam, I am new for R coding. Kindly help me out in sorting out the following problem. There are 50 rows with six coloumns(you could see in the attached .txt file). I wish to go for filtering this 50 rows for any one of the six coloumns satisfying the value >= 64. I need to have a final table with rows having >= 64 value in any one of the six coloumns and the rest could be

Dear all, I am just wondering if anybody has implemented a function that can give a piecewise linear interpolation in more than 2 dimensions? I have looked at the akima package, but I would rather like a piecewise linear interpolation rather than a spline and while it did the job quite satisfactory for 2 dimensions, I need to interpolate in at least three dimensions. If anybody has implemented

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

Why is it that when a vector contains no NAs, na.omit returns an empty vector rather than the object/vector that was passed to it? >na.omit(c(1,2,3,4)) numeric(0) in contrast ... >na.omit(c(1,2,3,NA)) [1] 1 2 3 attr(,"na.action") [1] 4 attr(,"na.action")attr(,"class") [1] "omit" Thanks, Alex --------------------------- platform

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

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

As of R 1.7.0, princomp no longer accept matrices with more coloumns than rows. I'm curious: Why was this decision made? I work a lot with data where more coloumns than rows is more of a rule than an exception (for instance spectroscopic data). To me, princomp have two advantages above prcomp: 1) It has a predict method, and 2) it has a biplot method. A biplot method shouldn't be too