An algorithm for the minimum sum of weighted absolute errors regression

We propose an algorithm to estimate the unknown constants in a multiple linear regression model under the minimum sum of weighted absolute errors (MSWAE). The proposed algorithm, a generalization of an earlier algorithm, is compared to a bounded variable algorithm. Some somputational experience is reported.     

An Alternative to Transfer Function Forecasting Based on Subspace Methods

In the time series literature, recent interest has focused on the so-called subspace methods. These techniques use canonical correlations and linear regressions to estimate the system matrices of an ARMAX model expressed in state space form. In this article, we use subspace methods to forecast two series with the help of some exogenous variables related to them. We compare the results with those obtained using traditional transfer function models and find that the forecasts obtained with both methods are similar. This result is very encouraging because, in contrast to transfer function models, subspace methods can be considered as almost automatic.     

Partially Accelerated Life Tests for the Weibull Distribution Under Multiply Censored Data

This article aims to estimate the parameters of the Weibull distribution in step-stress partially accelerated life tests under multiply censored data. The step partially acceleration life test is that all test units are first run simultaneously under normal conditions for a pre-specified time, and the surviving units are then run under accelerated conditions until a predetermined censoring time. The maximum likelihood estimates are used to obtaining the parameters of the Weibull distribution and the acceleration factor under multiply censored data. Additionally, the confidence intervals for the estimators are obtained. Simulation results show that the maximum likelihood estimates perform well in most cases in terms of the mean bias, errors in the root mean square and the coverage rate. An example is used to illustrate the performance of the proposed approach.     

Further Results on Order Statistics Generated by Two Simulation Methods

Uniform order statistics generated by two simulation methods are compared by means of Pitman’s measure of closeness. This measure, as a probability, is shown to be asymptotically 1/2. Some results are also established for fixed points of the cumulative distribution function (CDF) for a uniform order statistic. These fixed points are important for calculations involving the joint distribution of these order statistics.     

Exponentiated power Lindley power series class of distributions: Theory and applications

ABSTRACT

We introduce a new four-parameter generalization of the exponentiated power Lindley (EPL) distribution, called the exponentiated power Lindley power series (EPLPS) distribution. The new distribution arises on a latent complementary risks scenario, in which the lifetime associated with a particular risk is not observable; rather, we observe only the minimum lifetime value among all risks. The distribution exhibits a variety of bathtub-shaped hazard rate functions. It contains as particular cases several lifetime distributions. Various properties of the distribution are investigated including closed-form expressions for the density function, cumulative distribution function, survival function, hazard rate function, the rth raw moment, and also the moments of order statistics. Expressions for the Rényi and Shannon entropies are also given. Moreover, we discuss maximum likelihood estimation and provide formulas for the elements of the Fisher information matrix. Finally, two data applications are given showing flexibility and potentiality of the EPLPS distribution.     

Optimal allocation of time points for the random effects model

In this article the problem of the optimal selection and allocation of time points in repeated measures experiments is considered. D‐ optimal designs for linear regression models with a random intercept and first order auto‐regressive serial correlations are computed numerically and compared with designs having equally spaced time points. When the order of the polynomial is known and the serial correlations are not too small, the comparison shows that for any fixed number of repeated measures, a design with equally spaced time points is almost as efficient as the D‐ optimal design. When, however, there is no prior knowledge about the order of the underlying polynomial, the best choice in terms of efficiency is a D‐ optimal design for the highest possible relevant order of the polynomial. A design with equally‐spaced time points is the second best choice     

Optimal testing in a three way anova model

We consider likelihood ratio, score and Wald tests for a three-way random effects ANOVA model. Competitor tests are compared using criteria such as small sample power, asymptotic relative efficiency, and convenient null distribution. The final choice is between a new test and two tests long used in practice.     

On constructing sequences estimating the mixing distribution with applications

This paper shows that by minimizing a Chebychev norm a mixing distribution can be constructed which converges weakly to the true mixing distribution with probability one. Deely and Kruse (1968) established a similar result for the supremum norm. For both norms the constructed mixing distribution is computed by solving a linear programming problem, but this problem is considerably smaller when the Chebychev norm is used. Thus a suitable mixing distribution can be constructed from solving a linear programming problem with considerably less computational work than was previously known. To illustrate the application of this simpler procedure it is applied to derive nonparametric empirical Bayes estimates in a simulation study. Some density estimates are also illustrated.     

A monte carlo investigation of a statistic for a bivariate missing data problem

Testing the equal means hypothesis of a bivariate normal distribution with homoscedastic varlates when the data are incomplete is considered. If the correlational parameter, ρ, is known, the well-known theory of the general linear model is easily employed to construct the likelihood ratio test for the two sided alternative. A statistic, T, for the case of ρ unknown is proposed by direct analogy to the likelihood ratio statistic when ρ is known. The null and nonnull distribution of T is investigated by Monte Carlo techniques. It is concluded that T may be compared to the conventional t distribution for testing the null hypothesis and that this procedure results in a substantial increase in power-efficiency over the procedure based on the paired t test which ignores the incomplete data. A Monte Carlo comparison to two statistics proposed by Lin and Stivers (1974) suggests that the test based on T is more conservative than either of their statistics.