Source author record

Murad S. Taqqu

Murad S. Taqqu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

29works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2021arXiv

A Survey on Limit Theorems for Toeplitz Type Quadratic Functionals of Stationary Processes and Applications

This is a survey of recent results on central and non-central limit theorems for quadratic functionals of stationary processes. The underlying processes are Gaussian, linear or Lévy-driven linear processes with memory, and are defined either in discrete or continuous time. We focus on limit theorems for Toeplitz and tapered Toeplitz type quadratic functionals of stationary processes with applications in parametric and nonparametric statistical estimation theory. We discuss questions concerning Toeplitz matrices and operators, Fejér-type singular integrals, and Lévy-Itô-type and Stratonovich-type multiple stochastic integrals. These are the main tools for obtaining limit theorems.

preprint2020arXiv

Extensions of Rosenblatt's results on the asymptotic behavior of the prediction error for deterministic stationary sequences

One of the main problem in prediction theory of discrete-time second-order stationary processes $X(t)$ is to describe the asymptotic behavior of the best linear mean squared prediction error in predicting $X(0)$ given $ X(t),$ $-n\le t\le-1$, as $n$ goes to infinity. This behavior depends on the regularity (deterministic or non-deterministic) of the process $X(t)$. In his seminal paper {\it "Some purely deterministic processes" (J. of Math. and Mech.,} 6(6), 801-810, 1957), M. Rosenblatt has described the asymptotic behavior of the prediction error for discrete-time deterministic processes in the following two cases: (a) the spectral density $f(λ)$ of $X(t)$ is continuous and vanishes on an interval, (b) the spectral density $f(λ)$ has a very high order contact with zero. He showed that in the case (a) the prediction error variance behaves exponentially, while in the case (b), it behaves hyperbolically as $n\to\infty$. In this paper, using a new approach, we describe extensions of Rosenblatt's results to broader classes of spectral densities. Examples illustrate the obtained results.

preprint2016arXiv

A unified approach to self-normalized block sampling

The inference procedure for the mean of a stationary time series is usually quite different under various model assumptions because the partial sum process behaves differently depending on whether the time series is short or long-range dependent, or whether it has a light or heavy-tailed marginal distribution. In the current paper, we develop an asymptotic theory for the self-normalized block sampling, and prove that the corresponding block sampling method can provide a unified inference approach for the aforementioned different situations in the sense that it does not require the {\em a priori} estimation of auxiliary parameters. Monte Carlo simulations are presented to illustrate its finite-sample performance. The R function implementing the method is available from the authors.

preprint2016arXiv

Behavior of the generalized Rosenblatt process at extreme critical exponent values

The generalized Rosenblatt process is obtained by replacing the single critical exponent characterizing the Rosenblatt process by two different exponents living in the interior of a triangular region. What happens to that generalized Rosenblatt process as these critical exponents approach the boundaries of the triangle? We show by two different methods that on each of the two symmetric boundaries, the limit is non-Gaussian. On the third boundary, the limit is Brownian motion. The rates of convergence to these boundaries are also given. The situation is particularly delicate as one approaches the corners of the triangle, because the limit process will depend on how these corners are approached. All limits are in the sense of weak convergence in $C[0,1]$. These limits cannot be strengthened to convergence in $L^2(Ω)$.

preprint2016arXiv

On the validity of resampling methods under long memory

For long-memory time series, inference based on resampling is of crucial importance, since the asymptotic distribution can often be non-Gaussian and is difficult to determine statistically. However due to the strong dependence, establishing the asymptotic validity of resampling methods is nontrivial. In this paper, we derive an efficient bound for the canonical correlation between two finite blocks of a long-memory time series. We show how this bound can be applied to establish the asymptotic consistency of subsampling procedures for general statistics under long memory. It allows the subsample size $b$ to be $o(n)$, where $n$ is the sample size, irrespective of the strength of the memory. We are then able to improve many results found in the literature. We also consider applications of subsampling procedures under long memory to the sample covariance, M-estimation and empirical processes.

preprint2016arXiv

The impact of the diagonals of polynomial forms on limit theorems with long memory

We start with an i.i.d. sequence and consider the product of two polynomial-forms moving averages based on that sequence. The coefficients of the polynomial forms are asymptotically slowly decaying homogeneous functions so that these processes have long memory. The product of these two polynomial forms is a stationary nonlinear process. Our goal is to obtain limit theorems for the normalized sums of this product process in three cases: exclusion of the diagonal terms of the polynomial form, inclusion, or the mixed case (one polynomial form excludes the diagonals while the other one includes them). In any one of these cases, if the product has long memory, then the limits are given by Wiener chaos. But the limits in each of the cases are quite different. If the diagonals are excluded, then the limit is expressed as in the product formula of two Wiener-Itô integrals. When the diagonals are included, the limit stochastic integrals are typically due to a single factor of the product, namely the one with the strongest memory. In the mixed case, the limit stochastic integral is due to the polynomial form without the diagonals irrespective of the strength of the memory.

preprint2015arXiv

Convergence of long-memory discrete $k$-th order Volterra processes

We obtain limit theorems for a class of nonlinear discrete-time processes $X(n)$ called the $k$-th order Volterra processes of order $k$. These are moving average $k$-th order polynomial forms: \[ X(n)=\sum_{0<i_1,\ldots,i_k<\infty}a(i_1,\ldots,i_k)ε_{n-i_1}\ldotsε_{n-i_k}, \] where $\{ε_i\}$ is i.i.d.\ with $\mathbb{E} ε_i=0$, $\mathbb{E} ε_i^2=1$, where $a(\cdot)$ is a nonrandom coefficient, and where the diagonals are included in the summation. We specify conditions for $X(n)$ to be well-defined in $L^2(Ω)$, and focus on central and non-central limit theorems. We show that normalized partial sums of centered $X(n)$ obey the central limit theorem if $a(\cdot)$ decays fast enough so that $X(n)$ has short memory. We prove a non-central limit theorem if, on the other hand, $a(\cdot)$ is asymptotically some slowly decaying homogeneous function so that $X(n)$ has long memory. In the non-central case the limit is a linear combination of Hermite-type processes of different orders. This linear combination can be expressed as a centered multiple Wiener-Stratonovich integral.

preprint2015arXiv

The universality of homogeneous polynomial forms and critical limits

Nourdin et al. [9] established the following universality result: if a sequence of off-diagonal homogeneous polynomial forms in i.i.d. standard normal random variables converges in distribution to a normal, then the convergence also holds if one replaces these i.i.d. standard normal random variables in the polynomial forms by any independent standardized random variables with uniformly bounded third absolute moment. The result, which was stated for polynomial forms with a finite number of terms, can be extended to allow an infinite number of terms in the polynomial forms. Based on a contraction criterion derived from this extended universality result, we prove a central limit theorem for a strongly dependent nonlinear processes, whose memory parameter lies at the boundary between short and long memory.

preprint2014arXiv

Generalized Hermite processes, discrete chaos and limit theorems

We introduce a broad class of self-similar processes $\{Z(t),t\ge 0\}$ called generalized Hermite process. They have stationary increments, are defined on a Wiener chaos with Hurst index $H\in (1/2,1)$, and include Hermite processes as a special case. They are defined through a homogeneous kernel $g$, called "generalized Hermite kernel", which replaces the product of power functions in the definition of Hermite processes. The generalized Hermite kernels $g$ can also be used to generate long-range dependent stationary sequences forming a discrete chaos process $\{X(n)\}$. In addition, we consider a fractionally-filtered version $Z^β(t)$ of $Z(t)$, which allows $H\in (0,1/2)$. Corresponding non-central limit theorems are established. We also give a multivariate limit theorem which mixes central and non-central limit theorems.

preprint2014arXiv

Structure of the third moment of the generalized Rosenblatt distribution

The Rosenblatt distribution appears as limit in non-central limit theorems. The generalized Rosenblatt distribution is obtained by allowing different power exponents in the kernel that defines the usual Rosenblatt distribution. We derive an explicit formula for its third moment, correcting the one in \citet{maejima:tudor:2012:selfsimilar} and \citet{tudor:2013:analysis}. Evaluating this formula numerically, we are able to confirm that the class of generalized Hermite processes is strictly richer than the class of Hermite processes.

preprint2013arXiv

High order chaotic limits of wavelet scalograms under long--range dependence

Let $G$ be a non--linear function of a Gaussian process $\{X_t\}_{t\in\mathbb{Z}}$ with long--range dependence. The resulting process $\{G(X_t)\}_{t\in\mathbb{Z}}$ is not Gaussian when $G$ is not linear. We consider random wavelet coefficients associated with $\{G(X_t)\}_{t\in\mathbb{Z}}$ and the corresponding wavelet scalogram which is the average of squares of wavelet coefficients over locations. We obtain the asymptotic behavior of the scalogram as the number of observations and scales tend to infinity. It is known that when $G$ is a Hermite polynomial of any order, then the limit is either the Gaussian or the Rosenblatt distribution, that is, the limit can be represented by a multiple Wiener-Itô integral of order one or two. We show, however, that there are large classes of functions $G$ which yield a higher order Hermite distribution, that is, the limit can be represented by a a multiple Wiener-Itô integral of order greater than two.

preprint2013arXiv

Multivariate limit theorems in the context of long-range dependence

We study the limit law of a vector made up of normalized sums of functions of long-range dependent stationary Gaussian series. Depending on the memory parameter of the Gaussian series and on the Hermite ranks of the functions, the resulting limit law may be (a) a multivariate Gaussian process involving dependent Brownian motion marginals, or (b) a multivariate process involving dependent Hermite processes as marginals, or (c) a combination. We treat cases (a), (b) in general and case (c) when the Hermite components involve ranks 1 and 2. We include a conjecture about case (c) when the Hermite ranks are arbitrary.

preprint2013arXiv

Multivariate limits of multilinear polynomial-form processes with long memory

We consider the multilinear polynomial-form process \[X(n)=\sum_{1\le i_1<\ldots<i_k<\infty}a_{i_1}\ldots a_{i_k}ε_{n-i_1}\ldotsε_{n-i_k},\] obtained by applying a multilinear polynomial-form filter to i.i.d.\ sequence $\{ε_i\}$ where $\{a_i\}$ is regularly varying. The resulting sequence $\{X(n)\}$ will then display either short or long memory. Now consider a vector of such X(n), whose components are defined through different $\{a_i\}$'s, that is, through different multilinear polynomial-form filters, but using the same $\{ε_i\}$. What is the limit of the normalized partial sums of the vector? We show that the resulting limit is either a) a multivariate Gaussian process with Brownian motion as marginals, or b) a multivariate Hermite process, or c) a mixture of the two. We also identify the independent components of the limit vectors.

preprint2013arXiv

Power of Change-Point Tests for Long-Range Dependent Data

We investigate the power of the CUSUM test and the Wilcoxon change-point test for a shift in the mean of a process with long-range dependent noise. We derive analytiv formulas for the power of these tests under local alternatives. These results enable us to calculate the asymptotic relative efficiency (ARE) of the CUSUM test and the Wilcoxon change point test. We obtain the surprising result that for Gaussian data, the ARE of these two tests equals 1, in contrast to the case of i.i.d. noise when the ARE is known to be $3/π$.

preprint2013arXiv

Properties and numerical evaluation of the Rosenblatt distribution

This paper studies various distributional properties of the Rosenblatt distribution. We begin by describing a technique for computing the cumulants. We then study the expansion of the Rosenblatt distribution in terms of shifted chi-squared distributions. We derive the coefficients of this expansion and use these to obtain the Lévy-Khintchine formula and derive asymptotic properties of the Lévy measure. This allows us to compute the cumulants, moments, coefficients in the chi-square expansion and the density and cumulative distribution functions of the Rosenblatt distribution with a high degree of precision. Tables are provided and software written to implement the methods described here is freely available by request from the authors.

preprint2013arXiv

The Trace Problem for Toeplitz Matrices and Operators and its Impact in Probability

The trace approximation problem for Toeplitz matrices and its applications to stationary processes dates back to the classic book by Grenander and Szegö, "Toeplitz forms and their applications". It has then been extensively studied in the literature. In this paper we provide a survey and unified treatment of the trace approximation problem both for Toeplitz matrices and for operators and describe applications to discrete- and continuous-time stationary processes. The trace approximation problem serves indeed as a tool to study many probabilistic and statistical topics for stationary models. These include central and non-central limit theorems and large deviations of Toeplitz type random quadratic functionals, parametric and nonparametric estimation, prediction of the future value based on the observed past of the process, etc. We review and summarize the known results concerning the trace approximation problem, prove some new results, and provide a number of applications to discrete- and continuous-time stationary time series models with various types of memory structures, such as long memory, anti-persistent and short memory.

preprint2013arXiv

Wavelet estimation of the long memory parameter for Hermite polynomial of Gaussian processes

We consider stationary processes with long memory which are non-Gaussian and represented as Hermite polynomials of a Gaussian process. We focus on the corresponding wavelet coefficients and study the asymptotic behavior of the sum of their squares since this sum is often used for estimating the long-memory parameter. We show that the limit is not Gaussian but can be expressed using the non-Gaussian Rosenblatt process defined as a Wiener Itô integral of order 2. This happens even if the original process is defined through a Hermite polynomial of order higher than 2.

preprint2010arXiv

Asymptotic properties of U-processes under long-range dependence

Let $(X_i)_{i\geq 1}$ be a stationary mean-zero Gaussian process with covariances $ρ(k)=\PE(X_{1}X_{k+1})$ satisfying: $ρ(0)=1$ and $ρ(k)=k^{-D} L(k)$ where $D$ is in $(0,1)$ and $L$ is slowly varying at infinity. Consider the $U$-process $\{U_n(r),\; r\in I\}$ defined as $$ U_n(r)=\frac{1}{n(n-1)}\sum_{1\leq i\neq j\leq n}\1_{\{G(X_i,X_j)\leq r\}}\; , $$ where $I$ is an interval included in $\rset$ and $G$ is a symmetric function. In this paper, we provide central and non-central limit theorems for $U_n$. They are used to derive the asymptotic behavior of the Hodges-Lehmann estimator, the Wilcoxon-signed rank statistic, the sample correlation integral and an associated scale estimator. The limiting distributions are expressed through multiple Wiener-Itô integrals.

preprint2010arXiv

Berry-Esseen and Edgeworth approximations for the tail of an infinite sum of weighted gamma random variables

Consider the sum $Z = \sum_{n=1}^\infty λ_n (η_n - \mathbb{E}η_n)$, where $η_n$ are i.i.d.~gamma random variables with shape parameter $r > 0$, and the $λ_n$'s are predetermined weights. We study the asymptotic behavior of the tail $\sum_{n=M}^\infty λ_n (η_n - \mathbb{E}η_n)$ which is asymptotically normal under certain conditions. We derive a Berry-Essen bound and Edgeworth expansions for its distribution function. We illustrate the effectiveness of these expansions on an infinite sum of weighted chi-squared distributions.

preprint2010arXiv

Distribution functions of Poisson random integrals: Analysis and computation

We want to compute the cumulative distribution function of a one-dimensional Poisson stochastic integral $I(\krnl) = \displaystyle \int_0^T \krnl(s) N(ds)$, where $N$ is a Poisson random measure with control measure $n$ and $\krnl$ is a suitable kernel function. We do so by combining a Kolmogorov-Feller equation with a finite-difference scheme. We provide the rate of convergence of our numerical scheme and illustrate our method on a number of examples. The software used to implement the procedure is available on demand and we demonstrate its use in the paper.

preprint2010arXiv

Large scale behavior of wavelet coefficients of non-linear subordinated processes with long memory

We study the asymptotic behavior of wavelet coefficients of random processes with long memory. These processes may be stationary or not and are obtained as the output of non--linear filter with Gaussian input. The wavelet coefficients that appear in the limit are random, typically non--Gaussian and belong to a Wiener chaos. They can be interpreted as wavelet coefficients of a generalized self-similar process.

preprint2010arXiv

Max-stable sketches: estimation of Lp-norms, dominance norms and point queries for non-negative signals

Max-stable random sketches can be computed efficiently on fast streaming positive data sets by using only sequential access to the data. They can be used to answer point and Lp-norm queries for the signal. There is an intriguing connection between the so-called p-stable (or sum-stable) and the max-stable sketches. Rigorous performance guarantees through error-probability estimates are derived and the algorithmic implementation is discussed.

preprint2010arXiv

Technique for computing the PDFs and CDFs of non-negative infinitely divisible random variables

We present a method for computing the PDF and CDF of a non-negative infinitely divisible random variable $X$. Our method uses the Lévy-Khintchine representation of the Laplace transform $\mathbb{E} e^{-λX} = e^{-ϕ(λ)}$, where $ϕ$ is the Laplace exponent. We apply the Post-Widder method for Laplace transform inversion combined with a sequence convergence accelerator to obtain accurate results. We demonstrate this technique on several examples including the stable distribution, mixtures thereof, and integrals with respect to non-negative Lévy processes. Software to implement this method is available from the authors and we illustrate its use at the end of the paper.

preprint2007arXiv

Bounds for the covariance of functions of infinite variance stable random variables with applications to central limit theorems and wavelet-based estimation

We establish bounds for the covariance of a large class of functions of infinite variance stable random variables, including unbounded functions such as the power function and the logarithm. These bounds involve measures of dependence between the stable variables, some of which are new. The bounds are also used to deduce the central limit theorem for unbounded functions of stable moving average time series. This result extends the earlier results of Tailen Hsing and the authors on central limit theorems for bounded functions of stable moving averages. It can be used to show asymptotic normality of wavelet-based estimators of the self-similarity parameter in fractional stable motions.

preprint2004arXiv

Integral representations of periodic and cyclic fractional stable motions

Stable non-Gaussian self-similar mixed moving averages can be decomposed into several components. Two of these are the periodic and cyclic fractional stable motions which are the subject of this study. We focus on the structure of their integral representations and show that the periodic fractional stable motions have, in fact, a canonical representation. We study several examples and discuss questions of uniqueness, namely how to determine whether two given integral representations of periodic or cyclic fractional stable motions give rise to the same process.

preprint2004arXiv

Stable stationary processes related to cyclic flows

We study stationary stable processes related to periodic and cyclic flows in the sense of Rosinski [Ann. Probab. 23 (1995) 1163-1187]. These processes are not ergodic. We provide their canonical representations, consider examples and show how to identify them among general stationary stable processes. We conclude with the unique decomposition in distribution of stationary stable processes into the sum of four major independent components: 1. A mixed moving average component. 2. A harmonizable (or ``trivial'') component. 3. A cyclic component 4. A component which is different from these.