Small sample size comparisons of tests for homogeneity of variances by Monte-Carlo

A number of tests are available for testing the equality of several population variances. Some are claimed to be robust. We compared six of those claimed robust procedures by Monte Carlo simulated experiments, particularly for cases of small and unequal sample sizes. Our results show that the jack-knife test compares favorably with the other tests.     

Conditional specifications of multivariate pareto and student distributions

A random vector has a multivariate Pareto distribution if one of its univariate conditional distribution is Pareto and some of its marginals are identically distributed.A general method developed in the course of the proof of this result is applied also to characterize the multivariate Student (Cauchy) measure by one univariate Student conditional distribution.     

An Improvement of the Nonparametric Bootstrap Test for the Comparison of the Coefficient of Variations

In this article, we propose a new test for examining the equality of the coefficient of variation between two different populations. The proposed test is based on the nonparametric bootstrap method. It appears to yield several appreciable advantages over the current tests. The quick and easy implementation of the test can be considered as advantages of the proposed test. The test is examined by the Monte Carlo simulations, and also evaluated using various numerical studies.     

Confidence intervals for product of powers of the generalized variances of k multivariate normal populations

In this paper, some confidence intervals (CIs) for the product of powers of the generalized variances of k multivariate normal populations with possibly different dimensions are proposed. The performance of these CIs in terms of the coverage probabilities and average lengths were evaluated via a Monte Carlo simulation study. The results were found to be satisfactory. To demonstrate utility of the proposed CIs, applications on three real data sets were provided.     

A MATLAB package for multivariate normality test

A MATLAB package testing for multivariate normality (TMVN) is implemented as an interactive and graphical tool to examine multivariate normality (MVN). Monte Carlo simulation studies have failed to find a uniformly most powerful MVN test, which requires a rather extensive statistical inference procedure. TMVN contains several competitive MVN tests and provides a flexible and extensive testing environment for univariate or multivariate data analyses. Simulated results provide information of which test may possess more power for the selected non-MVN alternatives. Fisher's Iris data are used to show how TMVN can be used in practice.     

Nonparametric estimation of the bivariate survival function for one modified form of current-status data

In this article, the estimation of the bivariate survival function for one modified form of current-status data is considered. Two types of estimators, which are generalizations of the estimators by Campbell and Földes [G. Campbell and A. Földes, Large sample properties of nonparametric statistical inference, in Nonparametric Statistical Inference, B.V. Gnredenko, M.L. Puri, and I. Vineze, eds., North-Holland, Amsterdam, 1982, pp. 103–122] and Dabrowska [D.M. Dabrowska, Kaplan-Meier estimate on the plane, Ann. Stat. 18 (1988), pp. 1475–1489; D.M. Dabrowska, Kaplan-Meier estimate on the plane: weak convergence, LIL, and the bootstrap, J. Multivariate Anal. 29 (1989), pp. 308–325], are proposed. The consistency of the proposed estimators is established. A simulation study is conducted to investigate the performance of the proposed estimators.     

A Comparison of Frailty and Other Models for Bivariate Survival Data

Multivariate survival data arise when eachstudy subject may experience multiple events or when study subjectsare clustered into groups. Statistical analyses of such dataneed to account for the intra-cluster dependence through appropriatemodeling. Frailty models are the most popular for such failuretime data. However, there are other approaches which model thedependence structure directly. In this article, we compare thefrailty models for bivariate data with the models based on bivariateexponential and Weibull distributions. Bayesian methods providea convenient paradigm for comparing the two sets of models weconsider. Our techniques are illustrated using two examples.One simulated example demonstrates model choice methods developedin this paper and the other example, based on a practical dataset of onset of blindness among patients with diabetic Retinopathy,considers Bayesian inference using different models.     

Simultaneous Confidence Intervals for the Parameters of Pareto Distribution under Progressive Censoring

In this article, we consider the progressive Type II right censored sample from Pareto distribution. We introduce a new approach for constructing the simultaneous confidence interval of the unknown parameters of this distribution under progressive censoring. A Monte Carlo study is also presented for illustration. It is shown that this confidence region has a smaller area than that introduced by Ku? and Kaya (2007 Ku? , C. , Kaya , M. F. ( 2007 ). Estimation for the parameters of the Pareto distribution under progressive censoring . Commun. Statist. Theor. Meth. 36 : 1359 – 1365 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).     

Sharp bounds for the mean of the total time on test for distributions with increasing generalized failure rate

The results of this paper are the continuation of the research presented by Bieniek [Optimal bounds for the mean of the total time on test for distributions with decreasing generalized failure rate. Statistics. 2016;50:1206–1220]. We consider the remaining total time on test after a given failure in a life test experiment. We derive sharp upper bounds on the mean of the total time on test optimal in the class of distributions with increasing generalized failure rate with respect to generalized Pareto distributions. Specific results are obtained for distributions with increasing density and increasing failure rate. We also provide exemplary numerical values of the obtained bounds and we compare them with the corresponding bounds for distributions with decreasing generalized failure rate.     

Critical values improvement for the standard normal homogeneity test by combining Monte Carlo and regression approaches

The distribution of the test statistics of homogeneity tests is often unknown, requiring the estimation of the critical values through Monte Carlo (MC) simulations. The computation of the critical values at low α, especially when the distribution of the statistics changes with the series length (sample cardinality), requires a considerable number of simulations to achieve a reasonable precision of the estimates (i.e. 10⁶ simulations or more for each series length). If, in addition, the test requires a noteworthy computational effort, the estimation of the critical values may need unacceptably long runtimes.

To overcome the problem, the paper proposes a regression-based refinement of an initial MC estimate of the critical values, also allowing an approximation of the achieved improvement. Moreover, the paper presents an application of the method to two tests: SNHT (standard normal homogeneity test, widely used in climatology), and SNH2T (a version of SNHT showing a squared numerical complexity). For both, the paper reports the critical values for α ranging between 0.1 and 0.0001 (useful for the p-value estimation), and the series length ranging from 10 (widely adopted size in climatological change-point detection literature) to 70,000 elements (nearly the length of a daily data time series 200 years long), estimated with coefficients of variation within 0.22%. For SNHT, a comparison of our results with approximated, theoretically derived, critical values is also performed; we suggest adopting those values for the series exceeding 70,000 elements.     

A test for homogeneity when the sample is a positive lagrangian poisson distribution

In this note, we obtain, based on the sample sum, a statistic to test the homogeneity of a random sample from a positive (zero truncated) Lagrangian Poisson distribution given in Consul and Jain (1973). This test statistic conforms, in a special case, to Singh (1978). A goodness-of-fit test statistic for the Borel-Tanner distribution is obtained as a particular case cf our results.     

Construction of Simultaneous Confidence Intervals for Ratios of Means of Lognormal Distributions

For constructing simultaneous confidence intervals for ratios of means for lognormal distributions, two approaches using a two-step method of variance estimates recovery are proposed. The first approach proposes fiducial generalized confidence intervals (FGCIs) in the first step followed by the method of variance estimates recovery (MOVER) in the second step (FGCIs–MOVER). The second approach uses MOVER in the first and second steps (MOVER–MOVER). Performance of proposed approaches is compared with simultaneous fiducial generalized confidence intervals (SFGCIs). Monte Carlo simulation is used to evaluate the performance of these approaches in terms of coverage probability, average interval width, and time consumption.     

Bootstrap approach to test the homogeneity of order restricted mean vectors when the covariance matrices are unknown

Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are unknown, arbitrary positive definite and unequal are considered. This problem of testing has been studied to some extent, for example, by Kulatunga and Sasabuchi (1984 Kulatunga, D. D. S., Sasabuchi, S. (1984). A test of homogeneity of mean vectors against multivariate isotonic alternatives. Mem Fac Sci, Kyushu Univ Ser A Mathemat 38:151–161. [Google Scholar]) when the covariance matrices are known and also Sasabuchi et al. (2003 Sasabuchi, S., Tanaka, K., Tsukamodo, T. (2003). Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are common but unknown. Annals of Statistics. 31(5):1517–1536.[Web of Science ®] , [Google Scholar]) and Sasabuchi (2007 Sasabuchi, S. (2007). More powerful tests for homogeneity of multivariate normal mean vectors under an order restriction. Sankhya 69(4):700–716. [Google Scholar]) when the covariance matrices are unknown but common. In this paper, a test statistic is proposed and because of the main advantage of the bootstrap test is that it avoids the derivation of the complex null distribution analytically, a bootstrap test statistic is derived and since the proposed test statistic is location invariance the bootstrap p-value defined logical and some steps are presented to estimate it. Our numerical studies via Monte Carlo simulation show that the proposed bootstrap test can correctly control the type I error rates. The power of the test for some of the p-dimensional normal distributions is computed by Monte Carlo simulation. Also, the null distribution of test statistic is estimated using kernel density. Finally, the bootstrap test is illustrated using a real data.     

A nonparametric test for Type I censored data: r samples

This paper is concerned with testing that r random samples are all from the same population when the data are left (or right) censored. A statistic is developed which has elements of a goodness-of-fit statistic and of a modified Kruskal-Wallis statistic. The efficiency of this statistic relative to some other commonly used statistics is calculated. Some Monte Carlo comparisons are given.     

Estimation for the generalized half-normal distribution based on progressive type-II censoring

Based on progressively type-II censored data, the problem of estimating unknown parameters and reliability function of a two-parameter generalized half-normal distribution is considered. Maximum likelihood estimates are obtained by applying expectation-maximization algorithm. Since they do not have closed forms, approximate maximum likelihood estimators are proposed. Several Bayesian estimates with respect to different symmetric and asymmetric loss functions such as squared error, LINEX and general entropy are calculated. The Lindley approximation method is applied to determine Bayesian estimates. Monte Carlo simulations are performed to compare the performances of the different methods. Finally, one real data set is analysed.     

A fully nonparametric diagnostic test for homogeneity of variances

The authors propose a new nonparametric diagnostic test for checking the constancy of the conditional variance function σ²(x) in the regression model Y_i = m(x_i) + σ(x_i)?_i, i = 1,…, m. Their test, which does not assume a known parametric form for the conditional mean function m(x), is inspired by a recent asymptotic theory in the analysis of variance when the number of factor levels is large. The authors demonstrate through simulations the good finite‐sample properties of the test and illustrate its use in a study on the effect of drug utilization on health care costs.     

Nonparametric analysis of aggregate loss models

This paper describes a nonparametric approach to make inferences for aggregate loss models in the insurance framework. We assume that an insurance company provides a historical sample of claims given by claim occurrence times and claim sizes. Furthermore, information may be incomplete as claims may be censored and/or truncated. In this context, the main goal of this work consists of fitting a probability model for the total amount that will be paid on all claims during a fixed future time period. In order to solve this prediction problem, we propose a new methodology based on nonparametric estimators for the density functions with censored and truncated data, the use of Monte Carlo simulation methods and bootstrap resampling. The developed methodology is useful to compare alternative pricing strategies in different insurance decision problems. The proposed procedure is illustrated with a real dataset provided by the insurance department of an international commercial company.     

Optimal bounds for the mean of the total time on test for distributions with decreasing generalized failure rate

During any life test experiment it is of interest to study potential costs (or profits) of performing the test to the very end. Assuming that these costs are proportional to lifetimes of analysed items the experimenter needs to know the remaining total time on test, having just observed successive failure in the test. We derive sharp upper bounds on the expectation of the remaining total time on test statistic when the underlying distributions have decreasing generalized failure rate with respect to generalized Pareto distributions. In particular we obtain the bounds valid for distributions with decreasing density or failure rate. The results are illustrated with numerical examples.     

Goodness of fit test for the generalized Rayleigh distribution with unknown parameters

The use of goodness-of-fit test based on Anderson–Darling (AD) statistic is discussed, with reference to the composite hypothesis that a sample of observations comes from a generalized Rayleigh distribution whose parameters are unspecified. Monte Carlo simulation studies were performed to calculate the critical values for AD test. These critical values are then used for testing whether a set of observations follows a generalized Rayleigh distribution when the scale and shape parameters are unspecified and are estimated from the sample. Functional relationship between the critical values of AD is also examined for each shape parameter (α), sample size (n) and significance level (γ). The power study is performed with the hypothesized generalized Rayleigh against alternate distributions.