Testing for Stochastic Monotonicity

We propose a test of the hypothesis of stochastic monotonicity. This hypothesis is of interest in many applications in economics. Our test is based on the supremum of a rescaled U‐statistic. We show that its asymptotic distribution is Gumbel. The proof is difficult because the approximating Gaussian stochastic process contains both a stationary and a nonstationary part, and so we have to extend existing results that only apply to either one or the other case. We also propose a refinement to the asymptotic approximation that we show works much better in finite samples. We apply our test to the study of intergenerational income mobility.     

Mean functions based on meta-mixtures in nonhomogeneous Poisson processes

This paper deals with the software reliability model based on a nonhomogeneous Poisson process. We introduce new types of mean functions which can be either NHPP-I or NHPP-II according to the choice of the distribution function. The proposed mean function is motivated by the fact that a strictly monotone increasing function can be modeled by a distribution function and an unknown distribution function approximated by a mixture of beta distributions. Some existing mean functions can be regarded as special cases of the proposed mean functions. The EM algorithm is used to obtain maximum likelihood estimates of the parameters in the proposed model.     

A Methodology for More Efficient Tail Area Sampling with Discrete Probability Distributions

Monte Carlo Method is commonly used to observe the overall distribution and to determine the lower or upper bound value in statistical approach when direct analytical calculation is unavailable. However, this method would not be efficient if the tail area of a distribution is concerned. A new method, entitled Two-Step Tail Area Sampling, is developed, which uses the assumption of discrete probability distribution and samples only the tail area without distorting the overall distribution. This method uses a two-step sampling procedure. First, sampling at points separated by large intervals is done and second, sampling at points separated by small intervals is done with some check points determined at first-step sampling. Comparison with Monte Carlo Method shows that the results obtained from the new method converge to analytic value faster than Monte Carlo Method if the numbers of calculation of both methods are the same. This new method is applied to DNBR (Departure from Nucleate Boiling Ratio) prediction problem in design of the pressurized light water nuclear reactor.     

Robust precision matrix estimation via weighted median regression with regularization

A precision matrix is an important parameter of interests because its elements describe useful association information among multiple variables, which has a wide variety of applications. For example, it is used for inferring gene regulation networks in genomic studies and stock association networks in financial studies. However, in many cases, the precision matrix needs to be robustly estimated due to the presence of outliers. We propose estimating a sparse scaled precision matrix via weighted median regression with regularization. Our weighted median regression approach is consistent under various distributional assumptions including multivariate t‐ or contaminated Gaussian distributions. This fact is illustrated with simulation studies and a real data analysis with monthly stock return data. The Canadian Journal of Statistics 46: 265–278; 2018 © 2018 Statistical Society of Canada     

Generalized autoregressive and moving average models: multicollinearity,interpretation and a new modified model

In this paper, we call attention of two observed features in practical applications of the Generalized Autoregressive Moving Average (GARMA) model due to the structure of its linear predictor. One is the multicollinearity which may lead to a non-convergence of the maximum likelihood, using iteratively reweighted least squares, and the inflation of the estimator's variance. The second is that the inclusion of the same lagged observations into the autoregressive and moving average components confounds the interpretation of the parameters. A modified model, GAR-M, is presented to reduce the multicollinearity and to improve the interpretation of the parameters. The expectation and variance under stationarity conditions are presented for the identity and logarithm link function. In a general sense, simulation studies show that the maximum likelihood estimators based on the GARMA and GAR-M models are equivalent but the GAR-M estimators presented a little lower standard errors and some restrictions in the parametric space are imposed to guarantee the stationarity of the process. Also, a real data analysis illustrates the GAR-M fit for daily hospitalization rates of elderly people due to respiratory diseases from October 2012 to April 2015 in São Paulo city, Brazil.     

A note on the differentiation formula in Stratonovich type for fractional Brownian sheet

We consider a new proof on the differentiation formula in Stratonovich type for fractional Brownian sheet. Our proof is based on the repeated applications of differentiation formulas in Stratonovich form for a one-parameter Gaussian process.     

Graphical methods for evaluating ridge regression estimator in mixture experiments

When the component proportions in mixture experiments are restricted by lower and upper bounds, multicollinearity appears all too frequently. Thus, we can suggest the use of ridge regression as a mean for stabilizing the coefficient estimates in the fitted model. We propose graphical methods for evaluating the effect of ridge regression estimator with respect to the predicted response value and the prediction variance.     

Hierarchical Bayesian Approach to a Multi-Site Hidden Markov Model

Multivariate data with a sequential or temporal structure occur in various fields of study. The hidden Markov model (HMM) provides an attractive framework for modeling long-term persistence in areas of pattern recognition through the extension of independent and identically distributed mixture models. Unlike in typical mixture models, the heterogeneity of data is represented by hidden Markov states. This article extends the HMM to a multi-site or multivariate case by taking a hierarchical Bayesian approach. This extension has many advantages over a single-site HMM. For example, it can provide more information for identifying the structure of the HMM than a single-site analysis. We evaluate the proposed approach by exploiting a spatial correlation that depends on the distance between sites.     

Simultaneous prediction intervals for all patrwise differences among several future sample means

Hahn (1977) suggested a procedure for constructing prediction intervals for the difference between the means of two future samples from normal populations having equal variance, based on past samples selected from both populations. In this paper, we extend Hahn's work by constructing simultaneous prediction intervals for all pairwise differences among the means of k ≥ 2 future samples from normal populations with equal variances, using past samples taken from each of the k populations. For K = 2, this generalization reduces to Hahn's special case. These prediction intervals may be used when one has sampled the performance of several products and wishes to simultaneously as- sess the differences in future sample mean performance of these products with a predetermined overall coverage probability. The use of the new procedure is demonstrated with a numerical example.     

Some rememarks on locally most powerful unbiased tests for the detection of mixtures of poisson or binomial distributions

In t h i s note mixture models are used to represent overdispersion relative to Poisson or binomial distributions. We flnd a sufflclent condition on the mixing distribution underich the detection of mixture departures from the Poisson or binomial adrnits a locally most powerful unbiased test. The conditions specify plynoria: relations between the variance and mean of Le glxing distribution.