similar to: rpois hangs (PR#530)

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

Full_Name: Kjetil Kjernsmo Version: 0.65.1 OS: Digital UNIX 4.0 Submission from: (NULL) (129.240.28.172) I'm experiencing problems with rpois. Occasionally, it draws a very high number. Yeah, I know, this is statistics, things like that does happen, but this really strange because a poisson distribution with a parameter of 3 shouldn't see the number 1932 very often, but the same,

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

The R poisson random generator rpois appears to have a bug for theta 10 or larger. The sample mean of the pseudo variates is too small: sample mean approx theta - 0.5. I use Version 1.1.1 (August 15, 2000) Of R on a Dell OptiPlex computer with the Windows XP Professional operating system. Has this bug been fixed in later versions? (I found another reported rpois bug, but it appears to be

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

Hello, All: ????? Consider: Browse[2]> set.seed(1) Browse[2]> rpois(9, 1e10) NAs produced[1] NA NA NA NA NA NA NA NA NA ????? Should this happen? ????? I think that for, say, lambda>1e6, rpois should return rnorm(., lambda, sqrt(lambda)). ????? For my particular Monte Carlo, I have replaced my call to rpois with a call to the following: ?rpois. <- function(n,

Hi all, It's been a while since i've used R and I can't remember how to do the following: i have a = rpois (100000, x) b = rpois (100000, y) what is the code to show that a>b, b>a and a=b to show just the number of occurances? at the moment when I type a>b I get a nice long list of true or false. so i'm hoping to have the following a>b > 35000 b>a > 25000

On 2020-01-19 09:34, Benjamin Tyner wrote: >> ------------------------------------------------------------------------ >> Hello, All: >> >> >> ? ????? Consider: >> >> >> Browse[2]> set.seed(1) >> Browse[2]> rpois(9, 1e10) >> NAs produced[1] NA NA NA NA NA NA NA NA NA >> >> >> ? ????? Should this happen? >>

On 1/20/20 12:33 PM, Martin Maechler wrote: > > It's really something that should be discussed (possibly not > here, .. but then I've started it here ...). > > The NEWS for R 3.0.0 contain (in NEW FEATURES) : > > * Functions rbinom(), rgeom(), rhyper(), rpois(), rnbinom(), > rsignrank() and rwilcox() now return integer (not double) > vectors.

So imagine rpois is changed, such that the storage mode of its return value is sometimes integer and sometimes numeric. Then imagine the case where lambda is itself a realization of a random variable. Do we really want the storage mode to inherit that randomness? On 1/19/20 10:47 AM, Avraham Adler wrote: > Maybe there should be code for 64 bit R to use long long or the like? > > On

Dear all, Uwe e-mailed me yesterday about the problems some time ago. I had hoped someone would have a look at our computers, but due to some unfortunate circumstances, it hasn't happened. However, it seems unlikely that it is connected to a specific piece of hardware. It could be an architecture-dependent problem though. I also had a word with those who did the installation of R, and they

????? This issue arose for me in simulations to estimate confidence, prediction, and tolerance intervals from glm(., family=poisson) fits embedded in a BMA::bic.glm fit using a simulate.bic.glm function I added to the development version of Ecfun, available at "https://github.com/sbgraves237/Ecfun".? This is part of a vignette I'm developing, available at

On 2020-01-19 13:01, Avraham Adler wrote: > Crazy thought, but being that a sum of Poissons is Poisson in the sum, > can you break your ?big? simulation into the sum of a few smaller > ones? Or is the order of magnitude difference just too great? ????? I don't perceive that as feasible.? Once I found what was generating NAs, it was easy to code a function to return pseudo-random

Thanks to Luke and Avi for their comments.? I wrapped "round" around the call to "rnorm" inside my "rpois.".? For "lambda" really big, that "round" won't do anything.? However, it appears to give integers in floating point representation that are larger than .Machine$integer.max.? That sounds very much like what someone would want.?

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

There is a problem with `rpois'. It does seem to take care about the order of the arguments. This is an example: > rpois(n=1,lambda=2) [1] 3 > rpois(lambda=2,n=1) [1] 2 0 It obviously uses the first argument as the number of samples to be drawn, which is wrong. I used Version 0.49 Beta (April 23, 1997). Fredrik

On my Mac: str(.Machine) ... $ integer.max????????? : int 2147483647 ?$ sizeof.long????????? : int 8 ?$ sizeof.longlong????? : int 8 ?$ sizeof.longdouble??? : int 16 ?$ sizeof.pointer?????? : int 8 ????? On a Windows 10 machine I have, $ sizeof.long : int 4; otherwise the same as on my Mac. ????? Am I correct that $ sizeof.long = 4 means 4 bytes = 32 bits? log2(.Machine$integer.max)

Kjetil Kjernsmo <kjetil.kjernsmo@astro.uio.no> writes: > On 8 Nov 2000, Peter Dalgaard BSA wrote: > > >I'm not at all happy with this: > > > >Solaris : > >> range(sapply(1:2000, function(n) mean(rpois(10000, 15.0)))) > >[1] 15.0524 15.3403 > > Hm, OK, so it isn't just me.... I guess it is time to file a bug report, > should I do it,

Hi there, I need some helpo with one simulation. I have a matrix of 4x21 I want ro generate 100 random Chi2 distribution but inly from the first row of my matrix. This is what I try: Option (a): >lambda=1:100 >f=function(x)chisq(4,0,x) >lapply(lambda,f) It doesn't work! option(b): >m=mean(matrix[1,]) >x<-rpois(100,lambda=m) Works, but only generates a row of 25 values.

On 1/20/20 4:26 AM, Martin Maechler wrote: > Coming late here -- after enjoying a proper weekend ;-) -- > I have been agreeing (with Spencer, IIUC) on this for a long > time (~ 3 yrs, or more?), namely that I've come to see it as a > "design bug" that rpois() {and similar} must return return typeof() "integer". > > More strongly, I'm actually pretty