Spectral estimates for high‐frequency sampled continuous‐time autoregressive moving average processes

In this article, we consider a continuous‐time autoregressive moving average (CARMA) process driven by either a symmetric α‐stable Lévy process with α ∈ (0,2) or a symmetric Lévy process with finite second moments. In the asymptotic framework of high‐frequency data within a long time interval, we establish a consistent estimate for the normalized power transfer function by applying a smoothing filter to the periodogram of the CARMA process. We use this result to propose an estimator for the parameters of the CARMA process and exemplify the estimation procedure by a simulation study.     

Local Whittle estimation of multi‐variate fractionally integrated processes

This article derives a semi‐parametric estimator of multi‐variate fractionally integrated processes covering both stationary and non‐stationary values of d. We utilize the notion of the extended discrete Fourier transform and periodogram to extend the multi‐variate local Whittle estimator of Shimotsu (2007) to cover non‐stationary values of d. Consistency and asymptotic normality is shown for d ∈ (?1/2,∞). A simulation study illustrates the performance of the proposed estimator for relevant sample sizes. Empirical justification of the proposed estimator is shown through an empirical analysis of log spot exchange rates. We find that the log spot exchange rates of Germany, United Kingdom, Japan, Canada, France, Italy and Switzerland against the US Dollar for the period January 1974 until December 2001 are well decribed as I(1) processes.     

Testing non‐parametric hypotheses for stationary processes by estimating minimal distances

In this article, new tests for non‐parametric hypotheses in stationary processes are proposed. Our approach is based on an estimate of the L²‐distance between the spectral density matrix and its best approximation under the null hypothesis. We explain the main idea in the problem of testing for a constant spectral density matrix and in the problem of comparing the spectral densities of several correlated stationary time series. The method is based on direct estimation of integrals of the spectral density matrix and does not require the specification of smoothing parameters. We show that the limit distribution of the proposed test statistic is normal and investigate the finite sample properties of the resulting tests by means of a small simulation study.     

On the properties of the periodogram of a stationary long‐memory process over different epochs with applications

This article studies the asymptotic properties of the discrete Fourier transforms (DFT) and the periodogram of a stationary long‐memory time series over different epochs. The main theoretical result is a novel bound for the covariance of the DFT ordinates evaluated on two distinct epochs, which depends explicitly on the Fourier frequencies and the gap between the epochs. This result is then applied to obtain the limiting distribution of some nonlinear functions of the periodogram over different epochs, under the additional assumption of gaussianity. We then apply this result to construct an estimator of the memory parameter based on the regression in a neighbourhood of the zero‐frequency of the logarithm of the averaged periodogram, obtained by computing the empirical mean of the periodogram over adjacent epochs. It is shown that replacing the periodogram by its average has an effect similar to the frequency domain pooling to reduce the variance of the estimate. We also propose a simple procedure to test the stationarity of the memory coefficient. A limited Monte Carlo experiment is presented to support our findings.     

Bootstrapping the Local Periodogram of Locally Stationary Processes

Abstract. Locally stationary processes are non‐stationary stochastic processes the second‐order structure of which varies smoothly over time. In this paper, we develop a method to bootstrap the local periodogram of a locally stationary process. Our method generates pseudo local periodogram ordinates by combining a parametric time and non‐parametric frequency domain bootstrap approach. We first fit locally a time varying autoregressive model so as to capture the essential characteristics of the underlying process. A locally calculated non‐parametric correction in the frequency domain is then used so as to improve upon the locally parametric autoregressive fit. As an application, we investigate theoretically the asymptotic properties of the bootstrap method proposed applied to the class of local spectral means, local ratio statistics and local spectral density estimators. Some simulations demonstrate the ability of our method to give accurate estimates of the quantities of interest in finite sample situations and an application to a real‐life data‐set is presented.     

Consistent Estimation of the Value at Risk When the Error Distribution of the Volatility Model is Misspecified

A two‐step approach for conditional value at risk estimation is considered. First, a generalized quasi‐maximum likelihood estimator is employed to estimate the volatility parameter, then the empirical quantile of the residuals serves to estimate the theoretical quantile of the innovations. When the instrumental density h of the generalized quasi‐maximum likelihood estimator is not the Gaussian density, both the estimations of the volatility and of the quantile are generally asymptotically biased. However, the two errors counterbalance and lead to a consistent estimator of the value at risk. We obtain the asymptotic behavior of this estimator and show how to choose optimally h.     

NON‐PARAMETRIC ESTIMATION OF HIGH‐FREQUENCY SPOT VOLATILITY FOR BROWNIAN SEMIMARTINGALE WITH JUMPS

The availability of high‐frequency financial data has led to substantial improvements in our understanding of financial volatility. Most existing literature focuses on estimating the integrated volatility over a fixed period. This article proposes a non‐parametric threshold kernel method to estimate the time‐dependent spot volatility and jumps when the underlying price process is governed by Brownian semimartingale with finite activity jumps. The threshold kernel estimator combines the threshold estimation for integrated volatility and the kernel filtering approach for spot volatility when the price process is driven only by diffusions without jumps. The estimator proposed is consistent and asymptotically normal and has the same rate of convergence as the estimator studied by Kristensen (2010) in a setting without jumps. The Monte Carlo simulation study shows that the proposed estimator exhibits excellent performance over a wide range of jump sizes and for different sampling frequencies. An empirical example is given to illustrate the potential applications of the proposed method.     

A Modified Nonparametric Prewhitened Covariance Estimator

Abstract. This paper proposes a fully modified version of the spectral matrix estimator (and the long‐run variance estimator as a special case) proposed originally by Xiao and Linton [Journal of Time Series Analysis (2002) Vol. 23, pp. 215–250], and derives its asymptotic results. A striking feature of the modified spectral matrix estimator is to achieve the convergence rate of O(T ^?8/9) in the mean squared error (MSE), which is usually achieved under the fourth‐order spectral window. However, this estimator does not sacrifice the positive definiteness of the resulting estimate for the rate improvement; it is Hermitian and positive definite in finite samples by construction. The faster convergence rate is established by a multiplicative bias correction of the crude spectral estimator under the second‐order spectral window. The approximations to some sensible definitions of the MSE of the estimator and the bandwidths that minimize the asymptotic MSEs are also derived. Monte Carlo results indicate that for a wide variety of processes the modified spectral matrix estimator reduces the bias without inflating the variance and thus improves the MSE, compared with the crude, bias‐uncorrected estimator.     

UNIFIED INTERVAL ESTIMATION FOR RANDOM COEFFICIENT AUTOREGRESSIVE MODELS

The consistency of the quasi‐maximum likelihood estimator for random coefficient autoregressive models requires that the coefficient be a non‐degenerate random variable. In this article, we propose empirical likelihood methods based on weighted‐score equations to construct a confidence interval for the coefficient. We do not need to distinguish whether the coefficient is random or deterministic and whether the process is stationary or non‐stationary, and we present two classes of equations depending on whether a constant trend is included in the model. A simulation study confirms the good finite‐sample behaviour of our resulting empirical likelihood‐based confidence intervals. We also apply our methods to study US macroeconomic data.     

Inference for the Fourth‐Order Innovation Cumulant in Linear Time Series

The rescaled fourth‐order cumulant of the unobserved innovations of linear time series is an important parameter in statistical inference. This article deals with the problem of estimating this parameter. An existing nonparametric estimator is first discussed, and its asymptotic properties are derived. It is shown how the autocorrelation structure of the underlying process affects the behaviour of the estimator. Based on our findings and on an important invariance property of the parameter of interest with respect to linear filtering, a pre‐whitening‐based nonparametric estimator of the same parameter is proposed. The estimator is obtained using the filtered time series only; that is, an inversion of the pre‐whitening procedure is not required. The asymptotic properties of the new estimator are investigated, and its superiority is established for large classes of stochastic processes. It is shown that for the particular estimation problem considered, pre‐whitening can reduce the variance and the bias of the estimator. The finite sample performance of both estimators is investigated by means of simulations. The new estimator allows for a simple modification of the multiplicative frequency domain bootstrap, which extends its considerable range of validity. Furthermore, the problem of testing hypotheses about the rescaled fourth‐order cumulant of the unobserved innovations is also considered. In this context, a simple test for Gaussianity is proposed. Some real‐life data applications are presented.     

Inference on Multivariate Heteroscedastic Time Varying Random Coefficient Models

In this article, we introduce the general setting of a multivariate time series autoregressive model with stochastic time‐varying coefficients and time‐varying conditional variance of the error process. This allows modelling VAR dynamics for non‐stationary time series and estimation of time‐varying parameter processes by the well‐known rolling regression estimation techniques. We establish consistency, convergence rates, and asymptotic normality for kernel estimators of the paths of coefficient processes and provide pointwise valid standard errors. The method is applied to a popular seven‐variable dataset to analyse evidence of time variation in empirical objects of interest for the DSGE (dynamic stochastic general equilibrium) literature.     

Estimation of the long‐memory stochastic volatility model parameters that is robust to level shifts and deterministic trends

I provide conditions under which the trimmed FDQML estimator, advanced by McCloskey (2010) in the context of fully parametric short‐memory models, can be used to estimate the long‐memory stochastic volatility model parameters in the presence of additive low‐frequency contamination in log‐squared returns. The types of low‐frequency contamination covered include level shifts as well as deterministic trends. I establish consistency and asymptotic normality in the presence or absence of such low‐frequency contamination under certain conditions on the growth rate of the trimming parameter. I also provide theoretical guidance on the choice of trimming parameter by heuristically obtaining its asymptotic MSE‐optimal rate under certain types of low‐frequency contamination. A simulation study examines the finite sample properties of the robust estimator, showing substantial gains from its use in the presence of level shifts. The finite sample analysis also explores how different levels of trimming affect the parameter estimates in the presence and absence of low‐frequency contamination and long‐memory.     

On the Spectral Density of the Wavelet Coefficients of Long‐Memory Time Series with Application to the Log‐Regression Estimation of the Memory Parameter

Abstract. In recent years, methods to estimate the memory parameter using wavelet analysis have gained popularity in many areas of science. Despite its widespread use, a rigorous semi‐parametric asymptotic theory, comparable with the one developed for Fourier methods, is still lacking. In this article, we adapt to the wavelet setting, the classical semi‐parametric framework introduced by Robinson and his co‐authors for estimating the memory parameter of a (possibly) non‐stationary process. Our results apply to a class of wavelets with bounded supports, which include but are not limited to Daubechies wavelets. We derive an explicit expression of the spectral density of the wavelet coefficients and show that it can be approximated, at large scales, by the spectral density of the continuous‐time wavelet coefficients of fractional Brownian motion. We derive an explicit bound for the difference between the spectral densities. As an application, we obtain minimax upper bounds for the log‐scale regression estimator of the memory parameter for a Gaussian process and we derive an explicit expression of its asymptotic variance.     

Robust estimation of the scale and of the autocovariance function of Gaussian short‐ and long‐range dependent processes

A desirable property of an autocovariance estimator is to be robust to the presence of additive outliers. It is well known that the sample autocovariance, being based on moments, does not have this property. Hence, the use of an autocovariance estimator which is robust to additive outliers can be very useful for time‐series modelling. In this article, the asymptotic properties of the robust scale and autocovariance estimators proposed by Rousseeuw and Croux (1993) and Ma and Genton (2000) are established for Gaussian processes, with either short‐ or long‐range dependence. It is shown in the short‐range dependence setting that this robust estimator is asymptotically normal at the rate , where n is the number of observations. An explicit expression of the asymptotic variance is also given and compared with the asymptotic variance of the classical autocovariance estimator. In the long‐range dependence setting, the limiting distribution displays the same behaviour as that of the classical autocovariance estimator, with a Gaussian limit and rate when the Hurst parameter H is less than 3/4 and with a non‐Gaussian limit (belonging to the second Wiener chaos) with rate depending on the Hurst parameter when H ∈ (3/4,1). Some Monte Carlo experiments are presented to illustrate our claims and the Nile River data are analysed as an application. The theoretical results and the empirical evidence strongly suggest using the robust estimators as an alternative to estimate the dependence structure of Gaussian processes.     

Large‐scale volatility models: theoretical properties of professionals’ practice

Abstract. This article examines the way in which GARCH models are estimated and used for forecasting by practitioners in particular using the highly popular Riskmetrics^TM approach. Although it permits sizable computational gains and provide a simple way to impose positive semi‐definitiveness of multivariate version of the model, we show that this approach delivers non‐consistent parameter’ estimates. The novel theoretical result is corroborated by a set of Monte Carlo exercises. A set of empirical applications suggest that this could cause, in general, unreliable forecasts of conditional volatilities and correlations.     

Power Requirement for Mixing Shear‐Thinning Fluids with a Planetary Mixer

The optimal control of processes dealing with non‐Newtonian liquids requires the knowledge and control of the power demand of the mixing equipment. In this context, an extension of the Metzner and Otto concept to planetary mixers is proposed to adapt this concept to planetary mixers. The theoretical part of this work defines modified expressions of Reynolds and power numbers. These definitions introduce a characteristic velocity u_ch that is used to define the parameter K_s. A planetary mixer is employed to experimentally ascertain this guideline. Power consumption measurements carried out by mixing shear‐thinning fluids permit to determine the K_s factor. This factor varies only slightly with the flow behavior index and may be regarded as a defined constant for this geometry. Finally, experiments with an additional shear‐thickening fluid confirm the validity of this approach.     

TWO‐STEP ESTIMATION OF A MULTI‐VARIATE LÉVY PROCESS

Based on the concept of a Lévy copula to describe the dependence structure of a multi‐variate Lévy process, we present a new estimation procedure. We consider a parametric model for the marginal Lévy processes as well as for the Lévy copula and estimate the parameters by a two‐step procedure. We first estimate the parameters of the marginal processes and then estimate in a second step only the dependence structure parameter. For infinite Lévy measures, we truncate the small jumps and base our statistical analysis on the large jumps of the model. Prominent example will be a bivariate stable Lévy process, which allows for analytic calculations and, hence, for a comparison of different methods. We prove asymptotic normality of the parameter estimates from the two‐step procedure, and in particular, we derive the Godambe information matrix, whose inverse is the covariance matrix of the normal limit law. A simulation study investigates the loss of efficiency because of the two‐step procedure and the truncation.     

Square‐Root LASSO for High‐Dimensional Sparse Linear Systems with Weakly Dependent Errors

We study the square‐root LASSO method for high‐dimensional sparse linear models with weakly dependent errors. The asymptotic and non‐asymptotic bounds for the estimation errors are derived. Our results cover a wide range of weakly dependent errors, including α‐mixing, ρ‐mixing, ?‐mixing, and m‐dependent types. Numerical simulations are conducted to show the consistency property of square‐root LASSO. An empirical application to financial data highlights the importance of the results and method.     

A Robbins–Monro Algorithm for Non‐Parametric Estimation of NAR Process with Markov Switching: Consistency

We approach the problem of non‐parametric estimation for autoregressive Markov switching processes. In this context, the Nadaraya–Watson‐type regression functions estimator is interpreted as a solution of a local weighted least‐square problem, which does not admit a closed‐form solution in the case of hidden Markov switching. We introduce a non‐parametric recursive algorithm to approximate the estimator. Our algorithm restores the missing data by means of a Monte Carlo step and estimates the regression function via a Robbins–Monro step. We prove that non‐parametric autoregressive models with Markov switching are identifiable when the hidden Markov process has a finite state space. Consistency of the estimator is proved using the strong α‐mixing property of the model. Finally, we present some simulations illustrating the performances of our non‐parametric estimation procedure.     

Identification of the multiscale fractional Brownian motion with biomechanical applications

Abstract. In certain applications, for instance, biomechanics, turbulence, finance or internet traffic, it seems suitable to model the data by a generalization of a fractional Brownian motion (FBM) for which the Hurst parameter H depends on the frequency as a piece‐wise constant function. These processes are called multiscale fractional Brownian motions. In this article, we provide a statistical study of the multiscale fractional Brownian motions. We developed a method based on wavelet analysis. By using this method, we calculated the frequency changes, estimated the different parameters, tested the goodness‐of‐fit and gave the numerical algorithm. Biomechanical data are then studied with these new tools.