[R] strange behaviour of negative binomial

2005-01-05 Thread dax42
Dear list, I ran into a strange behaviour of the pnbinom function - or maybe I just made a stupid mistake. First thing is that pnbinom seems to be very slow. The other - more interesting one - is that I get two different curves when I plot the estimated density and the density given by pnbinom.

Re: [R] strange behaviour of negative binomial

2005-01-05 Thread Thomas Lumley
On Wed, 5 Jan 2005, dax42 wrote: Dear list, I ran into a strange behaviour of the pnbinom function - or maybe I just made a stupid mistake. First thing is that pnbinom seems to be very slow. The other - more interesting one - is that I get two different curves when I plot the estimated density

Re: [R] strange behaviour of negative binomial

2005-01-05 Thread Peter Dalgaard
dax42 [EMAIL PROTECTED] writes: Dear list, I ran into a strange behaviour of the pnbinom function - or maybe I just made a stupid mistake. First thing is that pnbinom seems to be very slow. The other - more interesting one - is that I get two different curves when I plot the estimated

Re: [R] strange behaviour of negative binomial

2005-01-05 Thread Gabor Grothendieck
dax42 Dax42 at web.de writes: : : Dear list, : : I ran into a strange behaviour of the pnbinom function - or maybe I : just made a stupid mistake. : First thing is that pnbinom seems to be very slow. The other - more : interesting one - is that I get two different curves when I plot the :

Re: [R] strange behaviour of negative binomial

2005-01-05 Thread dax42
Hi, I am sorry for the incorrect terms. Of course I mean the cumulative distribution functions. As I just got two answers telling me my problem does not exist, I tried again on a different setup - and indeed it works. So, my first run was done on a Mac OS X 10.3.7 running R 2.0.0, GUI version.

Re: [R] strange behaviour of negative binomial

2005-01-05 Thread Gabor Grothendieck
I tried it on Windows XP under both 2.0.1 and 2.1.0 and it seemed to work both times. I issued the command: set.seed(1) before both runs to ensure that I was actually using the same random numbers. I suggest you try that in all your runs too just in case its a problem with some sets of