非线性计量经济模型的非参数估计方法研究

本文提出使用核估计的方法构造平滑转移模型(STR)的非参数模拟最大似然估计(NPSML),给出了NPSML估计量的构造方法、渐近性质以及相应的核函数和窗宽的选择准则,并利用滑动窗宽算法对估计量的构造过程进行了改进。通过Monte Carlo实验证明,该方法是可靠的,并且当误差项存在序列相关时,此种估计量是稳健的。     

Combining kernel estimators in the uniform deconvolution problem

We construct a density estimator and an estimator of the distribution function in the uniform deconvolution model. The estimators are based on inversion formulas and kernel estimators of the density of the observations and its derivative. Initially the inversions yield two different estimators of the density and two estimators of the distribution function. We construct asymptotically optimal convex combinations of these two estimators. We also derive pointwise asymptotic normality of the resulting estimators, the pointwise asymptotic biases and an expansion of the mean integrated squared error of the density estimator. It turns out that the pointwise limit distribution of the density estimator is the same as the pointwise limit distribution of the density estimator introduced by Groeneboom and Jongbloed (Neerlandica, 57, 2003, 136), a kernel smoothed nonparametric maximum likelihood estimator of the distribution function.     

MCMC maximum likelihood for latent state models

This paper develops a pure simulation-based approach for computing maximum likelihood estimates in latent state variable models using Markov Chain Monte Carlo methods (MCMC). Our MCMC algorithm simultaneously evaluates and optimizes the likelihood function without resorting to gradient methods. The approach relies on data augmentation, with insights similar to simulated annealing and evolutionary Monte Carlo algorithms. We prove a limit theorem in the degree of data augmentation and use this to provide standard errors and convergence diagnostics. The resulting estimator inherits the sampling asymptotic properties of maximum likelihood. We demonstrate the approach on two latent state models central to financial econometrics: a stochastic volatility and a multivariate jump-diffusion models. We find that convergence to the MLE is fast, requiring only a small degree of augmentation.     

On finite sample properties of alternative estimators of coefficients in a structural equation with many instruments

We compare four different estimation methods for the coefficients of a linear structural equation with instrumental variables. As the classical methods we consider the limited information maximum likelihood (LIML) estimator and the two-stage least squares (TSLS) estimator, and as the semi-parametric estimation methods we consider the maximum empirical likelihood (MEL) estimator and the generalized method of moments (GMM) (or the estimating equation) estimator. Tables and figures of the distribution functions of four estimators are given for enough values of the parameters to cover most linear models of interest and we include some heteroscedastic cases and nonlinear cases. We have found that the LIML estimator has good performance in terms of the bounded loss functions and probabilities when the number of instruments is large, that is, the micro-econometric models with “many instruments” in the terminology of recent econometric literature.     

Density estimation in the uniform deconvolution model

We consider the problem of estimating a probability density function based on data that are corrupted by noise from a uniform distribution. The (nonparametric) maximum likelihood estimator for the corresponding distribution function is well defined. For the density function this is not the case. We study two nonparametric estimators for this density. The first is a type of kernel density estimate based on the empirical distribution function of the observable data. The second is a kernel density estimate based on the MLE of the distribution function of the unobservable (uncorrupted) data.     

Strong Consistency of the Maximum Likelihood Estimator in Generalized Linear and Nonlinear Mixed-Effects Models

Generalized linear and nonlinear mixed-effects models are used extensively in biomedical, social, and agricultural sciences. The statistical analysis of these models is based on the asymptotic properties of the maximum likelihood estimator. However, it is usually assumed that the maximum likelihood estimator is consistent, without providing a proof. A rigorous proof of the consistency by verifying conditions from existing results can be very difficult due to the integrated likelihood. In this paper, we present some easily verifiable conditions for the strong consistency of the maximum likelihood estimator in generalized linear and nonlinear mixed-effects models. Based on this result, we prove that the maximum likelihood estimator is consistent for some frequently used models such as mixed-effects logistic regression models and growth curve models.     

Finite sample properties of maximum likelihood estimator in spatial models

We investigate the finite sample properties of the maximum likelihood estimator for the spatial autoregressive model. A stochastic expansion of the score function is used to develop the second-order bias and mean squared error of the maximum likelihood estimator. We show that the results can be expressed in terms of the expectations of cross products of quadratic forms, or ratios of quadratic forms in a normal vector which can be evaluated using the top order invariant polynomial. Our numerical calculations demonstrate that the second-order behaviors of the maximum likelihood estimator depend on the degree of sparseness of the weights matrix.     

Pseudo-maximum likelihood estimation in two classes of semiparametric diffusion models

Two classes of semiparametric diffusion models are considered, where either the drift or the diffusion term is parameterized, while the other term is left unspecified. We propose a pseudo-maximum likelihood estimator (PMLE) of the parametric component that maximizes the likelihood with a preliminary estimator of the unspecified term plugged in. It is demonstrated how models and estimators can be used in a two-step specification testing strategy of semiparametric and fully parametric models, and shown that approximate/simulated versions of the PMLE inherit the properties of the actual but infeasible estimator. A simulation study investigates the finite sample performance of the PMLE.     

On a two-step estimation of a multivariate logit model

In this article the author studies the properties of the two-step estimation method proposed by Domencich and McFadden (Urban Travel Demand, North-Holland, 1975) for a multivariate logit model and shows that it is consistent but asymptotically less efficient than the maximum likelihood estimator. Its computation, however, can be considerably simpler than that of the maximum likelihood estimator, especially in models involving several dependent variables.     

Bias-corrected smoothed score function for single-index models

We in this paper investigate smoothed score function based confidence regions for parameters in single-index models. Because a plug-in estimator of nonparametric link function causes the bias of smoothed score function to be non-negligible, the limit of the score function is asymptotically normal with a non-zero mean due to the slow convergence rate of nonparametric estimation. A bias-corrected smoothed score function is recommended for achieving centered normal limit without under-smoothing or high order kernel, and then the confidence region can be constructed by chi-square distribution. Simulation studies are carried out to assess the performance of bias-corrected local likelihood, and to compare with normal approximation approach.     

On minimum distance estimators for unimodal densities

It is shown that minimum distance estimators for families of unimodal densities are always consistent; the rate of convergence is indicated. An algorithm is proposed for computing the minimum distance estimator for the family of all unimodal densities. References are given to the maximum likelihood method and the kernel method.     

Estimation of affine term structure models with spanned or unspanned stochastic volatility

We develop new procedures for maximum likelihood estimation of affine term structure models with spanned or unspanned stochastic volatility. Our approach uses linear regression to reduce the dimension of the numerical optimization problem yet it produces the same estimator as maximizing the likelihood. It improves the numerical behavior of estimation by eliminating parameters from the objective function that cause problems for conventional methods. We find that spanned models capture the cross-section of yields well but not volatility while unspanned models fit volatility at the expense of fitting the cross-section.     

Comovements and asymmetric tail dependence in state housing prices in the USA: A nonparametric approach

We reexamine the methods used in estimating comovements among US regional home prices and find that there are insufficient moments to ensure a normal limit necessary for employing the quasi‐maximum likelihood estimator. Hence we propose applying the self‐weighted quasi‐maximum exponential likelihood estimator and a bootstrap method to test and account for the asymmetry of comovements as well as different magnitudes across state pairs. Our results reveal interstate asymmetric tail dependence based on observed house price indices rather than residuals from fitting autoregressive–generalized autoregressive conditional heteroskedasticity (AR‐GARCH) models.     

Semiparametric maximum likelihood estimation in Cox proportional hazards model with covariate measurement errors

This paper studies semiparametric maximum likelihood estimators in the Cox proportional hazards model with covariate error, assuming that the conditional distribution of the true covariate given the surrogate is known. We show that the estimator of the regression coefficient is asymptotically normal and efficient, its covariance matrix can be estimated consistently by differentiation of the profile likelihood, and the likelihood ratio test is asymptotically chi-squared. We also provide efficient algorithms for the computations of the semiparametric maximum likelihood estimate and the profile likelihood. The performance of this method is successfully demonstrated in simulation studies.     

Analysis of Divergence in Loglinear Models When Expected Frequencies are Subject to Linear Constraints

Consider the loglinear model for categorical data under the assumption of multinomial sampling. We are interested in testing between various hypotheses on the parameter space when we have some hypotheses relating to the parameters of the models that can be written in terms of constraints on the frequencies. The usual likelihood ratio test, with maximum likelihood estimator for the unspecified parameters, is generalized to tests based on -divergence statistics, using minimum -divergence estimator. These tests yield the classical likelihood ratio test as a special case. Asymptotic distributions for the new -divergence test statistics are derived under the null hypothesis.     

On the inefficiency of the restricted maximum likelihood

The restricted maximum likelihood is preferred by many to the full maximum likelihood for estimation with variance component and other random coefficient models, because the variance estimator is unbiased. It is shown that this unbiasedness is accompanied in some balanced designs by an inflation of the mean squared error. An estimator of the cluster‐level variance that is uniformly more efficient than the full maximum likelihood is derived. Estimators of the variance ratio are also studied.     

Semiparametric efficient adaptive estimation of asymmetric GARCH models

In this paper we derive a semiparametric efficient adaptive estimator of an asymmetric GARCH model. Applying some general results from Drost et al. [1997. The Annals of Statistics 25, 786–818], we first estimate the unknown density function of the disturbances by kernel methods, then apply a one-step Newton–Raphson method to obtain a more efficient estimator than the quasi-maximum likelihood estimator. The proposed semiparametric estimator is adaptive for parameters appearing in the conditional standard deviation model with respect to the unknown distribution of the disturbances.     

Practical Problems with Reduced‐rank ML Estimators for Cointegration Parameters and a Simple Alternative*

Johansen's reduced‐rank maximum likelihood (ML) estimator for cointegration parameters in vector error correction models is known to produce occasional extreme outliers. Using a small monetary system and German data we illustrate the practical importance of this problem. We also consider an alternative generalized least squares (GLS) system estimator which has better properties in this respect. The two estimators are compared in a small simulation study. It is found that the GLS estimator can indeed be an attractive alternative to ML estimation of cointegration parameters.     

The Estimation of Pareto Distribution by a Weighted Least Square Method

The two-parameter Pareto distribution provides reasonably good fit to the distributions of income and property value, and explains many empirical phenomena. For the censored data, the two parameters are regularly estimated by the maximum likelihood estimator, which is complicated in computation process. This investigation proposes a weighted least square estimator to estimate the parameters. Such a method is comparatively concise and easy to perceive, and could be applied to either complete or truncated data. Simulation studies are conducted in this investigation to show the feasibility of the proposed method. This report will demonstrate that the weighted least square estimator gives better performance than unweighted least square estimators with simulation cases. We also illustrate that the weighted least square estimator is very close to maximum likelihood estimator with simulation studies.     

Beta spatial linear mixed model with variable dispersion using Monte Carlo maximum likelihood

We propose a beta spatial linear mixed model with variable dispersion using Monte Carlo maximum likelihood. The proposed method is useful for those situations where the response variable is a rate or a proportion. An approach to the spatial generalized linear mixed models using the Box–Cox transformation in the precision model is presented. Thus, the parameter optimization process is developed for both the spatial mean model and the spatial variable dispersion model. All the parameters are estimated using Markov chain Monte Carlo maximum likelihood. Statistical inference over the parameters is performed using approximations obtained from the asymptotic normality of the maximum likelihood estimator. Diagnosis and prediction of a new observation are also developed. The method is illustrated with the analysis of one simulated case and two studies: clay and magnesium contents. In the clay study, 147 soil profile observations were taken from the research area of the Tropenbos Cameroon Programme, with explanatory variables: elevation in metres above sea level, agro‐ecological zone, reference soil group and land cover type. In the magnesium content, the soil samples were taken from 0‐ to 20‐cm‐depth layer at each of the 178 locations, and the response variable is related to the spatial locations, altitude and sub‐region.