Least Absolute Regression Revisited

We describe a method for fitting a least absolute residual (LAR) line through a set of two–dimensional points. The algorithm is based on a labeling technique derived from linear programming. It is suited for interactive data analysis and can be carried out with graph paper and a programmable hand calculator. Tests conducted with a Pascal program indicate that the algorithm is computationally efficient.     

A fourier integral density estimate: a Monte Carlo study

Davis (1977) proposed the use of a kernel density estimate which is the sample characteristic function integrated over (-A(n) , A(n)), where A(n) is chosen to minimize the mean integrated square error of the estimate. The scalar, A(n), is determined by the sample size and the population characteristic function. This paper investigates, in a Monte Carlo study, the mean integrated square error obtained under a procedure suggested by Davis (1977) for estimating A(n) when the population characteristic function is unknown.     

Polynomial Histograms for Multivariate Density and Mode Estimation

Abstract. We consider the problem of efficiently estimating multivariate densities and their modes for moderate dimensions and an abundance of data. We propose polynomial histograms to solve this estimation problem. We present first‐ and second‐order polynomial histogram estimators for a general d‐dimensional setting. Our theoretical results include pointwise bias and variance of these estimators, their asymptotic mean integrated square error (AMISE), and optimal binwidth. The asymptotic performance of the first‐order estimator matches that of the kernel density estimator, while the second order has the faster rate of O(n^?6/(d+6)). For a bivariate normal setting, we present explicit expressions for the AMISE constants which show the much larger binwidths of the second order estimator and hence also more efficient computations of multivariate densities. We apply polynomial histogram estimators to real data from biotechnology and find the number and location of modes in such data.     

Prior Distributions for Item Parameters in IRT Models

The focus of this article is on the choice of suitable prior distributions for item parameters within item response theory (IRT) models. In particular, the use of empirical prior distributions for item parameters is proposed. Firstly, regression trees are implemented in order to build informative empirical prior distributions. Secondly, model estimation is conducted within a fully Bayesian approach through the Gibbs sampler, which makes estimation feasible also with increasingly complex models. The main results show that item parameter recovery is improved with the introduction of empirical prior information about item parameters, also when only a small sample is available.     

A Simulation-Based Specification Test for Diffusion Processes

This article makes two contributions. First, we outline a simple simulation-based framework for constructing conditional distributions for multifactor and multidimensional diffusion processes, for the case where the functional form of the conditional density is unknown. The distributions can be used, for example, to form predictive confidence intervals for time period t + τ, given information up to period t. Second, we use the simulation-based approach to construct a test for the correct specification of a diffusion process. The suggested test is in the spirit of the conditional Kolmogorov test of Andrews. However, in the present context the null conditional distribution is unknown and is replaced by its simulated counterpart. The limiting distribution of the test statistic is not nuisance parameter-free. In light of this, asymptotically valid critical values are obtained via appropriate use of the block bootstrap. The suggested test has power against a larger class of alternatives than tests that are constructed using marginal distributions/densities. The findings of a small Monte Carlo experiment underscore the good finite sample properties of the proposed test, and an empirical illustration underscores the ease with which the proposed simulation and testing methodology can be applied.     

Maximum Likelihood Estimation of a Poisson Parameter in a Model of Equioverlapping Samples:A Simulation Study

In a model of equioverlapping samples maximum likelihood estimation of a Poisson parameter is examined and compared with two linear unbiased estimations by mean squared error. Since a likelihood estimator is not explicitly available in general, a simulation study has been performed and the results are illustrated     

Improved estimators of a location vector with unknown scale parameter

The problem of estimating, under arbitrary quadratic loss, the location vector parameter θ of a p-variate distribution (p ≥ 3) with unknown covari-ance matrix ∑ = α² D (where D is a known diagonal matrix) is considered. A large class of improved shrinkage estimators is developed for this problem. This work generalizes results of Berger and Brandwein and Strawderman for the case of a known scale parameter and extends the authors’ results for the class of scale mixtures of normal distributions.     

The Extended Skew-Exponential Power Distribution and Its Derivation

We consider an extended family of asymmetric univariate distributions generated using a symmetric density, f, and the cumulative distribution function, G, of a symmetric distribution, which depends on two real-valued parameters λ and β and is such that when β = 0 it includes the entire class of distributions with densities of the form g(z | λ) = 2 G(λ z) f(z). A key element in the construction of random variables distributed according to the family is that they can be represented stochastically as the product of two random variables. From this representation we can readily derive theoretical properties, easy-to-implement simulation schemes, as well as extensions to the multivariate case and an explicit procedure for obtaining the moments. We give special attention to the extended skew-exponential power distribution. We derive its information matrix in order to obtain the asymptotic covariance matrix of the maximum likelihood estimators. Finally, an application to a real data set is reported, which shows that the extended skew-exponential power model can provide a better fit than the skew-exponential power distribution.     

A Distribution-Free Median-Unbiased Quantile Estimator

Given λ∈(0-,l), let x_λ(F) denote the unique λ-quantile of the distribution F. A distribution-free median-unbiased estimator of x_λ(F) is explicitly constructed     

OPTIMAL IV ESTIMATION OF SYSTEMS WITH STOCHASTIC REGRESSORS AND VAR DISTURBANCES WITH APPLICATIONS TO DYNAMIC SYSTEMS

This paper considers the general problem of Feasible Generalized Least Squares Instrumental Variables (FGLS IV) estimation using optimal instruments. First we summarize the sufficient conditions for the FGLS IV estimator to be asymptotically equivalent to an optimal GLS IV estimator. Then we specialize to stationary dynamic systems with stationary VAR errors, and use the sufficient conditions to derive new moment conditions for these models. These moment conditions produce useful IVs from the lagged endogenous variables, despite the correlation between errors and endogenous variables. This use of the information contained in the lagged endogenous variables expands the class of IV estimators under consideration and thereby potentially improves both asymptotic and small-sample efficiency of the optimal IV estimator in the class. Some Monte Carlo experiments compare the new methods with those of Hatanaka (1976). For the DGP used in the Monte Carlo experiments, asymptotic efficiency is strictly improved by the new IVs, and experimental small-sample efficiency is improved as well.