similar to: (NA || x) and (NA && y) {another "bug" } (PR#749)

This came up on r-help but indicates a bug. dnbinom(x,n,p) calls dbinom_raw(n-1,...) which returns 0 for n<1. -thomas ---------- Forwarded message ---------- Date: Thu, 08 Feb 2001 17:10:23 +0000 From: Yudi Pawitan <yudi@stat.ucc.ie> To: Mark Myatt <mark@myatt.demon.co.uk> Cc: R-Help <r-help@stat.math.ethz.ch> Subject: Re: [R] Goodness of fit to Poisson / NegBinomial

In several places, we have the convention that col = NA means `` don't draw '', e.g. for points() or rect(). For col = <vector> this only works when the vector is integer (numeric), but not in other cases. I think this is a bug, even though not in the strict sense, since we don't seem to claim it should work.. Example: plot(1:10, col = c(1:4,NA,NA,4:1)) ## fine

Full_Name: Kjetil Kjernsmo Version: Version 1.0.0 OS: osf1 Submission from: (NULL) (129.240.28.227) Dear all, I have come across another really strange bug in rpois. If I do > rpois(1, 1000000 * rpois(1, 10000)) sometimes it hangs, and hangs... On my XP1000 alphaev6 workstation, it seems to happen more frequently than on other machines. I started a session on an older computer, and from

I don't know if this is a bug report or a feature request, but a colleague complained that source(,echo=T) doesn't echo the comments from the source file. Thomas Lumley Assistant Professor, Biostatistics University of Washington, Seattle -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

Short example: > str(d <- data.frame(a = I(letters[1:5]))) `data.frame': 5 obs. of 1 variable: $ a:Class 'AsIs' chr [1:5] "a" "b" "c" "d" ... > str(d[TRUE,,drop = FALSE]) `data.frame': 5 obs. of 1 variable: $ a: chr "a" "b" "c" "d" ... The real problem about this is, that as soon

>>>>> "MM" == Martin Maechler <maechler@stat.math.ethz.ch> writes: >>>>> "TL" == Thomas Lumley <thomas@biostat.washington.edu> writes: TL> On Thu, 24 Aug 2000, Troels Ring wrote: >>> Dear friends. >>> >>> Is this as expected ? Is alpha and beta too large simply ? >>>

>>>>> "MM" == Martin Maechler <maechler@stat.math.ethz.ch> writes: >>>>> "TL" == Thomas Lumley <thomas@biostat.washington.edu> writes: TL> On Thu, 24 Aug 2000, Troels Ring wrote: >>> Dear friends. >>> >>> Is this as expected ? Is alpha and beta too large simply ? >>>

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

A colleague has a dozen of excel sheets and also expects to get updates regularly. He could open these in M$-Excel and export as *.csv manually, "bring back to Unix" and then read into R. Of course there must be options to start programming this in something like visual basic, but we wouldn't to really want to... We also know that probably Gnumeric could do the job (since the xls

Some of you may have seen the pythag() part in the R API definition in "Writing R Extensions" (source = doc/manual/R-exts.texinfo). or followed the report and Prof. Brian Ripley's answer about pythag()'s availability from R's binary. As we say in above manual >> `pythag(A, B)' computes `sqrt(A^2 + B^2)' without overflow or >> destructive

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

>>>>> "TL" == Thomas Lumley <tlumley@u.washington.edu> writes: on R-help and the Subject "Log likelihood from glm" TL> On Mon, 9 Jul 2001, Phil wrote: >> >> Is there any way of getting log likelihood values for each iteration from >> poisson glm model? >> TL> Yes, the trace=TRUE option. TL>

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

>>>>> "TL" == Thomas Lumley <thomas at biostat.washington.edu> writes: TL> I'm receiving duplicates of a number of old r-help messages, all with TL> headers having the rly-mx2.maxis.net.my section in the message below. Is TL> this at your end? Do you know about it? Doug Bates has averted me; and I have since stopped these to be propagated