On Estimation of Hurst Parameter Under Noisy Observations

It is widely accepted that some financial data exhibit long memory or long dependence, and that the observed data usually possess noise. In the continuous time situation, the factional Brownian motion B^H and its extension are an important class of models to characterize the long memory or short memory of data, and Hurst parameter H is an index to describe the degree of dependence. In this article, we estimate the Hurst parameter of a discretely sampled fractional integral process corrupted by noise. We use the preaverage method to diminish the impact of noise, employ the filter method to exclude the strong dependence, and obtain the smoothed data, and estimate the Hurst parameter by the smoothed data. The asymptotic properties such as consistency and asymptotic normality of the estimator are established. Simulations for evaluating the performance of the estimator are conducted. Supplementary materials for this article are available online.     

A wavelet lifting approach to long-memory estimation

Reliable estimation of long-range dependence parameters is vital in time series. For example, in environmental and climate science such estimation is often key to understanding climate dynamics, variability and often prediction. The challenge of data collection in such disciplines means that, in practice, the sampling pattern is either irregular or blighted by missing observations. Unfortunately, virtually all existing Hurst parameter estimation methods assume regularly sampled time series and require modification to cope with irregularity or missing data. However, such interventions come at the price of inducing higher estimator bias and variation, often worryingly ignored. This article proposes a new Hurst exponent estimation method which naturally copes with data sampling irregularity. The new method is based on a multiscale lifting transform exploiting its ability to produce wavelet-like coefficients on irregular data and, simultaneously, to effect a necessary powerful decorrelation of those coefficients. Simulations show that our method is accurate and effective, performing well against competitors even in regular data settings. Armed with this evidence our method sheds new light on long-memory intensity results in environmental and climate science applications, sometimes suggesting that different scientific conclusions may need to be drawn.     

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.     

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.     

Estimation Methods of the Long Memory Parameter: Monte Carlo Analysis and Application

Since the seminal paper of Granger & Joyeux (1980), the concept of a long memory has focused the attention of many statisticians and econometricians trying to model and measure the persistence of stationary processes. Many methods for estimating d, the long-range dependence parameter, have been suggested since the work of Hurst (1951). They can be summarized in three classes: the heuristic methods, the semi-parametric methods and the maximum likelihood methods. In this paper, we try by simulation, to verify the two main properties of [dcirc]: the consistency and the asymptotic normality. Hence, it is very important for practitioners to compare the performance of the various classes of estimators. The results indicate that only the semi-parametric and the maximum likelihood methods can give good estimators. They also suggest that the AR component of the ARFIMA (1, d, 0) process has an important impact on the properties of the different estimators and that the Whittle method is the best one, since it has the small mean squared error. We finally carry out an empirical application using the monthly seasonally adjusted US Inflation series, in order to illustrate the usefulness of the different estimation methods in the context of using real data.     

Temporal Aggregation of Stationary and Non-stationary Continuous-Time Processes

Abstract.  We study the autocorrelation structure of aggregates from a continuous-time process. The underlying continuous-time process or some of its higher derivative is assumed to be a stationary continuous-time auto-regressive fractionally integrated moving-average (CARFIMA) process with Hurst parameter H . We derive closed-form expressions for the limiting autocorrelation function and the normalized spectral density of the aggregates, as the extent of aggregation increases to infinity. The limiting model of the aggregates, after appropriate number of differencing, is shown to be some functional of the standard fractional Brownian motion with the same Hurst parameter of the continuous-time process from which the aggregates are measured. These results are then used to assess the loss of forecasting efficiency due to aggregation.     

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.     

Discrete Gamma Distributions: Properties and Parameter Estimations

A two-parameter discrete gamma distribution is derived corresponding to the continuous two parameters gamma distribution using the general approach for discretization of continuous probability distributions. One parameter discrete gamma distribution is obtained as a particular case. A few important distributional and reliability properties of the proposed distribution are examined. Parameter estimation by different methods is discussed. Performance of different estimation methods are compared through simulation. Data fitting is carried out to investigate the suitability of the proposed distribution in modeling discrete failure time data and other count data.     

Long-range dependence and approximate Bayesian computation

In this work, we propose a method for estimating the Hurst index, or memory parameter, of a stationary process with long memory in a Bayesian fashion. Such approach provides an approximation for the posterior distribution for the memory parameter and it is based on a simple application of the so-called approximate Bayesian computation (ABC), also known as likelihood-free method. Some popular existing estimators are reviewed and compared to this method for the fractional Brownian motion, for a long-range binary process and for the Rosenblatt process. The performance of our proposal is remarkably efficient.     

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.     

Multipower variation from generalized difference for fractional integral processes with jumps

This article presents limit theorems of the multipower variation based on a generalized difference for the fractional integral process with jumps observed in high frequency. In particular, we obtain the large number laws for threshold multipower variation and multipower variation and the associated central limit theorems. The limit theorems are applied to estimate Hurst parameter, and the consistence and asymptotic distribution of the estimator are established. These results will provide some new statistical tools to analyze long-memory effect in high-frequency situation.     

Confidence Interval for Shrinkage Parameters in Ridge Regression

In ridge regression, the estimation of ridge parameter k is an important problem. There are several methods available in the literature to do this job some what efficiently. However, no attempts were made to suggest a confidence interval for the ridge parameter using the knwoledge from the data. In this article, we propose a data dependent confidence interval for the ridge parameter k. The method of obtaining the confidence interval is illustrated with the help of a data set. A simulation study indicates that the empirical coverage probability of the suggested confidence intervals are quite high.     

Spatiotemporal prediction for log-Gaussian Cox processes

Space–time point pattern data have become more widely available as a result of technological developments in areas such as geographic information systems. We describe a flexible class of space–time point processes. Our models are Cox processes whose stochastic intensity is a space–time Ornstein–Uhlenbeck process. We develop moment-based methods of parameter estimation, show how to predict the underlying intensity by using a Markov chain Monte Carlo approach and illustrate the performance of our methods on a synthetic data set.     

Estimation of the Hurst parameter in some fractional processes

We propose to estimate the Hurst parameter involved in fractional processes via a method based on the Karhunen–Loève expansion of a Gaussian process. We specifically investigate the cases of the fractional Brownian motion, the fractional Ornstein–Uhlenbeck family and the fractional Brownian bridge. We numerically compare our results with the ones obtained by the maximum-likelihood method, which show the validity of our proposal.     

A comparison of Bayesian model selection based on MCMC with an application to GARCH-type models

This paper presents a comprehensive review and comparison of five computational methods for Bayesian model selection, based on MCMC simulations from posterior model parameter distributions. We apply these methods to a well-known and important class of models in financial time series analysis, namely GARCH and GARCH-t models for conditional return distributions (assuming normal and t-distributions). We compare their performance with the more common maximum likelihood-based model selection for simulated and real market data. All five MCMC methods proved reliable in the simulation study, although differing in their computational demands. Results on simulated data also show that for large degrees of freedom (where the t-distribution becomes more similar to a normal one), Bayesian model selection results in better decisions in favor of the true model than maximum likelihood. Results on market data show the instability of the harmonic mean estimator and reliability of the advanced model selection methods.     

Outlier identification and robust parameter estimation in a zero-inflated Poisson model

The Zero-inflated Poisson distribution has been used in the modeling of count data in different contexts. This model tends to be influenced by outliers because of the excessive occurrence of zeroes, thus outlier identification and robust parameter estimation are important for such distribution. Some outlier identification methods are studied in this paper, and their applications and results are also presented with an example. To eliminate the effect of outliers, two robust parameter estimates are proposed based on the trimmed mean and the Winsorized mean. Simulation results show the robustness of our proposed parameter estimates.     

Semiparametric Global Cross-ratio Models for Bivariate Censored Data

Abstract.  Multivariate correlated failure time data arise in many medical and scientific settings. In the analysis of such data, it is important to use models where the parameters have simple interpretations. In this paper, we formulate a model for bivariate survival data based on the Plackett distribution. The model is an alternative to the Gamma frailty model proposed by Clayton and Oakes. The parameter in this distribution has a very appealing odds ratio interpretation for dependence between the two failure times; in addition, it allows for negative dependence. We develop novel semiparametric estimation and inference procedures for the model. The asymptotic results of the estimator are developed. The performance of the proposed techniques in finite samples is examined using simulation studies; in addition, the proposed methods are applied to data from an observational study in cancer.     

Estimating the Hurst parameter in financial time series via heuristic approaches

This research investigates long memory financial equity markets using three heuristic methodologies namely a proposed modified variance time-aggregated plot, modified rescaled-range plot and periodogram approaches. The intensity of the long memory process is quantified in terms of Hurst parameter (H). Five Malaysian equity market indices are selected in the empirical studies with the inclusion of pre- and post-drastic economic events. Our empirical results evidenced dissimilar long memory behaviours in the different regimes of significant economic events. It is also found that after the short-memory adjustment, all the equity markets exhibited substantial reductions in long memory estimations.     

Minimum Contrast Estimation for Fractional Diffusions

Abstract.  When the Hurst coefficient of a fractional Brownian motion     is greater than 1/2 it is possible to define a stochastic integral with respect to     , as the pathwise limit of Riemann sums, and thus to consider pathwise solutions to fractional diffusion equations. In this paper, we consider the vanishing drift case and assume that the solution X _t is parameterized by θ in a compact parameter space Θ . Our main interest is the estimation of θ based on discrete time, but with very frequent observations. It is shown that the estimation problem in this context is locally asymptotically mixed normal. The asymptotic behaviour of a certain class of minimum contrast estimators is then studied and asymptotic efficiency is discussed.     

Bayesian estimation of the Hurst parameter of fractional Brownian motion

The primary purpose of this study was to find Bayesian estimates for the Hurst dimension of a Fractional Brownian motion with a Beta prior when the process is observed at discrete times. Overestimation is observed though the overestimation is less severe as real H goes up. In addition, the estimated H decreases as Beta parameters go up given an Alpha value. In contrast, the estimated H increases as Alpha parameters go up given a Beta value. For the real-world data, the 2011 daily Taiwan stock index was used and the estimated Hurst index was 0.21.