Continuous Invertibility and Stable QML Estimation of the EGARCH(1,1) Model

I introduce the notion of continuous invertibility on a compact set for volatility models driven by a stochastic recurrence equation. I prove strong consistency of the quasi‐maximum likelihood estimator (QMLE) when the quasi‐likelihood criterion is maximized on a continuously invertible domain. This approach yields, for the first time, the asymptotic normality of the QMLE for the exponential general autoregressive conditional heteroskedastic (EGARCH(1,1)) model under explicit but non‐verifiable conditions. In practice, I propose to stabilize the QMLE by constraining the optimization procedure to an empirical continuously invertible domain. The new method, called stable QMLE, is asymptotically normal when the observations follow an invertible EGARCH(1,1) model.     

Non‐parametric Threshold Estimation for Models with Stochastic Diffusion Coefficient and Jumps

Abstract. We consider a stochastic process driven by diffusions and jumps. Given a discrete record of observations, we devise a technique for identifying the times when jumps larger than a suitably defined threshold occurred. This allows us to determine a consistent non‐parametric estimator of the integrated volatility when the infinite activity jump component is Lévy. Jump size estimation and central limit results are proved in the case of finite activity jumps. Some simulations illustrate the applicability of the methodology in finite samples and its superiority on the multipower variations especially when it is not possible to use high frequency data.     

Central Limit Theorems for the Non‐Parametric Estimation of Time‐Changed Lévy Models

Abstract. Let {Z_t}_{t 0} be a Lévy process with Lévy measure ν and let be a random clock, where g is a non‐negative function and is an ergodic diffusion independent of Z. Time‐changed Lévy models of the form are known to incorporate several important stylized features of asset prices, such as leptokurtic distributions and volatility clustering. In this article, we prove central limit theorems for a type of estimators of the integral parameter β(?):=∫?(x)ν(dx), valid when both the sampling frequency and the observation time‐horizon of the process get larger. Our results combine the long‐run ergodic properties of the diffusion process with the short‐term ergodic properties of the Lévy process Z via central limit theorems for martingale differences. The performance of the estimators are illustrated numerically for Normal Inverse Gaussian process Z and a Cox–Ingersoll–Ross process .     

Estimation of the Lognormal-Pareto Distribution Using Probability Weighted Moments and Maximum Likelihood

This article deals with the estimation of the lognormal-Pareto and the lognormal-generalized Pareto distributions, for which a general result concerning asymptotic optimality of maximum likelihood estimation cannot be proved. We develop a method based on probability weighted moments, showing that it can be applied straightforwardly to the first distribution only. In the lognormal-generalized Pareto case, we propose a mixed approach combining maximum likelihood and probability weighted moments. Extensive simulations analyze the relative efficiencies of the methods in various setups. Finally, the techniques are applied to two real datasets in the actuarial and operational risk management fields.     

An Orthogonality‐Based Estimation of Moments for Linear Mixed Models

Abstract. Estimating higher‐order moments, particularly fourth‐order moments in linear mixed models is an important, but difficult issue. In this article, an orthogonality‐based estimation of moments is proposed. Under only moment conditions, this method can easily be used to estimate the model parameters and moments, particularly those of higher order than the second order, and in the estimators the random effects and errors do not affect each other. The asymptotic normality of all the estimators is provided. Moreover, the method is readily extended to handle non‐linear, semiparametric and non‐linear models. A simulation study is carried out to examine the performance of the new method.     

Model Selection Criterion Based on the Multivariate Quasi‐Likelihood for Generalized Estimating Equations

The generalized estimating equations (GEE) approach has attracted considerable interest for the analysis of correlated response data. This paper considers the model selection criterion based on the multivariate quasi‐likelihood (MQL) in the GEE framework. The GEE approach is closely related to the MQL. We derive a necessary and sufficient condition for the uniqueness of the risk function based on the MQL by using properties of differential geometry. Furthermore, we establish a formal derivation of model selection criterion as an asymptotically unbiased estimator of the prediction risk under this condition, and we explicitly take into account the effect of estimating the correlation matrix used in the GEE procedure.     

A Spline‐Based Semiparametric Maximum Likelihood Estimation Method for the Cox Model with Interval‐Censored Data

Abstract. We propose a spline‐based semiparametric maximum likelihood approach to analysing the Cox model with interval‐censored data. With this approach, the baseline cumulative hazard function is approximated by a monotone B‐spline function. We extend the generalized Rosen algorithm to compute the maximum likelihood estimate. We show that the estimator of the regression parameter is asymptotically normal and semiparametrically efficient, although the estimator of the baseline cumulative hazard function converges at a rate slower than root‐n. We also develop an easy‐to‐implement method for consistently estimating the standard error of the estimated regression parameter, which facilitates the proposed inference procedure for the Cox model with interval‐censored data. The proposed method is evaluated by simulation studies regarding its finite sample performance and is illustrated using data from a breast cosmesis study.     

A Class of Convolution‐Based Models for Spatio‐Temporal Processes with Non‐Separable Covariance Structure

Abstract. In this article, we propose a new parametric family of models for real‐valued spatio‐temporal stochastic processes S ( x , t ) and show how low‐rank approximations can be used to overcome the computational problems that arise in fitting the proposed class of models to large datasets. Separable covariance models, in which the spatio‐temporal covariance function of S ( x , t ) factorizes into a product of purely spatial and purely temporal functions, are often used as a convenient working assumption but are too inflexible to cover the range of covariance structures encountered in applications. We define positive and negative non‐separability and show that in our proposed family we can capture positive, zero and negative non‐separability by varying the value of a single parameter.     

Autocovariance Estimation in Regression with a Discontinuous Signal and m‐Dependent Errors: A Difference‐Based Approach

We discuss a class of difference‐based estimators for the autocovariance in nonparametric regression when the signal is discontinuous and the errors form a stationary m‐dependent process. These estimators circumvent the particularly challenging task of pre‐estimating such an unknown regression function. We provide finite‐sample expressions of their mean squared errors for piecewise constant signals and Gaussian errors. Based on this, we derive biased‐optimized estimates that do not depend on the unknown autocovariance structure. Notably, for positively correlated errors, that part of the variance of our estimators that depend on the signal is minimal as well. Further, we provide sufficient conditions for ‐consistency; this result is extended to piecewise Hölder regression with non‐Gaussian errors. We combine our biased‐optimized autocovariance estimates with a projection‐based approach and derive covariance matrix estimates, a method that is of independent interest. An R package, several simulations and an application to biophysical measurements complement this paper.     

存在自相关时检验和估计方法的改进——基于自回归分布滞后模型的自相关检验和参数估计

存在自相关时的自相关检验和参数估计是基础计量经济学的一个重要内容,并且存在自相关时的原模型已转化为自回归分布滞后模型。讨论存在自相关时的自相关检验和参数估计问题,提出了一种基于自回归分布滞后模型的自相关检验法,并同时给出了相应的参数估计。     

Spatio‐temporal Ornstein–Uhlenbeck Processes: Theory,Simulation and Statistical Inference

Spatio‐temporal modelling is an increasingly popular topic in Statistics. Our paper contributes to this line of research by developing the theory, simulation and inference for a spatio‐temporal Ornstein–Uhlenbeck process. We conduct detailed simulation studies and demonstrate the practical relevance of these processes in an empirical study of radiation anomaly data. Finally, we describe how predictions can be carried out in the Gaussian setting.     

Moments and generating functions for the absortion distribution and its negative binomial analogue

The moments of the absorption are difficult to obtain. The generating functions are basic hypergeometric functions. This paper shows how to define two shift operators that allow elementary arguments to be used to develop recursions for the expected values of general functions. The exact moments of the distribution follow. The generating function for the negative binomial analogue gives the moments directly.     

Empirical Likelihood Inference for the Rao‐Hartley‐Cochran Sampling Design

The Hartley‐Rao‐Cochran sampling design is an unequal probability sampling design which can be used to select samples from finite populations. We propose to adjust the empirical likelihood approach for the Hartley‐Rao‐Cochran sampling design. The approach proposed intrinsically incorporates sampling weights, auxiliary information and allows for large sampling fractions. It can be used to construct confidence intervals. In a simulation study, we show that the coverage may be better for the empirical likelihood confidence interval than for standard confidence intervals based on variance estimates. The approach proposed is simple to implement and less computer intensive than bootstrap. The confidence interval proposed does not rely on re‐sampling, linearization, variance estimation, design‐effects or joint inclusion probabilities.     

An Investigation of Parameter Recovery in MCMC Estimation for the Additive IRT Model

The article aims at evaluating the parameter recovery for the multidimensional additive IRT model (Sheng, 2005 Sheng , Y. ( 2005 ). Bayesian Analysis of Hierarchical IRT Models: Comparing and Combining the Unidimensional & Multi-Unidimensional IRT Models. PhD thesis , Faculty of the Graduate SchoolUniversity of Missouri-Columbia . [Google Scholar]; Sheng and Wikle, 2009 Sheng , Y. , Wikle , C. K. ( 2009 ). Bayesian IRT models incorporating general and specific abilities . Behaviormetrika 36 : 27 – 48 .[Crossref] , [Google Scholar]). By estimating the model parameters via Gibbs sampler, a simulation study is conducted under different testing conditions, e.g., dimensionality, test and subtest lengths, correlation matrices, and different values of discrimination parameters. The results show that, especially when the test length is short and the abilities are highly correlated, the accuracy of the parameter estimates is reduced and more iterations are required to convergence. An application in educational testing is also described to show the effectiveness of the model in use.     

Score estimation in the monotone single‐index model

We consider estimation in the single‐index model where the link function is monotone. For this model, a profile least‐squares estimator has been proposed to estimate the unknown link function and index. Although it is natural to propose this procedure, it is still unknown whether it produces index estimates that converge at the parametric rate. We show that this holds if we solve a score equation corresponding to this least‐squares problem. Using a Lagrangian formulation, we show how one can solve this score equation without any reparametrization. This makes it easy to solve the score equations in high dimensions. We also compare our method with the effective dimension reduction and the penalized least‐squares estimator methods, both available on CRAN as R packages, and compare with link‐free methods, where the covariates are elliptically symmetric.     

A non‐stationary bivariate INAR(1) process with a simple cross‐dependence: Estimation with some properties

This paper considers modelling of a non‐stationary bivariate integer‐valued autoregressive process of order 1 (BINAR(1)) where the cross‐dependence between the counting series is formed through the relationship of the current series with the previous‐lagged count series observations while the pair of innovations is independent and marginally Poisson. In addition, this paper proposes a generalised quasi‐likelihood (GQL) estimating equation based on the exact specification of the mean score and the auto‐covariance structure. The proposed approach is also compared with other popular techniques such as conditional maximum likelihood (CML), generalised least squares (GLS) and generalised method of moment (GMM) based on simulated data from the proposed BINAR(1). Moreover, the model is applied to weekly series of day and night road accidents arising in some regions of Mauritius and is compared with other existing BINAR(1) models.     

Estimation of the parameters of symmetric stable ARMA and ARMA–GARCH models

In this article, we first propose the modified Hannan–Rissanen Method for estimating the parameters of autoregressive moving average (ARMA) process with symmetric stable noise and symmetric stable generalized autoregressive conditional heteroskedastic (GARCH) noise. Next, we propose the modified empirical characteristic function method for the estimation of GARCH parameters with symmetric stable noise. Further, we show the efficiency, accuracy and simplicity of our methods with Monte-Carlo simulation. Finally, we apply our proposed methods to model the financial data.     

Composite Estimating Equation Method for the Accelerated Failure Time Model with Length‐biased Sampling Data

Length‐biased sampling data are often encountered in the studies of economics, industrial reliability, epidemiology, genetics and cancer screening. The complication of this type of data is due to the fact that the observed lifetimes suffer from left truncation and right censoring, where the left truncation variable has a uniform distribution. In the Cox proportional hazards model, Huang & Qin (Journal of the American Statistical Association, 107, 2012, p. 107) proposed a composite partial likelihood method which not only has the simplicity of the popular partial likelihood estimator, but also can be easily performed by the standard statistical software. The accelerated failure time model has become a useful alternative to the Cox proportional hazards model. In this paper, by using the composite partial likelihood technique, we study this model with length‐biased sampling data. The proposed method has a very simple form and is robust when the assumption that the censoring time is independent of the covariate is violated. To ease the difficulty of calculations when solving the non‐smooth estimating equation, we use a kernel smoothed estimation method (Heller; Journal of the American Statistical Association, 102, 2007, p. 552). Large sample results and a re‐sampling method for the variance estimation are discussed. Some simulation studies are conducted to compare the performance of the proposed method with other existing methods. A real data set is used for illustration.     

Estimation in an empirical bayes model for longitudinal and cross‐sectionally clustered binary data

Some studies generate data that can be grouped into clusters in more than one way. Consider for instance a smoking prevention study in which responses on smoking status are collected over several years in a cohort of students from a number of different schools. This yields longitudinal data, also cross‐sectionaliy clustered in schools. The authors present a model for analyzing binary data of this type, combining generalized estimating equations and estimation of random effects to address the longitudinal and cross‐sectional dependence, respectively. The estimation procedure for this model is discussed, as are the results of a simulation study used to investigate the properties of its estimates. An illustration using data from a smoking prevention trial is given.