Some joint distributions for conditional random walks

Joint distributions concerning maxima, minima, and their indices are determined for certain conditional random walks called Bernoulli excursion and Bernoulli meander. The distribution of the local time of these processes is treated by generating function technique. Limiting distributions are also given, providing some partial results for Brownian excursion and meander.     

A month carlo approximation of the distributions of the maximum of various brownjan bridges

A Brownian bridge of order a is the weak limit of a residual partial sum obtained from regression fitting. When q=0, the process is teh usual Brownian bridge, and the distribution of the maximum is known analytically. For q>1, the reflection principle does not easily apply. The supremum distributions are approximated by a Monte Carlo technique. Tables of these distributions, as well as a finite sample size correction are given.     

A note on the residual median process

We study the residual median process, defined as the median of those observations which are greater than a number t. Using appropriate limit theorems, it is shown that the stochastic process converges in law to a Gaussian process defined in terms of a Brownian bridge.     

Computing the distribution of sequential Hölder norms of the Brownian motion

ABSTRACT

In this article, we give explicit formulas and study practical computations for the distribution function of sequential Hölder norms of a Brownian motion and of a Brownian bridge. We also discuss some statistical applications in the detection of some short “epidemic” changes in a sample.     

Efficient computation of the cdf of the maximal difference between a Brownian bridge and its concave majorant

In this paper, we describe two computational methods for calculating the cumulative distribution function and the upper quantiles of the maximal difference between a Brownian bridge and its concave majorant. The first method has two different variants that are both based on a Monte Carlo approach, whereas the second uses the Gaver–Stehfest (GS) algorithm for the numerical inversion of the Laplace transform. If the former method is straightforward to implement, it is very much outperformed by the GS algorithm, which provides a very accurate approximation of the cumulative distribution as well as its upper quantiles. Our numerical work has a direct application in statistics: the maximal difference between a Brownian bridge and its concave majorant arises in connection with a nonparametric test for monotonicity of a density or regression curve on [0,1]. Our results can be used to construct very accurate rejection region for this test at a given asymptotic level.     

Stochastically Curtailed Tests Under Fractional Brownian Motion

Stochastic curtailment has been considered for the interim monitoring of group sequential trials (Davis and Hardy, 1994). Statistical boundaries in Davis and Hardy (1994) were derived using theory of Brownian motion. In some clinical trials, the conditions of forming a Brownian motion may not be satisfied. In this paper, we extend the computations of Brownian motion based boundaries, expected stopping times, and type I and type II error rates to fractional Brownian motion (FBM). FBM includes Brownian motion as a special case. Designs under FBM are compared to those under Brownian motion and to those of O’Brien–Fleming type tests. One- and two-sided boundaries for efficacy and futility monitoring are also discussed. Results show that boundary values decrease and error rates deviate from design levels when the Hurst parameter increases from 0.1 to 0.9, these changes should be considered when designing a study under FBM.     

On some results for Bernoulli excursions

We give a short survey of some distribution results for Bernoulli excursion. The main emphasis is given to quantities expressible in terms of additive functionals, such as local times, area, etc.     

Heuristic approach to some laws for Brownian motion

Csàki and Vincze have shown that for an elementary tied-down random walk, the pair (maximum, instant of maximum) has the same law as (time spent in (0, 1/2), time spent above 1/2). Formal passage to the limit indicates that the former pair has for a Brownian bridge the same law as (local time at 0, duration of positivity). A quadrivariate density of Karatzas and Shreve and an equivalence for Brownian motion with drift follow.     

The Representation of Hypergeometric Random Variables using Independent Bernoulli Random Variables

In this paper, we show that a hypergeometric random variable can be represented as a sum of independent Bernoulli random variables that are, except in degenerate cases, not identically distributed. In the proof, we use the factorial moment generating function. An asymptotic result on the probabilities of the Bernoulli random variables in the sum is also presented. Numerical examples are used to illustrate the results.     

Brownian Bridge and Other Path-dependent Gaussian Processes Vectorial Simulation

The iterative simulation of the Brownian bridge is well known. In this article, we present a vectorial simulation alternative based on Gaussian processes for machine learning regression that is suitable for interpreted programming languages implementations. We extend the vectorial simulation of path-dependent trajectories to other Gaussian processes, namely, sequences of Brownian bridges, geometric Brownian motion, fractional Brownian motion, and Ornstein–Ulenbeck mean reversion process.     

On a Simultaneous Characterization of the Poisson and Bernoulli Distributions

Kimeldorf et al. (1981) established a simultaneous characterization of the Poisson and Bernoulli distributions. In this note two variants of the authors' characterizing condition are considered each of which is shown also to characterize simultaneously the Poisson and Bernoulli distributions.     

A Latent Process Model for Joint Modeling of Events and Marker

The paper formulates joint modeling of a counting process and a sequence of longitudinal measurements, governed by a common latent stochastic process. The latent process is modeled as a function of explanatory variables and a Brownian motion process. The conditional likelihood given values of the latent process at the measurement times, has been drawn using Brownian bridge properties; then integrating over all possible values of the latent process at the measurement times leads to the desired joint likelihood. An estimation procedure using joint likelihood and a numerical optimization is described. The method is applied to the study of cognitive decline and Alzheimer's disease.     

Limit theorems for dependent Bernoulli variables with statistical inference

We consider a class of dependent Bernoulli variables where the conditional success probability is a linear combination of the last few trials and the original success probability. We obtain its limit theorems including the strong law of large numbers, weak invariance principle, and law of the iterated logarithm. We also derive some statistical inference results which make the model applicable. Simulation results are exhibited as well to show that with small sample size the convergence rate is satisfying and the proposed estimators behave well.     

Invariance principles for jackknifing U-statistics for finite population sampling and some applications

For simple random sampling (without replacement) from a finite population, suitable stochastic processes are constructed from the entire sequence of jackknife estimators based on smooth functions of U-statistics and these are approximated (in distributions) by some Brownian bridge processes. Strong convergence of the Tukey estimator of the variance of a jackknife U-statistic has been interpreted suitably and established. Some applications of these results in sequential analysis relating to finite population sampling are also considered.     

On drift parameter estimation in models with fractional Brownian motion

We consider a stochastic differential equation involving standard and fractional Brownian motion with unknown drift parameter to be estimated. We investigate the standard maximum likelihood estimate of the drift parameter, two non-standard estimates and three estimates for the sequential estimation. Model strong consistency and some other properties are proved. The linear model and Ornstein–Uhlenbeck model are studied in detail. As an auxiliary result, an asymptotic behaviour of the fractional derivative of the fractional Brownian motion is established.     

James Bernoulli与《推测术》

一、JamesBernoulli与Bernoulli家族1713年 ,也就是JamesBernoulli死后的第 8年 ,他的伟大著作《推测术》终于在数学家们的期盼中出版了。如果说赌博是概率的母亲的话 ,那么经过几千年的孕育 ,概率终于以JamesBernoulli的《推测术》(Arsconjectandi)的出版为标志诞生了。对此著名统计史学家IanHacking的评价是“此书标志着概率概念漫长形成过程的终结与数学概率论的开始”。16 5 4年 12月 2 7日 ,JamesBernoulli出生在瑞士的一个小城Basel。他是这个伟大的数学家族中的第一个数学家 ,也是最著名的一位。他的祖先早年为了逃避宗教迫害 …     

Asymptotic normality of MRPP statistics from invariance principles of u-statistics

Multi-response permutation procedures (MRPP) were recently introduced to test differences between a priori classified groups of objects ( Mielke, Berry Johnson, 1976; Mielke, 1979 ). The null distributions of the MRPP statistics were initially conjectured to be asymptotically normal for some specified conditions within the setting of a sequence of finite populations due to Madow ( 1948 ).

Asymptotic normality of a class of MRPP statistics (under the null hypothesis) is shown in two cases: (i) the setting which considers the populations to be the samples resulting from sequential independent identically distributed (i.i.d.) sampling (sampling from infinite populations) and (ii) the setting of a sequence of increasingly large finite populations (sampling from finite populations). The results are direct applications of the weak convergence of a U-statistic process in the i.i.d. case to a Brownian motion (Bhattacharyya and Sen, 1977) and of the weak convergence of a U-statistic process in the finite populations case to a Brownian bridge (Sen, 1972). The conditions are milder for the i.i.d. case than for the finite populations case. However, neither case provides a restriction of a practical consequence in applications of MRPP. In either case, convergence is shown to depend on the asymptotic ratios of the group sizes to the population size.     

Kolmogorov–Smirnov test for life test data with hybrid censoring

This work considers goodness-of-fit for the life test data with hybrid censoring. An alternative representation of the Kolmogorov–Smirnov (KS) statistics is provided under Type-I censoring. The alternative representation leads us to approximate the limiting distributions of the KS statistic as a functional of the Brownian bridge for Type-II, Type-I hybrid, and Type-II hybrid censored data. The approximated distributions are used to obtain the critical values of the tests in this context. We found that the proposed KS test procedure for Type-II censoring has more power than the available one(s) in literature.     

Repeated Confidence Intervals Under Fractional Brownian Motion in Long-Term Clinical Trials

Repeated confidence interval (RCI) is an important tool for design and monitoring of group sequential trials according to which we do not need to stop the trial with planned statistical stopping rules. In this article, we derive RCIs when data from each stage of the trial are not independent thus it is no longer a Brownian motion (BM) process. Under this assumption, a larger class of stochastic processes fractional Brownian motion (FBM) is considered. Comparisons of RCI width and sample size requirement are made to those under Brownian motion for different analysis times, Type I error rates and number of interim analysis. Power family spending functions including Pocock, O'Brien-Fleming design types are considered for these simulations. Interim data from BHAT and oncology trials is used to illustrate how to derive RCIs under FBM for efficacy and futility monitoring.     

A Simple Distribution for the Sum of Correlated,Exchangeable Binary Data

This article describes a generalization of the binomial distribution. The closed form probability function for the probability of k successes out of n correlated, exchangeable Bernoulli trials depends on the number of trials and its two parameters: the common success probability and the common correlation. The distribution is derived under the assumption that the common correlation between all pairs of Bernoulli trials remains unchanged conditional on successes in all completed trials. The distribution was developed to model bond defaults but may be suited to biostatistical applications involving clusters of binary data encountered in repeated measurements or toxicity studies of families of organisms. Maximum likelihood estimates for the parameters of the distribution are found for a set of binary data from a developmental toxicity study on litters of mice.