search for: cloudintensity

...amp(runif(n), type)) 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. More details: I'm calling this function a number of times through a function intensityfromcloud <- function(ncloudsbin, cloudintensity, amptype, exptime, ...) { 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 writin...

..., the code isn't very involved, really, so if you want to try, here it is: ncloudsbin <- function(binno, ntotalclouds, numberofbins = 100, linewidth = numberofbins / 6) return(ntotalclouds * dnorm(binno, sd = linewidth)) intensityfromcloud <- function(ncloudsbin, cloudintensity, exptime, amptype) return(cloudintensity * exptime * ramp(ncloudsbin, amptype)) photonsdetectedbin <- function(ncloudsbin, intensityfromcloud, continuumintensity, exptime, amptype, sensitivity) return(sum(rpois(ncloudsbin, sensit...