similar to: Reading single precision floats from binary file

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

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

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

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 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 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 have been away for a month, but some time ago, I noted something that was very surprising to me. I think this example illustrates it well: > test <- data.frame(col1 = c(1,2), col2=c(3, 4)) > test col1 col2 1 1 3 2 2 4 > apply(test, 1, function(x) browser()) Called from: FUN(newX[, i], ...) Browse[1]> x$col1 NULL Browse[1]> x[["col1"]] [1] 1

Dear all, I'm generally not a big fan of GUIs, but I have just discovered one thing where a GUI would be a nice thing: creating color palettes. I'm creating an "image", for transparencies for a lecture I'm giving, so this is going to look great... :-) What I'm out for, is to make green tones and red tones, and a relatively sharp break between them. This break will

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

On 8 Nov 2000, Peter Dalgaard BSA wrote: >Done... Great! >(This kind of behaviour has been observed on various platforms and >various compilers, sometimes getting stuck around 15.15 and sometimes >around 14.5. Others report no problems. Very strange. I've seen one >case where it apparently shifted from one regime to the other.) Yes, it does that here too.... The following is

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

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

Hello! First, I must also congratulate the R core team with their accomplishment! I have gotten to like R a lot, and I have recommended it for inclusion in an "Astronomy for Linux"-distribution which is in use by many professional astronomers and observatories. I'm currently working on importing data from files created by a C program (that I have not written myself, I have the

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

Full_Name: Kjetil Kjernsmo Version: 1.0.1 OS: osf1 Submission from: (NULL) (129.240.28.227) Hello! I'm getting strange numbers from a random number generation function again. I was wondering if it could be the same problem as I reported in (PR#439). I'm using the functions (some comments deleted for brevity): qamp <- function(p, type=c("point", "nolens")) {

Dear all, I have some difficulty getting R-1.3.0 to run on the alpha, with osf4.0e (Tru64, or whatever they call it... :-) ). configure reports the following configuration: R is now configured for alphaev6-dec-osf4.0e Source directory: . Installation directory: /astro/local C compiler: gcc -mieee -g -O2 C++ compiler: c++ -g -O2 FORTRAN

Dear all, I have some difficulty getting R-1.3.0 to run on the alpha, with osf4.0e (Tru64, or whatever they call it... :-) ). configure reports the following configuration: R is now configured for alphaev6-dec-osf4.0e Source directory: . Installation directory: /astro/local C compiler: gcc -mieee -g -O2 C++ compiler: c++ -g -O2 FORTRAN

Dear all, I have an image that I need to filter, and so I'm looking for a method to convolve it with a matrix. If I understood the docs for convolve correctly, it only works in 1D (and I have tried to convolve, it didn't look good). So, I wondered if anybody have implemented 2D convolution in R, or have any good advices to share (beyond having a look at mvfft), before I go hacking?