Semiparametric deconvolution with unknown error variance

Deconvolution is a useful statistical technique for recovering an unknown density in the presence of measurement error. Typically, the method hinges on stringent assumptions about the nature of the measurement error, more specifically, that the distribution is entirely known. We relax this assumption in the context of a regression error component model and develop an estimator for the unknown density. We show semi-uniform consistency of the estimator and provide an application to the stochastic frontier model.     

On the Kuiper test for normality with mean and variance unknown

Summary A modified form of the Kuiper statistic V _n is developed for testing the composite hypothesis that a sample of size n comes from a normal population with unspecified mean and variance. Its distribution is derived using Monte Carlo methods. Power comparison with the adjusted Kuiper test proposed by L outer and K oerts [6] indicates that our test is superior with respect to certain alternatives.     

On the Kuiper test for normality with mean and variance unknown

Summary If one wants to test the hypothesis as to whether a set of observations comes from a completely specified continuous distribution or not, one can use the Kuiper test. But if one or more parameters have to be estimated, the standard tables for the Kuiper test are no longer valid. This paper presents a table to use with the Kuiper statistic for testing whether a sample comes from a normal distribution when the mean and variance are to be estimated from the sample. The critical points are obtained by means of Monte-Carlo calculation; the power of the test is estimated by simulation; and the results of the powers for several alternative distributions are compared with the estimated powers of the Kolmogorov-Smirnov test.     

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.     

The effects of autocorrelation among errors on the consistency property of OLS variance estimator

In a linear regression model the ordinary least squares (OLS) variance estimator (S²) converges in probability to E(S²) even when the errors are autocorrelated. Of interest, however, is the rate of convergence. In this paper we shed some light on this question for the case of a linear trend model. In particular the relation between the rate of convergence and the correlation property of the errors is explored. It is shown that the retardation of the rate of convergence is not appreciable if the correlation is moderate, but it can be severe for extreme correlations.     

Authorized variance analysis with consumption rate

On the base of budget or plan information, comparing with actual results, making variance analysis, finding real reasons behind variance, this is important way of control and also important function of budget. However, without consideration of changes of environment, there are some limitations for static variance analysis with benchmark of plan data. Adjusting for benchmark according to actual condition, then doing variance analysis, these will improve utilization of variance analysis. Adjusted benchmark is often known as authorized cost/profit. For different understanding for authority concept, with RCA＇s （the abbreviation for Resource Consumption Accounting） view and numerical examples, this paper brings out the concept of consumption rate variance to promoting deeper understanding and analyzing reasons behind variance.     

Linear and nonlinear regression with stable errors

In this paper we describe methods and evaluate programs for linear regression by maximum likelihood when the errors have a heavy tailed stable distribution. The asymptotic Fisher information matrix for both the regression coefficients and the error distribution parameters are derived, giving large sample confidence intervals for all parameters. Simulated examples are shown where the errors are stably distributed and also where the errors are heavy tailed but are not stable, as well as a real example using financial data. The results are then extended to nonlinear models and to non-homogeneous error terms.     

Posterior Inference in Bayesian Quantile Regression with Asymmetric Laplace Likelihood

The paper discusses the asymptotic validity of posterior inference of pseudo‐Bayesian quantile regression methods with complete or censored data when an asymmetric Laplace likelihood is used. The asymmetric Laplace likelihood has a special place in the Bayesian quantile regression framework because the usual quantile regression estimator can be derived as the maximum likelihood estimator under such a model, and this working likelihood enables highly efficient Markov chain Monte Carlo algorithms for posterior sampling. However, it seems to be under‐recognised that the stationary distribution for the resulting posterior does not provide valid posterior inference directly. We demonstrate that a simple adjustment to the covariance matrix of the posterior chain leads to asymptotically valid posterior inference. Our simulation results confirm that the posterior inference, when appropriately adjusted, is an attractive alternative to other asymptotic approximations in quantile regression, especially in the presence of censored data.     

Forecasting economic time series with unconditional time-varying variance

The classical forecasting theory of stationary time series exploits the second-order structure (variance, autocovariance, and spectral density) of an observed process in order to construct some prediction intervals. However, some economic time series show a time-varying unconditional second-order structure. This article focuses on a simple and meaningful model allowing this nonstationary behaviour. We show that this model satisfactorily explains the nonstationary behaviour of several economic data sets, among which are the U.S. stock returns and exchange rates. The question of how to forecast these processes is addressed and evaluated on the data sets.     

Cointegrating rank selection in models with time-varying variance

Reduced rank regression (RRR) models with time varying heterogeneity are considered. Standard information criteria for selecting cointegrating rank are shown to be weakly consistent in semiparametric RRR models in which the errors have general nonparametric short memory components and shifting volatility provided the penalty coefficient _Cn→∞$C_n\to \infty$ and _Cn/n→0$C_n/n\to 0$ as n→∞$n\to \infty$. The AIC criterion is inconsistent and its limit distribution is given. The results extend those in Cheng and Phillips (2009a) and are useful in empirical work where structural breaks or time evolution in the error variances is present. An empirical application to exchange rate data is provided.     

Identification in models with autoregressive errors

Consider the model

$\text{A(L)x}{}_{\mathrm{t}}\text{=B(L)y}{}_{\mathrm{t}}\text{+C(L)z}{}_{\mathrm{t}}\text{=u}{}_{\mathrm{t}}\text{, t=1,…,T}$

, where

$\text{A(L)=(B(L):C(L))}$

is a matrix of polynomials in the lag operator so that $\text{L}{}^{\mathrm{r}}\text{x}{}_{\mathrm{t}}\text{=x}{}_{\mathrm{t?r}}$, and $\text{y}{}_{\mathrm{t}}$ is a vector of n endogenous variables,

$\text{B(L)=}\text{∑}\text{s=0}\text{k}\text{B}{}_{\mathrm{s}}\text{L}{}^{\mathrm{s}}$

$\text{B}{}_0\text{I}{}_{\mathrm{n}}$, and the remaining $\text{B}{}_{\mathrm{s}}$ are $\text{n × n}$ square matrices,

$\text{C(L)=}\text{∑}\text{s=0}\text{k}\text{C}{}_{\mathrm{s}}\text{L}{}^{\mathrm{s}}$

, and $\text{C}{}_{\mathrm{s}}$ is $\text{n × m}$.Suppose that $\text{u}{}_{\mathrm{t}}$ satisfies

$\text{R(L)u}{}_{\mathrm{t}}\text{=e}{}_{\mathrm{t}}$

, where

$\text{R(L)=}\text{∑}\text{s=0}\text{r}\text{R}{}_{\mathrm{s}}\text{L}{}^{\mathrm{s}}$

, $\text{R}{}_0\text{=I}{}_{\mathrm{n}}$, and $\text{R}{}_{\mathrm{s}}$ is a $\text{n × n}$ square matrix. $\text{e}{}_{\mathrm{t}}$ may be white noise, or generated by a vector moving average stochastic process.Now writing

$\text{Ψ(L)=R(L)A(L)}$

, it is assumed that ignoring the implicit restrictions which follow from eq. (1), $\text{Ψ(L)}$ can be consistently estimated, so that if the equation

$\text{Ψ(L)x}{}_{\mathrm{t}}\text{=e}{}_{\mathrm{t}}$

has a moving average error stochastic process, suitable conditions [see E.J. Hannan] for the identification of the unconstrained model are satisfied, and that the appropriate conditions (lack of multicollinearity) on the data second moments matrices discussed by Hannan are also satisfied. Then the essential conditions for identification of the A(L) and R(L) can be considered by requiring that for the true $\text{Ψ(L)}$ eq. (1) has a unique solution for $\text{A(L)}$ and $\text{R(L)}$.There are three types of lack of identification to be distinguished. In the first there are a finite number of alternative factorisations. Apart from a factorisation condition which will be satisfied with probability one a necessary and sufficient condition for lack of identification is that $\text{A(L)}$ has a latent root λ in the sense that for some non-zero vector β,

$\text{β′A(λ)=0}$

.The second concept of lack of identification corresponds to the Fisher conditions for local identifiability on the derivatives of the constraints. It is shown that a necessary and sufficient condition that the model is locally unidentified in this sense is that R(L) and A(L) have a common latent root, i.e., that for some vectors δ and β,

$\text{R(λ)δ=0 and β′A(λ)=0}$

.Firstly it is shown that only if further conditions are satisfied will this lead to local unidentifiability in the sense that there are solutions of the equation

$\text{Ψ(z)=R(z)A(z)}$

in any neighbourhood of the true values.     

Robust trend inference with series variance estimator and testing-optimal smoothing parameter

The paper develops a novel testing procedure for hypotheses on deterministic trends in a multivariate trend stationary model. The trends are estimated by the OLS estimator and the long run variance (LRV) matrix is estimated by a series type estimator with carefully selected basis functions. Regardless of whether the number of basis functions K is fixed or grows with the sample size, the Wald statistic converges to a standard distribution. It is shown that critical values from the fixed-K asymptotics are second-order correct under the large-K asymptotics. A new practical approach is proposed to select K that addresses the central concern of hypothesis testing: the selected smoothing parameter is testing-optimal in that it minimizes the type II error while controlling for the type I error. Simulations indicate that the new test is as accurate in size as the nonstandard test of Vogelsang and Franses (2005) and as powerful as the corresponding Wald test based on the large-K asymptotics. The new test therefore combines the advantages of the nonstandard test and the standard Wald test while avoiding their main disadvantages (power loss and size distortion, respectively).     

The analysis of unbalanced linear models with variance components

Statistical inference for fixed effects, random effects and components of variance in an unbalanced linear model with variance components will be discussed. Variance components will be estimated by Restricted Maximum Likelihood. Iterative procedures for computing the estimates, such as Fisher scoring and the EM-algorithm, are described.     

z-Transform and identification of linear econometric models with autocorrelated errors

The paper consists of two main parts. In the first part we derive the solution of systems of linear stochastic difference equations by means of thez-transform. In the second part thisz-transform is used to treat the problem of identification of linear econometric systems (the term econometric is used to stress the special aspects of the identification problem dealt with in econometrics). It is shown, that under suitable restrictions observationally equivalent structures are related by unimodular matrices. Using this result, we state (rank-) conditions which ensure, that the unimodular matrices are constant, such that the classical econometric identification theorems can be applied. These conditions are given for stationary errors in the general case as well as in the MA, AR and ARMA case.     

GMM estimation of spatial autoregressive models with unknown heteroskedasticity

In the presence of heteroskedastic disturbances, the MLE for the SAR models without taking into account the heteroskedasticity is generally inconsistent. The 2SLS estimates can have large variances and biases for cases where regressors do not have strong effects. In contrast, GMM estimators obtained from certain moment conditions can be robust. Asymptotically valid inferences can be drawn with consistently estimated covariance matrices. Efficiency can be improved by constructing the optimal weighted estimation.     

On new variance approximations for linear models with inequality constraints

In this paper, we examine the estimation of linear models subject to inequality constraints with a special focus on new variance approximations for the estimated parameters. For models with one inequality restriction, the proposed variance formulas are exact. The variance approximations proposed in this paper can be used in regression analysis, Kalman filtering, and balancing national accounts, when inequality constraints are to be incorporated in the estimation procedure.     

Fitting nonlinear time-series models with applications to stochastic variance models

New strategies for the implementation of maximum likelihood estimation of nonlinear time series models are suggested. They make use of recent work on the EM algorithm and iterative simulation techniques. The estimation procedures are applied to the problem of fitting stochastic variance models to exchange rate data.     

Copula-based multivariate GARCH model with uncorrelated dependent errors

Multivariate GARCH (MGARCH) models are usually estimated under multivariate normality. In this paper, for non-elliptically distributed financial returns, we propose copula-based multivariate GARCH (C-MGARCH) model with uncorrelated dependent errors, which are generated through a linear combination of dependent random variables. The dependence structure is controlled by a copula function. Our new C-MGARCH model nests a conventional MGARCH model as a special case. The aim of this paper is to model MGARCH for non-normal multivariate distributions using copulas. We model the conditional correlation (by MGARCH) and the remaining dependence (by a copula) separately and simultaneously. We apply this idea to three MGARCH models, namely, the dynamic conditional correlation (DCC) model of Engle [Engle, R.F., 2002. Dynamic conditional correlation: A simple class of multivariate generalized autoregressive conditional heteroskedasticity models. Journal of Business and Economic Statistics 20, 339–350], the varying correlation (VC) model of Tse and Tsui [Tse, Y.K., Tsui, A.K., 2002. A multivariate generalized autoregressive conditional heteroscedasticity model with time-varying correlations. Journal of Business and Economic Statistics 20, 351–362], and the BEKK model of Engle and Kroner [Engle, R.F., Kroner, K.F., 1995. Multivariate simultaneous generalized ARCH. Econometric Theory 11, 122–150]. Empirical analysis with three foreign exchange rates indicates that the C-MGARCH models outperform DCC, VC, and BEKK in terms of in-sample model selection and out-of-sample multivariate density forecast, and in terms of these criteria the choice of copula functions is more important than the choice of the volatility models.