On Blest’s Measure of Kurtosis Adjusted for Skewness

We reconsider the derivation of Blest’s (2003) skewness adjusted version of the classical moment-based coefficient of kurtosis and propose an adaptation of it which generally eliminates the effects of asymmetry a little more successfully. Lower bounds are provided for the two skewness adjusted kurtosis moment measures as functions of the classical coefficient of skewness. The results from a Monte Carlo experiment designed to investigate the sampling properties of numerous moment-based estimators of the two skewness adjusted kurtosis measures are used to identify those estimators with lowest mean squared error for small to medium sized samples drawn from distributions with varying levels of asymmetry and tailweight.     

The alpha–beta skew normal distribution: properties and applications

This paper introduces a new class of skew distributions by extending the alpha skew normal distribution proposed by Elal-Olivero [Elal-Olivero, D. Alpha-skew-normal distribution. Proyecciones. 2010;29:224–240]. Statistical properties of the new family are studied in details. In particular, explicit expressions for the moments and the shape parameters including the skewness and the kurtosis coefficients and the moment generating function are derived. The problem of estimating parameters on the basis of a random sample coming from the new class of distribution is considered. To examine the performance of the obtained estimators, a Monte Carlo simulation study is conducted. Flexibility and usefulness of the proposed family of distributions are illustrated by analysing three real data sets.     

Skewness and mixtures of normal distributions

The effect of skewness on hypothesis tests for the existence of a mixture of univariate and bivariate normal distributions is examined through a Monte Carlo study. A likelihood ratio test based on results of the simultaneous estimation of skewness parameters, derived from power transformations, with mixture parameters is proposed. This procedure detects the difference between inherent distributional skewness and the apparent skewness which is a manifestation of the mixture of several distributions. The properties of this test are explored through a simulation study.     

The skew generalized t distribution as the scale mixture of a skew exponential power distribution and its applications in robust estimation

In this paper, we consider the family of skew generalized t (SGT) distributions originally introduced by Theodossiou [P. Theodossiou, Financial data and the skewed generalized t distribution, Manage. Sci. Part 1 44 (12) ( 1998), pp. 1650–1661] as a skew extension of the generalized t (GT) distribution. The SGT distribution family warrants special attention, because it encompasses distributions having both heavy tails and skewness, and many of the widely used distributions such as Student's t, normal, Hansen's skew t, exponential power, and skew exponential power (SEP) distributions are included as limiting or special cases in the SGT family. We show that the SGT distribution can be obtained as the scale mixture of the SEP and generalized gamma distributions. We investigate several properties of the SGT distribution and consider the maximum likelihood estimation of the location, scale, and skewness parameters under the assumption that the shape parameters are known. We show that if the shape parameters are estimated along with the location, scale, and skewness parameters, the influence function for the maximum likelihood estimators becomes unbounded. We obtain the necessary conditions to ensure the uniqueness of the maximum likelihood estimators for the location, scale, and skewness parameters, with known shape parameters. We provide a simple iterative re-weighting algorithm to compute the maximum likelihood estimates for the location, scale, and skewness parameters and show that this simple algorithm can be identified as an EM-type algorithm. We finally present two applications of the SGT distributions in robust estimation.     

The odd power cauchy family of distributions: properties,regression models and applications

We study some mathematical properties of a new generator of continuous distributions with one extra parameter called the odd power Cauchy family including asymptotics, linear representation, moments, quantile and generating functions, entropies, order statistics and extreme values. We introduce two bivariate extensions of the new family. The maximum likelihood method is discussed to estimate the model parameters by means of a Monte Carlo simulation study. We define a new log-odd power Cauchy–Weibull regression model. The usefulness of the proposed models is proved empirically by means of three real data sets.     

Nonparametric test for the homogeneity of the overall variability

In this paper, we propose a nonparametric test for homogeneity of overall variabilities for two multi-dimensional populations. Comparisons between the proposed nonparametric procedure and the asymptotic parametric procedure and a permutation test based on standardized generalized variances are made when the underlying populations are multivariate normal. We also study the performance of these test procedures when the underlying populations are non-normal. We observe that the nonparametric procedure and the permutation test based on standardized generalized variances are not as powerful as the asymptotic parametric test under normality. However, they are reliable and powerful tests for comparing overall variability under other multivariate distributions such as the multivariate Cauchy, the multivariate Pareto and the multivariate exponential distributions, even with small sample sizes. A Monte Carlo simulation study is used to evaluate the performance of the proposed procedures. An example from an educational study is used to illustrate the proposed nonparametric test.     

Testing the equality means of several log-normal distributions

Researches propose various methods for comparing the means of two log-normal distributions. Some of these methods have been recently extended to test the equality means of several log-normal populations. Investigations show that none of the established methods is satisfactory. In this article, we provide three methods based on the computational approach test, which is a parametric bootstrap approach, for testing the means of several log-normal distributions. Further, we compare our methods with the existing methods through Monte Carlo simulation. The numerical results show that the Type I errors of these procedures are satisfactory regardless of the sample size, number of populations, and the true parameters. Finally, we explain the considered methods by real examples.     

Modified Tukey's Control Chart

Phase I of control analysis requires large amount of data to fit a distribution and estimate the corresponding parameters of the process under study. However, when only individual observations are available, and no a priori knowledge exists, the presence of outliers can bias the analysis. A relatively recent and successful approach to address this situation is Tukey's Control Chart (TCC), a charting method that applies the Box Plot technique to estimate the control limits. This procedure has proven to be effective for symmetric distributions. However, when skewness is present the average run length performance diminishes significantly. This article proposes a modified version of TCC to consider skewness with minimum assumptions on the underlying distribution of observations. Using theoretical results and Monte Carlo simulation, the modified TCC is tested over several distributions proving a better representation of skewed populations, even in cases when only a limited number of observations are available.     

A Monte Carlo comparison of Jarque–Bera type tests and Henze–Zirkler test of multivariate normality

In the paper, tests for multivariate normality (MVN) of Jarque-Bera type, based on skewness and kurtosis, have been considered. Tests proposed by Mardia and Srivastava, and the combined tests based on skewness and kurtosis defined by Jarque and Bera have been taken into account. In the Monte Carlo simulations, for each combination of p = 2, 3, 4, 5 number of traits and n = 10(5)50(10)100 sample sizes 10,000 runs have been done to calculate empirical Type I errors of tests under consideration, and empirical power against different alternative distributions. Simulation results have been compared to the Henze–Zirkler’s test. It should be stressed that no test yet proposed is uniformly better than all the others in every combination of conditions examined.     

A comparative study of tests for paired lifetime data

In this paper, we investigate different procedures for testing the equality of two mean survival times in paired lifetime studies. We consider Owen’s M-test and Q-test, a likelihood ratio test, the paired t-test, the Wilcoxon signed rank test and a permutation test based on log-transformed survival times in the comparative study. We also consider the paired t-test, the Wilcoxon signed rank test and a permutation test based on original survival times for the sake of comparison. The size and power characteristics of these tests are studied by means of Monte Carlo simulations under a frailty Weibull model. For less skewed marginal distributions, the Wilcoxon signed rank test based on original survival times is found to be desirable. Otherwise, the M-test and the likelihood ratio test are the best choices in terms of power. In general, one can choose a test procedure based on information about the correlation between the two survival times and the skewness of the marginal survival distributions.     

Robustness of subset selection procedures

The robustness (and the number of non-best populations selected) of 11 subset selection procedures is investigated by means of simulation experiments. If the underlying distributions differ only in their location parameter, the subset selection procedures are robust for symmetric distributions or distributions with negative skewness. With increasing positive skewness and increasing number of populations the considered parametric procedures fail in robustness slightly. This non-robustness is more serious in the case of unequal variances. Non-parametric subset selection rules show then an increasing non-robustness with increasing sample size.     

Bayesian multivariate GARCH models with dynamic correlations and asymmetric error distributions

The main goal in this paper is to develop and apply stochastic simulation techniques for GARCH models with multivariate skewed distributions using the Bayesian approach. Both parameter estimation and model comparison are not trivial tasks and several approximate and computationally intensive methods (Markov chain Monte Carlo) will be used to this end. We consider a flexible class of multivariate distributions which can model both skewness and heavy tails. Also, we do not fix tail behaviour when dealing with fat tail distributions but leave it subject to inference.     

Interval Estimations for the Mean of the Generalized Poisson and Generalized Negative Binomial Distributions

In this article, we extend the Wald, score, skewness-corrected score, likelihood ratio, and mid-P intervals for the means of the generalized Poisson and generalized negative binomial distributions. These distributions are the members of the discrete version of the natural exponential family (NEF) with cubic variance function (CVF). Also, the coverage probabilities, the distal and mesial noncoverage probabilities, and the lengths of the proposed confidence intervals are estimated by means of a Monte Carlo simulation study. Finally, some practical examples are provided to show the applicability of the proposed intervals in applied studies.     

Inference on proportional hazard rate model parameter under Type-I progressively hybrid censoring scheme

ABSTRACT

In this paper, under Type-I progressive hybrid censoring sample, we obtain maximum likelihood estimator of unknown parameter when the parent distribution belongs to proportional hazard rate family. We derive the conditional probability density function of the maximum likelihood estimator using moment-generating function technique. The exact confidence interval is obtained and compared by conducting a Monte Carlo simulation study for burr Type XII distribution. Finally, we obtain the Bayes and posterior regret gamma minimax estimates of the parameter under a precautionary loss function with precautionary index k = 2 and compare their behavior via a Monte Carlo simulation study.     

Bayesian analysis of circular distributions based on non-negative trigonometric sums

ABSTRACT

Fernández-Durán [Circular distributions based on nonnegative trigonometric sums. Biometrics. 2004;60:499–503] developed a new family of circular distributions based on non-negative trigonometric sums that is suitable for modelling data sets that present skewness and/or multimodality. In this paper, a Bayesian approach to deriving estimates of the unknown parameters of this family of distributions is presented. Because the parameter space is the surface of a hypersphere and the dimension of the hypersphere is an unknown parameter of the distribution, the Bayesian inference must be based on transdimensional Markov Chain Monte Carlo (MCMC) algorithms to obtain samples from the high-dimensional posterior distribution. The MCMC algorithm explores the parameter space by moving along great circles on the surface of the hypersphere. The methodology is illustrated with real and simulated data sets.     

Omnibus tests of multinormality based on skewness and kurtosis

Measures of univariate skewness and kurtosis have long been used as a test of univariate normality, several omnibus test procedures based on a combination of the measures having been proposed, see Pearson, D’Agestino and Bowman (1977) and Mardia (1979). Mardia (1970) proposed measures of multivariate skewness and kurtosis, and constructed a test of multinormality based on these measures. we obtain the correlation between these measures and propose several omnibus tests using the two measures. The performances of these tests are compared by means of a Monte Carlo study.     

Inference for Weibull distribution under adaptive Type-I progressive hybrid censored competing risks data

In this paper, a competing risks model is considered under adaptive type-I progressive hybrid censoring scheme (AT-I PHCS). The lifetimes of the latent failure times have Weibull distributions with the same shape parameter. We investigate the maximum likelihood estimation of the parameters. Bayes estimates of the parameters are obtained based on squared error and LINEX loss functions under the assumption of independent gamma priors. We propose to apply Markov Chain Monte Carlo (MCMC) techniques to carry out a Bayesian estimation procedure and in turn calculate the credible intervals. To evaluate the performance of the estimators, a simulation study is carried out.     

An Extension of the Normal/Independent Family for Studies of Robust Estimators

The family of normal scale mixture distributions, also called the Normal/Independent family, has been used for efficient Monte Carlo studies of robust estimators. The distributions in this family are unimodal. The Normal/Independent family is extended by introducing a location mixing in addition to the scale mixing. Distributions in this extension may be nonunimodal. The asymptotic variances of robust estimators of location are compared using the distributions from the extension. A Monte Carlo swindle similar to the one used in the Princeton study is given for the extended family. A small simulation study demonstrates the efficiency of the swindle. The swindle is compared with other swindle methods based on Fisher's score function and regression.     

The generalized odd log-logistic family of distributions: properties,regression models and applications

We propose a new class of continuous distributions with two extra shape parameters named the generalized odd log-logistic family of distributions. The proposed family contains as special cases the proportional reversed hazard rate and odd log-logistic classes. Its density function can be expressed as a linear combination of exponentiated densities based on the same baseline distribution. Some of its mathematical properties including ordinary moments, quantile and generating functions, two entropy measures and order statistics are obtained. We derive a power series for the quantile function. We discuss the method of maximum likelihood to estimate the model parameters. We study the behaviour of the estimators by means of Monte Carlo simulations. We introduce the log-odd log-logistic Weibull regression model with censored data based on the odd log-logistic-Weibull distribution. The importance of the new family is illustrated using three real data sets. These applications indicate that this family can provide better fits than other well-known classes of distributions. The beauty and importance of the proposed family lies in its ability to model different types of real data.     

Efficient posterior integration in stable paretian models

The paper proposes a Markov Chain Monte Carlo method for Bayesian analysis of general regression models with disturbances from the family of stable distributions with arbitrary characteristic exponent and skewness parameter. The method does not require data augmentation and is based on combining fast Fourier transforms of the characteristic function to get the likelihood function and a Metropolis random walk chain to perform posterior analysis. Both a validation nonlinear regression and a nonlinear model for the Standard and Poor’s composite price index illustrate the method.