Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Optimal Design of Step Stress Partially Accelerated Life Test under Progressive Type-II Censored Data with Random Removal for Gompertz Distribution

Dina H. Abdel Hady

Department of Statistics, Mathematics and Insurance, Faculty of Commerce, Tanta University

Abstract: This paper deals with random removal of progressively type II censored data. The removal of the data is assumed to follow a binomial or a uniform distribution, and the life time testing is assumed to follow a Gompertz distribution. Parameters of these distributions are estimated using the Maximum Likelihood estimation procedure. Fisher information matrix is used to estimate the asymptotic mean squared error and to construct confidence intervals of model parameters. The optimal partially accelerated lifetime testing (PALT) is estimated by minimizing the Generalized Asymptotic Variance (GAV). Simulation study is performed to Clarification the statistical properties of the parameters. A simulation results reveal that for the fixed values of the parameters, the error and optimal time decrease with increasing sample size n; estimates of binomial are more stable with a relatively small error with increasing sample size; and the test design is robust and works well for binomial removal.

Keywords: partially accelerated life test, PALT, SS-PALT, progressively type II censored data, Gompertz distribution optimal design, D-optimality

Received November 12, 2018; Revised December 19, 2018; Accepted January 01, 2019

Citation: Dina H. Abdel Hady, Optimal Design of Step Stress Partially Accelerated Life Test under Progressive Type-II Censored Data with Random Removal for Gompertz Distribution. American Journal of Applied Mathematics and Statistics, 7 (1). 37-42. January 2019. doi:10.12691/ajams-7-1-6

1. Introduction

Traditional life testing may show no failure or few failures of some highly reliable units because of some severe conditions (stresses) that may occur in the form of pressure, voltage, temperature, vibration, load, cycling rate, etc. In such situations life testing has to be performed at higher than usual conditions in order to obtain failures as fast as possible. The data collected under stresses or severe conditions should not be used to estimate the life distribution at normal use assumptions. There are three types of stress, a) step stress; b) progressive stress, and c) constant stress; Ismail, [1] when testing is conducted under stress, it is called either “Accelerated Life Testing (ALT)” or “Partially Accelerated Life Testing (PALT)”. In ‘Alt” testing, units are put under stress to reach more failure in short time, and the mathematical stress model is known; while in “PALT”, data can be extrapolated to normal use conditions, and the statistical model parameters could be estimated and step-stress could be applied. Nelson [2] proposed the step stress PALT (SS-PALT) where units are tested at normal use conditions within a specified time interval. Censored schemes arise when not all lifetime testing units are observed. When a unit does not fail within that time interval, it is put under stress, until failure or until the termination of the time interval (censored scheme). Rao [3] and Balakrishnan and Aggarwala [4] have mentioned that step stress saves money, effort and time. There are two censoring scheme types. Type I arises when the experiment continues up to a pre-specified time T and Type-II censoring scheme requires the experiment to continue until a pre-specified number of failures m ≤ n occur.

These two types do not allow for the removal of any unit until life testing is terminated. Nevertheless, this allowance may be desirable to compromise between at least one observation is sought and a reduced time interval for the experimentation. At this situation, progressive censoring is needed.

In the progressive Type II censoring scheme, the experimenter selects n (iid) units for life testing; and observes the occurrence of the first failure at time and this is removed from testing. When the second failure occurs at time one of the surviving units are randomly selected and removed from the test, the experiment terminates when the mth failure occurs at time and the remaining surviving units are all removed from the Test Progressive type-II censoring scheme with fixed number of removals is considered by Cohen [5] and Cohen and Norgaard [6]. In some reliability experiments and when the number of removals is not fixed, a progressive censoring with random removals model could be considered [7,8,9]. The random removal is assumed to follow the binomial distribution, and used to estimate the parameters of the reliability model.

Ismail [1] discusses step-stress partially-accelerated model. The model differs than the progressive type II censoring scheme in that the lifetime of the test follows the exponential distribution rather than the binomial.

This paper considers the step-stress partially-accelerated model, where the removals are assumed to follow the binomial and the life time-testing follows the Gompertz distribution. The optimal stress change time is also being determined by minimizing the generalized asymptotic variance of the MLE parameters. In Section 2, assumptions of the Gompretz distribution are given; Section 3 presents the assumptions of the partially accelerated model. Estimation of model parameters is given in Section 3; simulation study results are given in Section 4.

2. The Gompretz Distribution

The probability density function of the Gompretz distribution takes the form:

[picture] | (1)

And the cumulative distribution function is:

[picture] | (2)

The following assumptions are used:

1. n identical and independent units are selected for life testing.

2. The lifetime of each unit has Gompertz distribution.

The following steps are to be followed:

a) Each of the n units is first run under normal use condition. If it does not fail or remove from the test by a pre-specified time. it is put under accelerated condition (stress).

b) At the ith failure a random number of the surviving units, are randomly selected and removed from the test.

c) Finally, at the mth failure the remaining surviving units are all removed from the test and the test is terminated.

The lifetime, say X, of a unit under SS-PALT can be rewritten as

Where, the pdf of X is given by

[picture] | (3)

In addition, the survival functions under normal and accelerate use conditions respectively is given by

[picture] | (4)

And

[picture] | (5)

3. Estimation of Parameters

Let , denote the observation obtained form a progressively type-II censored sample with random removals in a step-stress PALT. Here

Given a pre-determined number of removals the conditional likelihood function of the observations takes the following form

[picture] | (6)

Equations (1) and (3) are inserted in (6) and simplify, we get

[picture] | (7)

3.1. Parameter Estimation with Binomial Removals

Given that the number of units removed from the test at each failure time follows a binomial distribution and for Thus, the number of units removed at each failure time follows a binomial distribution such that

And for

[picture] | (8)

The likelihood of the sample of size n is given as follows

[picture] | (9)

Where

That is

[picture] | (10)

The log-likelihood function [Eq. 10] can be written as follows:

[picture] | (11)

[picture] | (12)

To obtain the estimate of the the first partial derivatives with respect to and of Equation (12) are obtained and equated to zero as follows:

[picture] | (13)

[picture] | (14)

[picture] | (15)

[picture] | (16)

From Equation (16), is estimated as follows

[picture] | (17)

There is no closed-from solution to this system of equations 13, 14 and 15. using the Newton-Raphson iteration method, for root finding. is Iteratively as follows:

[picture] | (18)

where g is the vector of normal equations,

and

and G [Equation 17] is the matrix of second derivatives

[picture] | (19)

Where

[picture] | (20)

[picture] | (21)

[picture] | (22)

[picture] | (23)

[picture] | (24)

[picture] | (25)

Convergence of the Newton-Raphson algorithm for the estimates of and depends on the tolerance limit change with each successive iteration, to and

Numerically inverting the above G matrix above, we easily obtain Fisher Information matrix, i.e.

The approximate two sided confidence intervals for and can be, constructed as follows:

3.2. Parameter Estimation with the Uniform Removals

The model assumes that removed units are independent and that the probability of each removed unit is the same, such that:

[picture] | (26)

And for

The joint probability distribution of

is given by

[picture] | (27)

where

The maximum likelihood estimators can be derived directly by maximizing the equations (12) and then solving for equations (13), (14) and (15).

3.3. D-optimality

Fisher's information matrix is used to determine t in the SS-PALT type II progressive censoring scheme proposed criterion is based on the determinant of Fisher's information matrix. The generalized asymptotic variance (GAV) is the reciprocal of the determinant of Fisher's information matrix F [10], thus

[picture] | (28)

The D-optimality criterion is the optimal value of that maximizes the determinant of the Fisher's information matrix F and minimizes the GAV.

Table 1. Simulation study results with Binomial Removals for θ=2, α=2.1, β=2.3, τ=3

Table 2. Simulation study results with uniform Removals for θ=2, α=2.1, β=2.3, τ=3

4. Simulation Study

A simulation study is performed to study the properties of the estimators using ML method; the study involves the computation of Mean squared errors (MSEs), the construction of confidence intervals for different sample sizes; and the determination of, the optimal stress change time. The following steps were followed:

a)Value of n and m to be specified.

b)Value of the parameters were set as: .

c) A random sample with size n and censoring size m were generated, with random removals, from the random variable X given by (3).

d) Generate a group value and also where and

e) ML estimates were computed, for n= 20, 50, 80 and 100,

f) The mean squared error (MSE), the 95% confidence interval of parameters and the bias associated with the MLE of the parameters, optimal value of the Optimal GAV of the MLEs of the model parameters are obtained numerically for each sample size.

5. Conclusions

The (SS-PALT) under progressive type-II censored data with binomial and uniform random removal assuming Gompertz distribution was studied. The Newton-Raphson method is applied to obtain the optimal stress-change time which minimizes the Generalized Asymptotic Variance (GAV).

The maximum Likelihood estimation procedure was used the estimation of model parameters, mean squared error and the optimum plan for the binomial and uniform removals for different sample size were computed and shown. Simulation results show that the error and the optimal time decrease as sample size increase, when parameters are fixed. that both the average value of and the average value of GAV for type-II progressive censoring are getting close to those of complete sample with the bigger m and close faster for bigger n. Hence from the numerical result we can conclude that estimates of binomial are more stable with relatively small error with increasing sample size. Therefore, the test design obtained here is robust design and work well for binomial removal.

References

[1]   | A. Ismail, (2009). Optimal Design of Step-Stress Life Test with Progressively type-II Censored Exponential Data, International Mathematical Forum, 4 , no 40, 1963-1976.
[2]   | W. Nelson. (1990). Accelerated Life Testing: Statistical Models, Data Analysis and Test Plans, John Wiley and sons, New York.
[3]   | R. Rao, (1992), Equivalence of the tampered random variables and tampered failure rate models in ALT for a class of life distribution having the setting the clock back to zero property, Communication in Statistics – Theory and Methods, Vol. 21, No. 3, 647-664.
[4]   | N. Balakrishnan, and R. Aggarwala. (2000). Progressive Censoring: Theory, Methods, and Applications, Birkhauser, Boston.
[5]   | C. Cohen. (1963). Progressively Censored Samples in the Life Testing, Technometerics, 5, 327-339.
[6]   | C. Cohen and N. J. Norgaard. (1977). Progressively Censored Sampling in the Three Parameter Gamma Distribution, Technometerics, 19 , 333-340.
[7]   | H. K. Yuen. and S. K. (1996). Tse, Parameters Estimation for Weibull distributed Lifetime under progressive Censoring with random removals, Journal of Statistical Computation and Simulation, 55, 57-71.
[8]   | C. Wu, S. F. Wu, and, H. Y. Chan. (2004). MLE and the Estimated Expected Test Time for Pareto Distribution under Progressive Censoring Data, International Journal of Information Management Sciences, 15(3), 29-42.
[9]   | Tse, S. K., Yang, C. and. Yuen, H. K., 2000. Statistical Analysis of Weibull Distributed Lifetime Data under Type-II Progressive Censoring with Binomial Removals. Journal of Applied Statistics, vol. 27, no. 8, pp. 1033-1043.
[10]   | S. Bai, J. G. Kim and Y. R. Chun, (1993). Design of Failure-Censored Accelerated Life-Test Sampling Plans for Lognormal and Weibull Distributions, Eng. Opt. Vol. 21, 197-212.