similar to: Garbage Collecting

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

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

Dear all, I'm a bit concerned about underflow problems. Such problems typically occur when two almost equal numbers are subtracted. I was wondering if R has any mechanisms for warning users about potential problems, or if any policy has been developed? To check, I just did: > 1.0000000000000001 - 1 [1] 0 > 1.000000000000001 - 1 [1] 1.110223e-15 > 1.00000000000001 - 1 [1]

Dear all, A looong time ago, Witold Eryk Wolski asked here why there wasn't a log="xy" parameter to the hist() function <URL:http://www.R-project.org/nocvs/mail/r-help/2001/0267.html>, and Prof. Ripley responded that a loglog histogram does not make much sense, and that one should use a better density estimate if one seeks to plot log density. I understand the point and I

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

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

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 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''m working with the geometric distribution for the time being, and I''m confused. This may have more to do with statistics than R itself, but since I''m getting results from R I find counterintuitive (well, yeah, my statistical intuition has not been properly sharpened), I feel like asking. The point first: If I do > rgeom(1,prob=1) I get: [1] NaN Warning

Dear all, After looking at several packages I decided I would use R for my thesis work (we've even got a site license for S-plus), so I'm learning. I have run into a problem, and nobody on the house has a solution, so I'm wondering if anybody out there can help. I've got two functions of four variables, one outputs a scalar, the other an array, and I want to be able to give an

Hello! I have two matrices of equal dimension ll and lu that I want to do a ks.test on corresponding rows in these matrices (dim(ll) is [1] 101 100). If I do e.g. > ks.test(ll[50,],lu[50,]) just for testing, it displays a lot of numbers, and some more info: [116] -3.000000e-02 -4.000000e-02 -4.000000e-02 -3.000000e-02 -4.000000e-02 [121] -3.000000e-02 -4.000000e-02 -3.000000e-02 -2.000000e-02

Dear all, I'm making a system to create a star-chart, and I have encountered one of the intuitive and clearly thoughtful conventions that are so common in astronomy :-J... Well, the point being, I have to reverse the x-axis. I figured, I might just set xlim=rev(range(x)). Do e.g. > x <- 0:5 > y <- (0:5)^3 > plot(x, y) > plot(x, y, xlim=rev(range(x))) The problem is that here

Full_Name: Kjetil Kjernsmo Version: 1.0.0 OS: osf1 Submission from: (NULL) (129.240.28.227) Hello! I'm afraid I have another obscure bug. R has just crashed with a segmentation fault. This time, I have even failed to reproduce it, but fortunately, I have a core dump. It's at <URL:http://www.astro.uio.no/~kjetikj/tmp/core> The call I did was log10(quantile(chisqtests5M11001,

Dear all, I started a large run this afternoon, and when I came back after dinner, three of five processes had stopped with the error message: Error in matrix(nullprofile, length(nullprofile), nprofiles) : Unimplemented feature in copyVector I have seen anything about this anywhere (and I don't know where to look...). As you might guess, I'm trying to make a matrix where the