Estimating stage occupation probabilities in non-Markov models

We study non-Markov multistage models under dependent censoring regarding estimation of stage occupation probabilities. The individual transition and censoring mechanisms are linked together through covariate processes that affect both the transition intensities and the censoring hazard for the corresponding subjects. In order to adjust for the dependent censoring, an additive hazard regression model is applied to the censoring times, and all observed counting and “at risk” processes are subsequently given an inverse probability of censoring weighted form. We examine the bias of the Datta–Satten and Aalen–Johansen estimators of stage occupation probability, and also consider the variability of these estimators by studying their estimated standard errors and mean squared errors. Results from different simulation studies of frailty models indicate that the Datta–Satten estimator is approximately unbiased, whereas the Aalen–Johansen estimator either under- or overestimates the stage occupation probability due to the dependent nature of the censoring process. However, in our simulations, the mean squared error of the latter estimator tends to be slightly smaller than that of the former estimator. Studies on development of nephropathy among diabetics and on blood platelet recovery among bone marrow transplant patients are used as demonstrations on how the two estimation methods work in practice. Our analyses show that the Datta–Satten estimator performs well in estimating stage occupation probability, but that the censoring mechanism has to be quite selective before a deviation from the Aalen-Johansen estimator is of practical importance. N. Gunnes—Supported by a grant from the Norwegian Cancer Society.     

The effect of tapering on the semiparametric estimators for nonstationary long memory processes

In this paper, we study, by a Monte Carlo simulation, the effect of the order p of “Zhurbenko-Kolmogorov” taper on the asymptotic properties of semiparametric estimators. We show that p  =  [d + 1/2] + 1 gives the smallest variances and mean squared errors. These properties depend also on the truncation parameter m. Moreover, we study the impact of the short-memory components on the bias and variances of these estimators. We finally carry out an empirical application by using four monthly seasonally adjusted logarithm Consumer Price Index series.      

A note on conditional Wilcoxon tests with natural and mid-ranks

Consider the nonparametric 2-sample problem with ties. It is shown that the conditional variance (given the vector of lengths of ties) of the Wilcoxon statistic is with “natural” ranks at most as large as with “mid-” (also “mean”, “average”) ranks.     

On prediction and mean squared error

Practical questions motivate the search for predictors either of an as yet unobserved random vector, or of a random function of a parameter. An extension of the classical UMVUE theory is presented to cover such situations. In includes a Rao-Blackwell-type theorem, a Cramer-Rao-type inequality, and necessary and sufficient conditions for a predictor to minimize the mean squared error uniformly in the parameter. Applications are considered to the problem of selected means, the species problem, and the examination of some u-v estimates of Robbins (1988).     

Forecasting Contemporaneously Aggregated Vector ARMA Processes

Given a multiple time series that is generated by a multivariate ARMA process and assuming the objective is to forecast a weighted sum of the individual variables, then under a mean squared error measure of forecasting precision, it is preferable to forecast the disaggregated multiple time series and aggregate the forecasts, rather than forecast the aggregated series directly, if the involved processes are known. This result fails to hold if the processes used for forecasting are estimated from a given set of time series data. The implications of these results for empirical research are investigated using different sets of economic data.     

On the natural restrictions in the singular Gauss–Markov model

A Gauss–Markov model is said to be singular if the covariance matrix of the observable random vector in the model is singular. In such a case, there exist some natural restrictions associated with the observable random vector and the unknown parameter vector in the model. In this paper, we derive through the matrix rank method a necessary and sufficient condition for a vector of parametric functions to be estimable, and necessary and sufficient conditions for a linear estimator to be unbiased in the singular Gauss–Markov model. In addition, we give some necessary and sufficient conditions for the ordinary least-square estimator (OLSE) and the best linear unbiased estimator (BLUE) under the model to satisfy the natural restrictions.      

A comparison of indicators for evaluating x-11-arima seasonal adjustment

Summary The evaluation of the performance of seasonal adjustment procedures is an issue of practical importance in view of the unobservable nature of the components. Looking at just one indicator when judging the overall quality of a procedure may be misleading, even though this is common practice when many series are involved. The main purpose of this paper is to compare the information content of different synthetic indicators with reference to the X-11-ARIMA procedure. Sixty-six different types of monthly seasonal series are generated and the seasonal component then extracted by carrying out X-11-ARIMA with standard options. The correlation between the pseudo-true error for each series and various synthetic indicators allows us to compare the latter's reliability, under both the hypotheses of minimum and maximum variance of the pseudo-true seasonal component. We show that the overall quality indexQ-the indicator most commonly adopted by users of the X-11-ARIMA-is always outperformed by the simpler diagnostics based on the stability of the estimates. In particular, the “sliding-spans” indicator, proposed by Findley et al. (1990) and included in the diagnostics of the new X-12 procedure, shows a much stronger correlation with the pseudo-true error in the seasonal adjustment. We also show that the total forecasting errors in the one-year-ahead extrapolation of the seasonal component have a good informative power and perform almost as well as the “sliding-spans” indicator.     

An alternative stochastic restricted Liu estimator in linear regression

In this paper, we introduce an alternative stochastic restricted Liu estimator for the vector of parameters in a linear regression model when additional stochastic linear restrictions on the parameter vector are assumed to hold. The new estimator is a generalization of the ordinary mixed estimator (OME) (Durbin in J Am Stat Assoc 48:799–808, 1953; Theil and Goldberger in Int Econ Rev 2:65–78, 1961; Theil in J Am Stat Assoc 58:401–414, 1963) and Liu estimator proposed by Liu (Commun Stat Theory Methods 22:393–402, 1993). Necessary and sufficient conditions for the superiority of the new stochastic restricted Liu estimator over the OME, the Liu estimator and the estimator proposed by Hubert and Wijekoon (Stat Pap 47:471–479, 2006) in the mean squared error matrix (MSEM) sense are derived. Furthermore, a numerical example based on the widely analysed dataset on Portland cement (Woods et al. in Ind Eng Chem 24:1207–1241, 1932) and a Monte Carlo evaluation of the estimators are also given to illustrate some of the theoretical results.     

Bayesian inference and forecasting in the stationary bilinear model

A stationary bilinear (SB) model can be used to describe processes with a time-varying degree of persistence that depends on past shocks. This study develops methods for Bayesian inference, model comparison, and forecasting in the SB model. Using monthly U.K. inflation data, we find that the SB model outperforms the random walk, first-order autoregressive AR(1), and autoregressive moving average ARMA(1,1) models in terms of root mean squared forecast errors. In addition, the SB model is superior to these three models in terms of predictive likelihood for the majority of forecast observations.     

Criteria for the equality between ordinary least squares and best linear unbiased estimators under certain linear models

A new necessary and sufficient condition is derived for the equality between the ordinary least-squares estimator and the best linear unbiased estimator of the expectation vector in linear models with certain specific design matrices. This condition is then applied to special cases involving one-way and two-way classification models.     

Some results on random design regression with long memory errors and predictors

This paper studies nonparametric regression with long memory (LRD) errors and predictors. First, we formulate general conditions which guarantee the standard rate of convergence for a nonparametric kernel estimator. Second, we calculate the mean integrated squared error (MISE). In particular, we show that LRD of errors may influence MISE. On the other hand, an estimator for a shape function is typically not influenced by LRD in errors. Finally, we investigate properties of a data-driven bandwidth choice. We show that averaged squared error (ASE) is a good approximation of MISE; however, this is not the case for a cross-validation criterion.     

De gustibus non est disputandum? A new approach to the estimation of equivalence scales

We argue that when the household composition changes, consumption patterns vary not only because of the cost effect that equivalence scales try to measure, but also because of a “taste” or “style” effect. This effect can be identified and measured, under a few assumptions, with the use of a new methodology, calledDM ² (Decomposition Model of the effects of Demographic Metamorphosis), that can be viewed as a generalisation of Ray's (1983) price-scaling approach to the construction of equivalent scales. An empirical application to data drawn from the Istat 1995 Italian Household Budget Survey suggests that the proposed method improves our understanding of households' consumption patterns and the reliability of the equivalence scales that we derive. We gratefully acknowledge helpful comments from Giorgio Calzolari, Franco Polverini, Ugo Trivellato and an anonymous referee, but retain full responsibility for all errors, and for the processing of Istat (Italian Institute of Statistics) microdata on Household Budgets. Financial support for this research was provided by the italian MURST (Research project on “Equivalence scales” directed by Prof. Guido Ferrari, University of Firenze, Ref. No. 9913105354; and Research project on “Low fertility in Italy: between economic constraints and value changes”, directed by Prof. Massimo Livi Bacci, University of Firenze, Ref. No. MM13107238). Preliminary findings on this research topic have been presented in a few seminars and conferences: cf., e.g., De Santis and Maltagliati (2000 and 2001).     

Improvement of the Liu estimator in linear regression model

In the presence of stochastic prior information, in addition to the sample, Theil and Goldberger (1961) introduced a Mixed Estimator for the parameter vector β in the standard multiple linear regression model (T,Xβ,σ² I). Recently, the Liu estimator which is an alternative biased estimator for β has been proposed by Liu (1993). In this paper we introduce another new Liu type biased estimator called Stochastic restricted Liu estimator for β, and discuss its efficiency. The necessary and sufficient conditions for mean squared error matrix of the Stochastic restricted Liu estimator to exceed the mean squared error matrix of the mixed estimator will be derived for the two cases in which the parametric restrictions are correct and are not correct. In particular we show that this new biased estimator is superior in the mean squared error matrix sense to both the Mixed estimator and to the biased estimator introduced by Liu (1993).     

Multiple Time Scales and the Lifetime Coefficient of Variation: Engineering Applications

We consider linear combinations of “natural” timescales and choose the “best” one which provides the minimum coefficient of variation of the lifetime. Our time scale is in fact a generalized Miner time scale because the latter is based on an appropriate weighting of the times spent on low and high level loadings. The suggested modus operandi for finding the“best” time scale has many features in common with the approach suggested by Farewell and Cox (1979) and Oakes (1995) which is devoted to multiple time scales in survival analysis. This revised version was published online in July 2006 with corrections to the Cover Date.     

Improved estimation of the population parameters when some additional information is available

Estimation of population parameters is considered by several statisticians when additional information such as coefficient of variation, kurtosis or skewness is known. Recently Wencheko and Wijekoon (Stat Papers 46:101–115, 2005) have derived minimum mean square error estimators for the population mean in one parameter exponential families when coefficient of variation is known. In this paper the results presented by Gleser and Healy (J Am Stat Assoc 71:977–981, 1976) and Arnholt and Hebert (, 2001) were generalized by considering T (X) as a minimal sufficient estimator of the parametric function g(θ) when the ratio t²=[ g(q) ]^-2Var[ T(X ) ]{\tau^{2}=[ {g(\theta )} ]^{-2}{\rm Var}[ {T(\boldsymbol{X} )} ]} is independent of θ. Using these results the minimum mean square error estimator in a certain class for both population mean and variance can be obtained. When T (X) is complete and minimal sufficient, the ratio τ² is called “WIJLA” ratio, and a uniformly minimum mean square error estimator can be derived for the population mean and variance. Finally by applying these results, the improved estimators for the population mean and variance of some distributions are obtained.     

Forecasting with serially correlated regression models

In this article we investigate the asymptotic and finite-sample properties of predictors of regression models with autocorrelated errors. We prove new theorems associated with the predictive efficiency of generalized least squares (GLS) and incorrectly structured GLS predictors. We also establish the form associated with their predictive mean squared errors as well as the magnitude of these errors relative to each other and to those generated from the ordinary least squares (OLS) predictor. A large simulation study is used to evaluate the finite-sample performance of forecasts generated from models using different corrections for the serial correlation.     

Insights into the appropriate level of disaggregation for efficient time series model forecasting

This paper provides a potentially valuable insight on how to assess if the forecasts from an autoregressive moving average model based on aggregated data could be substantially improved through disaggregation. It is argued that, theoretically, the absence of moving average (MA) terms indicates that no forecasting efficiency improvements can be achieved through disaggregation. In practice, it is found that there is a strong correlation between the statistical significance of the MA component in the aggregate model and the magnitude of the forecast mean square error (MSE) decreases that can be achieved through disaggregation. That is, if a model includes significant MA terms, the forecast MSE improvements that may be gained from disaggregation could be substantial. Otherwise, they are more likely to be relatively small or non-existent.     

Fuzzy prior information and minimax estimation in the linear regression model

We consider the linear regression modely=Xβ+u with prior information on the unknown parameter vector β. The additional information on β is given by a fuzzy set. Using the mean squared error criterion we derive linear estimators that optimally combine the data with the fuzzy prior information. Our approach generalizes the classical minimax procedure firstly proposed by Kuks and Olman.