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.     

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.     

Ruin problem of a two-dimensional fractional Brownian motion risk process

This paper investigates ruin probability and ruin time of a two-dimensional fractional Brownian motion risk process. The net loss process of an insurance company is modeled by a fractional Brownian motion. The two-dimensional fractional Brownian motion risk process models the surplus processes of an insurance and a reinsurance company, where the net loss is divided between them in some specified proportions. The ruin problem considered is that of the two-dimensional risk process first entering the negative quadrant, that is, the simultaneous ruin problem. We derive both asymptotics of the ruin probability and approximations of the scaled conditional ruin time as the initial capital tends to infinity.     

Approximation to two independent Gaussian processes from a unique Lévy process and applications

Abstract

In this article, we construct two families of processes, from a unique Lévy process, the finite dimensional distributions of which converge in law towards the finite dimensional distributions of the two independent Gaussian processes. As applications of this result, we obtain families of processes that converge in law towards fractional Brownian motion, sub-fractional Brownian motion and bifractional Brownian motion, respectively.     

Operating Characteristics for Group Sequential Trials Monitored Under Fractional Brownian Motion

Brownian motion has been used to derive stopping boundaries for group sequential trials, however, when we observe dependent increment in the data, fractional Brownian motion is an alternative to be considered to model such data. In this article we compared expected sample sizes and stopping times for different stopping boundaries based on the power family alpha spending function under various values of Hurst coefficient. Results showed that the expected sample sizes and stopping times will decrease and power increases when the Hurst coefficient increases. With same Hurst coefficient, the closer the boundaries are to that of O'Brien-Fleming, the higher the expected sample sizes and stopping times are; however, power has a decreasing trend for values start from H = 0.6 (early analysis), 0.7 (equal space), 0.8 (late analysis). We also illustrate study design changes using results from the BHAT study.     

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.     

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.     

Innovative methods for modeling of scale invariant processes

This paper presents some innovative methods for modeling discrete scale invariant (DSI) processes and evaluation of corresponding parameters. For the case where the absolute values of the increments of DSI processes are in general increasing, we consider some moving sample variance of the increments and present some heuristic algorithm to characterize successive scale intervals. This enables us to estimate scale parameter of such DSI processes. To present some superior structure for the modeling of DSI processes, we consider the possibility that the variations inside the prescribed scale intervals show some further self-similar behavior. Such consideration enables us to provide more efficient estimators for Hurst parameters. We also present two competitive estimation methods for the Hurst parameters of self-similar processes with stationary increments and prove their efficiency. Using simulated samples of some simple fractional Brownian motion, we show that our estimators of Hurst parameter are more efficient as compared with the celebrated methods of convex rearrangement and quadratic variation. Finally we apply the proposed methods to evaluate DSI behavior of the S&P500 indices in some period.     

A strong convergence to the tempered fractional Brownian motion

In this paper, we give a strong approximation of the tempered fractional Brownian motion via transport processes and derive the rate of convergence.     

Statistical inference on the drift parameter in fractional Brownian motion with a deterministic drift

In statistical inference on the drift parameter a in the fractional Brownian motion W^H_t with the Hurst parameter H ∈ (0, 1) with a constant drift Y^H_t = at + W^H_t, there is a large number of options how to do it. We may, for example, base this inference on the properties of the standard normal distribution applied to the differences between the observed values of the process at discrete times. Although such methods are very simple, it turns out that more appropriate is to use inverse methods. Such methods can be generalized to non constant drift. For the hypotheses testing about the drift parameter a, it is more proper to standardize the observed process, and to use inverse methods based on the first exit time of the observed process of a pre-specified interval until some given time. These procedures are illustrated, and their times of decision are compared against the direct approach. Other generalizations are possible when the random part is a symmetric stochastic integral of a known, deterministic function with respect to fractional Brownian motion.     

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.     

Alternative forms of fractional Brownian motion

It is pointed out that two contradictory definitions of fractional Brownian motion are well-established, one prevailing in the probabilistic literature, the other in the econometric literature. Each is associated with a different definition of nonstationary fractional time series, arising in functional limit theorems based on such series. These various definitions have occasionally led to some confusion. The paper discusses the definitions and attempts a clarification.     

Group sequential tests under fractional Brownian motion in monitoring clinical trials

Group sequential tests have been effective tools in monitoring long term clinical trials. There have been several popular discrete sequential boundaries proposed for modeling interim analysis of clinical trials under the assumption of Brownian motion for the stochastic processes generated from test statistics. In this paper, we study the five sequential boundaries in Lan and DeMets (Biometrika 70:659–663, 1983) under the fractional Brownian motion. The fractional Brownian includes the classic Brownian motion as a special case. An example from a real data set is used to illustrate the applications of the boundaries.     

Fractional Brownian motion and long term clinical trial recruitment

Prediction of recruitment in clinical trials has been a challenging task. Many methods have been studied, including models based on Poisson process and its large sample approximation by Brownian motion (BM); however, when the independent incremental structure is violated for BM model, we could use fractional Brownian motion to model and approximate the underlying Poisson processes with random rates. In this paper, fractional Brownian motion (FBM) is considered for such conditions and compared to BM model with illustrated examples from different trials and simulations.     

The local fractional bootstrap

We introduce a bootstrap procedure for high‐frequency statistics of Brownian semistationary processes. More specifically, we focus on a hypothesis test on the roughness of sample paths of Brownian semistationary processes, which uses an estimator based on a ratio of realized power variations. Our new resampling method, the local fractional bootstrap, relies on simulating an auxiliary fractional Brownian motion that mimics the fine properties of high‐frequency differences of the Brownian semistationary process under the null hypothesis. We prove the first‐order validity of the bootstrap method, and in simulations, we observe that the bootstrap‐based hypothesis test provides considerable finite‐sample improvements over an existing test that is based on a central limit theorem. This is important when studying the roughness properties of time series data. We illustrate this by applying the bootstrap method to two empirical data sets: We assess the roughness of a time series of high‐frequency asset prices and we test the validity of Kolmogorov's scaling law in atmospheric turbulence data.     

Models for residual time to AIDS

The distributions of the time from Human Immunodeficiency Virus (HIV) infection to the onset of Acquired Immune Deficiency Syndrome (AIDS) and of the residual time to AIDS diagnosis are important for modeling the growth of the AIDS epidemic and for predicting onset of the disease in an individual. Markers such as CD4 counts carry valuable information about disease progression and therefore about the two survival distributions. Building on the framework set out by Jewell and Kalbfleisch (1992), we study these two survival distributions based on stochastic models for the marker process (X(t)) and a marker-dependent hazard (h()). We examine various plausible CD4 marker processes and marker-dependent hazard functions for AIDS proposed in recent literature. For a random effects plus Brownian motion marker process X(t)=(a+bt+BM(t))⁴, where a has a normal distribution, b<0 is an unknown parameter and BM(t) is Brownian motion, and marker-dependent hazard h(X(t)), we prove that, given CD4 cell count X(t), the residual time to AIDS distribution does not depend on the time since infection t. Using simulation and numerical integration, we find the marginal incubation period distribution, the marginal hazard and the residual time distribution for several combinations of marker processes and marker-dependent hazards. An example using data from the Multicenter AIDS Cohort Study is given. A simple regression model relating the cube root of residual time to AIDS to CD4 count is suggested.     

Fractional Brownian motion and clinical trials

The purpose of this paper is to extend the widely used classical Brownian motion technique for monitoring clinical trial data to a larger class of stochastic processes, i.e. fractional Brownian motion, and compare these results. The beta-blocker heart attack trial is presented as an example to illustrate both methods.     

Some Asymptotic Formulas for a Brownian Motion With a Regular Variation From a Parabolic Domain

Consider a Brownian motion with a regular variation starting at an interior point of a domain D in R^{d + 1}, d ? 1 and let τ_D denote the first time the Brownian motion exits from D. Estimates with exact constants for the asymptotics of log?P(τ_D > T) are given for T → ∞, depending on the shape of the domain D and the order of the regular variation. Furthermore, the asymptotically equivalence are obtained. The problem is motivated by the early results of Lifshits and Shi, Li in the first exit time, and Karamata in the regular variation. The methods of proof are based on their results and the calculus of variations.     

Variance Estimation for Fractional Brownian Motions with Fixed Hurst Parameters

Some real-world phenomena in geo-science, micro-economy, and turbulence, to name a few, can be effectively modeled by a fractional Brownian motion indexed by a Hurst parameter, a regularity level, and a scaling parameter σ², an energy level. This article discusses estimation of a scaling parameter σ² when a Hurst parameter is known. To estimate σ², we propose three approaches based on maximum likelihood estimation, moment-matching, and concentration inequalities, respectively, and discuss the theoretical characteristics of the estimators and optimal-filtering guidelines. We also justify the improvement of the estimation of σ² when a Hurst parameter is known. Using the three approaches and a parametric bootstrap methodology in a simulation study, we compare the confidence intervals of σ² in terms of their lengths, coverage rates, and computational complexity and discuss empirical attributes of the tested approaches. We found that the approach based on maximum likelihood estimation was optimal in terms of efficiency and accuracy, but computationally expensive. The moment-matching approach was found to be not only comparably efficient and accurate but also computationally fast and robust to deviations from the fractional Brownian motion model.     

The Asymptotic Behavior of a Brownian Motion with a Drift from a Random Domain

Consider a Brownian motion with drift starting at an interior point of a random domain D in R ^d+1, d ≥ 1, let τ _D denote the first time the Brownian motion exits from D. Estimates with exact constants for the asymptotics of log P(τ _D  > T) are given for T → ∞, depending on the shape of the domain D and the order of the drift. The problem is motivated by the model in insurance and early works of Lifshits and Shi. The methods of proof are based on the calculus of variations and early works of Li, Lifshits and Shi in the drift free case.