A STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS
|
|
- Denis Watson
- 6 years ago
- Views:
Transcription
1 A STATISTICAL TEST FOR MONOTONIC AND NON-MONOTONIC TREND IN REPAIRABLE SYSTEMS Jan Terje Kvaløy Department of Mathematics and Science, Stavanger University College, P.O. Box 2557 Ullandhaug, N-491 Stavanger, Norway, Bo Lindqvist, Håkon Malmedal Department of Mathematical Sciences, Norwegian University of Science and Technology, N-7491 Trondheim, Norway, ABSTRACT A statistical test for trend in recurrent event data, for instance failure data from a repairable system, based on the general null hypothesis of a renewal process and with power against both monotonic and nonmonotonic trends is proposed. Simulations show that the test has favorable qualities as a general test for trend by having good power properties against both monotonic and nonmonotonic alternatives. 1 INTRODUCTION In recurrent event data it is often of interest to detect possible systematic alterations in the pattern of events. An example is a repairable system for which it is important to detect possible systematic changes in the pattern of failures. Such changes can for instance be caused by various aging effects or reliability growth. We say that there is a trend in the pattern of failures if the inter-failure times tend to alter in some systematic way, which means that the inter-failure times are not identically distributed. By using statistical trend tests it is possible to decide whether such an alteration is statistically significant or not. In this paper we shall focus on the repairable system example, but the methods presented are relevant for any kind of recurrent event data which can be modeled by the same statistical models as used for repairable systems in this paper. A trend in the pattern of failures of a repairable system can be either monotonic, corresponding to an improving or deteriorating system, or non-monotonic, corresponding to for instance cyclic variations or a so called bathtub trend characterizing a system going through the three phases infant illness, useful life and wear out. The purpose of the paper is to present a general test for trend with power against both monotonic and nonmonotonic trends. Common statistical models for systems without trend are the homogeneous Poisson process (HPP) and the renewal process (RP). The HPP, assuming independent identically exponentially distributed interarrival times, is a special case of the far more general RP, assuming independent identically distributed interarrival times. For both models there exist a number of tests with power against monotonic trends but few tests with power against nonmonotonic trends. A test with power both against monotonic and nonmonotonic trends based on an HPP model, called the Anderson-Darling test for trend, was studied by Kvaløy & Lindqvist [1]. The problem, however, with this test and many other tests based on the null hypothesis of an HPP is that even if the null hypothesis is rejected, it does not necessarily indicate a trend, just that the process is not an HPP. An RP may for instance cause rejection (Lindqvist et al. [2]).
2 Elvebakk [3] has studied how the Anderson-Darling test for trend could be generalized to a test of the null hypothesis RP by use of resampling techniques. In the present paper we discuss how the Anderson-Darling test for trend can be directly generalized to a test of the null hypothesis RP, without using resampling techniques. This yields a straightforwardly calculated test, which is a generally applicable test for trend with good power properties both against monotonic and non-monotonic trends. We also show that the well known Lewis-Robinson test (Lewis & Robinson [4]) follows from a similar derivation taking the same starting point, and we compare the generalized Anderson-Darling test for trend to the Lewis-Robinson test and the test ( [5]) in a simulation study. 2 MODELS AND TESTS 2.1 Notation and terminology In this paper attention is restricted to one system which is observed from time t =. The successive failure times are denoted T 1,T 2,... and the corresponding interfailure or interarrival times are denoted X 1,X 2,... where X i = T i T i 1, i =1, 2,... (1) Another way to represent the same information is by the counting process representation N(t) = number of failures in (,t]. (2) A system is said to be failure censored if it is observed until a given number n of failures has occurred, and time censored if it is observed until a given time τ. We say that a system inhibits no trend if the marginal distributions of all interarrival times are identical, otherwise there is a trend. If the expected length of the interarrival times is monotonically increasing or decreasing with time, corresponding to an improving or a deteriorating system, there is a monotone trend (a decreasing or an increasing trend), otherwise the trend is nonmonotone. 2.2 Renewal process The stochastic process T 1,T 2,... is an RP if the interarrival times X 1,X 2,... are independent and identically distributed. The conditional intensity of an RP given the history F t up to, but not including time t, can be written γ(t F t )=z(t T N(t ) ) (3) where z(t) is the hazard rate of the distribution of the interarrival times and T N(t ) is the last failure time strictly before time t. 2.3 Trend-renewal process The trend-renewal process (TRP) is presented by Lindqvist et al. [6], and is a generalization of the RP to a process with trend. The RP, HPP and nonhomogenous Poisson process (NHPP) are all special cases of this process. Let λ(t) be a nonnegative function, called the trend function, defined for t andlet Λ(t) = t λ(u)du. Then the process T 1,T 2,... is a TRP if Λ(T 1 ), Λ(T 2 ),... is an RP. The conditional intensity of a TRP may be written (Lindqvist et al. [6]) γ(t F t )=z(λ(t) Λ(T N(t ) ))λ(t) (4)
3 where z(t) is the hazard rate of the distribution of the underlying RP. In other words the intensity depends on both the age of the system and the time since last failure. The NHPP follows as the special case when the RP part of the TRP is an HPP with intensity 1. In this case the hazard rate z(t) 1 and thus γ(t F t )=λ(t). The special case of an RP follows when λ(t) is constant, and the HPP follows when both λ(t) and z(t) are constant. 2.4 Generalized Anderson-Darling test for trend The construction of the generalized Anderson-Darling test for trend is based on the independent increment property of the RP and Donsker s theorem. Let ξ 1,...,ξ n be n independent identically distributed random variables with E(ξ i )= and Var(ξ i )=σ 2,let[u] denote the integer part of u, letz j = ξ ξ j and define Q n (t) = Z [nt] σ nt [nt] + n σ n ξ [nt]+1 (5) for t 1. Then according to Donsker s theorem (e.g Billingsley [7]) Q d n W,where d denotes convergence in distribution and W is a Brownian motion. Now consider a failure censored RP with interarrival times X 1,...,X n. Define ξ i = X i µ where µ = E(X i )andletz j = ξ ξ j = j i=1 X i jµ = T j jµ. It follows from Donsker s theorem that V d n,µ,σ W where V n,µ,σ (t) = T [nt] [nt]µ σ n + nt [nt] σ n (X [nt]+1 µ) (6) and σ 2 = Var(X i ). Further defining W (t) =W (t) tw (1), which is a Brownian bridge, we will have that Vn,σ d W,where n X Vn,σ(t) =V n,µ,σ (t) tv n,µ,σ (1) = σ (T [nt] +(nt [nt]) X [nt]+1 t) 2 (7) T n T n where X = T n /n. Notice that µ cancels in (7). Replacing σ by a consistent estimator ˆσ, we conclude that V n,ˆσ d W. It follows from the above that under the null hypothesis of no trend, Vn,ˆσ (t) will be approximately a Brownian bridge. On the other hand, if there is a trend in the data, V n,ˆσ (t) is likely to deviate from the Brownian bridge. Thus tests for trend can be based on measures of deviation from a Brownian bridge of Vn,ˆσ (t). A suggestion is to use the Anderson-Darling type statistic GAD = 1 (Vn,ˆσ (t))2 1 t(1 t) dt d (W (t)) 2 dt. (8) t(1 t) The test based on this statistic will be called the generalized Anderson-Darling test for trend, and the limit distribution is the usual limit distribution of the Anderson-Darling statistic (Anderson & Darling [8], [9]). It is easily realized that the statistic will have sensitivity against deviations from the Brownian bridge caused both by monotonic and nonmonotonic trends. Straightforward but tedious calculations show that GAD = n X 2 [ n q 2 ˆσ 2 i ln( i i 1 )+(q i + r i ) 2 ln( n i +1 ] r2 i n i n ) (9) i=1
4 where q i =(T i ix i )/T n and r i = nx i /T n 1 Various consistent estimators of σ 2 can be used. Rather than using the usual estimator, S 2 = n i=1 (X i X) 2 /(n 1), we recommend using the estimator ˆσ 2 = n 1 1 (X i+1 X i ) 2. (1) 2(n 1) i=1 This estimator is consistent under the null hypothesis of independent and identically distributed interarrival times, but tends to be smaller than S 2 in cases with positive dependence between subsequent interarrival times. Thus using this variance estimator yields a test statistic with larger power against alternatives like monotonic and bathtub shaped trends than using the estimator S 2. Finally, empirical explorations have shown that the small sample level properties of the test can be improved by multiplying the test statistic by an (asymptotically negligible) factor (1 4/n) (Malmedal [1]). Thus the final form of the proposed test statistic is 2 (n 4) X GAD = ˆσ 2 [ n qi 2 ln( i i=1 i 1 )+(q i + r i ) 2 ln( n i +1 n i r2 i n ) where q i =(T i ix i )/T n, r i = nx i /T n 1andˆσ 2 is calculated from (1). Using the limit distribution of the test statistic (Anderson & Darling [8], [9]), the null hypothesis is rejected on the 5% level if GAD > For time censored data we propose conditioning on the observed number of failures, and treating the data as if they were data from a failure censored system. 2.5 The Lewis-Robinson test Other trend tests can be constructed by following the same approach but using other measures of deviation from the Brownian bridge of Vn,ˆσ (t). A possibility is for instance to consider the statistic 1 1 Vn,ˆσ (t)dt d W (t)dt (12) where 1 W (t)dt is normally distributed with expectation and variance 1/12. This statistic will primarily have power against deviations from a Brownian bridge caused by monotonic trends. By scaling the test statistic to be asymptotically standard normally distributed, straightforward calculations yield the test statistic 1 LR Vn,ˆσ = (t) dt = X n 1 i=1 T i n 1T 2 n 12 ˆσ T n. (13) n 12 This test statistic is a variant of the classical Lewis-Robinson test statistic (Lewis & Robinson [4]), often written LR = X n 1 i=1 T i n 1T 2 n (14) S T n 1 n 12 For time censored data n 1 is replaced by n and T n by τ. The original derivation by Lewis & Robinson [4] differs from the derivation presented here by taking a different starting point and using other arguments to prove the asymptotic normality. ] (11)
5 2.6 The test A classical trend test is the rank test developed by [5]. This is a test of the null hypothesis RP versus the alternative of a monotonic trend. The test statistic is computed by counting the number of reverse arrangements, M, among the interarrival times X 1,...,X n. We have a reverse arrangement if X i <X j for i<j,andm is thus M = n 1 n i=1 j=i+1 I(X i <X j ) (15) where I(A) = 1 or according to whether the event A occurs or not. Under the null hypothesis M is approximately normally distributed with expectation n(n 1)/4 and variance (2n 3 +3n 2 5n)/72 for n 1. For n<1 there exist tables. 3 SIMULATION STUDY In this section the generalized Anderson-Darling test for trend (11), the Lewis-Robinson test (14) and the test (15) are compared in a simulation study. Rejection probabilities are estimated by simulating 1 data sets for each choice of model and parameter values, and recording the relative number of rejections of each test. In the curves of (Fig. 2) and (Fig. 4) the rejection probabilities are simulated and plotted in a number of points, with straight lines drawn between the points. Let ˆp denote the estimated rejection probability. Then the standard deviation of ˆp is ˆp(1 ˆp)/1.16. All simulations are done in C++. The nominal significance level has been set to 5%, and only simulations of failure censored processes are reported. Generating data T 1,T 2,... from a TRP is straightforward. Recall that the process T 1,T 2,... is a TRP if Λ(T 1 ), Λ(T 2 ),... is an RP. Then failure times T 1,T 2,... from a TRP are generated by first generating arrival times θ 1,θ 2,... from the relevant RP and transforming these to T 1 =Λ 1 (θ 1 ),T 2 =Λ 1 (θ 2 ),... Arrival times of the RP are generated by simulating interarrival times θ i θ i 1 from the RP-distribution. 3.1 Level properties First the level properties of the tests are studied by generating datasets from an HPP with intensity 1 (i.e. the interarrival times are independent and identically exponentially distributed with expectation 1). In (Fig. 1) the simulated level of the tests for samples of size 4 to 8 are reported. The plots in (Fig. 1) show that all tests have satisfying level properties for samples of sizes larger that around 15. For very small sample sizes all.9 Generalized Anderson Darling.9 Lewis Robinson Sample size Sample size Sample size Figure 1: Simulations of HPPs with intensity 1 and sample sizes from 4 to 8.
6 tests tend to be conservative and the generalized Anderson-Darling test for trend is the most conservative test. For the test the normal approximation has been used for all sample sizes. The level properties of this test for samples smaller than 1 are better if critical values from tables are used. 3.2 Power properties Monotonic trend Datasets with a monotonic trend are generated by simulating data from TRPs with the underlying RP distribution being Weibull and the trend function λ(t), governing the trend of the process, being on the power law form λ(t) =bt b 1. The intensity of the process can then be written γ(t F t )=αβ(t b TN(t ) b )β 1 bt b 1. (16) For samples of size 25 the rejection probability as a function of b is simulated, where b<1 corresponds to a decreasing trend, b = 1 corresponds to no trend and b>1 corresponds to an increasing trend. This is repeated for two different values of the shape parameter β of 1 Overdispersed TRP 1 Underdispersed TRP GAD LR.2.1 GAD LR b Figure 2: Simulations of TRPs with underlying Weibull RPs with shape parameters respectively.75 (overdispersed TRP) and 1.5 (underdispersed TRP), trend function λ(t) =bt b 1 and samples of size 25. b the underlying Weibull distribution, β =.75 and β = 1.5, corresponding respectively to a process which is overdispersed or underdispersed relative to an NHPP. The scale parameter α is set to 1. The results are displayed in (Fig. 2). We see in (Fig. 2) that the generalized Anderson-Darling test for trend and the Lewis- Robinson test have fairly similar power properties, with the generalized Anderson-Darling test being a bit more powerful than the Lewis-Robinson against increasing trend and opposite against decreasing trend. The test is clearly less powerful than the other tests against increasing trend but is the most powerful test against decreasing trend in the overdispersed case Nonmonotonic trend Datasets with a bathtub trend are generated by simulating data from TRPs with trend function on the form displayed in (Fig. 3). Here d represents the average of the function λ(t) over the interval [,τ]. The degree of bathtub shape can be expressed by the parameter c, withc = corresponding to a horizontal line (no trend).
7 (t) a d b c e fi e fi t Figure 3: Bathtub-shaped trend function. For samples of size 25 and 75, the rejection probability as a function of c is simulated with τ set to a value which implies that the expected number of failures in [,τ] becomes respectively approximately 25 and approximately 75. The parameter e is adjusted according to c to make the expected number of failures in each phase (decreasing, no, increasing trend) equal. The shape parameter of the underlying Weibull distribution is set to β = 1.5 (underdispersed). The results are displayed in (Fig. 4) failures GAD LR failures GAD LR c Figure 4: Simulations of TRPs with underlying Weibull RPs with shape parameter 1.5 (underdispersed TRP), trend function displayed in (Fig. 3) and samples of size respectively 25 and 75. We see in (Fig. 4) that the generalized Anderson-Darling test for trend has clearly better power properties than the other tests in these cases. Also note that the generalized Anderson-Darling test for trend improves substantially with increasing sample size while the other tests do not improve much. Similar patterns are seen if we consider overdispersed TRPs (Malmedal [1]). One difference seen for overdispersed TRPs is that none of the tests are particular powerful for small sample sizes due to the larger variance of the underlying RP (Malmedal [1]). 3.3 Discussion In the simulations presented the generalized Anderson-Darling test for trend is the best test against increasing and bathtub shaped trend, while the other tests are slightly better against decreasing trend. It is not surprising that the generalized Anderson-Darling test is better than the other tests against bathtub shaped trend since this test, contrary to c
8 the other tests, is constructed to have power against both nonmonotonic alternatives as well as monotonic ones. The fact that the generalized Anderson-Darling test also performs at least as good as the other tests against monotonic alternatives makes this test a very attractive choice as a trend test for general use. The only drawback with the generalized Anderson-Darling test for trend (11) is that it is too conservative for very small sample sizes. In such cases other level correcting factors than 1 4/n could be used. Another alternative is to use a bootstrap version of the test (Elvebakk [3]; Malmedal [1]). The bootstrap test has good level properties even for small samples and it also turns out that bootstrapping in many cases improves the power properties of the test, in some cases substantially (Malmedal [1]). 4 CONCLUSION A test for trend based on the general null hypothesis of an RP which has power against both monotonic and nonmonotonic trends has been presented. This test has good overall properties compared to existing trend tests and is recommended for general use. For very small sample sizes it is recommended to use a resample version of the test. ACKNOWLEDGEMENT We would like to thank associate professor Michael Kosorok for introducing us to how the Anderson-Darling test could be generalized from the null hypothesis HPP to RP by using the independent increment property of the problem. REFERENCES [1] Kvaløy, J. T. & Lindqvist, B. H., TTT-based Tests for Trend in Repairable Systems Data, Reliability Engineering and System Safety, 6, (1998), [2] Lindqvist, B. H., Kjønstad, G. A. & Meland, N., Testing for trend in repairable systems data, Proceedings of ESREL 94, La Baule, France, May 3 - June 3, [3] Elvebakk, G., Analysis of Repairable Systems Data: Statistical Inference for a Class of Models Involving Renewals, Heterogeneity and Time Trends, PhD thesis, Department of Mathematical Sciences, Norwegian University of Science and Technology, [4] Lewis, P. A. & Robinson, D. W., Testing for a monotone trend in a modulated renewal process, in F. Proschan & R.J Serfling (eds), Reliability and Biometry, , SIAM, Philadelphia, [5], H. B., Nonparametric tests against trend, Econometrica, 13, (1945), [6] Lindqvist, B. H., Elvebakk, G. & Heggland, K., The Trend-Renewal Process for Statistical Analysis of Repairable Systems, submitted for publication, [7] Billingsley, P., Convergence of Probability Measures, John Wiley, New York, [8] Anderson, T. W. & Darling, D. A., Asymptotic theory of certain goodness of fit criteria based on stochastic processes, Annals of Mathematical Statistics, 23, (1952), [9] Anderson, T. W. & Darling, D. A., A test of goodness of fit, Journal of the American Statistical Association, 49, (1954), [1] Malmedal, H., Trend testing, Diploma Thesis, (In Norwegian), Department of Mathematical Sciences, Norwegian University of Science and Technology, 21.
Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim
Tests for trend in more than one repairable system. Jan Terje Kvaly Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim ABSTRACT: If failure time data from several
More informationarxiv: v1 [stat.me] 22 Feb 2018
A Class of Tests for Trend in Time Censored Recurrent Event Data Jan Terje Kvaløy Department of Mathematics and Physics University of Stavanger, Norway arxiv:182.8339v1 [stat.me] 22 Feb 218 Bo Henry Lindqvist
More informationInvestigation of some tests for homogeneity of intensity with applications to insurance data
U.U.D.M. Project Report 2011:26 Investigation of some tests for homogeneity of intensity with applications to insurance data Sara Gustin Examensarbete i matematik, 15 hp Handledare och examinator: Jesper
More informationIntroduction to repairable systems STK4400 Spring 2011
Introduction to repairable systems STK4400 Spring 2011 Bo Lindqvist http://www.math.ntnu.no/ bo/ bo@math.ntnu.no Bo Lindqvist Introduction to repairable systems Definition of repairable system Ascher and
More informationLink to published article: (Access to content may be restricted)
Kvaløy, J.T. (2002) Covariate Order Tests for Covariate Effect. Lifetime Data Analysis, 8(1), pp. 35-51 Link to published article: http://link.springer.com/article/10.1023/a:1013518815447 (Access to content
More informationCHAPTER 1 A MAINTENANCE MODEL FOR COMPONENTS EXPOSED TO SEVERAL FAILURE MECHANISMS AND IMPERFECT REPAIR
CHAPTER 1 A MAINTENANCE MODEL FOR COMPONENTS EXPOSED TO SEVERAL FAILURE MECHANISMS AND IMPERFECT REPAIR Helge Langseth and Bo Henry Lindqvist Department of Mathematical Sciences Norwegian University of
More informationA nonparametric monotone maximum likelihood estimator of time trend for repairable systems data
A nonparametric monotone maximum likelihood estimator of time trend for repairable systems data Knut Heggland 1 and Bo H Lindqvist Department of Mathematical Sciences Norwegian University of Science and
More informationRepairable Systems Reliability Trend Tests and Evaluation
Repairable Systems Reliability Trend Tests and Evaluation Peng Wang, Ph.D., United Technologies Research Center David W. Coit, Ph.D., Rutgers University Keywords: repairable system, reliability trend test,
More informationA novel repair model for imperfect maintenance
IMA Journal of Management Mathematics (6) 7, 35 43 doi:.93/imaman/dpi36 Advance Access publication on July 4, 5 A novel repair model for imperfect maintenance SHAOMIN WU AND DEREK CLEMENTS-CROOME School
More informationTESTING A HOMOGENEITY OF STOCHASTIC PROCESSES
K Y B E R N E I K A V O L U M E 4 3 0 0 7, N U M B E R 4, P A G E S 4 5 4 3 0 ESING A HOMOGENEIY OF SOCHASIC PROCESSES Jaromír Antoch and Daniela Jarušová he paper concentrates on modeling the data that
More informationExamination paper for TMA4275 Lifetime Analysis
Department of Mathematical Sciences Examination paper for TMA4275 Lifetime Analysis Academic contact during examination: Ioannis Vardaxis Phone: 95 36 00 26 Examination date: Saturday May 30 2015 Examination
More informationAn Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems
An Integral Measure of Aging/Rejuvenation for Repairable and Non-repairable Systems M.P. Kaminskiy and V.V. Krivtsov Abstract This paper introduces a simple index that helps to assess the degree of aging
More informationAN INTEGRAL MEASURE OF AGING/REJUVENATION FOR REPAIRABLE AND NON REPAIRABLE SYSTEMS
R&RAA # 1 (Vol.1) 8, March AN INEGRAL MEASURE OF AGING/REJUVENAION FOR REPAIRABLE AND NON REPAIRABLE SYSEMS M.P. Kaminskiy and V.V. Krivtsov Abstract his paper introduces a simple index that helps to assess
More informationSoftware Reliability Growth Modelling using a Weighted Laplace Test Statistic
Software Reliability Growth Modelling using a Weighted Laplace Test Statistic Yan Luo Torsten Bergander A. Ben Hamza Concordia Institute for Information Systems Engineering Concordia University, Montréal,
More informationModelling the risk process
Modelling the risk process Krzysztof Burnecki Hugo Steinhaus Center Wroc law University of Technology www.im.pwr.wroc.pl/ hugo Modelling the risk process 1 Risk process If (Ω, F, P) is a probability space
More informationTime-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation
Time-varying failure rate for system reliability analysis in large-scale railway risk assessment simulation H. Zhang, E. Cutright & T. Giras Center of Rail Safety-Critical Excellence, University of Virginia,
More informationComputer Simulation of Repairable Processes
SEMATECH 1996 Applied Reliability Tools Workshop (ARTWORK IX) Santa Fe Computer Simulation of Repairable Processes Dave Trindade, Ph.D. Senior AMD Fellow Applied Statistics Introduction Computer simulation!
More informationStochastic Renewal Processes in Structural Reliability Analysis:
Stochastic Renewal Processes in Structural Reliability Analysis: An Overview of Models and Applications Professor and Industrial Research Chair Department of Civil and Environmental Engineering University
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued
Introduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations Research
More informationReliability Growth in JMP 10
Reliability Growth in JMP 10 Presented at Discovery Summit 2012 September 13, 2012 Marie Gaudard and Leo Wright Purpose of Talk The goal of this talk is to provide a brief introduction to: The area of
More informationManual for SOA Exam MLC.
Chapter 10. Poisson processes. Section 10.5. Nonhomogenous Poisson processes Extract from: Arcones Fall 2009 Edition, available at http://www.actexmadriver.com/ 1/14 Nonhomogenous Poisson processes Definition
More informationA comparison study of the nonparametric tests based on the empirical distributions
통계연구 (2015), 제 20 권제 3 호, 1-12 A comparison study of the nonparametric tests based on the empirical distributions Hyo-Il Park 1) Abstract In this study, we propose a nonparametric test based on the empirical
More informationMahdi karbasian* & Zoubi Ibrahim
International Journal of Industrial Engineering & Production Research (010) pp. 105-110 September 010, Volume 1, Number International Journal of Industrial Engineering & Production Research ISSN: 008-4889
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 29, 2014 Introduction Introduction The world of the model-builder
More informationMaintenance of Repairable Systems
1 Maintenance of Repairable Systems Bo Henry Lindqvist 1.1 Introduction A commonly used definition of a repairable system (Ascher and Feingold [3]) states that this is a system which, after failing to
More informationA hidden semi-markov model for the occurrences of water pipe bursts
A hidden semi-markov model for the occurrences of water pipe bursts T. Economou 1, T.C. Bailey 1 and Z. Kapelan 1 1 School of Engineering, Computer Science and Mathematics, University of Exeter, Harrison
More informationNon-Parametric Tests for Imperfect Repair
Non-Parametric Tests for Imperfect Repair GUSTAVO GILARDONI Departamento de Estatística, UnB Prédio CIC/EST, Campus Darcy Ribeiro Brasília - DF, 70910-900 gilardon@gmail.com ENRICO COLOSIMO Departamento
More informationI I FINAL, 01 Jun 8.4 to 31 May TITLE AND SUBTITLE 5 * _- N, '. ', -;
R AD-A237 850 E........ I N 11111IIIII U 1 1I!til II II... 1. AGENCY USE ONLY Leave 'VanK) I2. REPORT DATE 3 REPORT TYPE AND " - - I I FINAL, 01 Jun 8.4 to 31 May 88 4. TITLE AND SUBTITLE 5 * _- N, '.
More informationSystem Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models
System Simulation Part II: Mathematical and Statistical Models Chapter 5: Statistical Models Fatih Cavdur fatihcavdur@uludag.edu.tr March 20, 2012 Introduction Introduction The world of the model-builder
More informationRELIABILITY MODELING AND EVALUATION IN AGING POWER SYSTEMS. A Thesis HAG-KWEN KIM
RELIABILITY MODELING AND EVALUATION IN AGING POWER SYSTEMS A Thesis by HAG-KWEN KIM Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment of the requirements for the
More informationA Simulation Study on Confidence Interval Procedures of Some Mean Cumulative Function Estimators
Statistics Preprints Statistics -00 A Simulation Study on Confidence Interval Procedures of Some Mean Cumulative Function Estimators Jianying Zuo Iowa State University, jiyizu@iastate.edu William Q. Meeker
More informationReliability of Technical Systems
Main Topics 1. Introduction, Key Terms, Framing the Problem 2. Reliability Parameters: Failure Rate, Failure Probability, etc. 3. Some Important Reliability Distributions 4. Component Reliability 5. Software
More informationSolution: The process is a compound Poisson Process with E[N (t)] = λt/p by Wald's equation.
Solutions Stochastic Processes and Simulation II, May 18, 217 Problem 1: Poisson Processes Let {N(t), t } be a homogeneous Poisson Process on (, ) with rate λ. Let {S i, i = 1, 2, } be the points of the
More informationA Bivariate Weibull Regression Model
c Heldermann Verlag Economic Quality Control ISSN 0940-5151 Vol 20 (2005), No. 1, 1 A Bivariate Weibull Regression Model David D. Hanagal Abstract: In this paper, we propose a new bivariate Weibull regression
More informationChapter 5. Statistical Models in Simulations 5.1. Prof. Dr. Mesut Güneş Ch. 5 Statistical Models in Simulations
Chapter 5 Statistical Models in Simulations 5.1 Contents Basic Probability Theory Concepts Discrete Distributions Continuous Distributions Poisson Process Empirical Distributions Useful Statistical Models
More informationReliability analysis of power systems EI2452. Lifetime analysis 7 May 2015
Reliability analysis of power systems EI2452 Lifetime analysis 7 May 2015 Agenda Summary of content: Introduction nonreparable/reparable General information about statistical surveys Lifetime data Nonparametric
More informationModelling of Indian Stock Prices using Nonhomogeneous Poisson Processes with Time Trends
IOSR Journal of Business and Management (IOSR-JBM) e-issn: 2278-487X, p-issn: 2319-7668. Volume 13, Issue 4 (Sep. - Oct. 2013), PP 73-86 Modelling of Indian Stock Prices using Nonhomogeneous Poisson Processes
More informationPart III. Hypothesis Testing. III.1. Log-rank Test for Right-censored Failure Time Data
1 Part III. Hypothesis Testing III.1. Log-rank Test for Right-censored Failure Time Data Consider a survival study consisting of n independent subjects from p different populations with survival functions
More informationIntroduction to Reliability Theory (part 2)
Introduction to Reliability Theory (part 2) Frank Coolen UTOPIAE Training School II, Durham University 3 July 2018 (UTOPIAE) Introduction to Reliability Theory 1 / 21 Outline Statistical issues Software
More informationUNIT 5:Random number generation And Variation Generation
UNIT 5:Random number generation And Variation Generation RANDOM-NUMBER GENERATION Random numbers are a necessary basic ingredient in the simulation of almost all discrete systems. Most computer languages
More informationOne-Sample Numerical Data
One-Sample Numerical Data quantiles, boxplot, histogram, bootstrap confidence intervals, goodness-of-fit tests University of California, San Diego Instructor: Ery Arias-Castro http://math.ucsd.edu/~eariasca/teaching.html
More informationReliability Engineering I
Happiness is taking the reliability final exam. Reliability Engineering I ENM/MSC 565 Review for the Final Exam Vital Statistics What R&M concepts covered in the course When Monday April 29 from 4:30 6:00
More information(b) What is the variance of the time until the second customer arrives, starting empty, assuming that we measure time in minutes?
IEOR 3106: Introduction to Operations Research: Stochastic Models Fall 2006, Professor Whitt SOLUTIONS to Final Exam Chapters 4-7 and 10 in Ross, Tuesday, December 19, 4:10pm-7:00pm Open Book: but only
More informationComputer Science, Informatik 4 Communication and Distributed Systems. Simulation. Discrete-Event System Simulation. Dr.
Simulation Discrete-Event System Simulation Chapter 4 Statistical Models in Simulation Purpose & Overview The world the model-builder sees is probabilistic rather than deterministic. Some statistical model
More informationNotes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes.
Unit 2: Models, Censoring, and Likelihood for Failure-Time Data Notes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes. Ramón
More informationOther Survival Models. (1) Non-PH models. We briefly discussed the non-proportional hazards (non-ph) model
Other Survival Models (1) Non-PH models We briefly discussed the non-proportional hazards (non-ph) model λ(t Z) = λ 0 (t) exp{β(t) Z}, where β(t) can be estimated by: piecewise constants (recall how);
More informationInvestigation of goodness-of-fit test statistic distributions by random censored samples
d samples Investigation of goodness-of-fit test statistic distributions by random censored samples Novosibirsk State Technical University November 22, 2010 d samples Outline 1 Nonparametric goodness-of-fit
More informationCHAPTER 3 ANALYSIS OF RELIABILITY AND PROBABILITY MEASURES
27 CHAPTER 3 ANALYSIS OF RELIABILITY AND PROBABILITY MEASURES 3.1 INTRODUCTION The express purpose of this research is to assimilate reliability and its associated probabilistic variables into the Unit
More information10 Introduction to Reliability
0 Introduction to Reliability 10 Introduction to Reliability The following notes are based on Volume 6: How to Analyze Reliability Data, by Wayne Nelson (1993), ASQC Press. When considering the reliability
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models
Introduction to Empirical Processes and Semiparametric Inference Lecture 25: Semiparametric Models Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations
More informationConcepts and Tests for Trend in Recurrent Event Processes
JIRSS (2013) Vol. 12, No. 1, pp 35-69 Concepts and Tests for Trend in Recurrent Event Processes R. J. Cook, J. F. Lawless Department of Statistics and Actuarial Science, University of Waterloo, Ontario,
More informationThe exponential distribution and the Poisson process
The exponential distribution and the Poisson process 1-1 Exponential Distribution: Basic Facts PDF f(t) = { λe λt, t 0 0, t < 0 CDF Pr{T t) = 0 t λe λu du = 1 e λt (t 0) Mean E[T] = 1 λ Variance Var[T]
More informationCompeting risks for repairable systems: A data study
Competing risks for repairable systems: A data study Helge Langseth and Bo Henry Lindqvist Department of Mathematical Sciences, Norwegian University of Science and Technology, Norway Abstract We analyze
More informationTest for Discontinuities in Nonparametric Regression
Communications of the Korean Statistical Society Vol. 15, No. 5, 2008, pp. 709 717 Test for Discontinuities in Nonparametric Regression Dongryeon Park 1) Abstract The difference of two one-sided kernel
More information1 Delayed Renewal Processes: Exploiting Laplace Transforms
IEOR 6711: Stochastic Models I Professor Whitt, Tuesday, October 22, 213 Renewal Theory: Proof of Blackwell s theorem 1 Delayed Renewal Processes: Exploiting Laplace Transforms The proof of Blackwell s
More informationBootstrap tests. Patrick Breheny. October 11. Bootstrap vs. permutation tests Testing for equality of location
Bootstrap tests Patrick Breheny October 11 Patrick Breheny STA 621: Nonparametric Statistics 1/14 Introduction Conditioning on the observed data to obtain permutation tests is certainly an important idea
More informationLinear rank statistics
Linear rank statistics Comparison of two groups. Consider the failure time T ij of j-th subject in the i-th group for i = 1 or ; the first group is often called control, and the second treatment. Let n
More informationDESIGN OF PREVENTIVE MAINTENANCE SCHEDULING MODEL FOR DETERIORATING SYSTEMS
DESIGN OF PREVENTIVE MAINTENANCE SCHEDULING MODEL FOR DETERIORATING SYSTEMS P.A. Ozor a, S.O. Onyegegbu Department of Mechanical Engineering, University of Nigeria, Nsukka, NIGERIA. a Email: paul.ozor@unn.edu.ng
More informationCHAPTER 3 MATHEMATICAL AND SIMULATION TOOLS FOR MANET ANALYSIS
44 CHAPTER 3 MATHEMATICAL AND SIMULATION TOOLS FOR MANET ANALYSIS 3.1 INTRODUCTION MANET analysis is a multidimensional affair. Many tools of mathematics are used in the analysis. Among them, the prime
More informationLecturer: Olga Galinina
Renewal models Lecturer: Olga Galinina E-mail: olga.galinina@tut.fi Outline Reminder. Exponential models definition of renewal processes exponential interval distribution Erlang distribution hyperexponential
More informationPractice Problems Section Problems
Practice Problems Section 4-4-3 4-4 4-5 4-6 4-7 4-8 4-10 Supplemental Problems 4-1 to 4-9 4-13, 14, 15, 17, 19, 0 4-3, 34, 36, 38 4-47, 49, 5, 54, 55 4-59, 60, 63 4-66, 68, 69, 70, 74 4-79, 81, 84 4-85,
More informationResampling-Based Control of the FDR
Resampling-Based Control of the FDR Joseph P. Romano 1 Azeem S. Shaikh 2 and Michael Wolf 3 1 Departments of Economics and Statistics Stanford University 2 Department of Economics University of Chicago
More informationISyE 6644 Fall 2014 Test #2 Solutions (revised 11/7/16)
1 NAME ISyE 6644 Fall 2014 Test #2 Solutions (revised 11/7/16) This test is 85 minutes. You are allowed two cheat sheets. Good luck! 1. Some short answer questions to get things going. (a) Consider the
More informationSlides 8: Statistical Models in Simulation
Slides 8: Statistical Models in Simulation Purpose and Overview The world the model-builder sees is probabilistic rather than deterministic: Some statistical model might well describe the variations. An
More informationStochastic process. X, a series of random variables indexed by t
Stochastic process X, a series of random variables indexed by t X={X(t), t 0} is a continuous time stochastic process X={X(t), t=0,1, } is a discrete time stochastic process X(t) is the state at time t,
More informationOnline Supplement to Are Call Center and Hospital Arrivals Well Modeled by Nonhomogeneous Poisson Processes?
Online Supplement to Are Call Center and Hospital Arrivals Well Modeled by Nonhomogeneous Poisson Processes? Song-Hee Kim and Ward Whitt Industrial Engineering and Operations Research Columbia University
More informationLecture 3. Truncation, length-bias and prevalence sampling
Lecture 3. Truncation, length-bias and prevalence sampling 3.1 Prevalent sampling Statistical techniques for truncated data have been integrated into survival analysis in last two decades. Truncation in
More informationChap 4. Software Reliability
Chap 4. Software Reliability 4.2 Reliability Growth 1. Introduction 2. Reliability Growth Models 3. The Basic Execution Model 4. Calendar Time Computation 5. Reliability Demonstration Testing 1. Introduction
More informationMultiple Testing of One-Sided Hypotheses: Combining Bonferroni and the Bootstrap
University of Zurich Department of Economics Working Paper Series ISSN 1664-7041 (print) ISSN 1664-705X (online) Working Paper No. 254 Multiple Testing of One-Sided Hypotheses: Combining Bonferroni and
More informationOptional Stopping Theorem Let X be a martingale and T be a stopping time such
Plan Counting, Renewal, and Point Processes 0. Finish FDR Example 1. The Basic Renewal Process 2. The Poisson Process Revisited 3. Variants and Extensions 4. Point Processes Reading: G&S: 7.1 7.3, 7.10
More informationExercises. (a) Prove that m(t) =
Exercises 1. Lack of memory. Verify that the exponential distribution has the lack of memory property, that is, if T is exponentially distributed with parameter λ > then so is T t given that T > t for
More informationORDER RESTRICTED STATISTICAL INFERENCE ON LORENZ CURVES OF PARETO DISTRIBUTIONS. Myongsik Oh. 1. Introduction
J. Appl. Math & Computing Vol. 13(2003), No. 1-2, pp. 457-470 ORDER RESTRICTED STATISTICAL INFERENCE ON LORENZ CURVES OF PARETO DISTRIBUTIONS Myongsik Oh Abstract. The comparison of two or more Lorenz
More informationConfidence Intervals for Reliability Growth Models with Small Sample Sizes. Reliability growth models, Order statistics, Confidence intervals
Confidence Intervals for Reliability Growth Models with Small Sample Sizes John Quigley, Lesley Walls University of Strathclyde, Glasgow, Scotland Key Words Reliability growth models, Order statistics,
More informationAnalyzing Municipal Blockage Failure Datasets for Sewer Systems Yongliang Jin, 1 and Amlan Mukherjee 2
Analyzing Municipal Blockage Failure Datasets for Sewer Systems Yongliang Jin, 1 and Amlan Mukherjee 2 1 Graduate Research Assistant, Department of Civil and Environmental Engineering, Michigan Tech, Houghton,
More informationMath 181B Homework 1 Solution
Math 181B Homework 1 Solution 1. Write down the likelihood: L(λ = n λ X i e λ X i! (a One-sided test: H 0 : λ = 1 vs H 1 : λ = 0.1 The likelihood ratio: where LR = L(1 L(0.1 = 1 X i e n 1 = λ n X i e nλ
More informationSemiparametric Regression
Semiparametric Regression Patrick Breheny October 22 Patrick Breheny Survival Data Analysis (BIOS 7210) 1/23 Introduction Over the past few weeks, we ve introduced a variety of regression models under
More informationExam C Solutions Spring 2005
Exam C Solutions Spring 005 Question # The CDF is F( x) = 4 ( + x) Observation (x) F(x) compare to: Maximum difference 0. 0.58 0, 0. 0.58 0.7 0.880 0., 0.4 0.680 0.9 0.93 0.4, 0.6 0.53. 0.949 0.6, 0.8
More informationBIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY
BIAS OF MAXIMUM-LIKELIHOOD ESTIMATES IN LOGISTIC AND COX REGRESSION MODELS: A COMPARATIVE SIMULATION STUDY Ingo Langner 1, Ralf Bender 2, Rebecca Lenz-Tönjes 1, Helmut Küchenhoff 2, Maria Blettner 2 1
More informationSPRING 2007 EXAM C SOLUTIONS
SPRING 007 EXAM C SOLUTIONS Question #1 The data are already shifted (have had the policy limit and the deductible of 50 applied). The two 350 payments are censored. Thus the likelihood function is L =
More informationModeling of Dependence Between Critical Failure and Preventive Maintenance: The Repair Alert Model
Modeling of Dependence Between Critical Failure and Preventive Maintenance: The Repair Alert Model Bo Lindqvist Helge Langseth Bård Støve bo@math.ntnu.no Department of Mathematical Sciences Norwegian University
More informationPoisson Processes. Stochastic Processes. Feb UC3M
Poisson Processes Stochastic Processes UC3M Feb. 2012 Exponential random variables A random variable T has exponential distribution with rate λ > 0 if its probability density function can been written
More informationChapter 9 Part II Maintainability
Chapter 9 Part II Maintainability 9.4 System Repair Time 9.5 Reliability Under Preventive Maintenance 9.6 State-Dependent Systems with Repair C. Ebeling, Intro to Reliability & Maintainability Chapter
More informationDistribution-Free Procedures (Devore Chapter Fifteen)
Distribution-Free Procedures (Devore Chapter Fifteen) MATH-5-01: Probability and Statistics II Spring 018 Contents 1 Nonparametric Hypothesis Tests 1 1.1 The Wilcoxon Rank Sum Test........... 1 1. Normal
More informationAsymptotic Statistics-III. Changliang Zou
Asymptotic Statistics-III Changliang Zou The multivariate central limit theorem Theorem (Multivariate CLT for iid case) Let X i be iid random p-vectors with mean µ and and covariance matrix Σ. Then n (
More informationSummary: the confidence interval for the mean (σ 2 known) with gaussian assumption
Summary: the confidence interval for the mean (σ known) with gaussian assumption on X Let X be a Gaussian r.v. with mean µ and variance σ. If X 1, X,..., X n is a random sample drawn from X then the confidence
More informationPackage homtest. February 20, 2015
Version 1.0-5 Date 2009-03-26 Package homtest February 20, 2015 Title Homogeneity tests for Regional Frequency Analysis Author Alberto Viglione Maintainer Alberto Viglione
More informationA conceptual interpretation of the renewal theorem with applications
Risk, Reliability and Societal Safety Aven & Vinnem (eds) 2007 Taylor & Francis Group, London, ISBN 978-0-415-44786-7 A conceptual interpretation of the renewal theorem with applications J.A.M. van der
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationLecture 7. Poisson and lifetime processes in risk analysis
Lecture 7. Poisson and lifetime processes in risk analysis Jesper Rydén Department of Mathematics, Uppsala University jesper.ryden@math.uu.se Statistical Risk Analysis Spring 2014 Example: Life times of
More informationLet (Ω, F) be a measureable space. A filtration in discrete time is a sequence of. F s F t
2.2 Filtrations Let (Ω, F) be a measureable space. A filtration in discrete time is a sequence of σ algebras {F t } such that F t F and F t F t+1 for all t = 0, 1,.... In continuous time, the second condition
More informationSurvival Analysis. Lu Tian and Richard Olshen Stanford University
1 Survival Analysis Lu Tian and Richard Olshen Stanford University 2 Survival Time/ Failure Time/Event Time We will introduce various statistical methods for analyzing survival outcomes What is the survival
More informationBootstrapping Heteroskedasticity Consistent Covariance Matrix Estimator
Bootstrapping Heteroskedasticity Consistent Covariance Matrix Estimator by Emmanuel Flachaire Eurequa, University Paris I Panthéon-Sorbonne December 2001 Abstract Recent results of Cribari-Neto and Zarkos
More informatione 4β e 4β + e β ˆβ =0.765
SIMPLE EXAMPLE COX-REGRESSION i Y i x i δ i 1 5 12 0 2 10 10 1 3 40 3 0 4 80 5 0 5 120 3 1 6 400 4 1 7 600 1 0 Model: z(t x) =z 0 (t) exp{βx} Partial likelihood: L(β) = e 10β e 10β + e 3β + e 5β + e 3β
More information3 Modeling Process Quality
3 Modeling Process Quality 3.1 Introduction Section 3.1 contains basic numerical and graphical methods. familiar with these methods. It is assumed the student is Goal: Review several discrete and continuous
More informationHeteroskedasticity-Robust Inference in Finite Samples
Heteroskedasticity-Robust Inference in Finite Samples Jerry Hausman and Christopher Palmer Massachusetts Institute of Technology December 011 Abstract Since the advent of heteroskedasticity-robust standard
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationExact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring
Exact Inference for the Two-Parameter Exponential Distribution Under Type-II Hybrid Censoring A. Ganguly, S. Mitra, D. Samanta, D. Kundu,2 Abstract Epstein [9] introduced the Type-I hybrid censoring scheme
More informationApplication of Homogeneity Tests: Problems and Solution
Application of Homogeneity Tests: Problems and Solution Boris Yu. Lemeshko (B), Irina V. Veretelnikova, Stanislav B. Lemeshko, and Alena Yu. Novikova Novosibirsk State Technical University, Novosibirsk,
More information