Search

Sequential change-point analysis, which identifies a change of probability distribution in an infinite sequence of random observations, has important applications in many fields. A good method should detect a change point as soon as possible, and keep a low amount of false alarms. As one of the most popular methods, Shiryaev-Roberts (SR) procedure holds many optimalities. However, its implementation requires the pre-change and post-change distributions to be known, which is not achievable in practice. In this paper, we construct a nonparametric version of the SR procedure by embedding different versions of empirical likelihood, assuming two training samples, before and after change, are available for parameter estimations. Simulations are conducted to compare the performance of the proposed method with existing methods. The results show that when the underlying distribution is unknown, and training sample sizes are small, the proposed modified procedure shows advantage by giving a smaller delay of detection. [ABSTRACT FROM AUTHOR]

Nonparametric tests for equality of multivariate distributions are frequently desired in research. It is commonly required that test-procedures based on relatively small samples of vectors accurately control the corresponding Type I Error (TIE) rates. Often, in the multivariate testing, extensions of null-distribution-free univariate methods, e.g., Kolmogorov-Smirnov and Cramér-von Mises type schemes, are not exact, since their null distributions depend on underlying data distributions. The present paper extends the density-based empirical likelihood technique in order to nonparametrically approximate the most powerful test for the multivariate two-sample (MTS) problem, yielding an exact finite-sample test statistic. We rigorously apply one-to-one-mapping between the equality of vectors' distributions and the equality of distributions of relevant univariate linear projections. We establish a general algorithm that simplifies the use of projection pursuit, employing only a few of the infinitely many linear combinations of observed vectors' components. The displayed distribution-free strategy is employed in retrospective and group sequential manners. A novel MTS nonparametric procedure in the group sequential manner is proposed. The asymptotic consistency of the proposed technique is shown. Monte Carlo studies demonstrate that the proposed procedures exhibit extremely high and stable power characteristics across a variety of settings. Supplementary materials for this article are available online. [ABSTRACT FROM AUTHOR]

In this paper, we introduce a new distribution defined on $ \mathbb {Z} $ Z , called the $ {{ZPMOIEP}} $ ZPMOIEP distribution, which can be viewed as a natural extension of the zero-and-one-inflated Poisson ( $ {ZOIP} $ ZOIP ) distribution. It is designed to fit the count data with potentially excess zeros and/or ones, and/or minus ones. We explore its various properties and investigate the estimation of the unknown parameters. Moreover, simulation experiments are carried out to attest to the performance of the estimation. Through the use of a useful data set on football scores, the applicability of the proposed distribution is examined. [ABSTRACT FROM AUTHOR]

It is of great interest to compute optimal exact confidence intervals for the success probability (p) in a binomial distribution, the number of subjects with a certain attribute (M) or the total number of subjects (N) in a hypergeometric distribution, and the mean λ of a Poisson distribution. In this paper, efficient algorithms are proposed to compute an admissible exact interval for each of the four parameters when the sample size (n) or the random observation X is large. The algorithms are utilized in four practical examples: evaluating the relationship between two diseases, certifying companies, estimating the proportion of drug users, and analyzing earthquake frequency. The intervals computed by the algorithms are shorter, and the calculations are faster, demonstrating the accuracy of the results and the time efficiency of the proposed algorithms. [ABSTRACT FROM AUTHOR]

The exponential distribution can provide a simple and appealing survival-time model in reliability analysis and life tests. Due to certain experimental designs, time limitations, budgetary requirements, and other reasons, however, only censored data can be obtained. In this paper, we propose a novel estimation procedure for the mean survival time of an exponential distribution under a sequential censoring scheme, which can be treated as a combination of type I censoring and type II censoring. The procedure makes a trade-off between estimation error and sampling cost, using the minimum number of observations. An extensive set of Monte Carlo simulations is conducted to further validate its remarkable performance. To demonstrate the practical applicability, we then implement this newly proposed procedure to assess the reliability of various Backblaze's hard disk models. [ABSTRACT FROM AUTHOR]

In this paper we introduce a discrete-time and continuous state-space Markov stationary process { X n ; n = 1 , 2 , ... } , where X n has a two-parameter Weibull distribution, X n 's are dependent and there is a positive probability that X n = X n + 1 . The motivation came from the gold price data where there are several instances for which X n = X n + 1 . Hence, the existing methods cannot be used to analyze this data. We derive different properties of the proposed Weibull process. It is observed that the joint cumulative distribution function of X n and X n + 1 has a very convenient copula structure. Hence, different dependence properties and dependence measures can be obtained. The maximum likelihood estimators cannot be obtained in explicit forms, we have proposed a simple profile likelihood method to compute these estimators. We have used this model to analyze two synthetic data sets and one gold price data set of the Indian market, and it is observed that the proposed model fits quite well with the data set. [ABSTRACT FROM AUTHOR]

Comparative lifetime experiments are of particular importance in production processes when one wishes to determine the relative merits of several competing products with regard to their reliability. This paper confines itself to the data obtained by running a joint progressive Type-II censoring plan on samples in a combined manner. The problem of Bayesian predicting failure times of surviving units is discussed in details when parent populations are exponential. Two real data sets are analyzed in order to illustrate all the inferential procedures developed here. When destructive experiments under a censoring scheme finished, the researchers are usually interested to estimate remaining lifetimes of surviving units for sequel experiments. Findings of this paper are useful for these purposes specially when samples are non-homogeneous such as those taken from industrial storages. [ABSTRACT FROM AUTHOR]

In various scenarios where products and services are accompanied by warranties to ensure their reliability over a specified time, the two-parameter (shifted) exponential distribution serves as a fundamental model for time-to-event data. In modern production process, the products often come with warranties, and their quality can be manifested by the changes in the scale and origin parameters of a shifted exponential (SE) distribution. This paper introduces the Max-EWMA chart, employing maximum likelihood estimators and exponentially weighted moving average (EWMA) statistics, to jointly monitor SE distribution parameters. Additionally, we extend two additional charts, namely the Max-DEWMA and Max-TEWMA charts to enhance early-stage shift detection. Performance evaluations under zero-state and steady-state conditions compare these charts with the existing Max-CUSUM chart in terms of expected value and standard deviation of the run length (RL) distribution. Our findings reveal that among the Max-EWMA schemes, the Max-EWMA SE chart outperforms the others in terms of steady-state performance, while the Max-TEWMA chart surpasses the Max-EWMA and Max-DEWMA SE charts in respect to zero-state performance. Moreover, the proposed Max-EWMA schemes demonstrate advantages over Max-CUSUM, especially for small to moderate smoothing constants. We also provide an illustrative example to demonstrate the implementation of the proposed schemes. [ABSTRACT FROM AUTHOR]

Proportional data arise frequently in a wide variety of fields of study. Such data often exhibit extra variation such as over/under dispersion, sparseness and zero inflation. For example, the hepatitis data present both sparseness and zero inflation with 19 contributing non-zero denominators of 5 or less and with 36 having zero seropositive out of 83 annual age groups. The whitefly data consists of 640 observations with 339 zeros (53%), which demonstrates extra zero inflation. The catheter management data involve excessive zeros with over 60% zeros averagely for outcomes of 193 urinary tract infections, 194 outcomes of catheter blockages and 193 outcomes of catheter displacements. However, the existing models cannot always address such features appropriately. In this paper, a new two-parameter probability distribution called Lindley–binomial (LB) distribution is proposed to analyze the proportional data with such features. The probabilistic properties of the distribution such as moment, moment generating function are derived. The Fisher scoring algorithm and EM algorithm are presented for the computation of estimates of parameters in the proposed LB regression model. The issues on goodness of fit for the LB model are discussed. A limited simulation study is also performed to evaluate the performance of derived EM algorithms for the estimation of parameters in the model with/without covariates. The proposed model is illustrated through three aforementioned proportional datasets. [ABSTRACT FROM AUTHOR]

Comparative lifetime experiments are remarkable when the study is to ascertain the relative merits of two competing products regarding the duration of their service life. This paper considers the comparative lifetime experiments of two Gompertz populations under a balanced joint progressive Type-II censoring scheme. The lifetime distributions of the units are assumed to follow the Gompertz distribution with a common shape but different scale parameters. The maximum likelihood estimates of the unknown parameters are derived. The existence of the maximum likelihood estimates is proved. Expectation-maximization and stochastic expectation-maximization algorithms are provided to calculate the estimates. The bootstrap-p, bootstrap-t, and approximate confidence intervals are established. To obtain the Bayesian estimates, it is assumed that the prior of scale parameters is a Beta-Gamma distribution and the prior of the common shape parameter is an independent Gamma distribution. Under squared error loss and LINEX loss functions, the Metropolis-Hastings algorithm is provided to compute the Bayes estimates and the credible intervals. Further, the statistical inferences with order restriction are studied when it is known a priori that the expectation of the lifespan of one population is shorter than that of the other population. A wide range of simulation experiments is conducted to evaluate the performance of the proposed methods. Finally, the lifetimes of white organic light-emitting diodes and the breaking strengths of jute fiber of gauge lengths are analyzed to illustrate the practical application of the proposed model and methods. [ABSTRACT FROM AUTHOR]

The modality is an important topic for modelling. Using parametric models is an efficient way when real data set shows trimodality. In this paper, we propose a new class of trimodal probability distributions, that is, probability distributions that have up to three modes. Trimodality itself is achieved by applying a proper transformation to density function of certain continuous probability distributions. At first, we obtain preliminary results for an arbitrary density function $ g(x) $ g (x) and, next, we focus on the Gaussian case, studying trimodal Gaussian model more deeply. The Gaussian distribution is applied to produce the trimodal form of Gaussian known as normal distribution. The tractability of analytical expression of normal distribution and properties of the trimodal normal distribution are important reasons why we choose normal distribution. Furthermore, the existing distributions should be improved to be capable of modelling efficiently when there exists a trimodal form in a data set. After new density function is proposed, estimating its parameters is important. Since Mathematica 12.0 software has optimization tools and important modelling techniques, computational steps are performed using this software. The bootstrapped form of real data sets are applied to show the modelling ability of the proposed distribution when real data sets show trimodality. [ABSTRACT FROM AUTHOR]

This paper introduces a new family of quantile regression models whose response variable follows a reparameterized Marshall-Olkin distribution indexed by quantile, scale, and asymmetry parameters. The family has arisen by applying the Marshall-Olkin approach to distributions belonging to the location-scale family. Models of higher flexibility and whose structure is similar to generalized linear models were generated by quantile reparameterization. The maximum likelihood (ML) method is presented for the estimation of the model parameters, and simulation studies evaluated the performance of the ML estimators. The advantages of the family are illustrated through an application to a set of nutritional data, whose results indicate it is a good alternative for modeling slightly asymmetric response variables with support on the real line. [ABSTRACT FROM AUTHOR]

In this paper, we propose two multimodal circular distributions which are suitable for modeling circular data sets with two or more modes. Both distributions belong to the regular exponential family of distributions and are considered as extensions of the von Mises distribution. Hence, they possess the highly desirable properties, such as the existence of non-trivial sufficient statistics and optimal inferences for their parameters. Fine particulates (PM2.5) are generally emitted from activities such as industrial and residential combustion and from vehicle exhaust. We illustrate the utility of our proposed models using a real data set consisting of fine particulates (PM2.5) pollutant levels in Houston region during Fall season in 2019. Our results provide a strong evidence that its diurnal pattern exhibits four modes; two peaks during morning and evening rush hours and two peaks in between. [ABSTRACT FROM AUTHOR]

The extreme value distribution was developed for modeling extreme-order statistics or extreme events. In this study, we discuss the distribution of the largest extreme. The main objective of this paper is to determine the best estimators of the unknown parameters of the extreme value distribution. Thus, both classical and Bayesian methods are used. The classical estimation methods under consideration are maximum likelihood estimators, moment's estimators, least squares estimators, and weighted least squares estimators, percentile estimators, the ordinary least squares estimators, best linear unbiased estimators, L-moments estimators, trimmed L-moments estimators, and Bain and Engelhardt estimators. We also propose new estimators for the unknown parameters. Bayesian estimators of the parameters are derived by using Lindley's approximation and Markov Chain Monte Carlo methods. The asymptotic confidence intervals are considered by using maximum likelihood estimators. The Bayesian credible intervals are also obtained by using Gibbs sampling. The performances of these estimation methods are compared with respect to their biases and mean square errors through a simulation study. The maximum daily flood discharge (annual) data sets of the Meriç River and Feather River are analyzed at the end of the study for a better understanding of the methods presented in this paper. [ABSTRACT FROM AUTHOR]

This paper focuses on contingency tables (CTs) used to accept or reject hypothesis that two categorical features are independent. This paper promotes a new view of the essence of CTs. Particular CTs are created as a result of phenomena governed by random factors. We put forward four simple scenarios under which CTs may be created. Further in the paper, a measure of untruthfulness of , which also may serve as the test statistics, is defined. The paper recalls several CT-related test statistics. Each test statistics produces its own power-of-the-test (PoT) function. The PoTs were determined with the Monte Carlo method and then compared to each other. In order to ease such a PoT-to-PoT comparison, the measure named 'average outgoing untruthfulness level' is put forward. Derivation, properties and example how to apply the measure in practice are given in details. [ABSTRACT FROM AUTHOR]

Directional data appears in several branches of research. In some cases, those directional variables are only defined in subsets of the K-dimensional unit sphere. For example, in some applications, angles as measured responses are limited on the positive orthant. Analysis on subsets of the K-dimensional unit sphere is challenging and nowadays there are not many proposals that discuss this topic. Thus, from a methodological point of view, it is important to have probability distributions defined on bounded subsets of the K-dimensional unit sphere. Specifically, in this paper, we introduce a nonparametric Bayesian model to describe directional variables restricted to the first orthant. This model is based on a Dirichlet process mixture model with multivariate projected Gamma densities as kernel distributions. We show how to carry out inference for the proposed model based on a slice sampling scheme. The proposed methodology is illustrated using simulated data sets as well as a real data set. [ABSTRACT FROM AUTHOR]

Gompertz distribution is a significant and commonly used lifetime distribution, which plays an important role in reliability engineering. In this paper, we study the statistical inference of Gompertz distribution based on adaptive Type-II hybrid progressive censored schemes. From the perspective of frequentist, we derive the point estimations through the method of maximum likelihood estimation (MLE) and the existence of MLE is proved. Besides MLE, we propose the stochastic EM algorithm to reduce complexity and simplify computing. We also apply the method of Bootstraps (Bootstrap-p and Bootstrap-t) to construct confidence intervals. From Bayesian aspect, the Bayes estimates of the unknown parameters are evaluated by applying the MCMC method, the average length and coverage rate of credible intervals are also carried out. The Bayes inference is based on the squared error loss function and LINEX loss function. Furthermore, a numerical simulation is conducted to assess the performance of the proposed methods. Finally, a real-life example is considered to illustrate the application and development of the inference methods. In summary, the Bayesian method seems to perform the best among all approaches, while other approaches also present different advantages. [ABSTRACT FROM AUTHOR]

Statistical modelling of sports data is indispensable to analyse the sports behaviour and apprehend significant inferences that are helpful to adopt decisive strategies before or during the sports events. This paper introduces a stochastic model as the distribution of difference derived from the Bivariate Affine-Linear Exponential distribution. The distribution of difference is first ever used to model the margin of victory that provides an adequate fitting on the observed data. A simulation study is carried out to observe the stability of the model parameters through their average estimated values, biases, standard errors, root mean square errors and confidence intervals. The performance of the proposed model is examined by applying it on the real data of the National Football League and comparing the results with those of the existing models. Finally, the quantile function of the proposed distribution is used to assess the possible range of point spreads for winning the bet in a particular game. [ABSTRACT FROM AUTHOR]

In this paper, we introduce a new lifetime distribution as an alternative to generalized Gompertz, Gompertz distribution and its modified ones. This new distribution is a special case of the family of distributions introduced by Granzotto et al. [D.C.T. Granzotto, F. Louzada and N. Balakrishnan, Cubic rank transmuted distributions: inferential issues and applications., J. Stat. Comput. Simul. 87 (2017), pp. 2760–2778]. We obtain some characteristic properties of suggested distribution such as hazard function, ordinary moments, coefficient of skewness, coefficient of kurtosis, moment generating function, quantile function and median. We discuss three different methods of estimation to estimate the parameters of proposed distribution. A comprehensive Monte Carlo simulation study is performed in order to compare the performances of estimators according to mean square errors and biases. Finally, three real data applications are performed to illustrate usefulness of suggested distribution. [ABSTRACT FROM AUTHOR]

In this paper, a new bivariate discrete generalized exponential distribution, whose marginals are discrete generalized exponential distributions, is studied. It is observed that the proposed bivariate distribution is a flexible distribution whose cumulative distribution function has an analytical structure. In addition, a new bivariate geometric distribution can be obtained as a special case. We study different properties of this distribution and propose estimation of its parameters. We will see that the maximum of the variables involved in the proposed bivariate distribution defines some new classes of univariate discrete distributions, which are interesting in their own sake, and can be used to analyze some Reliability systems whose components are positive dependent. Some important futures of this new univariate family of discrete distributions are also studied in details. In addition, a general class of bivariate discrete distributions, whose marginals are exponentiated discrete distributions, is introduced. Moreover, the analysis of two real bivariate data sets is performed to indicate the effectiveness of the proposed models. Finally, we conclude the paper. [ABSTRACT FROM AUTHOR]

The widely used Fellegi–Sunter model for probabilistic record linkage does not leverage information contained in field values and consequently leads to identical classification of match status regardless of whether records agree on rare or common values. Since agreement on rare values is less likely to occur by chance than agreement on common values, records agreeing on rare values are more likely to be matches. Existing frequency-based methods typically rely on knowledge of error probabilities associated with field values and frequencies of agreed field values among matches, often derived using prior studies or training data. When such information is unavailable, applications of these methods are challenging. In this paper, we propose a simple two-step procedure for frequency-based matching using the Fellegi–Sunter framework to overcome these challenges. Matching weights are adjusted based on frequency distributions of the agreed field values among matches and non-matches, estimated by the Fellegi–Sunter model without relying on prior studies or training data. Through a real-world application and simulation, our method is found to produce comparable or better performance than the unadjusted method. Furthermore, frequency-based matching provides greater improvement in matching accuracy when using poorly discriminating fields with diminished benefit as the discriminating power of matching fields increases. [ABSTRACT FROM AUTHOR]

Although there have been a number of available tests of bivariate exchangeability, i.e. bivariate symmetry for bivariate distributions, the literature is void of tests whether a multivariate distribution with more than two dimensions is exchangeable or not. In this paper, multivariate permutation tests of exchangeability of multivariate distributions are proposed, which are based on the non-parametric combination methodology, i.e. on combining non-parametric bivariate exchangeability tests. Numerical experiments on real as well as simulated multivariate data with more than two dimensions are presented here. The multivariate permutation test turns out to be typically more powerful than a bivariate exchangeability test performed only over a single pair of variables, and also more suitable compared to tests exploiting the approaches of Benjamini–Yekutieli or Bonferroni. [ABSTRACT FROM AUTHOR]

In this paper, we introduced a new probability distribution named as inverse power Maxwell distribution. The proposal distribution can be seen as an extension of the Maxwell distribution with more flexibility in modeling upside-down lifetime data. Some statistical properties of this distribution are derived. In estimation viewpoint, five methods are used for estimating the unknown parameters of the distribution and these methods are performed through the simulation study. Finally, two real data sets were analyzed to illustrate the applicability of the proposed distribution, proving that it fits each real data set much better than some other existing distributions. [ABSTRACT FROM AUTHOR]

The aim of this paper is to compare passenger (pax) demand between airports based on the arithmetic mean (MPD) and the median pax demand (MePD). A three phases approach is applied. First phase, we use bootstrap procedures to estimate the distribution of the arithmetic MPD and the MePD for each block of routes distance; second phase, we use percentile, standard, bias corrected, and bias corrected accelerated methods to calculate bootstrap confidence bands for the MPD and the MePD; and third phase, we implement Monte Carlo (MC) experiments to analyse the finite sample performance of the applied bootstrap. Our results conclude that it is more meaningful to use the estimation of MePD rather than the estimation of MPD in the air transport industry. By carrying out MC experiments, we demonstrate that the bootstrap methods produce coverages close to the nominal for the MPD and the MePD. [ABSTRACT FROM AUTHOR]

In this paper, we first present a characterization of exponential distribution based on the extropy of record values and next introduce a goodness-of-fit test for exponentiality. Monte Carlo simulation is used to compute the critical values of our proposed test for different sample sizes and significance levels. To show the advantage of the proposed test, we adopt 58 competitor tests and compute the adjusted power against different alternatives with distinct types of hazard function. The power results show that our proposed test has superior adjusted power if the alternatives have increasing failure rates or bathtub decreasing-increasing failure rates, especially when the sample size is small. Finally, three real examples are used to illustrate the applicability and robustness of our proposed test by monitoring the p-values of the tests. [ABSTRACT FROM AUTHOR]

Parametric methods for the calculation of reference intervals in clinical studies often rely on the identification of a suitable transformation so that the transformed data can be assumed to be drawn from a Gaussian distribution. In this paper, the two-stage transformation recommended by the International Federation for Clinical Chemistry is compared with a novel generalised Box–Cox family of transformations. Investigation is also made of sample sizes needed to achieve certain criteria of reliability in the calculated reference interval. Simulations are used to show that the generalised Box–Cox family achieves a lower bias than the two-stage transformation. It was found that there is a possibility that the two-stage transformation will result in percentile estimates that cannot be back-transformed to obtain the required reference intervals, a difficulty not observed when using the generalised Box–Cox family introduced in this paper. [ABSTRACT FROM PUBLISHER]

The aim of this paper is to provide some practical aspects of point and interval estimates of the global maximum of a function using extreme value theory. Consider a real-valued function f:D→ defined on a bounded interval D such that f is either not known analytically or is known analytically but has rather a complicated analytic form. We assume that f possesses a global maximum attained, say, at u*∈D with maximal value x*=max u f(u)≐f(u*). The problem of seeking the optimum of a function which is more or less unknown to the observer has resulted in the development of a large variety of search techniques. In this paper we use the extreme-value approach as appears in Dekkers et al. [A moment estimator for the index of an extreme-value distribution, Ann. Statist. 17 (1989), pp. 1833-1855] and de Haan [Estimation of the minimum of a function using order statistics, J. Amer. Statist. Assoc. 76 (1981), pp. 467-469]. We impose some Lipschitz conditions on the functions being investigated and through repeated simulation-based samplings, we provide various practical interpretations of the parameters involved as well as point and interval estimates for x*. [ABSTRACT FROM AUTHOR]

This paper examines a number of statistics that have been proposed to test the normality assumption in the tobit (censored regression) model. It argues that a number of commonly proposed statistics can be interpreted as different versions of the Lagrange multiplier, or score, test for a common null hypothesis. This observation is useful in examining the Monte Carlo results presented in the paper. The Monte Carlo results suggest that the computational convenience of a number of statistics is obtained at the cost of poor finite sample performance under the null hypothesis. [ABSTRACT FROM AUTHOR]

The coefficients of regression are usually estimated for minimization problems with asymmetric loss functions. In this paper, we rather correct predictions so that the prediction error follows a generalized Gaussian distribution. In our method, we not only minimize the expected value of the asymmetric loss, but also lower the variance of the loss. Predictions usually have errors. Therefore, it is necessary to use predictions in consideration of these errors. Our approach takes into account prediction errors. Furthermore, even if we do not understand the prediction method, which is a possible circumstance in, e.g. deep learning, we can use our method if we know the prediction error distribution and asymmetric loss function. Our method can be applied to procurement of electricity from electricity markets. [ABSTRACT FROM AUTHOR]

In this paper, we consider two well-known parametric long-term survival models, namely, the Bernoulli cure rate model and the promotion time (or Poisson) cure rate model. Assuming the long-term survival probability to depend on a set of risk factors, the main contribution is in the development of the stochastic expectation maximization (SEM) algorithm to determine the maximum likelihood estimates of the model parameters. We carry out a detailed simulation study to demonstrate the performance of the proposed SEM algorithm. For this purpose, we assume the lifetimes due to each competing cause to follow a two-parameter generalized exponential distribution. We also compare the results obtained from the SEM algorithm with those obtained from the well-known expectation maximization (EM) algorithm. Furthermore, we investigate a simplified estimation procedure for both SEM and EM algorithms that allow the objective function to be maximized to split into simpler functions with lower dimensions with respect to model parameters. Moreover, we present examples where the EM algorithm fails to converge but the SEM algorithm still works. For illustrative purposes, we analyze a breast cancer survival data. Finally, we use a graphical method to assess the goodness-of-fit of the model with generalized exponential lifetimes. [ABSTRACT FROM AUTHOR]

It is well-known that classical Tobit estimator of the parameters of the censored regression (CR) model is inefficient in case of non-normal error terms. In this paper, we propose to use the modified maximum likelihood (MML) estimator under the Jones and Faddy's skew t-error distribution, which covers a wide range of skew and symmetric distributions, for the CR model. The MML estimators, providing an alternative to the Tobit estimator, are explicitly expressed and they are asymptotically equivalent to the maximum likelihood estimator. A simulation study is conducted to compare the efficiencies of the MML estimators with the classical estimators such as the ordinary least squares, Tobit, censored least absolute deviations and symmetrically trimmed least squares estimators. The results of the simulation study show that the MML estimators work well among the others with respect to the root mean square error criterion for the CR model. A real life example is also provided to show the suitability of the MML methodology. [ABSTRACT FROM AUTHOR]

In this paper, a new flexible generator of distributions is proposed. Some of its fundamental properties including quantile, skewness, kurtosis, hazard rate function, moments, mean deviations, mean time to failure, mean time between failure, availability and reliability function of consecutive linear and circular systems are studied. The hazard rate function can be increasing, decreasing, unimodal-bathtub, unimodal, bathtub, J and inverse J-shaped depending on its parameters values. After introducing the general class, two special models of the new family are discussed in detail. Maximum likelihood and Bayesian methods are used to estimate the model parameters. A detailed simulation study is carried out to examine the bias and mean square error of maximum likelihood and Bayesian estimators. We also illustrate the importance of the new family by means of two distinctive real data sets. It can serve as an alternative model to other lifetime distributions in the existing statistical literature for modeling positive and negative real data in many areas. [ABSTRACT FROM AUTHOR]

The first aim of this paper is to introduce a modular test for the three-way contingency table (TT). The second aim is to describe the procedure of generating TT using the bar method. The third aim is on the one hand to suggest the measure of untruthfulness of H0 and on the other hand to compare the quality of independence tests by using their power. Critical values for analyzed statistics were determined by simulating the Monte Carlo method. [ABSTRACT FROM AUTHOR]

In this paper we introduce a new class of multivariate unimodal distributions, motivated by Khintchine's representation for unimodal densities on the real line. We start by introducing a new class of unimodal distributions which can then be naturally extended to higher dimensions, using the multivariate Gaussian copula. Under both univariate and multivariate settings, we provide MCMC algorithms to perform inference about the model parameters and predictive densities. The methodology is illustrated with univariate and bivariate examples, and with variables taken from a real data set. [ABSTRACT FROM AUTHOR]

In this paper, the problem of estimating the mean vector under non-negative constraints on location vector of the multivariate normal distribution is investigated. The value of the wavelet threshold based on Stein's unbiased risk estimators is calculated for the shrinkage estimator in restricted parameter space. We suppose that covariance matrix is unknown and we find the dominant class of shrinkage estimators under Balance loss function. The performance evaluation of the proposed class of estimators is checked through a simulation study by using risk and average mean square error values. [ABSTRACT FROM AUTHOR]

In this paper, we study the problem of network discovery and influence propagation, and propose an integrated approach for the analysis of lead-lag synchronization in multiple choices. Network models for the processes by which decisions propagate through social interaction have been studied before, but only a few consider unknown structures of interacting agents. In fact, while individual choices are typically observed, inferring individual influences – who influences who – from sequences of dynamic choices requires strong modeling assumptions on the cross-section dependencies of the observed panels. We propose a class of parametric models which extends the vector autoregression to the case of pairwise influences between individual choices over multiple items and supports the analysis of influence propagation. After uncovering a collection of theoretical properties (conditional moments, parameter sensitivity, identifiability and estimation), we provide an economic application to music broadcasting, where a set of songs are diffused over radio stations; we infer station-to-station influences based on the proposed methodology and assess the propagation effect of initial launching stations to maximize songs diffusion. Both on the theoretical and empirical sides, the proposed approach connects fields which are traditionally treated as separated areas: the problem of network discovery and the one of influence propagation. [ABSTRACT FROM AUTHOR]

Count time series data with excess zeros are observed in several applied disciplines. When these zero-inflated counts are sequentially recorded, they might result in serial dependence. Ignoring the zero-inflation and the serial dependence might produce inaccurate results. In this paper, Markov zero-inflated count time series models based on a joint distribution on consecutive observations are proposed. The joint distribution function of the consecutive observations is constructed through copula functions. First- and second-order Markov chains are considered with the univariate margins of zero-inflated Poisson (ZIP), zero-inflated negative binomial (ZINB), or zero-inflated Conway–Maxwell–Poisson (ZICMP) distributions. Under the Markov models, bivariate copula functions such as the bivariate Gaussian, Frank, and Gumbel are chosen to construct a bivariate distribution of two consecutive observations. Moreover, the trivariate Gaussian and max-infinitely divisible copula functions are considered to build the joint distribution of three consecutive observations. Likelihood-based inference is performed and asymptotic properties are studied. To evaluate the estimation method and the asymptotic results, simulated examples are studied. The proposed class of models are applied to sandstorm counts example. The results suggest that the proposed models have some advantages over some of the models in the literature for modeling zero-inflated count time series data. [ABSTRACT FROM AUTHOR]

This paper considers the vector autoregressive model of order p, VAR(p), with multivariate t error distributions, the latter being more prevalent in real life than the usual multivariate normal distribution. It is believed that the maximum-likelihood equations for the multivariate t distribution have convergence problem, hence we develop estimation procedures for VAR(p) model using the normal mean–variance mixture representation of multivariate t distribution. The procedure relies on the computational ease available in Expectation Maximization-based algorithms. The estimators obtained are explicit functions of sample observations and therefore are easy to compute. Extensive simulation experiments show that the estimators have negligible bias and are considerably more efficient than an existing method that uses the least-squares error approach. It is shown that the proposed estimators are robust to plausible deviations from an assumed distribution and hence are more advantageous when compared with the other estimator. One real-life example is given for illustration purposes. [ABSTRACT FROM AUTHOR]

Many control charts have been developed for the simultaneous monitoring of the time interval T between successive occurrences of an event E and its magnitude X. All these TBEA (Time Between Events and Amplitude) control charts assume a known distribution for the random variables T and X. But, in practice, as it is rather difficult to know their actual distributions, proposing a distribution free approach could be a way to overcome this 'distribution choice' dilemma. For this reason, we propose in this paper a distribution free upper-sided EWMA (Exponentially Weighted Moving Average) type control chart, for simultaneously monitoring the time interval T and the magnitude X of an event. In order to investigate the performance of this control chart and obtain its run length properties, we also develop a specific method called 'continuousify' which, coupled with a classical Markov chain technique, allows to obtain reliable and replicable results. A numerical comparison shows that our distribution-free EWMA TBEA chart performs as the parametric Shewhart TBEA chart, but without the need to pre-specify any distribution. An illustrative example obtained from a French forest fire database is also provided to show the implementation of the proposed EWMA TBEA control chart. [ABSTRACT FROM AUTHOR]

In this paper, a new two-parameter discrete distribution is introduced. It belongs to the family of the weighted geometric distribution (GD), with the feature of using a particular trigonometric weight. This configuration adds an oscillating property to the former GD which can be helpful in analyzing the data with over-dispersion, as developed in this study. First, we present the basic statistical properties of the new distribution, including the cumulative distribution function, hazard rate function and moment generating function. Estimation of the related model parameters is investigated using the maximum likelihood method. A simulation study is performed to illustrate the convergence of the estimators. Applications to two practical datasets are given to show that the new model performs at least as well as some competitors. [ABSTRACT FROM AUTHOR]

The concept of ranked set sampling (RSS) is applicable whenever ranking on a set of sampling units can be done easily using a judgment method or based on an auxiliary variable. In this paper, we consider a study variable Y correlated with the auxiliary variable X and use it to rank the sampling units. Further (X , Y) is assumed to have Cambanis-type bivariate uniform (CTBU) distribution. We obtain an unbiased estimator of a scale parameter associated with the study variable Y based on different RSS schemes. We perform the efficiency comparison of the proposed estimators numerically. We present the trends in the efficiency performance of estimators under various RSS schemes with respect to parameters through line and surface plots. Further, we develop a Matlab function to simulate data from CTBU distribution and present the performance of proposed estimators through a simulation study. The results developed are implemented to real-life data also. [ABSTRACT FROM AUTHOR]

In this paper, we present a repetitive sampling method to construct control charts using exponentially weighted moving averages (EWMA) and double exponentially weighted moving averages (DEWMA) to monitor shift in the process. For non-normal processes, t-distribution with various degrees of freedom (i.e. df = 4 , 10 , 20 , 40 , 50) is used as symmetric distribution, gamma distribution with unit scale parameter and various shape parameters (i.e. 0.5 , 1 , 2 , 3 , 4) is used as positively skewed distribution and Weibull distribution with unit scale parameter and various shape parameters (i.e. 10 and 20) is used as negatively skewed distribution. We use Monte Carlo simulations to check whether the process is out of control. We use average run length as a tool to find the ability of proposed control charts to identify a shift earlier in a process, as compared to other control charts currently used to monitor the same type of process. The proposed control charts are applied to two real datasets. [ABSTRACT FROM AUTHOR]

The Birnbaum–Saunders distribution is a widely used distribution in reliability applications to model failure times. For several samples from possible different Birnbaum–Saunders distributions, if their means can be considered as the same, it is of importance to make inference for the common mean. This paper presents procedures for interval estimation and hypothesis testing for the common mean of several Birnbaum–Saunders populations. The proposed approaches are hybrids between the generalized inference method and the large sample theory. Some simulation results are conducted to present the performance of the proposed approaches. The simulation results indicate that our proposed approaches perform well. Finally, the proposed approaches are applied to analyze a real example on the fatigue life of 6061-T6 aluminum coupons for illustration. [ABSTRACT FROM AUTHOR]

Spatial sign and rank-based methods have been studied in the recent literature, especially when the dimension is smaller than the sample size. In this paper, a classification method based on the distribution of rank functions for high-dimensional data is considered with extension to functional data. The method is fully nonparametric in nature. The performance of the classification method is illustrated in comparison with some other classifiers using simulated and real data sets. Supporting code in R are provided for computational implementation of the classification method that will be of use to others. [ABSTRACT FROM AUTHOR]

Application of the exact statistical inference frequently leads to non-standard probability distributions of the considered estimators or test statistics. The exact distributions of many estimators and test statistics can be specified by their characteristic functions, as is the case for the null distribution of the Bartlett's test statistic. However, analytical inversion of the characteristic function, if possible, frequently leads to complicated expressions for computing the distribution function and the corresponding quantiles. An efficient alternative is the well-known method based on numerical inversion of the characteristic functions, which is, however, ignored in popular statistical software packages. In this paper, we present the explicit characteristic function of the corrected Bartlett's test statistic together with the computationally fast and efficient implementation of the approach based on numerical inversion of this characteristic function, suggested for evaluating the exact null distribution used for testing homogeneity of variances in several normal populations, with possibly unequal sample sizes. [ABSTRACT FROM AUTHOR]

This paper analyses the statistical properties of the threshold model of collective actions. We establish the conditions of the threshold model under which the sequence of decisions is exchangeable. Then, if the probability of participation is a continuous and strictly increasing function of the willingness to participate, then the unique urn representation of the threshold model is the Pólya urn. In such a case, the sequence of participation rates converges almost surely to a random variable Z that has a beta distribution. Also, we show the critical mass and self-sustaining participation rates are fixed points of the distribution function. Further, we discuss the non-exchangeable model of collective actions, and we define the feedback leadership condition. We apply these results to characterise the process of the adoption of hybrid corn in Iowa, and the diffusion of the Internet for different levels of income. [ABSTRACT FROM AUTHOR]

This paper introduces a new two-parameter exponentiated discrete Lindley distribution. A wide range of its structural properties are investigated. This includes the shape of the probability mass function, hazard rate function, moments, skewness, kurtosis, stress–strength reliability, mean residual lifetime, mean past lifetime, order statistics and L-moment statistics. The hazard rate function can be increasing, decreasing, decreasing–increasing–decreasing, increasing–decreasing–increasing, unimodal, bathtub, and J-shaped depending on its parameters values. Two methods are used herein to estimate the model parameters, namely, the maximum likelihood, and the proportion. A detailed simulation study is carried out to examine the bias and mean square error of maximum likelihood and proportion estimators. The flexibility of the proposed model is explained by using four distinctive data sets. It can serve as an alternative model to other lifetime distributions in the existing statistical literature for modeling positive real data in many areas. [ABSTRACT FROM AUTHOR]

In this paper, we study empirical distribution function (EDF) estimators based on ranked set sampling (RSS). Ranked set samples are constructed by using three different sampling designs which are level-0, level-1 and level-2. These designs are different in terms of their replacement policies. The sample is constructed by sampling with replacement in level-0 sampling design. In level-1 sampling design, the sample is constructed by sampling with replacement of the unmeasured units in each set. Also, we construct the sample without replacement on the entire set in level-2 sampling design. In simulation study, we compare EDF estimators when using sampling designs with EDF estimator when using simple random sampling (SRS) in finite populations having different distribution functions. The effects of both perfect and imperfect ranking on the estimators of the sampling designs are investigated. Relative efficiencies of the EDF estimators are obtained by using their integrated mean squared errors (IMSEs), numerically. For illustrative purposes, these EDF estimators based on sampling designs are applied to a real data set. By using both the simulation data and real data, we show that these EDF estimators have higher efficiencies against EDF based on SRS. [ABSTRACT FROM AUTHOR]

This paper focuses on estimating the number of species and the number of abundant species in a specific geographic region and, consequently, draw inferences on the number of rare species. The word 'species' is generic referring to any objects in a population that can be categorized. In the areas of biology, ecology, literature, etc, the species frequency distributions are usually severely skewed, in which case the population contains a few very abundant species and many rare ones. To model a such situation, we develop an asymmetric multinomial-Dirichlet probability model using species frequency data. Posterior distributions on the number of species and the number of abundant species are obtained and posterior inferences are induced using MCMC simulations. Simulations are used to demonstrate and evaluate the developed methodology. We apply the method to a DNA segment data set and a butterfly data set. Comparisons among different approaches to inferring the number of species are also discussed in this paper. [ABSTRACT FROM AUTHOR]

In this paper, we are employing the generalized linear model (GLM) in the form lij=χλ to decompose the symmetry model into the class of models discussed in Tomizawa (1992). In this formulation, the random component would be the observed counts fij with an underlying Poisson distribution. This approach utilizes the non-standard log-linear model and our focus in this paper therefore relates to models that are decompositions of the complete symmetry model. That is, models that are implied by the symmetry models. We develop factor and regression variables required for the implementation of these models in SAS PROC GENMOD and SPSS PROC GENLOG. We apply this methodology to analyse the three 4 × 4 contingency table, one of which is the Japanese Unaided distance vision data. Results obtained in this study are consistent with those from the numerous literature on the subject. We further extend our applications to the 6 × 6 Brazilian social mobility data. We found that both the quasi linear diagonal-parameters symmetry (QLDPS) and the quasi 2-ratios parameter symmetry (Q2RPS) models fit the Brazilian data very well. Parsimonious models being the QLDPS and the quasi-conditional symmetry (QCS) models. The SAS and SPSS programs for implementing the models discussed in this paper are presented in Appendices A, B and C. [ABSTRACT FROM AUTHOR]