rstream: Streams of Random Numbers for Stochastic Simulation

The package rstream provides a unified interface to streams of random numbers for the R statistical computing language. Features are: * independent streams of random numbers * substreams * easy handling of streams (initialize, reset) * antithetic random variates The paper describes this packages and demonstrates an simple example the usefulness of this approach.Series: Preprint Series / Department of Applied Statistics and Data Processin

Computing the Two-Sided Kolmogorov-Smirnov Distribution

We propose an algorithm to compute the cumulative distribution function of the two-sided Kolmogorov-Smirnov test statistic D_n and its complementary distribution in a fast and reliable way. Different approximations are used in different regions of n, x. Java and C programs are available.

Random number generation with multiple streams for sequential and parallel computing

International audienceWe provide a review of the state of the art on the design and implementation of random number generators (RNGs) for simulation, on both sequential and parallel computing environments. We focus on the need for multiple streams and substreams of random numbers, explain how they can be constructed and managed, review software libraries that offer them, and illustrate their usefulness via examples. We also review the basic quality criteria for good random number generators and their theoretical and empirical testing

Computing the Two-Sided Kolmogorov-Smirnov Distribution

We propose an algorithm to compute the cumulative distribution function of the two-sided Kolmogorov-Smirnov test statistic Dn and its complementary distribution in a fast and reliable way. Different approximations are used in different regions of n, x. Java and C programs are available

Estimating the Probability of a Rare Event Over a Finite Time Horizon

We study an approximation for the zero-variance change of measure to estimate the probability of a rare event in a continuous-time Markov chain. The rare event occurs when the chain reaches a given set of states before some fixed time limit. The jump rates of the chain are expressed as functions of a rarity parameter in a way that the probability of the rare event goes to zero when the rarity parameter goes to zero, and the behavior of our estimators is studied in this asymptotic regime. After giving a general expression for the zero-variance change of measure in this situation, we develop an approximation of it via a power series and show that this approximation provides a bounded relative error when the rarity parameter goes to zero. We illustrate the performance of our approximation on small numerical examples of highly reliableMarkovian systems. We compare it to a previously proposed heuristic that combines forcing with balanced failure biaising. We also exhibit the exact zero-variance change of measure for these examples and compare it with these two approximations

Randomized Quasi-Monte Carlo: An Introduction for Practitioners

International audienceWe survey basic ideas and results on randomized quasi-Monte Carlo (RQMC) methods, discuss their practical aspects, and give numerical illustrations. RQMC can improve accuracy compared with standard Monte Carlo (MC) when estimating an integral interpreted as a mathematical expectation. RQMC estimators are unbiased and their variance converges at a faster rate (under certain conditions) than MC estimators, as a function of the sample size. Variants of RQMC also work for the simulation of Markov chains, for function approximation and optimization, for solving partial differential equations, etc. In this introductory survey, we look at how RQMC point sets and sequences are constructed, how we measure their uniformity, why they can work for high-dimensional integrals, and how can they work when simulating Markov chains over a large number of steps