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

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

逐次定数截尾下Weibull分布参数的逆矩估计

基于逐次定数截尾样本下,讨论了两参数Weibull分布的参数估计,得到了两参数的逆矩估计.并利用模拟方法与极大似然估计作比较,模拟结果表明逆矩估计优于极大似然估计。     

非参数核密度估计在异方差模型中的应用

针对非参数核密度估计中最优窗宽的选择在实际建模中的不足,提出了一个新的最优窗宽选择的迭代方法,克服了使用传统的经验法则所带来的局限性。并在此基础上用一种新的非参数核密度估计ML方法应用到了中国股票市场,通过与极大似然估计对比论证了此方法的有效性和可行性。实证分析表明,通过与实际值的模拟对比,运用非参数估计技术得到上证指数日收益率的拟合值要优于极大似然估计的拟合值。     

逐步增加的Ⅱ型截尾下双参数指数分布的可靠性分析

基于逐步增加的Ⅱ型截尾模型,讨论了双参数指数分布的可靠性指标的估计。导出了形状参数、尺度参数及可靠度函数的极大似然估计(MLE)和Bayes估计,最后运用Monte Carlo方法对Bayes估计和极大似然估计的MSE,进行了模拟比较。     

浅谈煤层瓦斯吸附常数分布参数的极大似然估计

极大似然估计是估计的另一种计算方法,最早由高斯先生提出,后来由英国的统计学家费歇先生进行了命名和定义。此方法得到了广泛的应用,当前实验室研究的重要课题内容是在极大似然原理的基础上,采取随机试验的方法对结果进行相关的数据统计和分析。文章利用极大似然这种估计方法,对煤层瓦斯吸附常数分布参数进行了研究。     

Logistic回归的ArctanLASSO惩罚似然估计及应用

Logistic回归是计量经济学中应用最广的离散选择模型。当变量个数较多时,极大似然估计解释性较差,为此本文基于新的惩罚函数ArctanLASSO,给出Logistic回归的一种非凸惩罚似然估计进行参数估计和变量选取,并证明了估计量的n1/2相合性和Oracle性质。本文结合二阶近似处理、LLA方法和梯度下降法给出估计算法,并通过最小化BIC准则对正则化参数进行选取。模拟数据分析显示,当样本量较大时,该方法在参数估计和变量选取两个方面都优于传统的LASSO、SCAD和MCP方法,样本量较小时,该方法同样具有很大优势。实际数据分析表明,该方法很好地权衡了拟合程度和非零系数的选择,是最优的备选模型,具有重要的实际意义。     

半参数变系数空间滞后模型的截面极大似然估计

在线性参数空间滞后模型中,解释变量的系数一般假设为固定常数,本文首先放松了这种假设,将解释变量的系数设定为某一变量的未知函数,提出一类全新的半参数变系数空间滞后模型;其次导出了该模型的截面极大似然估计,并证明了该估计的一致性;最后用蒙特卡洛数值模拟方法考察了该估计在小样本条件下的性质,数值模拟结果显示我们提出的估计方法在小样本条件下依然有优良的表现。     

空间动态面板模型拟极大似然估计的渐近效率改进

Lee和Yu(2008)研究了一类同时带个体与时间固定效应的空间动态面板模型的拟极大似然估计量的大样本性质.本文说明当扰动项非正态时,拟极大似然估计量的渐近效率可以被进一步提高.为此,我们构造了一组合待定矩阵且形式一般的矩条件以用来包含对数似然函数一阶条件的特殊形式.从无冗余矩条件的角度,选取最优待定矩阵得到了最佳广义矩估计量.本文证明了当扰动项正态分布时,最佳广义矩估计量和拟极大似然估计量渐近等价;当扰动项非正态分布时,广义矩估计量具有比拟极大似然估计量更高的渐近效率.Monte Carlo实验结果与本文的理论预期一致.     

利用HMM算法确定重组自交系基因型

重组自交系是杂交一代经过连续自交而获得的自交系群体。隐形马尔科夫模型(HMM)是一种极大似然估计算法,在很多生物信息研究中取得了理想的结果。文章提供了利用一阶隐型马尔科夫模型来确定重组自交系个体基因型的方法,并论述了其准确性和时效性。该算法在模拟数据及小鼠基因型数据上取得了理想的结果。     

估计MS-DSGE模型脉冲响应函数的EM算法及其应用

研究目标：考察不同区制下外生冲击对中国宏观经济的非对称性效应。研究方法：引入两状态的Markov区制转换过程建立MS-DSGE模型,并基于MS-DSGE模型的Markov区制转换动态因子模型的表示提出了估计MS-DSGE模型脉冲响应函数的极大似然估计EM算法。研究发现：本文提出的估计方法具有良好的有限样本性质和收敛性,参数估计量具有渐近正态分布。实证分析发现,应持续施行扩张性政策以刺激经济稳定增长,对冲挤占效应以及稳定物价水平。尤其,当经济处于“衰退”区制时,政府应实施及时有效的调控政策刺激经济运行区制的转移。研究创新：与Bayesian分析方法比较,本文提出的估计方法避免了对数线性化MS-DSGE模型的随机奇异性以及对先验分布的设定和观测变量选取的非稳健性。研究价值：提出了一种估计MS-DSGE模型脉冲响应函数的方法。     

Mean squared error of ridge estimators in logistic regression

It is well known that the maximum likelihood estimator (MLE) is inadmissible when estimating the multidimensional Gaussian location parameter. We show that the verdict is much more subtle for the binary location parameter. We consider this problem in a regression framework by considering a ridge logistic regression (RR) with three alternative ways of shrinking the estimates of the event probabilities. While it is shown that all three variants reduce the mean squared error (MSE) of the MLE, there is at the same time, for every amount of shrinkage, a true value of the location parameter for which we are overshrinking, thus implying the minimaxity of the MLE in this family of estimators. Little shrinkage also always reduces the MSE of individual predictions for all three RR estimators; however, only the naive estimator that shrinks toward 1/2 retains this property for any generalized MSE (GMSE). In contrast, for the two RR estimators that shrink toward the common mean probability, there is always a GMSE for which even a minute amount of shrinkage increases the error. These theoretical results are illustrated on a numerical example. The estimators are also applied to a real data set, and practical implications of our results are discussed.     

Asymptotic distribution of the cointegrating vector estimator in error correction models with conditional heteroskedasticity

This paper explores the asymptotic distribution of the cointegrating vector estimator in error correction models with conditionally heteroskedastic errors. Asymptotic properties of the maximum likelihood estimator (MLE) of the cointegrating vector, which estimates the cointegrating vector and the multivariate GARCH process jointly, are provided. The MLE of the cointegrating vector follows mixture normal, and its asymptotic distribution depends on the conditional heteroskedasticity and the kurtosis of standardized innovations. The reduced rank regression (RRR) estimator and the regression-based cointegrating vector estimators do not consider conditional heteroskedasticity, and thus the efficiency gain of the MLE emerges as the magnitude of conditional heteroskedasticity increases. The simulation results indicate that the relative power of the t-statistics based on the MLE improves significantly as the GARCH effect increases.     

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.     

Prediction from ARFIMA models: Comparisons between MLE and semiparametric estimation procedures

This paper considers the effects on multi-step prediction of using semiparametric local Whittle estimators rather than MLE for long memory ARFIMA models. We consider various representations of the minimum MSE predictor with known parameters. We then conduct a detailed simulation study for when the true parameters are replaced with estimates. The predictor based on MLE is found to be superior, in the MSE sense, to the predictor based on the two-step local Whittle estimation. The “optimal” bandwidth local Whittle estimator produces worse predictions than the local Whittle using an agnostic bandwidth of the square root of the sample size.     

Consistent likelihood-based estimation of a star-shaped distribution

We consider the classical problem of nonparametrically estimating a star-shaped distribution, i.e., a distribution function F on [0,∞) with the property that F(u)/u is nondecreasing on the set {u : F(u) < 1}. This problem is intriguing because of the fact that a well defined maximum likelihood estimator (MLE) exists, but this MLE is inconsistent. In this paper, we argue that the likelihood that is commonly used in this context is somewhat unnatural and propose another, so called ‘smoothed likelihood’. However, also the resulting MLE turns out to be inconsistent. We show that more serious smoothing of the likelihood yields consistent estimators in this model.     

Some model-restricted shrinkage estimators for contingency tables with missing data

For contingency tables with extensive missing data, the unrestricted MLE under the saturated model, computed by the EM algorithm, is generally unsatisfactory. In this case, it may be better to fit a simpler model by imposing some restrictions on the parameter space. Perlman and Wu (1999) propose lattice conditional independence (LCI) models for contingency tables with arbitrary missing data patterns. When this LCI model fits well, the restricted MLE under the LCI model is more accurate than the unrestricted MLE under the saturated model, but not in general. Here we propose certain empirical Bayes (EB) estimators that adaptively combine the best features of the restricted and unrestricted MLEs. These EB estimators appear to be especially useful when the observed data is sparse, even in cases where the suitability of the LCI model is uncertain. We also study a restricted EM algorithm (called the ER algorithm) with similar desirable features. Received: July 1999     

GMM estimation of spatial autoregressive models with moving average disturbances

In this paper, we introduce the one-step generalized method of moments (GMM) estimation methods considered in Lee (2007a) and Liu, Lee, and Bollinger (2010) to spatial models that impose a spatial moving average process for the disturbance term. First, we determine the set of best linear and quadratic moment functions for GMM estimation. Second, we show that the optimal GMM estimator (GMME) formulated from this set is the most efficient estimator within the class of GMMEs formulated from the set of linear and quadratic moment functions. Our analytical results show that the one-step GMME can be more efficient than the quasi maximum likelihood (QMLE), when the disturbance term is simply i.i.d. With an extensive Monte Carlo study, we compare its finite sample properties against the MLE, the QMLE and the estimators suggested in Fingleton (2008a).     

Estimation and forecasting in vector autoregressive moving average models for rich datasets

We address the issue of modelling and forecasting macroeconomic variables using rich datasets by adopting the class of Vector Autoregressive Moving Average (VARMA) models. We overcome the estimation issue that arises with this class of models by implementing an iterative ordinary least squares (IOLS) estimator. We establish the consistency and asymptotic distribution of the estimator for weak and strong VARMA(p,q) models. Monte Carlo results show that IOLS is consistent and feasible for large systems, outperforming the MLE and other linear regression based efficient estimators under alternative scenarios. Our empirical application shows that VARMA models are feasible alternatives when forecasting with many predictors. We show that VARMA models outperform the AR(1), ARMA(1,1), Bayesian VAR, and factor models, considering different model dimensions.     

Distribution free estimation of heteroskedastic binary response models using Probit/Logit criterion functions

In this paper estimators for distribution free heteroskedastic binary response models are proposed. The estimation procedures are based on relationships between distribution free models with a conditional median restriction and parametric models (such as Probit/Logit) exhibiting (multiplicative) heteroskedasticity. The first proposed estimator is based on the observational equivalence between the two models, and is a semiparametric sieve estimator (see, e.g. Gallant and Nychka (1987), Ai and Chen (2003) and Chen et al. (2005)) for the regression coefficients, based on maximizing standard Logit/Probit criterion functions, such as NLLS and MLE. This procedure has the advantage that choice probabilities and regression coefficients are estimated simultaneously. The second proposed procedure is based on the equivalence between existing semiparametric estimators for the conditional median model (,  and ) and the standard parametric (Probit/Logit) NLLS estimator. This estimator has the advantage of being implementable with standard software packages such as Stata. Distribution theory is developed for both estimators and a Monte Carlo study indicates they both perform well in finite samples.     

Small sample considerations in estimation from panel data

Recent interest in statistical inference for panel data has focused on the problem of unobservable, individual-specific, random effects and the inconsistencies they introduce in estimation when they are correlated with other exogenous variables. Analysis of this problem has always assumed the variance components to be known. In this paper, we re-examine some of these questions in finite samples when the variance components must be estimated. In particular, when the effects are uncorrelated with other explanatory variables, we show that (i) the feasible Gauss-Markov estimator is more efficient than the within groups estimator for all but the fewest degrees of freedom and its variance is never more than 17% above the Cramer-Rao bound, (ii) the asymptotic approximation to the variance of the feasible Gauss-Markov estimator is similarly within 17% of the true variance but remains significantly smaller for moderately large samples sizes, and (iii) more efficient estimators for the variance components do not necessarily yield more efficient feasible Gauss-Markov estimators.