Analytical quasi maximum likelihood inference in multivariate volatility models

Quasi maximum likelihood estimation and inference in multivariate volatility models remains a challenging computational task if, for example, the dimension of the parameter space is high. One of the reasons is that typically numerical procedures are used to compute the score and the Hessian, and often they are numerically unstable. We provide analytical formulae for the score and the Hessian for a variety of multivariate GARCH models including the Vec and BEKK specifications as well as the recent dynamic conditional correlation model. By means of a Monte Carlo investigation of the BEKK–GARCH model we illustrate that employing analytical derivatives for inference is clearly preferable to numerical methods.     

Pseudo maximum likelihood estimation of spatial autoregressive models with increasing dimension

Pseudo maximum likelihood estimates are developed for higher-order spatial autoregressive models with increasingly many parameters, including models with spatial lags in the dependent variables both with and without a linear or nonlinear regression component, and regression models with spatial autoregressive disturbances. Consistency and asymptotic normality of the estimates are established. Monte Carlo experiments examine finite-sample behaviour.     

A local maximum likelihood estimator for Poisson regression

A local maximum likelihood estimator based on Poisson regression is presented as well as its bias, variance and asymptotic distribution. This semiparametric estimator is intended to be an alternative to the Poisson, negative binomial and zero-inflated Poisson regression models that does not depend on regularity conditions and model specification accuracy. Some simulation results are presented. The use of the local maximum likelihood procedure is illustrated on one example from the literature. This procedure is found to perform well. This research was partially supported by Calouste Gulbenkian Foundation and PRODEP III.     

On maximum likelihood prediction based on Type II doubly censored exponential data

In this paper, the maximum likelihood predictor (MLP) of the kth ordered observation, t _k, in a sample of size n from a two-parameter exponential distribution as well as the predictive maximum likelihood estimators (PMLE's) of the location and scale parameters, θ and β, based on the observed values t _r, …, t _s (1≤r≤s<k≤n), are obtained in closed forms, contrary to the belief they cannot be so expressed. When θ is known, however, the PMLE of β and MLP of t _k do not admit explicit expressions. It is shown here that they exist and are unique; sharp lower and upper bounds are also provided. The derived predictors and estimators are reasonable and also have good asymptotic properties. As applications, the total duration time in a life test and the failure time of a k-out-of-n system may be predicted. Finally, an illustrative example is included. Received: August 1999     

The efficiency of Spanish arable crop organic farms,a local maximum likelihood approach

The methodologies that have been used in existing research to assess the efficiency with which organic farms are operating are generally based either on the stochastic frontier methodology or on a deterministic non-parametric approach. Recently, Kumbhakar et al. (J Econom 137:1–27, 2007) proposed a new nonparametric, stochastic method based on the local maximum likelihood principle. We use this methodology to compare the efficiency ratings of organic and conventional arable crop farms in the Spanish region of Andalucía. Nonparametrically encompassing the stochastic frontier model is especially useful when comparing the performance of two groups that are likely to be characterized by different production technologies.

Estimation of multiple poisson means: A compromise between maximum likelihood and empirical bayes estimators

For estimatingp(⩾ 2) independent Poisson means, the paper considers a compromise between maximum likelihood and empirical Bayes estimators. Such compromise estimators enjoy both good componentwise as well as ensemble properties. Research supported by the NSF Grant Number MCS-8218091.     

Maximum likelihood methods for fitting the burr type XII distribution to multiply (progressively) censored life test data

Summary This paper develops mathematical and computational methodology for fitting, by the method of maximum likelihood (ML), the Burr Type XII distribution to multiply (or progressively) censored life test data. Mathematical expressions are given for approximating the asymptotic variances and covariances of the ML estimates (MLEs) of the parameters of the Burr Type XII distribution. A rigorous mathematical analysis is undertaken to investigate the existence and uniqueness of the MLEs for arbitrary sample data. The methodology of this paper is applied to progressively censored sample data arising in a life test experiment.     

Nonparametric maximum likelihood estimation in a non locally compact setting

Maximum likelihood estimation is considered in the context of infinite dimensional parameter spaces. It is shown that in some locally convex parameter spaces sequential compactness of the bounded sets ensures the existence of minimizers of objective functions and the consistency of maximum likelihood estimators in an appropriate topology. The theory is applied to revisit some classical problems of nonparametric maximum likelihood estimation, to study location parameters in Banach spaces, and finally to obtain Varadarajan’s theorem on the convergence of empirical measures in the form of a consistency result for a sequence of maximum likelihood estimators. Several parameter spaces sharing the crucial compactness property are identified. This research was supported by grants from the National Sciences and Engineering Research Council of Canada and the Fonds FCAR de la Province de Québec.     

Robust Likelihood Methods Based on the Skew-t and Related Distributions

The robustness problem is tackled by adopting a parametric class of distributions flexible enough to match the behaviour of the observed data. In a variety of practical cases, one reasonable option is to consider distributions which include parameters to regulate their skewness and kurtosis. As a specific representative of this approach, the skew‐t distribution is explored in more detail and reasons are given to adopt this option as a sensible general‐purpose compromise between robustness and simplicity, both of treatment and of interpretation of the outcome. Some theoretical arguments, outcomes of a few simulation experiments and various wide‐ranging examples with real data are provided in support of the claim.     

The empirical relevance of the competitive storage model

The empirical relevance of models of competitive storage arbitrage in explaining commodity price behavior has been seriously challenged in a series of pathbreaking papers by [Deaton and Laroque, 1992], [Deaton and Laroque, 1995] and [Deaton and Laroque, 1996]. Here we address their major criticism, that the model is in general unable to explain the degree of serial correlation observed in the prices of twelve major commodities. First, we present a simple numerical version of their model which, contrary to Deaton and Laroque (1992), can generate the high levels of serial correlation observed in commodity prices, if it is parameterized to generate realistic levels of price variation. Then, after estimating the [Deaton and Laroque, 1995] and [Deaton and Laroque, 1996] model using their data set, model specification and econometric approach, we show that the use of a much finer grid to approximate the equilibrium price function yields quite different estimates for most commodities. Results are obtained for coffee, copper, jute, maize, palm oil, sugar and tin that support the specifications of the storage model with positive constant marginal storage cost and no deterioration as in Gustafson (1958a). Consumption demand has a low response to price and, except for sugar, stockouts are infrequent. The observed magnitudes of serial correlation of price match those implied by the estimated model.     

A note on multinomial maximum likelihood estimation under ordered restrictions and the EM algorithm

A frequently occurring problem is to find the maximum likelihood estimation (MLE) of p subject to p∈C (C⊂ P the probability vectors in R ^k ). The problem has been discussed by many authors and they mainly focused when p is restricted by linear constraints or log-linear constraints. In this paper, we construct the relationship between the the maximum likelihood estimation of p restricted by p∈C and EM algorithm and demonstrate that the maximum likelihood estimator can be computed through the EM algorithm (Dempster et al. in J R Stat Soc Ser B 39:1–38, 1997). Several examples are analyzed by the proposed method.     

Pseudo empirical likelihood method in the presence of missing data

In this paper, we propose an estimator for the population mean when some observations on the study and auxiliary variables are missing from the sample. The proposed estimator is valid for any unequal probability sampling design, and is based upon the pseudo empirical likelihood method. The proposed estimator is compared with other estimators in a simulation study.     

THE BIVARIATE EXPONENTIAL DISTRIBUTION; A REVIEW AND SOME NEW RESULTS

Various models have been proposed as bivariate forms of the exponential distribution. A brief but comprehensive review is presented which classifies, interrelates and contrasts the different models and outlines what is known about distributional properties, applicability and estimation and testing of parameters (particularly the association parameter). Some new results are presented for one particular model. Maximum likelihood, and moment–type, estimators of the association parameter are examined. Asymptotic variances are derived and attention is given to the relative efficiency of the estimators and to problems of their evaluation.     

Smooth estimators for estimating order restricted scale parameters of two gamma distributions

We consider the problem of component-wise estimation of ordered scale parameters of two gamma populations, when it is known apriori which population corresponds to each ordered parameter. Under the scale equivariant squared error loss function, smooth estimators that improve upon the best scale equivariant estimators are derived. These smooth estimators are shown to be generalized Bayes with respect to a non-informative prior. Finally, using Monte Carlo simulations, these improved smooth estimators are compared with the best scale equivariant estimators, their non-smooth improvements obtained in Vijayasree, Misra & Singh (1995), and the restricted maximum likelihood estimators. Acknowledgments. Authors are thankful to a referee for suggestions leading to improved presentation.     

基于定数截尾样本的指数分布的环境因子估计

讨论了定数截尾样本下,指数分布环境因子的极大似然估计和区间估计,为研究估计的精度,运用随机模拟方法,对环境因子的置信区间的精度进行了讨论。     

Small sample properties of maximum likelihood versus generalized method of moments based tests for spatially autocorrelated errors

Many applied researchers have to deal with spatially autocorrelated residuals (SAR). Available tests that identify spatial spillovers as captured by a significant SAR parameter, are either based on maximum likelihood (MLE) or generalized method of moments (GMM) estimates. This paper illustrates the properties of various tests for the null hypothesis of a zero SAR parameter in a comprehensive Monte Carlo study. The main finding is that Wald tests generally perform well regarding both size and power even in small samples. The GMM-based Wald test is correctly sized even for non-normally distributed disturbances and small samples, and it exhibits a similar power as its MLE-based counterpart. Hence, for the applied researcher the GMM Wald test can be recommended, because it is easy to implement.     

Analysis of Categorical Data for Complex Surveys

This paper reviews methods for handling complex sampling schemes when analysing categorical survey data. It is generally assumed that the complex sampling scheme does not affect the specification of the parameters of interest, only the methodology for making inference about these parameters. The organisation of the paper is loosely chronological. Contingency table data are emphasised first before moving on to the analysis of unit‐level data. Weighted least squares methods, introduced in the mid 1970s along with methods for two‐way tables, receive early attention. They are followed by more general methods based on maximum likelihood, particularly pseudo maximum likelihood estimation. Point estimation methods typically involve the use of survey weights in some way. Variance estimation methods are described in broad terms. There is a particular emphasis on methods of testing. The main modelling methods considered are log‐linear models, logit models, generalised linear models and latent variable models. There is no coverage of multilevel models.     

Multiple-dose design and bias-reducing methods for limiting dilution assays

This paper gives an overview of several (mostly recent) statistical contributions to the theory of Limiting and Serial Dilution Assays (LDA's, SDA's). A simple and useful method is presented for the setup of a design for an LDA or an SDA. This method is based on several user-supplied design parameters, consisting in the researcher's advance information and other parameters inherent to the particular problem. The commonly used Maximum Likelihood (ML) and Minimum Chi-square methods for the estimation of the unknown parameter in an LDA or an SDA are described and compared to several bias-reducing estimation methods, e.g. jackknife and bootstrap versions of the ML method. One particular jackknife version is recommended.     

Convenient links for the estimation of hedonic price indexes: the case of unique,infrequently traded assets

Hedonic methods are a prominent approach in the construction of quality‐adjusted price indexes. This paper shows that the process of computing such indexes is substantially simplified if arithmetic (geometric) price indexes are computed based on exponential (log‐linear) hedonic functions estimated by the Poisson pseudo‐maximum likelihood (ordinary least squares) method. A Monte Carlo simulation study based on housing data illustrates the convenience of the links identified and the very attractive properties of the Poisson estimator in the hedonic framework.