On hysteretic vector autoregressive model with applications

This paper proposes a new hysteretic vector autoregressive (HVAR) model in which the regime switching may be delayed when the hysteresis variable lies in a hysteresis zone. We integrate an adapted multivariate Student-t distribution from amending the scale mixtures of normal distributions. This HVAR model allows for a higher degree of flexibility in the degrees of freedom for each time series. We use the proposed model to test for a causal relationship between any two target time series. Using posterior odds ratios, we overcome the limitations of the classical approach to multiple testing. Both simulated and real examples herein help illustrate the suggested methods. We apply the proposed HVAR model to investigate the causal relationship between the quarterly growth rates of gross domestic product of United Kingdom and United States. Moreover, we check the pairwise lagged dependence of daily PM2.5 levels in three districts of Taipei.     

Accurate computation of single and product moments of order statistics under progressive censoring

An accurate procedure is proposed to calculate approximate moments of progressive order statistics in the context of statistical inference for lifetime models. The study analyses the performance of power series expansion to approximate the moments for location and scale distributions with high precision and smaller deviations with respect to the exact values. A comparative analysis between exact and approximate methods is shown using some tables and figures. The different approximations are applied in two situations. First, we consider the problem of computing the large sample variance–covariance matrix of maximum likelihood estimators. We also use the approximations to obtain progressively censored sampling plans for log-normal distributed data. These problems illustrate that the presented procedure is highly useful to compute the moments with precision for numerous censoring patterns and, in many cases, is the only valid method because the exact calculation may not be applicable.     

Progressive censoring with fixed censoring times

In this paper, we revisit the progressive Type-I censoring scheme as it has originally been introduced by Cohen [Progressively censored samples in life testing. Technometrics. 1963;5(3):327–339]. In fact, original progressive Type-I censoring proceeds as progressive Type-II censoring but with fixed censoring times instead of failure time based censoring times. Apparently, a time truncation has been added to this censoring scheme by interpreting the final censoring time as a termination time. Therefore, not much work has been done on Cohens's original progressive censoring scheme with fixed censoring times. Thus, we discuss distributional results for this scheme and establish exact distributional results in likelihood inference for exponentially distributed lifetimes. In particular, we obtain the exact distribution of the maximum likelihood estimator (MLE). Further, the stochastic monotonicity of the MLE is verified in order to construct exact confidence intervals for both the scale parameter and the reliability.     

Design selection for strong orthogonal arrays

Strong orthogonal arrays (SOAs) were recently introduced and studied as a class of space‐filling designs for computer experiments. An important problem that has not been addressed in the literature is that of design selection for such arrays. In this article, we conduct a systematic investigation into this problem, and we focus on the most useful SOA(n,m,4,2 + )s and SOA(n,m,4,2)s. This article first addresses the problem of design selection for SOAs of strength 2+ by examining their three‐dimensional projections. Both theoretical and computational results are presented. When SOAs of strength 2+ do not exist, we formulate a general framework for the selection of SOAs of strength 2 by looking at their two‐dimensional projections. The approach is fruitful, as it is applicable when SOAs of strength 2+ do not exist and it gives rise to them when they do. The Canadian Journal of Statistics 47: 302–314; 2019 © 2019 Statistical Society of Canada     

On the power of goodness-of-fit tests for the exponential distribution under progressive Type-II censoring

The aim of this article is to review existing goodness-of-fit tests for the exponential distribution under progressive Type-II censoring and to provide some new ideas and adjustments. In particular, we consider two-parameter exponentially distributed random variables and adapt the proposed test procedures to our scenario if necessary. Then, we compare their power by an extensive simulation study. Furthermore, we propose five new test procedures that provide reasonable alternatives to those already known.     

On MLEs of the parameters of a modified Weibull distribution based on record values

The modified Weibull distribution can be used quite effectively to model complex data from mechanical engineering or survival analysis studies that posses a monotonic or a bathtub-shape hazard rate. In this paper, we study the MLEs of the parameters of a modified Weibull distribution model in the presence of upper kth record values. The existence and uniqueness of the MLEs are proven in this case. Real data analysis is performed for illustrative purposes.     

Inference and optimal design of multiple constant-stress testing for generalized half-normal distribution under type-II progressive censoring

The generalized half-normal (GHN) distribution and progressive type-II censoring are considered in this article for studying some statistical inferences of constant-stress accelerated life testing. The EM algorithm is considered to calculate the maximum likelihood estimates. Fisher information matrix is formed depending on the missing information law and it is utilized for structuring the asymptomatic confidence intervals. Further, interval estimation is discussed through bootstrap intervals. The Tierney and Kadane method, importance sampling procedure and Metropolis-Hastings algorithm are utilized to compute Bayesian estimates. Furthermore, predictive estimates for censored data and the related prediction intervals are obtained. We consider three optimality criteria to find out the optimal stress level. A real data set is used to illustrate the importance of GHN distribution as an alternative lifetime model for well-known distributions. Finally, a simulation study is provided with discussion.     

Confidence regions for Pareto parameters from a single and independent samples

[Abstract] Based on a single and on two independent samples, joint confidence regions for parameters of Pareto distributions are proposed with minimum volume properties and without assigning the confidence level to dimensions. In the one-sample case, comparisons are made to former simultaneous confidence sets for Pareto parameters by means of simulation and a real data set. The two-sample case is studied in various set-ups and comprises simultaneous confidence regions for the shape parameters, the scale parameters, and higher-dimensional vectors of these parameters, where common shape and common scale models are also considered.     

Modeling spatial overdispersion via the generalized Waring point process

Modeling spatial overdispersion requires point process models with finite‐dimensional distributions that are overdisperse relative to the Poisson distribution. Fitting such models usually heavily relies on the properties of stationarity, ergodicity, and orderliness. In addition, although processes based on negative binomial finite‐dimensional distributions have been widely considered, they typically fail to simultaneously satisfy the three required properties for fitting. Indeed, it has been conjectured by Diggle and Milne that no negative binomial model can satisfy all three properties. In light of this, we change perspective and construct a new process based on a different overdisperse count model, namely, the generalized Waring (GW) distribution. While comparably tractable and flexible to negative binomial processes, the GW process is shown to possess all required properties and additionally span the negative binomial and Poisson processes as limiting cases. In this sense, the GW process provides an approximate resolution to the conundrum highlighted by Diggle and Milne.     

Tempered Mittag-Leffler Lévy processes

In this article, we introduce tempered Mittag-Leffler Lévy processes (TMLLP). TMLLP is represented as tempered stable subordinator delayed by a gamma process. Its probability density function and Lévy density are obtained in terms of infinite series and Mittag-Leffler function, respectively. Asymptotic forms of the tails and moments are given. A step-by-step procedure of the parameters estimation and simulation of sample paths is given. We also provide main results available for Mittag-Leffler Lévy processes (MLLP) and some extensions which are not available in a collective way in a single article. Our results generalize and complement the results available on Mittag-Leffler distribution and MLLP in several directions. Further, the asymptotic forms of the moments of the first-exit times of the TMLLP are also discussed.