On the autocorrelation structure of inversive congruential pseudorandom number sequences

The inversive congruential method for generating uniform pseudorandom numbers is a particularly attractive alternative to linear congruential generators with their well known deficiencies. In the present paper inversive congruential generators with power of two modulus are considered. Known favourable results on the period length and the statistical independence properties of the generated sequences are summarized. Besides that the autocorrelation structure of inversive congruential sequences is studied. The main result reveals a remarkable autocorrelation property of the inversive congruential method.     

Normal random numbers lying on spirals and clubs

Two methods for transforming uniformly distributed random numbers into normally distributed random numbers are considered in conjunction with linear congruential generators. The two-dimensional lattice structure of the uniform random numbers is transformed by the Box-Muller method into a spiral structure and by the polar method into a club-shaped structure. The approximation of the two-dimensional normal distribution and the independence of the associated random variables are discussed.     

An exhaustive analysis of prime modulus multiplicative congruential random number generators with modulus smaller than 215

This paper presents the results of an exhaustive analysis of all of the prime modulus multiplicative congruential random number (RN) generators with moduli smaller than 2¹⁵. In an amount of around 20 million multipliers which are able to produce a full period of RNs, 239 multipliers have a good lattice structure. Among which 52 multipliers further pass a comprehensive battery of empirical tests. These 52 multipliers thus possess good local and global statistical properties. It is worthwhile to note that some empirically tested multipliers recommended in some previous studies are not on this list. The conclusion is that both theoretical and empirical tests are mandatory to sieve out good multipliers. To generate RNs of very long period, many existing techniques can be applied without further effort.     

Long-range correlation analysis of the Wichmann-Hill random number generator

The distribution of points (r _n,r _n+s), n = 0, 1, 2,...whose coordinates are terms at distance s of the pseudorandom sequence generated by the Wichmann and Hill method is studied. It is known that for many congruential generators critical values of the distance s exist such that these points, far from being uniformly distributed, are concentrated on very few lines. An algorithm is described for computing the critical distances within the Wichmann-Hill sequence and the results obtained are compared with those of other linear congruential generators.     

Generalized Linear Models

The importance of random number generators has increased over the years. This follows from the fact that contemporary research methods rely more and more on simulation and the increased importance of encryption technology. The output of a random number generator is created by either an algorithm or a physical device. The most popular method for random number generation is through the use of an algorithm. This article presents a new category of physical random bit generator that is packaged by several manufacturers. A statistical analysis of the output from the generators is given.     

Regularities and irregularities in a random 0,1 sequence

A few properties of the random 0,1 sequences are studied. These properties can be utilized to test the randomness of a 0,1 sequence produced by a computer (random generator). This test seems to be appropriate for discovering the periodicity (or the length of the period) of a pseudorandom sequence. A possible application of the presented results in regression analysis is also discussed.     

On the Deng-Lin random number generators and related methods

We study the structure and point out weaknesses of recently proposed random number generators based on special types of linear recurrences with small coefficients, which allow fast implementations. Our theoretical analysis is complemented by the results of simple empirical statistical tests that the generators fail decisively. Directions for improvement and alternative generators are also pointed out.     

Modified explicit inversive congruential pseudorandom numbers with power of 2 modulus

A modified version of the explicit inversive congruential method with power of 2 modulus for generating uniform pseudorandom numbers is introduced. The statistical independence behaviour of the generated sequences is studied based on the distribution of all pairs of successive pseudorandom numbers over the entire period. Lower and upper bounds for the discrepancy of the corresponding two-dimensional point sets are established. These results certainly play only a minor part in studying the statistical independence behaviour of the generated sequences, but they show that modified explicit inversive congruential pseudorandom numbers have some attractive properties at least regarding their two-dimensional discrepancy. The method of proof relies heavily on a thorough analysis of certain exponential sums.     

Dependent radius marks of Laguerre tessellations: a case study

We study a particular marked three-dimensional point process sample that represents a Laguerre tessellation. It comes from a polycrystalline sample of aluminium alloy material. The ‘points’ are the cell generators while the ‘marks’ are radius marks that control the size and shape of the tessellation cells. Our statistical mark correlation analyses show that the marks of the sample are in clear and plausible spatial correlation: the marks of generators close together tend to be small and similar and the form of the correlation functions does not justify geostatistical marking. We show that a simplified modelling of tessellations by Laguerre tessellations with independent radius marks may lead to wrong results. When we started from the aluminium alloy data and generated random marks by random permutation we obtained tessellations with characteristics quite different from the original ones. We observed similar behaviour for simulated Laguerre tessellations. This fact, which seems to be natural for the given data type, makes fitting of models to empirical Laguerre tessellations quite difficult: the generator points and radius marks have to be modelled simultaneously. This may imply that the reconstruction methods are more efficient than point-process modelling if only samples of similar Laguerre tessellations are needed. We also found that literature recipes for bandwidth choice for estimating correlation functions should be used with care.     

BSDEs driven by time-changed Lévy noises with non-Lipschitz generators

This paper deals with a class of backward stochastic differential equations (BSDEs for short) driven by time-changed Lévy noises. The existence and uniqueness of L^p(p ? 2) solutions for this kind of BSDEs with non-Lipschitz generators are obtained, which extend the corresponding results of Di Nunno and Sjursen (2014) [Stochastic Process. Appl. 124(4):1679-1709]. Furthermore, representation theorem for generators as well as converse comparison theorem for this kind of BSDEs are also studied.     

Matrix generators for pseudo-random vector generation

In this paper a method for the direct generation of pseudo-random vectors is considered. Thereby the n-th pseudo-random vector is recursively generated from the (n?1)-th pseudo-random vector by multiplication with a matrix. The period lengths of this generator type are examined and characterized. Furthermore it is shown that some popular pseudo-random number generators can be regarded as special cases of this method.     

Quadratic reflected BSDEs and related obstacle problems for PDEs

Abstract

In this paper, we study a kind of reflected backward stochastic differential equations (BSDEs) whose generators are of quadratic growth in z and linear growth in y. We first give an estimate of solutions to such reflected BSDEs. Then under the condition that the generators are convex with respect to z, we can obtain a comparison theorem, which implies the uniqueness of solutions for this kind of reflected BSDEs. Besides, the assumption of convexity also leads to a stability property in the spirit of above estimate. We further establish the nonlinear Feynman-Kac formula of the related obstacle problems for partial differential equations (PDEs) in our framework. At last, a numerical example is given to illustrate the applications of our theoretical results, as well as its connection with an optimal stopping time problem.     

On tail index estimation based on multivariate data

This article is devoted to the study of tail index estimation based on i.i.d. multivariate observations, drawn from a standard heavy-tailed distribution, that is, of which Pareto-like marginals share the same tail index. A multivariate central limit theorem for a random vector, whose components correspond to (possibly dependent) Hill estimators of the common tail index α, is established under mild conditions. We introduce the concept of (standard) heavy-tailed random vector of tail index α and show how this limit result can be used in order to build an estimator of α with small asymptotic mean squared error, through a proper convex linear combination of the coordinates. Beyond asymptotic results, simulation experiments illustrating the relevance of the approach promoted are also presented.     

概述随机非线性动力系统在经济周期研究中的运用

非线性动力学为经济周期的动态分析提供了全新的思路和方法,打破了传统的均衡线性分析的范式。考虑到复杂经济系统中本质的表现为非线性,而且还不可避免地存在随机噪声。因此,为了深入地探究经济周期的动力学形成机理,将随机非线性动力系统引入到经济周期问题的研究中。通过研究随机模型的稳定性、分岔、混沌和随机最优控制,实现对宏观经济动态演化和运行的评价、监测与控制。这不仅拓宽了随机动力学在宏观经济领域中的应用,而且也为宏观经济运行的研究提供了一个全新的思路和方法。     

On the distribution of quadratic forms in normal random variables

This paper provides necessary and sufficient conditions for a quadratic form in singular normal random variables to be distributed as a given linear combination of independent noncentral chi-square variables. Using this result, an extension of Cochran's theorem to quadratic forms of noncentral chi-square variables is derived.     

Asymptotic properties of nonparametric regression for long memory random fields

The kernel estimator of spatial regression function is investigated for stationary long memory (long range dependent) random fields observed over a finite set of spatial points. A general result on the strong consistency of the kernel density estimator is first obtained for the long memory random fields, and then, under some mild regularity assumptions, the asymptotic behaviors of the regression estimator are established. For the linear long memory random fields, a weak convergence theorem is also obtained for kernel density estimator. Finally, some related issues on the inference of long memory random fields are discussed through a simulation example.     

Estimating Archimedean Copulas in High Dimensions

Abstract. This article presents a novel estimation procedure for high‐dimensional Archimedean copulas. In contrast to maximum likelihood estimation, the method presented here does not require derivatives of the Archimedean generator. This is computationally advantageous for high‐dimensional Archimedean copulas in which higher‐order derivatives are needed but are often difficult to obtain. Our procedure is based on a parameter‐dependent transformation of the underlying random variables to a one‐dimensional distribution where a minimum‐distance method is applied. We show strong consistency of the resulting minimum‐distance estimators to the case of known margins as well as to the case of unknown margins when pseudo‐observations are used. Moreover, we conduct a simulation comparing the performance of the proposed estimation procedure with the well‐known maximum likelihood approach according to bias and standard deviation.     

Random number generators for the generalized Birnbaum–Saunders distribution

The generalized Birnbaum–Saunders distribution pertains to a class of lifetime models including both lighter and heavier tailed distributions. This model adapts well to lifetime data, even when outliers exist, and has other good theoretical properties and application perspectives. However, statistical inference tools may not exist in closed form for this model. Hence, simulation and numerical studies are needed, which require a random number generator. Three different ways to generate observations from this model are considered here. These generators are compared by utilizing a goodness-of-fit procedure as well as their effectiveness in predicting the true parameter values by using Monte Carlo simulations. This goodness-of-fit procedure may also be used as an estimation method. The quality of this estimation method is studied here. Finally, through a real data set, the generalized and classical Birnbaum–Saunders models are compared by using this estimation method.     

Robust estimation and wavelet thresholding in partially linear models

This paper is concerned with a semiparametric partially linear regression model with unknown regression coefficients, an unknown nonparametric function for the non-linear component, and unobservable Gaussian distributed random errors. We present a wavelet thresholding based estimation procedure to estimate the components of the partial linear model by establishing a connection between an l ₁-penalty based wavelet estimator of the nonparametric component and Huber’s M-estimation of a standard linear model with outliers. Some general results on the large sample properties of the estimates of both the parametric and the nonparametric part of the model are established. Simulations are used to illustrate the general results and to compare the proposed methodology with other methods available in the recent literature.     

A semiparametric multilevel survival model

Summary.  We propose a semiparametric multilevel survival model for clustered duration data in which the effect of a continuous covariate is represented by an unspecified, possibly non-linear, function. This model makes no distributional assumption about the cluster level random effects. The performance of the method is assessed via Monte Carlo simulations. The model is applied in an analysis of first-birth intervals in Bangladesh to examine period effects in the timing of first births, while allowing for clustering within communities; the analysis reveals a non-linear trend in the first-birth interval over time.