Estimating Quantiles of Normal Populations with a Common Mean

Suppose independent random samples are available from two normal populations with a common mean and unequal variances. Estimation of a quantile of the first population is considered with respect to the quadratic loss. Some new estimators for the quantile are proposed using some previously known estimators of a common mean. Inadmissibility results are proved for estimators which are equivariant under affine and location groups of transformations. Risk values of various estimators of a quantile are compared numerically using a detailed simulation study.     

Lack-of-fit testing of a regression model with response missing at random

This paper proposes a class of lack-of-fit tests for fitting a linear regression model when some response variables are missing at random. These tests are based on a class of minimum integrated square distances between a kernel type estimator of a regression function and the parametric regression function being fitted. These tests are shown to be consistent against a large class of fixed alternatives. The corresponding test statistics are shown to have asymptotic normal distributions under null hypothesis and a class of nonparametric local alternatives. Some simulation results are also presented.     

Defining equations for two-level factorial designs

Defining equations are introduced in the context of two-level factorial designs and they are shown to provide a concise specification of both regular and nonregular designs. The equations are used to find orthogonal arrays of high strength and some optimal designs. The latter optimal designs are formed in a new way by augmenting notional orthogonal arrays which are allowed to have some runs with a negative number of replicates before augmentation. Defining equations are also shown to be useful when the factorial design is blocked.     

A Semiparametric Bayesian Model for Circular-Linear Regression

Circular data are observations that are represented as points on a unit circle. Times of day and directions of wind are two such examples. In this work, we present a Bayesian approach to regress a circular variable on a linear predictor. The regression coefficients are assumed to have a nonparametric distribution with a Dirichlet process prior. The semiparametric Bayesian approach gives added flexibility to the model and is useful especially when the likelihood surface is ill behaved. Markov chain Monte Carlo techniques are used to fit the proposed model and to generate predictions. The method is illustrated using an environmental data set.     

Mixed effects smoothing spline analysis of variance

We propose a general family of nonparametric mixed effects models. Smoothing splines are used to model the fixed effects and are estimated by maximizing the penalized likelihood function. The random effects are generic and are modelled parametrically by assuming that the covariance function depends on a parsimonious set of parameters. These parameters and the smoothing parameter are estimated simultaneously by the generalized maximum likelihood method. We derive a connection between a nonparametric mixed effects model and a linear mixed effects model. This connection suggests a way of fitting a nonparametric mixed effects model by using existing programs. The classical two-way mixed models and growth curve models are used as examples to demonstrate how to use smoothing spline analysis-of-variance decompositions to build nonparametric mixed effects models. Similarly to the classical analysis of variance, components of these nonparametric mixed effects models can be interpreted as main effects and interactions. The penalized likelihood estimates of the fixed effects in a two-way mixed model are extensions of James–Stein shrinkage estimates to correlated observations. In an example three nested nonparametric mixed effects models are fitted to a longitudinal data set.     

ESTIMATING POPULATION SIZE FROM CAPTURE-RECAPTURE EXPERIMENTS IN CONTINUOUS TIME

Estimates for the size of a closed population are given for multiple recapture studies in continuous time. The estimates are derived by a method of moments for martingales. An estimate and associated standard error of the population size are derived for a homogeneous population when the capture rates are permitted to depend on time in an unspecified manner. Corresponding results are obtained when the capture rates vary among individuals as well. Explicit expressions are given for these estimates and standard errors which involve only simple computation.     

Some Results on Modeling Biological Sequences and Web Navigation with a Semi Markov Chain

In the article, the entrance probabilities and the probability distribution of the number of transitions to a state are studied to provide some answers to questions related to state occupancies for the semi Markov model. Biological sequences and Web navigation are two cases that initially seem to be different but to a certain extent they do have similarities. Two main aspects of word occurrences in biological sequences are: (a) where do they occur and (b) how many times do they occur. In Web navigation the similar questions are (a) when a node is visited and (b) how many times a node is visited. So, the theoretical results of this study are applied to model these two cases and derive distributions of word location or node occurrence and frequency of occurrences. Rewards/costs are included in the Web navigation model and analytic forms for the means, variances, and moments of total interval rewards/costs are provided.     

SELECTION OF WEIGHTS FOR WEIGHTED MODEL AVERAGING

We address the task of choosing prior weights for models that are to be used for weighted model averaging. Models that are very similar should usually be given smaller weights than models that are quite distinct. Otherwise, the importance of a model in the weighted average could be increased by augmenting the set of models with duplicates of the model or virtual duplicates of it. Similarly, the importance of a particular model feature (a certain covariate, say) could be exaggerated by including many models with that feature. Ways of forming a correlation matrix that reflects the similarity between models are suggested. Then, weighting schemes are proposed that assign prior weights to models on the basis of this matrix. The weighting schemes give smaller weights to models that are more highly correlated. Other desirable properties of a weighting scheme are identified, and we examine the extent to which these properties are held by the proposed methods. The weighting schemes are applied to real data, and prior weights, posterior weights and Bayesian model averages are determined. For these data, empirical Bayes methods were used to form the correlation matrices that yield the prior weights. Predictive variances are examined, as empirical Bayes methods can result in unrealistically small variances.     

Testing for the Cointegrating Rank of a VAR Process With Structural Shifts

Tests for the cointegrating rank of a vector autoregressive process are considered that allow for possible exogenous shifts in the mean of the data-generation process. The break points are assumed to be known a priori. It is proposed to estimate and remove the deterministic terms such as mean, linear-trend term, and a shift in a first step. Then systems cointegration tests are applied to the adjusted series. The resulting tests are shown to have known limiting null distributions that are free of nuisance parameters and do not depend on the break point. The tests are applied for analyzing the number of cointegrating relations in two German money-demand systems.     

A functional approach to estimating finite population distribution functions

Existing estimators of a finite population distribution function that utilize auxiliary information are often constructed by a point wise argument. As a result, these estimators are not always monotone. We adopt a functional approach to the problem and propose two estimators based on compositions of functions. Asymptotic variance formulae are derived for the proposed es-timators. Comparisons are made with existing estimators in a simulation study using three natural populations.     

Weighing the Evidence for Hypotheses with Small Samples of Right-censored Exponential Data

The p-value evidence for an alternative to a null hypothesis regarding the mean lifetime can be unreliable if based on asymptotic approximations when there is only a small sample of right-censored exponential data. However, a guarded weight of evidence for the alternative can always be obtained without approximation, no matter how small the sample, and has some other advantages over p-values. Weights of evidence are defined as estimators of 0 when the null hypothesis is true and 1 when the alternative is true, and they are judged on the basis of the ensuing risks, where risk is mean squared error of estimation. The evidence is guarded in that a preassigned bound is placed on the risk under the hypothesis. Practical suggestions are given for choosing the bound and for interpreting the magnitude of the weight of evidence. Acceptability profiles are obtained by inversion of a family of guarded weights of evidence for two-sided alternatives to point hypotheses, just as confidence intervals are obtained from tests; these profiles are arguably more informative than confidence intervals, and are easily determined for any level and any sample size, however small. They can help understand the effects of different amounts of censoring. They are found for several small size data sets, including a sample of size 12 for post-operative cancer patients. Both singly Type I and Type II censored examples are included. An examination of the risk functions of these guarded weights of evidence suggests that if the censoring time is of the same magnitude as the mean lifetime, or larger, then the risks in using a guarded weight of evidence based on a likelihood ratio are not much larger than they would be if the parameter were known.     

Empirical and rank processes of observations and residuals

Empirical and rank processes of a sample are shown to converge under weak conditions. These results are then extended to analogous empirical and rank processes formed from the residuals of a linear model. Representations of the limiting processes are obtained by means of a special construction. The technicalities are honed to require few regularity conditions.     

Temporal pattern in number of staff on sick leave: the effect of downsizing

Summary.  The pattern of absenteeism in the downsizing process of companies is a topic in focus in economics and social science. A general question is whether employees who are frequently absent are more likely to be selected to be laid off or in contrast whether employees to be dismissed are more likely to be absent for the remaining time of their working contract. We pursue an empirical and microeconomic investigation of these theses. We analyse longitudinal data that were collected in a German company over several years. We fit a semiparametric transition model based on a mixture Poisson distribution for the days of absenteeism per month. Prediction intervals are considered and the primary focus is on the period of downsizing. The data reveal clear evidence for the hypothesis that employees who are to be laid off are more frequently absent before leaving the company. Interestingly, though, no clear evidence is seen that employees being selected to leave the company are those with a bad absenteeism profile.     

Unit root tests and dramatic shifts with infinite variance processes

A model which explains data that is subject to sudden structural changes of unspecified nature is presented. The structural shifts are generated by a random walk component whose innovations belong to the normal domain of attraction of a symmetric stable law. To test the model against the stationarity case, several non-parametric, and regression-based statistics are studied. The non-parametric tests are a generalization of the variance ratio test to innovations with heavy-tailed distributions. The tests are consistent and shown to have good finite sample size and power properties and are applied to a set of economic variables.     

C93. The relevance of imaginary alternatives

Some practical approaches to the problem of choosing parameters which control the smoothness of kernel-based density estimators are investigated. Fixed and variable kernels are considered, and particularly simple approaches are investigated in the latter case. The performances of a wide range of estimators are compared in a simulation study.     

A Note on Nonparametric Regression with β-Mixing Sequences

Statistical risks of least-squares estimates of a regression function (possibly nonlinear) are studied under a β-mixing setup. The assumptions are minimal and are virtually the same as those of the i.i.d. case. The bounds obtained are optimal, up to a term of order log n. Similar results are also obtained for complexity regularized regression estimates in the presence of penalty functions.     

A class of mean residual life regression models with censored survival data

When describing a failure time distribution, the mean residual life is sometimes preferred to the survival or hazard rate. Regression analysis making use of the mean residual life function has recently drawn a great deal of attention. In this paper, a class of mean residual life regression models are proposed for censored data, and estimation procedures and a goodness-of-fit test are developed. Both asymptotic and finite sample properties of the proposed estimators are established, and the proposed methods are applied to a cancer data set from a clinic trial.     

Optimal record-based statistical procedures for the two-parameter exponential distribution

There are a number of situations in which the experimental data observed are record statistics. In this paper, optimal confidence intervals as well as uniformly most powerful (MP) tests for one-sided alternatives are developed. Since a uniformly MP test for a two-sided alternative does not exist, generalized likelihood ratio and uniformly unbiased and invariant tests are derived for the two parameters of the exponential distribution based on record data. For illustrative purposes, a data set on the times between consecutive telephone calls to a company's switchboard is analysed using the proposed procedures. Finally, some open problems in this direction are pointed out.     

Mixture designs for product improvement studies

Many products are mixtures of several components (ingredient). Characteristics of the products such as the strength of steel, the efficacy of a chemical pesticide, or the viscosity of a liquid detergent, depend only on the relative proportions of the components in the mixture. Studying changes in a product' properties caused by varying the ingredient proportions is the objective of performing mixture experiments. The inherent restriction that the sum of the component proportions equal unity creates different design strategies than are usually employed with independent factors where factorial arrangements are quite common. Experimental designs for exploring the entire mixture simplex region as well as for exploring only a subregion of the simplex are presented. In those cases where four or more components are considered and a subregion is to be investigated, computer-aided designs are the rule rather than the exception. Design criterion based on the properties (variance and bias) of the prediction equation are mentioned briefly and some suggestions are made for future research in mixture experiments.     

Multivariate outlier detection in incomplete survey data: the epidemic algorithm and transformed rank correlations

Summary.  As a part of the EUREDIT project new methods to detect multivariate outliers in incomplete survey data have been developed. These methods are the first to work with sampling weights and to be able to cope with missing values. Two of these methods are presented here. The epidemic algorithm simulates the propagation of a disease through a population and uses extreme infection times to find outlying observations. Transformed rank correlations are robust estimates of the centre and the scatter of the data. They use a geometric transformation that is based on the rank correlation matrix. The estimates are used to define a Mahalanobis distance that reveals outliers. The two methods are applied to a small data set and to one of the evaluation data sets of the EUREDIT project.