[Rd] Random number generation problem (PR#554)
Tue, 23 May 2000 20:33:30 +0200 (MET DST)
Full_Name: Kjetil Kjernsmo
Submission from: (NULL) (220.127.116.11)
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"))
if(type == "point")
return(1 / (2 * sqrt(1 - p)))
ramp <- function(n, type=c("point", "nolens"))
if(type == "point")
else return(rep(1, n))
...to generate random numbers. I think there might be something wrong
with runif, possibly the same problem as rpois had.
I'm calling this function a number of times through a function
intensityfromcloud <- function(ncloudsbin, cloudintensity, amptype,
amps <- ramp(ncloudsbin, amptype)
write(c(max(amps), sum(amps), sum(amps - 1)),
"lastamps.dat", ncolumns=3, append=T)
return(cloudintensity * exptime * amps)
...that is called a lot. As you can see, I'm writing the maximum
(among other things) of the generated numbers to file
(this is a preliminary hack, for a (quite unrelated) sanity
check I'm running).
Now, it turns out that in 8080 maxima, the number 14654.29 appears
twice. Things like that happen, this _is_ statistics,
and my pdf doesn't fall that fast, but you know...
I thought I'd report this sooner the better...
I made a histogram of the logarithm of the maxima for the 8080 points
and put it at <URL:http://www.astro.uio.no/~kjetikj/tmp/ampsmax.ps>
BTW, update on my other reports:
The problem I had with match.arg() (PR#536) was in the same piece of code,
and unfortunately, it seems to happen with 1.0.1 as well.
About my crashes, PR#545 and PR#551, I had another crash today, and I
have a 250MB (Warning! :-)) core dump at
I was doing a summary of a rather strange data.frame, I was just curious
as to what summary would say to such a thing. I don't know if it would
be of any use. Uwe suggested that it might be a hardware problem, and I
have contacted my sysadm, and a check will hopefully be performed on Friday.
If it would be useful to see my actions prior to the crash, please let me
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: email@example.com