Search

Let X1,…,Xn be a sample from a N(μ,σ2) distribution. Consider testing H0:μ=μ0 vs. Ha:μ>μ0 The usual t-test is based on the statistic t=X¯−μ0S/n where X¯=1n∑i=1nXi and S2=1n−1∑i=1n(Xi−X¯)2 Kim (2019...

Ebaid et al. (2020) proposed the following new extended FGM copula (1) C ( u , v ) = u v [ 1 + a ( 1 − u ) ( 1 − v ) ( 1 − b u ) ( 1 − b v ) ] , 0 ≤ b < 1. Ebaid et al. (2020) claimed that, the adm...

Kageyama Mohan (1984) have presented three methods of constructing new incomplete block designs from balanced incomplete block designs, They raise questions about the designs which come from each of their methods, These questions are answered, Another series of group divisible designs is derived as a special case of their second method.

The problem of estimation of the total weight or objects using a spring balance weighing design has been deait with in this paper Based on a theorem by Dey and Gupta (1977) giving a lower bound for the variance of the estimated total weight, a necessary and sufficient condition for this lower bound to be attained is obtained. A few special cases where the lower bound is attained are enumerated.

Researchers faced with a sequence of candidate model specifications must often choose the best specification that does not violate a testable identification assumption. One option in this scenario is sequential specification tests: hypothesis tests of the identification assumption over the sequence. Borrowing an idea from the change-point literature, this paper shows how to use the distribution of p-values from sequential specification tests to estimate the point in the sequence where the identification assumption ceases to hold. Unlike current approaches, this method is robust to individual errant p-values and does not require choosing a test level or tuning parameter. This paper demonstrates the method's properties with a simulation study, and illustrates it by application to the problems of choosing a bandwidth in a regression discontinuity design while maintaining covariate balance and of choosing a lag order for a time series model., Forthcoming in Communications in Statistics -- Theory and Methods

In this paper, we propose to consider the dependence structure of the trade/no trade categorical sequence of individual illiquid stocks returns. The framework considered here is wide as constant and time-varying zero returns probability are allowed. The ability of our approach in highlighting illiquid stock's features is underlined for a variety of situations. More specifically, we show that long-run effects for the trade/no trade categorical sequence may be spuriously detected in presence of a non-constant zero returns probability. Monte Carlo experiments, and the analysis of stocks taken from the Chilean financial market, illustrate the usefulness of the tools developed in the paper.

Various objective priors have been defined for the common correlation coefficient concerning several bivariate normal populations. In this paper, the proposed approach relies on the asymptotic matching of coverage probabilities corresponding to Bayesian credible intervals considering the corresponding frequentist ones. In the present paper, we focus on several matching criteria including quantile matching, distribution function matching, highest posterior density matching, and matching via inversion of test statistics. In addition, we consider reference priors for different groups of ordering. The proposed methods are investigated and compared between each other in terms of a frequentist coverage probability and then, they are illustrated through a simulation study and two real data examples.

In general, the transition matrix of a Markov chain is a stochastic matrix. For a system that is modeled by a Markov chain, the transition matrix must reflect the characteristics of that system. The present paper introduces a particular class of transition matrices in order to model Markov systems for which, as the length of the time interval becomes greater, a transition from one state to another is more likely. We call these transition matrices stepwise increasing. Moreover in some contexts it is less desirable that the stocks fluctuate over time. In those situations, one is interested in monotonically proceeding stock vectors. This paper examines monotonically proceeding stock vectors and stepwise increasing transition matrices. We present conditions on the transition matrix such that all stock vectors are monotonically proceeding. In particular, the set of monotonically proceeding vectors is characterized for the two-state and three-state cases.

In this paper, we introduce the concept of model quality for two-level regular fractional factorial designs. Under the effect hierarchy principle, this paper raises the definition of model ...

In this paper we introduce a sparse Bayesian kernel multinomial probit regression model for multi-class cancer classification. The relationship between the cancer types and gene expression measurements is explained by an unknown function which belongs to an abstract functional space like the reproducing kernel Hilbert space. We assign a sparse prior for regression parameters and perform variable selection by indexing the covariates of the model with a binary vector. The correlation prior for the binary vector assigned in this paper is able to distinguish models with the same size. The proposed method is successfully tested on one simulated data set and two publicly available real life data sets.

Wavelet analysis has been proved to be a powerful statistical technique in the non parametric regression. In this paper, we propose non linear wavelet-based estimators for multivariable mean regression function with long-memory data. We also provide an asymptotic expansion for the mean integrated squared error (MISE) of the function estimators. This MISE expansion still works even when the underlying mean regression function is only piecewise smooth. This paper extends the corresponding results in the literature for single variable to multivariable case.

The accurate estimation of an individual's usual dietary intake is an important topic in nutritional epidemiology. This paper considers the best linear unbiased predictor (BLUP) computed from repeatedly measured dietary data and derives several nonparametric prediction intervals for true intake. However, the performance of the BLUP and the validity of prediction intervals depends on whether required model assumptions for the true intake estimation problem hold. To address this issue, the paper examines how the BLUP and prediction intervals behave in the case of a violation of model assumptions, and then proposes an analysis pipeline for checking them with data.

This paper proposes new symmetric and asymmetric distributions applying methods analogous as the ones in Kim (2005) and Arnold et al. (2009) to the exponentiated normal distribution studied in Durrans (1992), that we call the power-normal (PN) distribution. The proposed bimodal extension, the main focus of the paper, is called the bimodal power-normal model and is denoted by BPN(α) model, where α is the asymmetry parameter. The authors give some properties including moments and maximum likelihood estimation. Two important features of the model proposed is that its normalizing constant has closed and simple form and that the Fisher information matrix is nonsingular, guaranteeing large sample properties of the maximum likelihood estimators. Finally, simulation studies and real applications reveal that the proposed model can perform well in both situations.

Multistage ranked-set sampling (MRSS) is a generalization of ranked-set sampling in which multiple stages of ranking are used. It is known that for a fixed distribution under perfect rankings, each additional stage provides a gain in efficiency when estimating the population mean. However, the maximum possible efficiency for the MRSS sample mean relative to the simple random sampling sample mean has not previously been determined. In this paper, we provide a method for computing this maximum possible efficiency under perfect rankings for any choice of the set size and the number of stages. The maximum efficiency tends to infinity as the number of stages increases, and, for large numbers of stages, the efficiency-maximizing distributions are symmetric multi-modal distributions where the number of modes matches the set size. The results in this paper correct earlier assertions in the literature that the maximum efficiency is bounded and that it is achieved when the distribution is uniform.

The paper deals with the problem of bounded risk point estimation for a linear combination of location parameters of two negative exponential distributions. Isogai and Futschik considered the situation when the location and scale parameters are all unknown. They proposed purely sequential procedures and gave second order expansions of the average sample sizes and risks. In this paper we propose three-stage procedures and derive second order expansions of the average sample sizes and risks. Further, we compare the results with those from previous work.

This paper presents nonparametric predictive inference for future order statistics. Given data consisting of n real-valued observations, m future observations are considered and predictive probabilities are presented for the r-th ordered future observation. In addition, joint and conditional probabilities for events involving multiple future order statistics are presented. The paper further presents the use of such predictive probabilities for order statistics in statistical inference, in particular considering pairwise and multiple comparisons based on two or more independent groups of data.

In this paper, we introduce a new concept of Poisson Stepanov-like almost automorphy (or Poisson S2-almost automorphy). Under some suitable conditions on the coefficients, we establish the existence and uniqueness of Stepanov-like almost automorphic mild solution to a class of semilinear stochastic differential equations with infinite dimensional Levy noise. We further discuss the global asymptotic stability of these solution. Finally, we give an example to illustrate the theoretical results obtained in this paper.

In this paper, the classical statistical test based on intuitionistic fuzzy hypotheses in relation to the underlying population parametric is extended. In this approach, the type-I, type-II, power of test, and p-value are extended for intuitionistic fuzzy hypotheses. Throughout the paper, some applied examples are provided for both parametric and non parametric cases to clarify the discussions.

Copulas have proved to be very successful tools for the flexible modeling of dependence. Bivariate copulas have been deeply researched in recent years, while building higher-dimensional copulas is still recognized to be a difficult task. In this paper, we study the higher-dimensional dependent reliability systems using a type of decomposition called “vine,” by which a multivariate distribution can be decomposed into a cascade of bivariate copulas. Some equations of system reliability for parallel, series, and k-out-of-n systems are obtained and then decomposed based on C-vine and D-vine copulas. Finally, a shutdown system is considered to illustrate the results obtained in the paper.

In this paper, we investigate an M/M/1 queue with a two-stage vacation policy which comprises of single working vacation and single vacation. Using the matrix-analytic method, we obtain the distribution of stationary system size, and then the decomposition structures of the stationary system size and the sojourn time are demonstrated. Furthermore, we study the waiting time by first-passage time analysis. Meanwhile, the busy-cycle analysis is provided by the limiting theorem of alternative renewal process. Finally, several numerical examples are presented in the paper.

Elasticity (or elasticity function) is a new concept that allows us to characterize the probability distribution of any random variable in the same way as characteristic functions and hazard and reverse hazard functions do. Initially defined for continuous variables, it was necessary to extend the definition of elasticity and study its properties in the case of discrete variables. A first attempt to define discrete elasticity is seen in Veres-Ferrer and Pavia (2014a). This paper develops this definition and makes a comparative study of its properties, relating them to the properties shown by discrete hazard and reverse hazard, as both defined in Chechile (2011). Similar to continuous elasticity, one of the most interesting properties of discrete elasticity focuses on the rate of change that this undergoes throughout its support. This paper centers on the study of the rate of change and develops a set of properties that allows us to carry out a detailed analysis. Finally, it addresses the calculati...

In this paper, the class of Lamé Lorenz curves is studied. This family has the advantage of modeling inequality with a single parameter. The family has a double motivation: it can be obtain from an economic model and from simple transformations of classical Lorenz curves. The underlying cumulative distribution functions have a simple closed form, and correspond to the Singh-Maddala and Dagum distributions, which are well known in the economic literature. The Lorenz order is studied and several inequality and polarization measures are obtained, including Gini, Donaldson-Weymark-Kakwani, Pietra and Wolfson indices. Some extensions of the Lamé family are obtained. Fitting and estimation methods under two different data configuration are proposed. Empirical applications with real data are given. Finally, some relationships with other curves are included. The authors thank to Ministerio de Econom a y Competitividad, project ECO2010-15455, for partial support. The second author thanks to the Ministerio de Educaci on (FPU AP-2010-4907) for partial support. We are grateful for the constructive suggestions provided by the reviewers, which improved the paper.

Experience ratemaking plays a crucial role in general insurance in determining future premiums of individuals in a portfolio by assessing observed claims from the whole portfolio. This paper investigates this problem in which claims can be modeled by certain parametric family of distributions. The Dirichlet process mixtures are employed to model the distributions of the parameters so as to make two advantages: to produce exact Bayesian experience premiums for a class of premium principles generated from generic error functions and, at the same time, provide robust and flexible ways to avoid possible bias caused by traditionally used priors such as non informative priors or conjugate priors. In this paper, the Bayesian experience ratemaking under Dirichlet process mixture models are investigated and due to the lack of analytical forms of the conditional expectations of the quantities concerned, the Gibbs sampling schemes are designed for the purpose of approximations.

In this paper, we investigate the complete moment convergence and Lr convergence for maximal partial sums of asymptotically almost negatively associated random variables under some general conditions. The results obtained in the paper generalize some corresponding ones for negatively associated random variables.

International audience; The zero-inflated binomial (ZIB) regression model was proposed to account for excess zeros in binomial regression. Since then, the model has been applied in various fields, such as ecology and epidemiology. In these applications, maximum-likelihood estimation (MLE) is used to derive parameter estimates. However, theoretical properties of the MLE in ZIB regression have not yet been rigorously established. The current paper fills this gap and thus provides a rigorous basis for applying the model. Consistency and asymptotic normality of the MLE in ZIB regression are proved. A consistent estimator of the asymptotic variance-covariance matrix of the MLE is also provided. Finite-sample behavior of the estimator is assessed via simulations. Finally, an analysis of a data set in the field of health economics illustrates the paper.

The aim of this paper is obtaining the amount of information there exists in the Pareto distribution in the presence of outliers. For the sake of this purpose, Shannon entropy, ϵ-entropy, Fisher information, and Kullback–Leibler distance are computed. Furthermore, a section has been devoted to compare these quantities in these two cases of the Pareto distribution (with outliers and the homogenous case). At the end of this paper, two actual examples, which are related to insurance companies, are brought. A brief summary of which is done in this work is also reported.

When are inferences (whether Direct-Likelihood, Bayesian, or Frequentist) obtained from partial data valid? This paper answers this question by offering a new asymptotic theory about inference with missing data that is more general than existing theories. It proves that as the sample size increases and the extent of missingness decreases, the average-loglikelihood function generated by partial data and that ignores the missingness mechanism will converge in probability to that which would have been generated by complete data; and if the data are Missing at Random, this convergence depends only on sample size. Thus, inferences from partial data, such as posterior modes, confidence intervals, likelihood ratios, test statistics, and indeed, all quantities or features derived from the partial-data loglikelihood function, will be consistently estimated. Additionally, the missing data mechanism has asymptotically no effect on parameter estimation and hypothesis testing if the data are Missing at Random. This adds to previous research which has only proved the consistency and asymptotic normality of the posterior mode. Practical implications are discussed, and the theory is illustrated through simulation using a previous study of International Human Rights Law., Submitted

In this note, we make some comments about the paper of Alheety and Kibria (2014) and correct the wrongly proved Theorems in that paper. © 2017 Taylor & Francis Group, LLC.

Belzunce et al. (1995) define the elasticity for non negative random variables as the reversed proportional failure rate (RPFR). Veres-Ferrer and Pavia (2012, 2014b) interpret it in economic terms, extending its definition to variables that can also take negative values, and briefly present the role of elasticity in characterizing probability distributions. This paper highlights a set of properties demonstrated by elasticity, which shows many similar properties to the reverse hazard function. This paper pays particular attention to studying the increase/decrease and the speed of change of the elasticity function. These are important properties because of the characterizing role of elasticity, which makes it possible to introduce our hypotheses and knowledge about the random process in a more meaningful and intuitive way. As a general rule, it is observed the need for distinguishing between positive and negative areas of the support.

This paper suggests a new stratified randomized response model based on Kuk's [Biometrika (1990), 77, 2, pp.436–438] model that has Neyman allocation and considerable gain in precision. It has been identified that the stratified randomized response models due to Kim and Warde (2004), Kim and Elam's (2005), and Kim and Elam's (2007) are members of the proposed model. It is shown that the proposed model is more efficient than Kuk's (1990) model both theoretically and empirically. The results of this paper are also extended in the situation when trials are repeated.

In this paper, we provide necessary conditions for discrete-time symmetric α-stable processes to be linear 2-ple Markov. The aim of this paper is to extend the results given by Adler et al. (1990) ...

This paper considers posterior consistency in the context of high-dimensional variable selection using the Bayesian lasso algorithm. In a frequentist setting, consistency is perhaps the most basic property that we expect any reasonable estimator to achieve. However, in a Bayesian setting, consistency is often ignored or taken for granted, especially in more complex hierarchical Bayesian models. In this paper, we have derived sufficient conditions for posterior consistency in the Bayesian lasso model with the orthogonal design, where the number of parameters grows with the sample size.

The basic model in this paper is an AR(1) model with a structural break in the AR parameter β at an unknown time k0. That is, yt = β1yt − 1I{t ⩽ k0} + β2yt − 1I{t > k0} + ϵt, t = 1, 2, ⋅⋅⋅, T, where I{ · } denotes the indicator function. Suppose |β1| < 1, |β2| < 1, and {ϵt, t ⩾ 1} is a sequence of i.i.d. random variables which are in the domain of attraction of the normal law with zero mean and possibly infinite variance, then the limiting distributions for the least squares estimators of β1 and β2 are studied in the present paper, which extend some results in Chong (2001).

The paper reviews finite mixture models for binomial counts with concomitant variables. These models are well known in theory, but they are rarely applied. We use a binomial finite mixture to model the number of credits gained by freshmen during the first year at the School of Economics of the University of Florence. The finite mixture approach allows us to appropriately account for the large number of zeroes and the multimodality of the observed distribution. Moreover, we rely on a concomitant variable specification to investigate the role of student background characteristics and of a compulsory pre-enrollment test in predicting gained credits. In the paper, we deal with model selection, including the choice of the number of components, and we devise numerical and graphical summaries of the model results in order to exploit the information content of the concomitant variable specification. The main finding is that the introduction of the pre-enrollment test gives additional information for student tutor...

This paper investigates the ruin probabilities of a renewal risk model with stochastic investment returns and dependent claim sizes. The investment is described as a portfolio of one risk-free asset and one risky asset whose price process is an exponential Levy process. The claim sizes are assumed to follow a two-sided linear process with independent and identically distributed step sizes. When the step-size distribution is heavy tailed, the paper establishes some uniform asymptotic formulas of ruin probabilities.

There are n cards serially numbered from 1 to n. The cards are shuffled and placed in a line one after the other on top of a table with faces up. The numbers on the faces are read from left to right. If there are consecutive numbers in increasing order of magnitude the corresponding cards are merged into one. After the merger, the cards are numbered serially from one to whatever the number of cards we now have. The cards are shuffled and placed in a line one after another on top of the table with faces up. The process continues until we have only one card left. In this paper, we develop a probabilistic recurrence relation approach to obtain the mean, variance, and distribution of the number of shuffles needed. A Markov chain formulation and its properties are discussed in the paper as well.

Multiple comparisons of the effects of several treatments with a control (MCC) has been a central problem in medicine and other areas. Nearly all of existing papers are devoted to comparing means of the effects. To study medical problems more deeply, one needs more information than mean relationship from the given data. It can be expected to get more useful and deeper conclusion by comparing the probability distributions, i.e., by comparison under stochastic orders. This paper presents a likelihood ratio testing procedure to compare effects under stochastic order for MCC problems, controlling the false discovery rate (FDR). Setting a test controlling FDR under stochastic order faces several non trivial problems. These problems are analyzed and solved in this paper. To facilitate the test more easily, the asymptotic p values for the test are used and their distributions are derived. It is shown that controllability of FDR for this comparison procedure can be guaranteed. A real data example is used to illus...

The score statistic S2 is commonly used for general likelihood-based inference. Pearson’s Chi-squared statistic X2 = ∑(O − E)2/E is ubiquitous in contingency table inference. Because tests and confidence intervals based on S2 have been shown to work well in practice and theory and because X2 has such a simple and intuitively appealing form, it is of interest to know when S2 is identical to X2 and when X2 has an approximate Chi-squared distribution. Toward these ends, this paper gives a simple proof that S2 = X2 for the broad class of multinomial-Poisson distributions when the alternative hypothesis is unrestricted in a certain sense. This paper also gives a sufficient condition under which the null distribution of the Pearson score statistic is approximately Chi-squared. Several examples illustrate the utility of the results and counter-examples highlight the importance of the sufficient conditions of the results.

In this paper, we have presented a proportional hazard version of the Marshall–Olkin extended family of distributions. This family of distributions has been compared in terms of stochastic orderings with the Marshall-Olkin extended family of distributions. Considering the Gompertz distribution as the baseline, the monotonicity of the resulting failure rate is shown to be either increasing or bathtub, even though the Gompertz distribution has an increasing failure rate. The maximum likelihood estimation of the parameters has been studied and a data set, involving the serum–reversal times, has been analyzed and it has been shown that the model presented in this paper fit better than the Gompertz or even the Mrashall–Olkin Gompertz distribution. The extension presented in this paper can be used in other family of distributions as well.

The ratio of normal tail probabilities and the ratio of Student’s t tail probabilities have gained an increased attention in statistics and related areas. However, they are not well studied in the literature. In this paper, we systematically study the functional behaviors of these two ratios. Meanwhile, we explore their difference as well as their relationship. It is surprising that the two ratios behave very different to each other. Finally, we conclude the paper by conducting some lower and upper bounds for the two ratios.

This paper considers the problem of variable selection for a class of linear regression models with instrumental variables. We focus on the case that the covariates are endogenous variables, and some auxiliary instrumental variables are available. An instrumental variable based variable selection procedure is proposed by using a modified smooth-threshold estimating equations (SEE). The proposed method can attenuate the effect of endogeneity of covariates, and can avoid the convex optimization problem. Hence, it is flexible and easy to implement. Simulation results indicate that the proposed variable selection method is workable.

In this paper, we extend the AENP and the GMC criterion proposed by Zhang et al. (2008) to the case of nonregular orthogonal designs. A G-AENP and correspondingly a G-GMC criterion are proposed. The confounding frequency vector (CFV) and the generalized wordlength pattern (GWLP), as the base of MGA and GMA criteria, are shown to be functions of the G-AENP. Some optimal properties of G-GMC designs are obtained. At the last, we give an efficient algorithm for finding optimal designs and tabulate some G-GMC designs with 16- and 18-run for application and comparison with GMA and MGA designs.

In this paper we study a smooth estimator of the regression quantile function in the censorship model when the covariates take values in some abstract function space. The main goal of this paper is to establish the asymptotic normality of the kernel estimator of the regression quantile, under α-mixing condition and, on the concentration properties on small balls probability measure of the functional regressors. Some applications and particular cases are studied. This study can be applied in time series analysis to the prediction and building confidence bands. Some simulations are drawn to lend further support to our theoretical results and to compare the quality of behavior of the estimator for finite samples with different rates of censoring and sizes.

The paper considers a lognormal model for the survival times and obtains a Bayes solution by means of Gibbs sampler algorithm when the priors for the parameters are vague. The formulation given in the paper is mainly focused for censored data problems though it is equally well applicable for complete data scenarios as well. For the purpose of numerical illustration, we considered two real data sets on head and neck cancer patients when they have been treated using either radiotherapy or chemotherapy followed by radiotherapy. The paper not only compares the survival functions for the two therapies assuming a lognormal model but also provides a model compatibility study based on predictive simulation results so that the choice of lognormal model can be justified for the two data sets. The ease of our analysis as compared to an earlier approach is certainly an advantage.

The common parameter estimation methods of panel data linear model include least square dummy variable estimation, two-stage least square estimation, quasi-maximum likelihood estimation and generalized moment estimation. However, these estimation methods are not robust and are easily affected by outliers. Firstly, this paper extends support vector regression algorithm to fit several parallel super-plane simultaneously and provide a novel robust estimation of fixed-effect panel data linear model; then using the kernel trick, a robust estimation for fixed effect panel data nonlinear model is introduced. Finally, the proposed model (linear or nonlinear) is applied in forecasting air quality index of the cities of Jing-Jin-Ji district in China. Experiments shows that our proposed model are robust and have good generalization performance.

Linear regression models which account for skewed error distributions with fat tails have been previously studied. These two important features, skewness, and fat tails, are often observed in real data analyses. Covariates measured with an error also happen frequently in the observational data set-up. As a motivating example, wind speed as a covariate is usually used, among other covariates, to estimate the particulate matter (PM) which is one of the most critical air pollutants and has a major impact on human health and on the environment. However, the wind speed is measured with error and the distribution of PM is neither symmetric nor normally distributed (see Section “PM data application in Canada” for more details). Ignoring the issue of measurement error in covariates may produce bias in model parameters estimate and lead to wrong conclusions. In this paper, we propose an approach to study properly linear regression models where the covariates are measured with error and the error distribution is skewed with fat tails. We use a hierarchical Bayesian approach for inference, addressing also sensitivity of the results to priors. Performance of the proposed approach is evaluated through a simulation study and also by a real data application (PM in Canada).