A test statistic for weighted runs

A new test statistic based on runs of weighted deviations is introduced. Its use for observations sampled from independent normal distributions is worked out in detail. It supplements the classic χ² test which ignores the ordering of observations and provides additional sensitivity to local deviations from expectations. The exact distribution of the statistic in the non-parametric case is derived and an algorithm to compute p-values is presented. The computational complexity of the algorithm is derived employing a novel identity for integer partitions.     

Generalized choice models for categorical time series

The NDARMA models of Jacobs and Lewis (1983) allow the modeling of categorical processes with an ARMA-like serial dependence structure. These models can be represented through a backshift mechanism, and we analyze marginal and bivariate properties of the resulting backshift process. Motivated by this backshift mechanism, we define the new class of generalized choice (GC) models, which include the usual NDARMA models as a special case, and we derive results describing the marginal and bivariate distribution of the GC model. We discuss implications concerning DMA(∞) models and the serial dependence structure of NDARMA models. Examples show that the family of GC models allows creating sparsely parametrized models for categorical processes with different types of serial dependence structure.     

A note on Läuter-type robust design for polynomial regression models

In this paper, we consider the problem of model robust design for simultaneous parameter estimation among a class of polynomial regression models with degree up to k. A generalized D-optimality criterion, the Ψ_α‐optimality criterion, first introduced by Läuter (1974) is considered for this problem. By applying the theory of canonical moments and the technique of maximin principle, we derive a model robust optimal design in the sense of having highest minimum Ψ_α‐efficiency. Numerical comparison indicates that the proposed design has remarkable performance for parameter estimation in all of the considered rival models.     

Lower bounds of the wrap-around -discrepancy and relationships between MLHD and uniform design with a large size

The wrap-around (WD) L₂-discrepancy has been commonly used in experimental designs. In this paper, some lower bounds of the WD L₂-discrepancy for asymmetrical U-type designs are given and the expectation and variance of midpoint Latin hypercube designs (LHD) are also obtained. Relationships between midpoint LHD and uniform designs for symmetrical and asymmetrical cases are discussed in the sense of comparing the lower bound and the expectation of squared wrap-around L₂-discrepancy of U-type designs. Some comparisons between simple random sampling and the lower bounds of U-type designs are given.     

Sensitivity analysis of longitudinal data with intermittent missing values

Several models for longitudinal data with nonrandom missingness are available. The selection model of Diggle and Kenward is one of these models. It has been mentioned by many authors that this model depends on untested modelling assumptions, such as the response distribution, from the observed data. So, a sensitivity analysis of the study’s conclusions for such assumptions is needed. The stochastic EM algorithm is proposed and developed to handle continuous longitudinal data with nonrandom intermittent missing values when the responses have non-normal distribution. This is a step in investigating the sensitivity of the parameter estimates to the change of the response distribution. The proposed technique is applied to real data from the International Breast Cancer Study Group.     

Empirical Bayes vs. fully Bayes variable selection

For the problem of variable selection for the normal linear model, fixed penalty selection criteria such as AIC, _Cp${\mathrm{C}}_p$, BIC and RIC correspond to the posterior modes of a hierarchical Bayes model for various fixed hyperparameter settings. Adaptive selection criteria obtained by empirical Bayes estimation of the hyperparameters have been shown by George and Foster [2000. Calibration and Empirical Bayes variable selection. Biometrika 87(4), 731–747] to improve on these fixed selection criteria. In this paper, we study the potential of alternative fully Bayes methods, which instead margin out the hyperparameters with respect to prior distributions. Several structured prior formulations are considered for which fully Bayes selection and estimation methods are obtained. Analytical and simulation comparisons with empirical Bayes counterparts are studied.     

A comparative study on a permutation statistic

We studied asymptotic distribution and finite sample properties of a randomly weighted permutation statistic. The asymptotic normality and the finite sample simulations derived from our studies provided theoretical and numerical justifications for distributional assumption of many useful test statistics used in identifying spatial autocorrelations of mapped data. We compared a new method in computing the mean and the approximated variance of the randomly weighted D statistic, a special permutation statistic, with the Walter’s conditional method. In the numerical illustration of the method, we calculated the standardized values of the D statistic by subtracting the mean from the D statistic and dividing the difference by the standard deviation for the standardized mortality ratios (SMRs) and the life expectancies among the 48 states of the continental USA. Spatial autocorrelations of the SMRs and the life expectancies were found to be statistically significant.     

The two-sample -test with a known ratio of variances

We consider the two-sample t-test where error variances are unknown but with known relationships between them. This situation arises, for example, when two measuring instruments average different number of replicates to report the response. In particular we compare our procedure with the usual Satterthwaite approximation in the two sample t-test with variances unequal. Our procedure uses the knowledge of a known ratio of variances while the Satterthwaite approximation assumes only that the two variances are unequal. Simulations show that our procedure has both better size and better power than the Satterthwaite approximation. Finally, we consider an extension of our results to the General Linear Model.     

A note on the queue with a random set-up time under a restricted admissibility policy with a Bernoulli vacation schedule

This paper deals with the queue size distribution of an M^x/G/1 queue with a random set-up time and with a Bernoulli vacation schedule under a restricted admissibility policy. This generalizes the model studied by Madan and Choudhury [Sankhyá 66 (2004) 175–193].     

Least squares estimation of the coefficients of a partially explosive model,with polynomial regressions of same degree,generating a pair of related time series

This paper is concerned with a partially explosive linear model with polynomial regression components generating a pair of related time series. The least squares estimates of the coefficients are shown to be √N-consistent and asymptotically singular normal, when the degrees of polynomial regression components are same, thus generalising a result due to Venkataraman (1974).     

Competing risks simulation with the survsim R package

The goal of this work is to describe how the last version available of the survsim R package can be used to simulate a cohort in a competing risks context by means of a cause-specific hazards model following the ideas introduced by Beyersmann in 2009, and also allowing for individual heterogeneity through a random effect. An example of its application based on a real cohort will be discussed.     

Nonparametric control charts based on runs and Wilcoxon-type rank-sum statistics

In this article, we introduce three new distribution-free Shewhart-type control charts that exploit run and Wilcoxon-type rank-sum statistics to detect possible shifts of a monitored process. Exact formulae for the alarm rate, the run length distribution, and the average run length (ARL) are all derived. A key advantage of these charts is that, due to their nonparametric nature, the false alarm rate (FAR) and in-control run length distribution is the same for all continuous process distributions. Tables are provided for the implementation of the charts for some typical FAR values. Furthermore, a numerical study carried out reveals that the new charts are quite flexible and efficient in detecting shifts to Lehmann-type out-of-control situations.     

A two phases queueing system with Bernoulli vacation schedule under multiple vacation policy

This paper deals with a single server Poisson arrival queue with two phases of heterogeneous service along with a Bernoulli schedule vacation model, where after two successive phases service the server either goes for a vacation with probability p (0≤p≤1) or may continue to serve the next unit, if any, with probability q(=1−p). Further the concept of multiple vacation policy is also introduced here. We obtained the queue size distributions at a departure epoch and at a random epoch, Laplace Stieltjes Transform of the waiting time distribution and busy period distribution along with some mean performance measures. Finally we discuss some statistical inference related issues.     

Estimating the Lyapunov exponent from chaotic time series with dynamic noise

In this paper, we propose an estimator of the Lyapunov exponent of the skeleton for chaotic time series with dynamic noise and prove the consistency of the estimator under some assumptions.     

Inference for a simultaneous linear partially explosive model with polynomial regression components of different degrees

This paper is concerned with the estimation of the coefficients of simultaneous partially explosive model with polynomial regression components of different degrees in its equations. Since the least squares method breaks down in this case, a three stage estimation procedure is suggested for obtaining CAN estimates of the coefficients.     

Robust multivariate Bayesian forecasting under functional distortions in the -metric

This paper investigates robustness of multivariate forecasting in the Bayesian framework. The minimax approach is used to construct robust statistical procedures under deviations from hypothetical assumptions. The deviations are defined as functional distortions using the χ²-pseudo-metric. Two cases of deviations are considered: distortions of parameter distribution and distortions of joint distribution of observations and parameters. Explicit forms for the guaranteed upper risk functional are obtained and integral equations for robust prediction statistics are given for both cases.     

Quantile estimation in ultra-high frequency financial data: a comparison between parametric and semiparametric approach

In the context of ACD models for ultra-high frequency data different specifications are available to estimate the conditional mean of intertrade durations, while quantiles estimation has been completely neglected by literature, even if to trading extent it can be more informative. The main problem arising with quantiles estimation is the correct specification of durations probability law: the usual assumption of Exponentially distributed residuals, is very robust for the estimation of parameters of the conditional mean, but dramatically fails the distributional fit. In this paper a semiparametric approach is formalized, and compared with the parametric one, deriving from Exponential assumption. Empirical evidence for a stock of Italian financial market strongly supports the former approach.Paola Zuccolotto: The author wishes to thank Prof. A. Mazzali, Dott. G. De Luca, Dott. M. Sandri for valuable comments.