Re: [Rd] bug in modulus operator %% (PR#7852)

2005-05-11 Thread p . dalgaard
"McGehee, Robert" <[EMAIL PROTECTED]> writes: > Yes, but from ?"%%": > "It is guaranteed that 'x == (x %% y) + y * (x %/% y)' (up to rounding > error) ..." > > (R 2.1.0) > > x <- 1 > > y <- 0.2 > > x %% y > [1] 0.2 > > (x %% y) + y * (x %/% y) > [1] 1.2 > > Certainly 1 does not equal 1.2 as the

Re: [Rd] RAW types not restored from Rda files (PR#7812)

2005-04-25 Thread p . dalgaard
Byron Ellis <[EMAIL PROTECTED]> writes: > There's an endian difference, but other than that our results are > identical: > > 000 524458320a580a0002000202 > 020 010400040200010000100900 > 040 01

Re: [Rd] make check-all fails (PR#7784)

2005-04-10 Thread p . dalgaard
Peter Dalgaard <[EMAIL PROTECTED]> writes: > Aha! 100 times machine precision in not all that much when the numbers > themselves are in double digits. In fact, one is over 100. The case > that triggers the failure is #149 > > > 147 148 149 151 152

Re: [Rd] make check-all fails (PR#7784)

2005-04-09 Thread p . dalgaard
"M. Edward (Ed) Borasky" <[EMAIL PROTECTED]> writes: [snip] > >Hmm, could you replace the a1 == a2 with all.equal(a1, a2) instead? > >(inside reg-tests-1.R of course) > > > >Asking for identity up to machine precision does look a bit optimistic... > > > > > That worked ... it got through reg-tests

Re: [Rd] sys.on.exit not working (PR#7665)

2005-02-07 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > It is working as documented: there is no on.exit set for str, is there? > > 'sys.on.exit()' retrieves the expression stored for use by > 'on.exit' in the function currently being evaluated. (Note that > ^^

[Rd] Non-invisible returns from title and mtext (PR#7397)

2004-11-29 Thread p . dalgaard
A small cosmetic item: In current versions we have > plot(0) > title("Foobar") NULL in which the printing of NULL didn't happen in 1.9.1. Same story with mtext(), and both pretty apparent in the output from demo(graphics). I can't offhand see what caused the change. -- O__ Peter Dal

Re: [Rd] ecdf with lots of ties is inefficient (PR#7292)

2004-10-17 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > vals <- sort(unique(x)) > y <- tabulate(match(x, vals)) > rval <- approxfun(vals, cumsum(y)/n, method = "constant", yleft = 0, > yright = 1, f = 0, ties = "ordered") > > should work better for you and may be little

Re: [Rd] typo in R online manual (PR#7200)

2004-08-25 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > That is a mistake at your own department's site, but not in the current R > sources. It seems those docs are from R 1.2.2! Even R 1.6.2 has the > correct formula. > > Please ask whoever maintains that site to correct it, and don't expect the > R

Re: [Rd] Similar error as (PR#6958)

2004-06-24 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > > BTW, we just had two weeks of alpha and beta releases in which to > > find this sort of stuff, you know > > Ye...s! [I even might have found it myself, but I set LANG=C, LC_CTYPE=da_DK in my shells because I otherwise go nuts from the semi-tr

[Rd] Re: [R] Printing Lattice Graphs from Windows (PR#6950)

2004-06-06 Thread p . dalgaard
Dirk Eddelbuettel <[EMAIL PROTECTED]> writes: > On Sun, Jun 06, 2004 at 01:32:00PM -0700, Spencer Graves wrote: > > I agree: It sounds like a bug, as you said, Irk, in that > > That brilliant :) Many people, in Germany as well as abroad, managed to > chop Eddelbuettel quite well. Turning

Re: [Rd] error in make.names docs (PR#6912)

2004-05-24 Thread p . dalgaard
Peter Dalgaard <[EMAIL PROTECTED]> writes: > Prof Brian Ripley <[EMAIL PROTECTED]> writes: > > > Example output in comment is now corrected in R-patched. > > What has the rest to do with a bug report? > > He's hoping (with some reason) that the change of behaviour was the > bug. I seem to recall

Re: [Rd] error in make.names docs (PR#6912)

2004-05-24 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > Example output in comment is now corrected in R-patched. > What has the rest to do with a bug report? He's hoping (with some reason) that the change of behaviour was the bug. I seem to recall some internal discussion of the effect, but not whether i

Re: [Rd] optim-Bug (PR#6720)

2004-04-21 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > On 21 Apr 2004, Peter Dalgaard wrote: > ... > > I have this down to the use of R_alloc in the vect() function (line > > 40, optim.c). Replacing with S_alloc (which zeros memory) removes the > > issue for me on this machine and on the Opteron. Could

[Rd] Re: [R] Unexpected behaviour of identical (PR#6799)

2004-04-20 Thread p . dalgaard
"Swinton, Jonathan" <[EMAIL PROTECTED]> writes: > # works as expected > > ac <- c('A','B'); > > identical(ac,ac[1:2]) > [1] TRUE > > #but > > af <- factor(ac) > > identical(af,af[1:2]) > [1] FALSE > > Any opinions? Did a cross-check with Splus and it doesn't do that , so I think it qualifies

Re: [Rd] Incorrect handling of NA's in cor() (PR#6750)

2004-04-09 Thread p . dalgaard
Marek Ancukiewicz <[EMAIL PROTECTED]> writes: > Dear Thomas, > > The question becomes: how do we rank missing values? In > version 1.8.1 at least, cor () uses default handling of > missing values by rank() [by na.last parameter], that is > missing values are assigned the highest rank. However, i

Re: [Rd] Bug in points.formula (PR#6652)

2004-03-08 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > What I think we need is the number of rows of the model frame *before > subsetting*, so that any further arguments (and there are none such in > Martin's example) of the same length as the response get subsetted. > > Berwin's fix is ugly, and I thin

[Rd] Spelling (PR#6570)

2004-02-10 Thread p . dalgaard
I came across this in connection with an unrelated issue > beta[2] Error in beta[2] : object is not subsettable > beta[2] <- 5 Error in "[<-"(`*tmp*`, 2, value = 5) : object is not subsetable One of the messages must be wrong, but I need a native English speaker to tell me which one. -- O__

[Rd] Testing change to web interface (PR#6299)

2004-01-02 Thread p . dalgaard
Full_Name: Version: OS: Submission from: (NULL) (192.38.19.4) I'm playing with some modifications to work around the problem where people would inadvertently submit a report a second time by refreshing the main page in their browser. __ [EMAIL PROTE

[Rd] Bug repository problem (PR#6168)

2003-12-28 Thread p . dalgaard
Full_Name: Peter Dalgaard Version: -irrelevant- OS: -Linux/Mozilla 1.4- Submission from: (NULL) (80.199.16.96) It has been reported that after sending a report and returning to the main page, refreshing the main page (obviously a natural thing to want to do) resends the report as a duplicate. Thi

[Rd] Bug repository problem (PR#6167)

2003-12-28 Thread p . dalgaard
Full_Name: Peter Dalgaard Version: -irrelevant- OS: -Linux/Mozilla 1.4- Submission from: (NULL) (80.199.16.96) It has been reported that after sending a report and returning to the main page, refreshing the main page (obviously a natural thing to want to do) resends the report as a duplicate. Thi

[Rd] Segfault in foreign (PR#5049)

2003-11-12 Thread p . dalgaard
> library(foreign) > lookup.xport(T) Program received signal SIGSEGV, Segmentation fault. 0x4207a547 in strchr () from /lib/i686/libc.so.6 (1.8.0 on Linux) This actually did bite one of our users. Not that he expected to be able to pass a logical, but he accidentally defined all <-function(

[Rd] Server down and (hopefully) back up. (PR#4974)

2003-11-07 Thread p . dalgaard
The R-bugs service has been down for a while since a quick decision was taken to move the mailserver at pubhealth.ku.dk to a different machine yesterday. (The move was planned, just not at that time, and the old server had gottten into a disk full condition.) Since that implied a switch of machine

Re: [Rd] Something strange in cor.test in R-1.8.0 (PR#4718)

2003-10-23 Thread p . dalgaard
Kurt Hornik <[EMAIL PROTECTED]> writes: > > maechler writes: > > > "Duncan" == Duncan Murdoch <[EMAIL PROTECTED]> > > on Wed, 22 Oct 2003 08:17:32 -0400 writes: > > Duncan> I've now tried the code > > > set.seed(1) > > x <- rnorm(50) > > y <- rnorm(50) > > cor.test(x,y,met

[Rd] Re: [R] File Reading Problem (PR#4044)

2003-09-01 Thread p . dalgaard
Duncan Murdoch <[EMAIL PROTECTED]> writes: > On Mon, 1 Sep 2003 14:26:43 -0700, "Jiming Yu" > <[EMAIL PROTECTED]> wrote: > > >Dear all, > >I am trying to read characters byte by byte(in their ASCII codes) from a > >file > > I was going to suggest using readBin, but there seems to be a bug: >

Re: [Rd] Update on "Shift-Return causes segfault (PR#3293)"

2003-06-27 Thread p . dalgaard
Graeme Ambler <[EMAIL PROTECTED]> writes: > > What version of readline is linked to R 1.6.0 on that system? > Probably version 4.1 --- I remember in the past always having to install the > readline41-4.1-xx package to get R to work. Perhaps Martyn should go back to > linking against that in the

Re: [Rd] Shift-Return causes segfault (PR#3293)

2003-06-19 Thread p . dalgaard
Marc Schwartz <[EMAIL PROTECTED]> writes: > On Thu, 2003-06-19 at 05:58, [EMAIL PROTECTED] wrote: > > Full_Name: Graeme Ambler > > Version: 1.7.1 > > OS: Red Hat Linux 9 > > Submission from: (NULL) (137.222.80.161) > > > > > > Typing Shift-Return in R causes a segmentation fault. I am using the

Re: [Rd] Re: [R] Postscript query: plotting long vectors (PR#3132)

2003-05-30 Thread p . dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > Except that there is nothing reproducible in that report, not even the > claimed (by Don MacQueen) incorrect lines! Er, right, sorry. Here's a version of the effect, for inclusiont with the report: postscript() x <- seq(0,2*pi,,5) plot(x,sin(x

[Rd] Re: [R] Postscript query: plotting long vectors (PR#3132)

2003-05-30 Thread p . dalgaard
Don MacQueen <[EMAIL PROTECTED]> writes: > When I run the example in R 1.6.2, and view it with gs, I get a good plot. > When I run the example in R 1.7.0, and view it with gs, I get a bad plot. > (run on the same host) > > My "bad plot" is as described by Stephen. ... > (followed by ~20 lines

[Rd] unique turns ordered into factor (PR#2591)

2003-02-28 Thread p . dalgaard
[reported originally by Christoffer Tornøe] > f <- ordered(month.name) > f [1] January February March April May June July [8] AugustSeptember October November December 12 Levels: April < August < December < February < January < July < ... < September > unique(f) [1]