Maximum likelihood estimation in Dagum distribution with censored samples

In this work, we show that the Dagum distribution [3 Dagum, C. 1977. A new model of personal income distribution: Specification and estimation. Econ. Appl., XXX: 413–436.  [Google Scholar]] may be a competitive model for describing data which include censored observations in lifetime and reliability problems. Maximum likelihood estimates of the three parameters of the Dagum distribution are determined from samples with type I right and type II doubly censored data. We perform an empirical analysis using published censored data sets: in certain cases, the Dagum distribution fits the data better than other parametric distributions that are more commonly used in survival and reliability analysis. Graphical comparisons confirm that the Dagum model behaves better than a number of competitive distributions in describing the empirical hazard rate of the analyzed data. A probability plot to provide graphical check of the appropriateness of the Dagum model for right censored data is constructed, and the details are given in the appendix. Finally, a simulation study that shows the good performance of the maximum likelihood estimators of the Dagum shape parameters for finite type II doubly censored samples is carried out.     

A test for equality of variances with censored samples

A variance homogeneity test for type II right-censored samples is proposed. The test is based on Bartlett's statistic. The asymptotic distribution of the statistic is investigated. The limiting distribution is that of a linear combination of i.i.d. chi-square variables with 1 degree of freedom. By using simulation, the critical values of the null distribution of the modified Bartlett's statistic for testing the homogeneity of variances of two normal populations are obtained when the sample sizes and censoring levels are not equal. Also, we investigate the properties of the proposed test (size, power and robustness). Results show that the distribution of the test statistic depends on the censoring level. An example of the use of the new methodology in animal science involving reproduction in ewes is provided.     

Gap-ratio goodness of fit tests for weibull or extreme value distribution assumptions with left or right censored data

With data collection in environmental science and bioassay, left censoring because of nondetects is a problem. Similarly in reliability and life data analysis right censoring frequently occurs. There is a need for goodness of fit tests that can adapt to left or right censored data and be used to check important distributional assumptions without becoming too difficult to regularly implement in practice. A new test statistic is derived from a plot of the standardized spacings between the order statistics versus their ranks. Any linear or curvilinear pattern is evidence against the null distribution. When testing the Weibull or extreme value null hypothesis this statistic has a null distribution that is approximately F for most combinations of sample size and censoring of practical interest. Our statistic is compared to the Mann-Scheuer-Fertig statistic which also uses the standardized spacings between the order statistics. The results of a simulation study show the two tests are competitive in terms of power. Although the Mann-Scheuer-Fertig statistic is somewhat easier to compute, our test enjoys advantages in the accuracy of the F approximation and the availability of a graphical diagnostic.     

Asymptotic distributions for quadratic forms with applications to censored data tests of fit

We obtain the asymptotic distributions of certain forms in observations that are possible Type I and Type II censored. This result is directly applicable to the study of asympototic distributions for censored data versions of the Shapiro- wilk test for normality. Moreove, it applies more generally than just to the null hypothesis of normality.     

Generalized two sample kolmogorov-smirnov test involving a possibly censored sample

Let x be a sample from an unknown population X, and x' be another sample from which some of the smallest and/or the largest observations may have been removed . A nonparametric procedure, to test the null hypothesis stating that x' also comes from X, against the alternative hypothesis stating that this is false, is provided. Simulation results are included     

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.     

Asymptotocally nonparametric tests with censored paired data

A class of asymptotically nonparametric test with contains a test proposed by Wei(1980), is considered for testing the equality of two continuous distribution funcitons when paired observations are subject to arbitrary right censorship. It is shown that under the null hypothesis each test statistic converges in distribution to the standard normal random variable. Furthermore. the Monte Carlo simulation results indicate that some tests in this class are more powerful than Wei's test. A generalization to incomplete censored paired data is also included.     

Two-sample spatial rank test using projection

This paper explores in high-dimensional settings how to test the equality of two location vectors. We introduce a rank-based projection test under elliptical symmetry. Optimal projection direction is derived according to asymptotically and locally best power criteria. Data-splitting strategy is used to estimate optimal projection and construct test statistics. The limiting null distribution and power function of the proposed statistics are thoroughly investigated under some mild assumptions. The test is shown to keep type I error rates pretty well and outperforms several existing methods in a broad range of settings, especially in the presence of large correlation structures. Simulation studies are conducted to confirm the asymptotic results and a real data example is applied to demonstrate the advantage of the proposed procedure.     

Adjusted Wilcoxon signed rank test tables for ratio of percentiles

The ordinary Wilcoxon signed rank test table provides confidence intervals for the median of one population. Adjusted Wilcoxon signed rank test tables which can provide confidence intervals for the median and the 10th percentile of one population are created in this paper. Base-(n + 1) number system and theorems about property of symmetry of the adjusted Wilcoxon signed rank test statistic are derived for programming. Theorem 1 states that the adjusted Wilcoxon signed rank test statistic are symmetric around n(n + 1)/4. Theorem 2 states that the adjusted Wilcoxon signed rank test statistic with the same number of negative ranks m are symmetric around m(n+1)/2. 87.5% and 85% confidence intervals of the median are given in the table for n = 12, 13,…, 29 to create approximated 95% confidence intervals of the ratio of medians for two independent populations. 95% and 92.5% confidence intervals of the 10th percentile are given in the table for n = 26, 27, 28, 29 to create approximated 95% confidence regions of the ratio of the 10th percentiles for two independent populations. Finally two large datasets from wood industry will be partitioned to verify the correctness of adjusted Wilcoxon signed rank test tables for small samples.     

Linear signed rank tests for hultivariate location

Previously proposed linear signed rank tests for multivariate location are not invariant under linear transformations of the observations, The asymptotic relative efficiencies of the tests 2 with respect to Hotelling's T²test depend on the direction of shift and the covariance matrix of the alternative distributions. For distributions with highly correlated components, the efficiencies of some of these tests can be arbitrarily low; they approach zero for certain multivariate normal alternatives, This article proposes a transformation of the data to be performed prior to standard linear signed rank tests, The resulting procedures have attractive power and efficiency properties compared to the original tests, In particular, for elliptically symmetric contiguous alternafives, the efficiencies of the new tests equal those of corresponding univariate linear signed rank tests with respect to the t test.     

Two-sample prediction for progressively Type-II censored Weibull lifetimes

Prediction on the basis of censored data has an important role in many fields. This article develops a non-Bayesian two-sample prediction based on a progressive Type-II right censoring scheme. We obtain the maximum likelihood (ML) prediction in a general form for lifetime models including the Weibull distribution. The Weibull distribution is considered to obtain the ML predictor (MLP), the ML prediction estimate (MLPE), the asymptotic ML prediction interval (AMLPI), and the asymptotic predictive ML intervals of the sth-order statistic in a future random sample (Y_s) drawn independently from the parent population, for an arbitrary progressive censoring scheme. To reach this aim, we present three ML prediction methods namely the numerical solution, the EM algorithm, and the approximate ML prediction. We compare the performances of the different methods of ML prediction under asymptotic normality and bootstrap methods by Monte Carlo simulation with respect to biases and mean square prediction errors (MSPEs) of the MLPs of Y_s as well as coverage probabilities (CP) and average lengths (AL) of the AMLPIs. Finally, we give a numerical example and a real data sample to assess the computational comparison of these methods of the ML prediction.     

A bivariate signed rank test for two sample location problem

An affine-invariant signed rank test for the difference in location between two symmetric populations is proposed. The proposed test statistic is compared with Hotelling's T2 test statistic, Mardia's(1967)test statistic, Peters-Randles(1991) test statistic and Wilcoxon's rank sum test statistic using a Monte Carlo Study. It performs better than Mardia's test statistic under almost all populations considered. Under the bivariate normal distribution, it performs better than other test statistics compared for small differences in location between two populations except Hotelling's T2. It performs better than all statistics, including Hotelling's T , for sample size 15 when samples are drawn from Pearson type.     

Lifetime estimation from doubly censored data and dirichlet process prior knowledge on the observable random vector

The purpose of this paper is to present a nonparametric Bayesian procedure for estimating a survival curve in a double censoring situation. Assuming a proportional hazard rates model, we propose a consistent estimation of lifetime, based on a Dirichlet process prior knowledge on the observable random vector. Some large sample properties of this estimator are also derived, We prove strong consistency and asymptotic weak convergence to a Gaussian pro cess. Finally, a simulation study is presented in order to analyze the behavior of the proposed estimator, and establish some comparisons to other estimators.     

Linear rank tests for censored survival data

We describe a class of rank test procedures for the two-sample problem with right censored survival data. The class of tests is directly generalized from the linear rank tests by assigning each observation a rank according to its corresponding Wilcoxon scores. It allows a flexible choice of score functions, in particular, those powerful against scale differences between the two survival distributions. Monte Carlo simulations have shown that some members of this class have great power in detecting crossing-curve alternatives (alternatives where underlying survival curves cross over). The class also contains tests essentially equivalent to the Gehan-Wilcoxon and the logrank tests.     

Nonparametric analysis of interval censored and doubly truncated data

Doubly truncated data appear in a number of applications, including astronomy and survival analysis. For double-truncated data, the lifetime T is observable only when U≤T≤V, where U and V are the left-truncated and right-truncated time, respectively. In some situations, the lifetime T also suffers interval censoring. Using the EM algorithm of Turnbull [The empirical distribution function with arbitrarily grouped censored and truncated data, J. R. Stat. Soc. Ser. B 38 (1976), pp. 290–295] and iterative convex minorant algorithm [P. Groeneboom and J.A. Wellner, Information Bounds and Nonparametric Maximum Likelihood Estimation, Birkhäuser, Basel, 1992], we study the performance of the nonparametric maximum-likelihood estimates (NPMLEs) of the distribution function of T. Simulation results indicate that the NPMLE performs adequately for the finite sample.     

A new ranked set sampling protocol for the signed rank test

In this paper, we construct a new ranked set sampling protocol that maximizes the Pitman asymptotic efficiency of the signed rank test. The new sampling design is a function of the set size and independent order statistics. If the set size is odd and the underlying distribution is symmetric and unimodal, then the new sampling protocol quantifies only the middle observation. On the other hand, if the set size is even, the new sampling design quantifies the two middle observations. This data collection procedure for use in the signed rank test outperforms the data collection procedure in the standard ranked set sample. We show that the exact null distribution of the signed rank statistic W_RSS⁺ based on a data set generated by the new ranked set sample design for odd set sizes is the same as the null distribution of the simple random sample signed rank statistic W_SRS⁺ based on the same number of measured observations. For even set sizes, the exact null distribution of W_RSS⁺ is simulated.     

Nonparametric estimation of the bivariate cdf for arbitrarily censored data

Right, left or interval censored multivariate data can be represented by an intersection graph. Focussing on the bivariate case, the authors relate the structure of such an intersection graph to the support of the nonparametric maximum likelihood estimate (NPMLE) of the cumulative distribution function (CDF) for such data. They distinguish two types of non‐uniqueness of the NPMLE: representational, arising when the likelihood is unaffected by the distribution of the estimated probability mass within regions, and mixture, arising when the masses themselves are not unique. The authors provide a brief overview of estimation techniques and examine three data sets.     

A family of shrinkage estimators for Weibull shape parameter in censored sampling

In this paper a new class of shrinkage estimators has been introduced for the shape parameter in an independently identically distributed two-parameterWeibull model under censored sampling. The main idea is to incorporate the prior guessed value by correcting the standard estimator, which is essentially an unbiased estimator, with optimally weighted ratios of the guessed value and the standard estimator, instead of considering a convex combination of the standard estimator and the difference of the guessed value and the standard estimator. The resulting estimator dominates the standard estimator in a surprisingly large neighborhood of the guessed value. The suggested estimator has also been compared with the minimum mean squared error estimator and a class of estimators suggested by Singh and Shukla in IAPQR Trans 25(2), 107–118, 2000. It is found that the suggested class of estimators has lesser bias as well as lesser mean squared error than its competitors subject to certain conditions.      

Accuracy of approximate progressively censored reliability sampling plans for exponential models

Reliability sampling plans provide an efficient method to determine the acceptability of a product based upon the lifelengths of some test units. Usually, they depend on the producer and consumer’s quality requirements and do not admit closed form solutions. Acceptance sampling plans for one- and two-parameter exponential lifetime models, derived by approximating the operating characteristic curve, are presented in this paper. The accuracy of these approximate plans, which are explicitly expressible and valid for failure and progressive censoring, is assessed. The approximation proposed in the one-parameter case is found to be practically exact. Explicit lower and upper bounds on the smallest sample size are given in the two-parameter case. Some additional advantages are also pointed out.     

Interval estimation for the Pareto distribution based on the progressive Type II censored sample

For the complete sample and the right Type II censored sample, Chen [Joint confidence region for the parameters of Pareto distribution. Metrika 44 (1996), pp. 191–197] proposed the interval estimation of the parameter θ and the joint confidence region of the two parameters of Pareto distribution. This paper proposed two methods to construct the confidence region of the two parameters of the Pareto distribution for the progressive Type II censored sample. A simulation study comparing the performance of the two methods is done and concludes that Method 1 is superior to Method 2 by obtaining a smaller confidence area. The interval estimation of parameter ν is also given under progressive Type II censoring. In addition, the predictive intervals of the future observation and the ratio of the two future consecutive failure times based on the progressive Type II censored sample are also proposed. Finally, one example is given to illustrate all interval estimations in this paper.