Reliability Analysis of Repairable Systems with Dependent Component Failures. under Partially Perfect Repair

Size: px
Start display at page:

Download "Reliability Analysis of Repairable Systems with Dependent Component Failures. under Partially Perfect Repair"

Transcription

1 Reliability Analysis of Repairable Systems with Dependent Component Failures under Partially Perfect Repair Qingyu Yang & Nailong Zhang Department of Industrial and System Engineering Wayne State University Detroit, MI 48202, USA Yili Hong Department of Statistics Virginia Tech Blacksburg, VA 2406, USA Abstract: Existing reliability models for repairable systems with a single component can be well applied for a range of repair actions from perfect repair to minimal repair. Establishing reliability models for multi-component repairable systems, however, is still a challenge problem when considering the dependency of component failures. This paper focuses on a special repair assumption, called partially perfect repair, for repairable systems with dependent component failures, where only the failed component is repaired to as good as new condition. A parametric reliability model is proposed to capture the statistical dependency among different component failures, in which the oint distribution of the latent component failure time is established using copula functions. The model parameters are estimated by using the maximum likelihood method, and the maximum likelihood function is calculated based on the conditional probability. Based on the proposed reliability model, statistical hypothesis testing procedures are developed to determine the dependency of component failures. The developed methods are illustrated with an application in a cylinder head assembling cell that consists of multiple stations.

2 Keywords: Competing risks, recurrent events, Archimedean copula, Gaussian copula, lognormal, Weibull. ACRONYMS ML Maximum likelihood cdf Cumulative distribution function pdf AIC NLV Probability density function Akaike information criterion Negative log-likelihood value NOTATIONS N( t ) The number of system failures of all components occurring in the time interval (0, t ] r ( t ) The most recently observed failure time of component occurring before time t K τ The number of failure types The predetermined ending time for failure observation (i.e., last follow up time) 2

3 T i Random variables to denote the failure time of the th i failure t i Observed failure time of the th i failure ( ) a t, the age of component at time t since its last repair D The random vector of the latent ages to failure of all components D = [,,, ] T, where D k is the latent age to failure of the i th D D2 D K component µσ, µ = [ µ,, µ ] T K, and σ = [ σ,, σ ] T K, denoting the location, and scale parameters in the Weibull marginals, respectively M A lower triangular matrix with positive diagonal elements ω A matrix of angles with dimension K ( K ) C u u K (,, ) Copula function with K marginal distributions λ Lθ ( ) ˆθ Σ The association parameter in the Gumbal-Hougaard copula Likelihood function for unknown parameter vector θ ML estimator of θ Covariance matrix of multivariate lognormal distribution and Gaussian copula Λ Correlation matrix of multivariate lognormal distribution and Gaussian copula 3

4 . Introduction A repairable system is defined as a system that can be restored to fully satisfactory operation by parts replacements or adustments when failing to perform one or more of its functions satisfactorily []. Traditional study on repairable systems mainly focuses on reliability models for single-component repairable systems under different repair actions. A comprehensive review on statistical methods of single-component repairable systems is provided by Lindqvist [2]. In recent years, more attention has been put on the reliability modeling for multi-component systems under competing risks. The competing risks model assumes that each component s failure will cause the whole system s failure, thus a repair action should be taken whenever a failure occurs regardless of the failure type. A real example of a repairable competing risks system can be found in [3]. For reliability models under competing risks, most of the existing research assumes statistical independency of component failure [3]-[2]. Thus, the reliability analysis of the entire system subect to competing risks can be simplified by analyzing each component individually. When considering component failure dependency, most existing reliability models [3]-[20] assumes that, when a failure of one component occurs, it will result in a possible shock to the other components with a certain probability. In this paper, we establish reliability modeling of repairable complex systems subect to statistically dependent competing 4

5 risks due to degradation. Pan and Balakrishnan [2] studied a reliability model where the bivariate degradation of two statistically dependent performance characteristics are modeled by a gamma process. Other research [22]-[24] discussed the dependency between failures due to random shocks and degradation processes. However, these methods cannot deal with the dependency of failures due to different degradation processes. To deal with the dependency of competing risks in complex systems, copula functions have been applied in the literature [25]-[28]. These works, however, are only capable of dealing with competing systems that are either non-repairable or perfectly repaired. In this paper, we study repairable statistically dependent competing risks systems under a special repair action called partially perfect repair, i.e., only the failed component, rather than the whole system, is replaced whenever a failure happens [5], [7]. The partially perfect repair action is general for many complex multi-component engineering systems in real applications. For a simple example, if the bumper of a vehicle has a failure, only the failed bumper will be repaired, while the other parts of the vehicle will still be used. On the other hand, under most situations, it is more economical to replace the failed component (e.g., bumper), rather than repair it, due to high labor cost. For repairable competing risks systems under a perfect repair assumption, after each repair action, the system will have the same failure mechanism as the original one. In 5

6 contrast, under partially perfect repair actions, the replacement of a failed component affects the failure mechanism of all the other components that will be continually used in the system. Thus, the methods developed in the literature cannot be directly applied. In this research, a parametric statistical model is established to capture the statistical dependency of repairable multiple component failures under partially perfect repair assumptions. In the proposed model, the oint distribution of components latent ages to failure is established using copula functions. In particular, two types of copulas, i.e., Archimedean copula [29], and Gaussian copula [30], are applied to study the overall statistical dependency and pairwise statistical dependencies among different components, respectively. The framework of the paper is as follows. After the introduction, Section 2 proposes a general statistical model for multiple statistically dependent competing risks. When different marginal distributions of latent ages to failure are selected, different parametric forms of the proposed model are investigated in Sections 3. In Section 4, the model parameters are estimated by using the maximum likelihood (ML) method, and the ML function is calculated based on the conditional probability. Hypothesis testing procedures are developed in Section 5 to determine the failure dependency. Section 6 applies the proposed methods to an application of an automobile cylinder head assembly cell. Finally, Section 7 gives the summary. 6

7 2. Statistical multiple dependent competing risks model We consider a single repairable multi-component system, i.e., the failure data are collected from only one realization of the repairable system. The research results can be easily extended to multiple systems. Suppose the system consists of K different components or subsystems corresponding to K different types of failures. The system is under competing risk, and the repair action is assumed to be a partially perfect repair. In addition, it is assumed that two or more failures do not occur simultaneously. The observed failure events include failure time t, t 2,..., and associated failure types, which occur before a pre-determined time τ. The counting process, denoted by N( t ), records the number of failures regardless of failure types occurring in the time interval (0, t ]. Consider a system that consists of K new components starting to work at time 0. Because the components are under competing risks, the system fails if any component fails, while the failure time of all the other components cannot be observed. The latent age to failure of each component is denoted in this paper as a random variable D ; =, 2,..., K, with the oint survival function S(d,, d K ) = Pr(D > d,, D K > d K ). Let r ( t ) be the most recent failure time of component before time t. The running time of component at time t since its last replacement, which is defined as age and is denoted as a ( t ), can be calculated as a ( t) = t r ( t). Note that a ( t ) is a left-continuous function. The next failure time of the system on or after time t is 7

8 determined as the minimum value of { r ( t) + D, =,2,..., K } under the condition that D > a ( t), =,2,..., K. Fig. illustrates the first two failures in a competing risks T system. Let vector d,,..., d, K be the first realization of the latent ages to failures D. The first failure is determined by the minimal value of d ;, i i =,..., K. Suppose the first failure is due to component i, and occurs at time t (Fig. left). Under partially perfect repair assumptions, only component i is replaced, and the most recent failure times are updated as ri ( t) = t, and r ( t ) = 0, i, =,..., K. The second failure can be calculated as the minimum value of { d2, i + t, d2, ; i, =,..., K} (Fig. right), where T d2,,..., d 2, K is the second realization of D. d, d 2, r ( t ) = 0 r ( t 2) = 0 t 2 d,i d 2,i r ( t ) = 0 i t r ( t ) = t i 2 d,k d 2,K r ( ) 0 K t = r3 ( t 2) = 0 Fig.. Illustrations of the first (left), and the second (right) failures in a competing risks system. 3. Parametric forms for the proposed model The oint distribution of the random vector D describes the statistical failure mechanism of multiple components, and thus captures their statistical failure dependency. In this section, parametric models are proposed to characterize the oint distribution of D. 8

9 Lognormal and Weibull distributions are from the log-location-scale family, and are commonly-used as failure-time distributions. In this paper, they are chosen as the marginal distribution of random vector D. It is well known that, when the oint distribution of a random vector is multivariate lognormal, the marginal distributions are also lognormal. As for a multivariate Weibull distribution, we use a copula function to link the oint distributions to the marginal distributions. 3. Parametric forms for multivariate lognormal distribution When the oint distribution of random vector D is multivariate lognormal, the oint probability density function (pdf) is calculated as T p( d,, d K ; µσ, ) = exp /2 [ log( ) ] [ log( ) ] K /2 d µ Σ d µ () (2 π ) Σ 2 The model parameters θ include µ and Σ. K µ R, and K K Σ R are the mean, and covariance matrix of the multivariate lognormal, respectively. 3.2 Parametric forms for a multivariate Weibull distribution via Archimedean copula The Archimedean family of copulas are frequently used for the construction of multivariate distributions due to their simple forms [29] C( u,, uk ) = ψ ψ ( u) + + ψ ( uk ) (2) 9

10 where ψ is the generator of the Archimedean copula. Different generators will generate different Archimedean copulas. For example, ψ ( t) / λ = ( + t), and ψ = / λ ( t) exp( t ) are generators for Clayton, and Gumbel-Hougaard copulas, respectively. In this paper, the Gumbel-Hougaard copula is selected as an example to illustrate the application of the Archimedean copula family in the proposed reliability model. The K-dimensional Gumbel-Hougaard copula has the form { } λ λ C( u,, uk ) exp [ log( u )] [ log( uk )] λ = + + (3) where λ is called the association parameter. To construct the multivariate Weibull distribution, we choose the log-location-scale form [3] of the Weibull distribution as marginal distributions of D, i.e., log( d ) µ F ( d; µ, σ ) = Φ SEV σ (4) where Φ ( x SEV ) = exp[ exp( x )] is the cumulative distribution function (cdf) of the standard smallest extreme value distribution. Let F( d; µ, σ ) in (4) be u ; =,..., K, and substitute (4) into (3); then the oint cdf of random vector D following a multivariate Weibull distribution is obtained as λ λ K log( d ) µ F( d,, dk ; µ,σ, λ) = exp log Φ SEV (5) = σ 0

11 In this multivariate Weibull distribution constructed via the Gumbel-Hougaard copula, the model parameters θ = ( µ,σ, λ). 3.3 Parametric forms for multivariate Weibull distribution via the Gaussian copula The Gaussian copula is a special copula taking advantage of the pdf of the multivariate normal distribution. Specifically, a Gaussian copula has the form Gauss C ( u,,, u K ) = Φ Σ[ Φ ( u),, Φ ( u K )] (6) where Φ is the inverse of the cdf of the standard normal distribution, and Φ Σ is the cdf of a multivariate normal distribution with zero mean vector, and its covariance matrix equals its correlation matrix. The Gaussian copula density function is given as [32] T Φ ( u) Φ ( u ) Gauss c ( u,, uk ) = exp ( ) /2 2 Σ I (7) Σ Φ ( uk ) Φ ( uk ) where Σ is the correlation matrix, and I is the identity matrix. When applying a Gaussian copula to construct the oint Weibull distribution, the survival function of random vector D is obtained as = Gauss K µσσ K K d dk S( d,, d ;,, ) f ( x,, x ) dx dx (8) Gauss where f ( x,, x K ) denotes the pdf of the oint Weibull distribution obtained by the chain rule, i.e.,

12 K Gauss Gauss C du du f ( x,, xk ) = u u dx dx K Gauss = c ( u, u ) f ( x ) f ( x ) K K K K, (9) and f ( ) denotes the pdf of the marginal. In the multivariate Weibull distribution constructed via the Gaussian copula, the model parameters θ = ( µ,σσ, ). 4. Parameters estimation based on ML approach An ML method is developed in this paper to estimate model parameters. To implement the ML approach, we first calculate the likelihood function. Suppose the th i failure is due to component k, and occurs at time t i. The latent age to failure of the failed component k is equal to a ( t ), while the latent ages of all other components k i should be larger than a ( t ); k, =,..., K. Thus, the unconditional probability to i observe failure i is calculated as Pr( D = a ( t ), D > a ( t ); k, =,..., K) k k i i S[ a ( t ),, D,, a ( t )] (0) = i k K i D k Dk = ak ( ti ) Note that (0) gives the probability at a given time t for a continuous random variable. Although technically Pr( X = t) = 0 for a continuous random variable X with pdf f ( t ), Pr( X = t) can be interpreted as Pr( t X t + dt) = f ( t) dt, which is proportional to f ( t ). We ignore dt only for notational convenience in the calculation of the likelihood function for all continuous random variables in the rest of the paper. 2

13 Equation (0) only accounts for the probability of an observed failure at time t i regardless of previous failure data. The conditional probability of observing failure i given all the previous i- failures is solely determined by the ages of all component after the repair action of the (i-) th failure. Specifically, the likelihood of failure i, i =,2,, N( τ), conditioned on all the previous i- failures, can be calculated as L i [ ( ),,,, ( )] S a ti Dk ak ti D = S[ a ( t ),, a ( t )] k Dk = ak ( ti ) + + i K i () where, S( ) is the oint survival function of the latent ages to failures. For example, consider the first two failures illustrated in Fig.. After the repair action of the first failure at time t, the age of the failed component i is updated to a ( t + + ) = 0, while the ages of all the other components are given as a ( t ) = t, i i, =,..., K. As the second failure is due to component one, and occurs at time t 2, the likelihood of the second failure conditioned on the first failure can be calculated as Pr( D = a ( t ), D > a ( t );, =,..., K D > 0, D > t ; k i, k =,..., K ), where 2 2 i k ai ( t2) = t2 t, and ak ( t2 ) = t2; k i, k =,..., K. As there is no failure observed from t to the predetermined ending time τ, the N ( τ ) likelihood LN ( τ ) + can be calculated as L Pr{ D > [ τ r ( τ )],..., D > [ τ rk ( τ )]} = (2) [ ( ),..., ( )] K N ( τ ) S a tn ( τ ) ak tn ( τ ) 3

14 S[ a ( t ),..., a ( t )] =. S[ a ( t ),..., a ( t )] τ K τ + + N ( τ ) K N ( τ ) Combining the results in () and (2), the following Proposition can be used to calculate the likelihood function based on the observed failure data. Proposition. Given the observed failure data, the likelihood function can be calculated by N ( τ ) + L( θ ) = Li (3) i= where L can be calculated based on (), and (2) for i =, 2,, N( τ ), and i i = N( τ ) +, respectively. The estimated model parameters ˆθ are obtained by maximizing Lθ ( ). Based on the ML theory [33], [34], the estimated parameters ˆθ are asymptotically normally distributed with a large sample assumption. In the following sections, we will discuss how to apply the developed ML methods for different parametric models. 4. ML method applied for multivariate lognormal distribution When the oint distribution of D is selected as multivariate lognormal, the oint survival function S( ) can be calculated by integrating (), i.e., = d dk S ( dµ ;, Σ) p ( xµ ;, Σ) dx dx. (4) K 4

15 The likelihood function Lθ ( ) can be obtained by substituting (4) into () and (2). To calculate L i, the first order partial derivative of the survival function in (4) needs to be evaluated. We extend the first order partial derivative of the survival function as S[ a ( ti ),, Dk,, ak ( ti )] D k D = a ( t ) k k i = f ( a ( t )) Pr( D > a ( t ); k D = a ( t )) k k i i k k i (5) In (5), fk ( ak ( t i )) denotes the marginal lognormal probability density of the latent age to failure of component k. The detailed calculation of the conditional probability Pr( D > a ( t ); k D = a ( t )) is given in the Appendix. i k k i A challenge to implement the parameter estimation of the multivariate lognormal comes from the positive definite constraint of the covariance matrix Σ. We use Cholesky decomposition to achieve this constraint. Specifically, if M is a lower triangular matrix with positive diagonal elements, Σ = M T M is positive definite. Therefore, we directly estimate the lower triangular matrix ˆM, from which the covariance estimator ˆΣ can be calculated. 4.2 ML method applied for the multivariate Weibull distribution via the Archimedean copula When the oint Weibull function is constructed via a Gumbel-Hougaard copula, we need to substitute the survival function of random vector D into () and (2) to implement the ML approach. The cdf of random vector D is given in (5). Based on (5), 5

16 the survival function can be calculated using (25), and the details to obtain (25) are given in the Appendix. Thus, the likelihood function Lθ ( ) is calculated by substituting the survival function (25) into () and (2), where θ = ( µσ,, λ). A scale-invariant variable, called Kendall s tau, and calculated as τ = 4 E( C( U, V )), is generally used to measure the dependency in the bivariate C copula. For the Gumbal-Hougaard copula, Kendall s tau can be calculated based on association parameter λ [29], i.e., τ = λ. (6) Kendall s tau can only capture the overall statistical dependency of component failures. To capture the pairwise dependencies of component failures, the Gaussian copula function is applied to derive the oint distribution of random vector D, which will be discussed in the next section. 4.3 ML method applied to a multivariate Weibull distribution via the Gaussian copula Similar to the parameters estimation for the multivariate lognormal distribution, the positive definite constraint is also required for the multivariate normal distribution. The Cholesky decomposition, however, fails to construct the positive definite correlation matrix due to the additional constraint that requires a correlation matrix equal to the covariance matrix. To tackle this problem, we use a hypersphere decomposition proposed 6

17 by Rebonato [35], which provides an approach to construct the valid correlation matrix satisfying all the constraints. Specifically, based on an arbitrary matrix of angles ω of dimension K ( K ), matrix B of dimension K K is constructed by the following rules. Bi, = cosωi, sin ωi, k (if =,, K ) k= Bi, = sin ωi, k (if = K) k= (7) where the subscripts i and are used to represent the element s position in matrix B. Then T Σ = B B is a valid correlation matrix satisfying both constraints. Therefore, the likelihood function is a function of θ = ( µσω,, ). In addition, the correlation matrix is an immediate measure of failure dependency. 5. Hypothesis testing for dependency Based on the proposed reliability model for the competing risks systems, statistical hypothesis testing procedures are developed in this section to determine the component failure dependencies. In Section 5., a dependency test based on the multivariate lognormal distribution is proposed. Then in Section 5.2, a dependency test for the multivariate Weibull distribution derived by the Archimedean copula and the Gaussian copula are discussed, respectively. 7

18 5. A dependency test for the multivariate lognormal distribution In the multivariate lognormal distribution, the dependency information is captured by the correlation matrix Λ, which can be calculated based on covariance matrix Σ : Λ i, Σ i, = i = Σ Σ i, i, ;,,, K (8) where Λ, and i, Σ are the elements of the i, th i row, and the th column in Λ, and Σ, respectively. In addition, T Σ = M M, and M is the lower triangular matrix. Based on the correlation matrix Λ, a hypothesis testing procedure is developed to determine the statistical dependency among latent ages to failures of components. H 0 : component failures i, are statistically independent. H : component failures i, are statistically dependent. (9) When the asymptotic result is applied, we use a normal approximation to construct the test statistics. The test statistics W are used to test the dependency between i, component failures i and, which equals the estimate of correlation Λ ˆ i, divided by its estimated standard error var( Λ ˆ i, ), i.e., W =Λˆ / var( Λ ˆ ) (20) i, i, i, where the asymptotic variance var( Λ ˆ i, ) is calculated using the delta method, and the detailed calculation is given in the Appendix. In hypothesis testing (9), H 0 is reected if Wi, Z α / 2 >, or Wi, <, where α is the test significance level, and Z / 2 α /2 is the Z α 8

19 upper quantile of the standard normal distribution. Based on hypothesis testing (9), pairwise statistical dependencies between different component failures can be tested. 5.2 A dependency test for multivariate Weibull distributions In this section, a statistical dependency test for oint distribution D constructed via the Archimedean copula is first discussed. The Archimedean copula can capture the overall statistical dependency, which is determined by the global association parameter λ. We developed the following hypothesis testing procedure to test the overall failure dependency, i.e., to see whether Kendall s tau is equal to 0. H 0 : all failure types are statistically independent. H : not all failure types are statistically independent. (2) Here, the asymptotic test statistic W overall is constructed using the estimate of Kendall s tau ˆ τ divided by its estimated standard error var( ˆ τ ). H 0 is reected if Woverall > or Woverall < /2. As Kendall s tau is a function of the association Z α /2 Z α parameter λ, the variance of the estimate of Kendall s tau ˆ τ can be calculated by using the delta method. dτ T var( ˆ) ( ) var( ˆ dτ τ = λ)( ) dλ dλ var( ˆ λ) = ˆ4 λ λ = ˆ λ (22) 9

20 where var( ˆ λ ) denotes the asymptotic variance of ˆλ, which is calculated in the same way as the asymptotic covariance matrix of ˆM in the Appendix. Based on (22), the asymptotic test statistic of the overall dependency is given as W overall ˆ τ ˆ2 λ =. (23) var( ˆ λ) For the multivariate Weibull distribution constructed via the Gaussian copula, the correlation matrix determines the pairwise dependencies of latent ages to component failures. Thus, the test procedure is the same as that discussed for the multivariate lognormal distribution in Section 5.. The only difference here is that the correlation matrix Λ is a function of the angle matrix ω for hypersphere decomposition, rather than a function of lower triangular matrix M in Cholesky decomposition. Specifically, the asymptotic variance of the correlation coefficient Λ ˆ i, in (25) is calculated as ˆ Λ Λ var( Λ ) = ( ) Σˆ ( ) ω ω i, T i, i, ωˆ ω= ωˆ. (24) The test statistics for Multivariate Weibull via the Gaussian copula has the same form as (20). Thus, hypothesis testing (9) can be used here to test the pairwise statistical dependencies among different failure types. 6. An application in a cylinder head production process We consider an application from an automobile power-train plant in the United States. The plant has a manufacturing cell that consists of four stations, which are denoted by 20

21 stations, 2, 3, and 4, to protect proprietary sensitivity information. A cylinder head is assembled by going through all stations sequentially, and all stations move synchronized with each other. The statistical failure dependency among these 4 stations is not well understood, thus we will investigate the dependency based on the proposed methodology. Fig. shows the cumulative number of station failures collected from the cylinder head manufacturing process for about one and a half years. 70 Cumulative number of failures Station Station 2 Station 3 Station Time/Hours Fig. 2. Failure data from stations, 2, 3, and 4. Given the failure time data of the manufacturing cell, the statistical reliability model is established when the marginal distributions are selected as lognormal and Weibull, 2

22 respectively. The ML method is used to estimate the model paremeters. All the compuation is done in the R environment. For the numerical optimization methods, we applied both Nelder-Mead and simulated annealing methods. The results obtained from these two methods are very close. Based on the developed reliability model, the hypothesies testing methods can be used to investigate the dependency among stations, 2, 3, and A reliability model with a multivariate lognormal distribution When the oint distribution of D is selected as a multivariate lognormal distribution, the statistical reliability model is established, and the model parameters are estimated by maximizing the likelihood function in (3). The parameter estimates, θ ˆ = ( µ ˆ, M ˆ ), are [ ] µ ˆ = ˆ M = The estimate of the covariance matrix, calculated as Σˆ = Mˆ M ˆ T, is given as ˆ Σ = To test the statistical dependency of component failures, the ML estimate of covariance coefficient Λ is calculated using (8), and the asymptotic variance of i, Λ, T ˆ i 22

23 is obtained based on (26) in the Appendix. The test statistics W obtained by (20) are i, used to calculate the p-values. Table I lists calculated p-values and test results of pairwise statistical dependencies with a 95% confidence level, where indicates statistical dependency, and 0 indicates statistical independency. From Table I, it can be seen that stations (, 2) and stations (, 4) fail statistically dependently. Table I p-values, and statistical dependency tests for the multivariate lognormal distribution Station A dependency test for the multivariate Weibull distriubtions When the Gumbel-Hougaard copula is applied to derive the oint Weibull distribution, the ML estimates of the parameters in the marginal Weibull distribution [ µσ ˆ, ˆ ] T is calculated as T [ ] T [ ] µ ˆ = σˆ = The estimate of association parameter ˆ λ = Based on (23), Woverall is calculated as The p-value is 0.532, which indicates that the null hypothesis is not reected. Thus, the dependency test result shows that all failure types are statistically independent. 23

24 When applying the Gaussian copula to construct the oint distribution, the model parameters θ = ( µσω,, ), including the parameters from both Weibull marginals, and the multivariate normal. The ML estimates are obtained as [ ] [ ] µ ˆ = σˆ = ωˆ = Based on (7) of the hypersphere decomposition, the estimate of ˆB is calculated based on ˆω as ˆ B = Then the correlation matrix estimate ˆΛ in the Gaussian copula is calculated by Λˆ = Bˆ B ˆ T as ˆ Λ = To test the statistical dependency of component failures, Table II lists the p-values and test results of pairwise statistical dependency for the Gaussian copula based on (24), and (20). Table II T T 24

25 p-values, and statistical dependency test when using a Gaussian copula to construct the oint Weibull distribution. Station Compared with the results from the multivariate lognormal distribution, the statistical dependency result of stations (, 2) is different, while the dependency results of the other 5 pairs are consistent. The difference is due to the different multivariate distributions applied in the statistical reliability model. After the hypothesis testing, if a group of components fails statistically independently from other components, the corresponding observed failure data can be used to estimate their statistical failure distributions individually. The statistical reliability model can be simplified, and the statistical reliability model can be fitted with better precision. Thus, we fit the reliability model for a multivariate lognormal distribution with the restirction that station 3 is statistically independent from all the others. Similarly, the reliability model for the multivariate Weibull distribution is also fitted with the restirction that both station 2 and station 3 are statistically independent from all the other stations. In the literature, the likelihood values, and the Akaike information criterion (AIC) values that take the number of parameters into account are generally used to compare different models for model fitting [26]. Table III lists the negative log-likelihood values (NLV), and AIC values for different reliability models when the oin distributions of D 25

26 are selected as a mulitvariate lognormal distribution, a mulitvariate Weibull distribution constructed via the Gumbel-Hougaard copula, and a mulitvariate Weibull distribution constructed via the Gaussian copula, respectively. The NLV for the corresponding models are 46.7,.609, and 0.546, respectively; and the AIC values are , , and , respectively. As smaller NLV or AIC values indicate better model fitting, the performance of the multivariate Weibull distribution constructed via a Gaussian copula is better than the other two models. Thus, the results of the statistical dependency tests based on a reliabity model with Weiblull marginal distributions and a Gaussian copula is selected. We also compare the proposed statistical model that captures component statistical failure dependency to the model that assumes the statistical failure independency in Table III. Under the statistical failure independency assumption, the NLV are , and 3.42, and the AIC values are , and , when lognormal, and Weibull are selected as the failure time distributions for all stations, respectively. Similar to models developed in this research, the reliability model with a Weibull distribution has a better performance. By comparing NLV and AIC values of the fitted models, it can be seen that the proposed statistical model that captures failure dependency fits the observed data best. Table III NLV, and AIC values for different reliability models multivariate lognormal multivariate Weibull distribution 26

27 distribution Our model (station 3 fails statistically independently from others) Gumbel-Hougaard copula Gaussian copula (stations 2&3 fail statistically independently from others) NLV: 46.7 AIC: NLV:.609 AIC: NLV: AIC: Independent model NLV: AIC: NLV: 3.42 AIC: Summary In this paper, a general statistical reliability model is proposed for repairable multi-component systems considering statistical dependent competing risks under a partially perfect repair assumption. For the reliability analysis of repairable multi-component systems, most of the research in the literature assumes component failure statistical independence. The failure mechanism (marginal distribution) of each component can thus be estimated individually based on its failure data. This paper focuses on the reliability analysis of repairable systems with component failure dependency. In the developed model, copula functions are used to model the oint distribution of component failure times. Specifically, two types of copulas, i.e., the Archimedean copula, and the Gaussian copula, are applied to study the overall dependency, and pairwise dependencies among different components, respectively. 27

28 Although the copula function method is also applied in the literature to study non-repairable systems, or systems under perfect repair action (replace the whole system when a failure happens), the problem studied in this paper is much more complex than those in the literature. When the whole system is replaced after a failure, the system will have the same failure mechanism as the original one. In contrast, when only the failed component is replaced, replacement of failed component affects the failure mechanism of the other components when considering failure dependency. Thus, the methods in the literature cannot be directly applied. Under competing risks assumptions, only the failed component is recorded as the latent ages to failures of other components cannot be observed. After a repair action under the partially perfect repair assumption, the failure mechanism of the new system and components will be changed. Thus, for a single repairable system in which the failure data can only be collected from a single realization, model parameters estimation is challenging. To tackle this problem, an ML method is developed in this research, and the ML function is calculated based on conditional probability. The partially perfect repair action is useful for many complex multi-component engineering systems when only the failed component is repaired, and the repair action is a replacement due to high labor cost. In addition, the partially perfect repair can be used to approximate the general imperfect repair actions. But the accuracy of approximation is 28

29 difficult to be evaluated under most situations. A future research topic is to consider repairable dependent competing risks under general repair actions. Hypothesis testing is developed in this paper to test the statistical dependency of component failures. The obtained statistical failure dependence provides more accurate information for reliability predictions, which can be used for system maintenance. It is also interesting to compare the maintenance performance with different repairable actions, e.g., perfect repair, partially perfect repair, and minimal repair. This work is another research topic that will be studied in the future. Appendix The conditional probability of the multivariate normal distribution Let X with mean µ X denote the latent age to failure of component k, and X 2 with mean µ X 2 denote the latent ages to failure of components, k, =,..., K. The variance and covariance matrix of X, X 2 are denoted by 2 σ, and Σ X 2, respectively. In addition, denote the covariance matrix between X and X 2 by Σ. X,X 2 The conditional distribution of X 2 conditioned on X = x is multivariate lognormal distributed with mean µ + Σ ( x µ ) Σ, and covariance matrix T X 2 X,X2 X X2 T 2 X X / X σ 2,X 2,X2 Σ Σ Σ. Thus, the conditional likelihood in (5) can be calculated as a cumulative probability of a multivariate normal distribution, which is well noted in the literature. 29

30 The calculation of the survival function of a multivariate Weibull constructed via the copula function It is well known that the cdf and survival function of a bivariate distribution has the following relationship. S( d, d ) = F(, ) [ F(, d ) + F( d, )] + F( d, d ) Here, we extend this relationship from two dimensions to K dimension as S d µ,σ λ = F F + + F (25) ( ;, ) 0 ( ) K K where F k denotes the summation of cumulative probabilities of the vector d, in which there are exactly K k coordinates replaced by. For example, when K = 3, we can obtain F 2 F d d 2 F d 2 d 3 F d d 3 = (,, ) + (,, ) + (,, ). The calculation of the asymptotic covariance of the ML estimate by using the Delta method Based on ML theory, the ML estimate θ ˆ = ( µ ˆ, M ˆ ) is asymptotically normally distributed with a large-sample assumption [33]. Thus, the asymptotic covariance Σ ˆ θˆ for ˆθ is achievable from the local information matrix ( ˆ Ι θ ), i.e., Σˆ = Ι( θ ˆ). θˆ Ι( θ ˆ) is the negative of the Hessian matrix Hθ ( ) evaluated at θ = θ, ˆ which can be calculated as [33] 2 ɵ log( L( θ)) Iθ ( ) = T θ θ θ= θɵ 30

31 Because the asymptotic covariance matrix of ˆM is part of Σ ˆ ˆ, we use θ Σ ˆ to Mˆ denote it for notation convenience. Given the asymptotic covariance matrix of the estimate θ ˆ = ( µ ˆ, M ˆ ), the asymptotic covariance of the ML estimate Λ can be, calculated based on (8) by using the delta method as ˆ i ˆ Λ Λ var( Λ ) = ( ) Σˆ ( ) M M i, T i, i, Mˆ M= Mˆ (26) References [] H. Ascher and H. Feingold, Repairable systems reliability : modeling, inference, misconceptions and their causes. New York: M. Dekker, 984. [2] B. H. Lindqvist, "On the statistical modeling and analysis of repairable systems," Statistical Science, vol. 2, pp , Nov [3] H. Langseth and B. H. Lindqvist, "Competing risks for repairable systems: A data study," Journal of Statistical Planning and Inference, vol. 36, pp , [4] H. Pham and H. Wang, "Optimal (τ, T) opportunistic maintenance of a k outof n: G system with imperfect PM and partial failure," Naval Research Logistics, vol. 47, pp , [5] L. Wang, J. Chu, and W. Mao, "A condition-based replacement and spare provisioning policy for deteriorating systems with uncertain deterioration to failure," European Journal of Operational Research, vol. 94, pp , [6] Q. Yang, Y. Hong, Y. Chen, and J. Shi, "Failure Profile Analysis of Complex Repairable Systems With Multiple Failure Modes," IEEE Trans. Reliability, vol. 6, pp. 80-9,

32 [7] A. Lehmann, "Joint modeling of degradation and failure time data," Journal of Statistical Planning and Inference, vol. 39, pp , [8] W. J. Li and H. Pham, "Reliability modeling of multi-state degraded systems with multi-competing failures and random shocks," IEEE Trans. Reliability, vol. 54, pp , [9] Y. Zhu, E. A. Elsayed, H. Liao, and L. Y. Chan, "Availability optimization of systems subect to competing risk," European Journal of Operational Research, vol. 202, pp , 200. [0] Q. Yang and Y. Chen, "Sensor system reliability modeling and analysis for fault diagnosis in multistage manufacturing processes," IIE Transactions, vol. 4, pp , [] Q. Yang and Y. Chen, "Reliability of coordinate sensor systems under the risk of sensor precision degradations," IEEE Transactions on Automation Science and Engineering, vol. 7, pp , 200. [2] Q. Yang and Y. Chen, "Monte Carlo Methods for Reliability Evaluation of Linear Sensor Systems," IEEE Transactions on Reliability, vol. 60, pp , 20. [3] J. P. Jhang and S. H. Sheu, "Optimal age and block replacement policies for a multi-component system with failure interaction," International Journal of Systems Science, vol. 3, pp , [4] D. Murthy and D. Nguyen, "Study of a multi-component system with failure interaction," European Journal of Operational Research, vol. 2, pp , 985. [5] T. Nakagawa and D. Murthy, "Optimal replacement policies for a two-unit system with failure interactions," RAIRO Recherche Opérationnelle/Operations Research, vol. 27, pp ,

33 [6] T. Satow and S. Osaki, "Optimal replacement policies for a two-unit system with shock damage interaction," Computers & Mathematics with Applications, vol. 46, pp , [7] P. A. Scarf and M. Deara, "Block replacement policies for a two component system with failure dependence," Naval Research Logistics, vol. 50, pp , [8] R. I. Zequeira and C. Bérenguer, "On the inspection policy of a two-component parallel system with failure interaction," Reliability Engineering & System Safety, vol. 88, pp , [9] A. Barros, C. Berenguer, and A. Grall, "A maintenance policy for two-unit parallel systems based on imperfect monitoring information," Reliability Engineering & System Safety, vol. 9, pp. 3-36, [20] S. H. Sheu and C. T. Liou, "Optimal replacement of a k-out-of-n system subect to shocks," Microelectronics Reliability, vol. 32, pp , 992. [2] Z. Pan and N. Balakrishnan, "Reliability modeling of degradation of products with multiple performance characteristics based on gamma processes," Reliability Engineering & System Safety, vol. 96, pp , 20. [22] H. Peng, Q. M. Feng, and D. W. Coit, "Reliability and maintenance modeling for systems subect to multiple dependent competing failure processes," IIE Transactions, vol. 43, pp. 2-22, 20. [23] L. Jiang, Q. M. Feng, and D. W. Coit, "Reliability and Maintenance Modeling for Dependent Competing Failure Processes With Shifting Failure Thresholds," IEEE Trans. Reliability, vol. 6, pp , 202. [24] Y. P. Wang and H. Pham, "A Multi-Obective Optimization of Imperfect Preventive Maintenance Policy for Dependent Competing Risk Systems With Hidden Failure," IEEE Trans. Reliability, vol. 60, pp ,

34 [25] S. M. S. Lo and R. A. Wilke, "A copula model for dependent competing risks," Journal of the Royal Statistical Society: Series C (Applied Statistics), vol. 59, pp , 200. [26] Y. P. Wang and H. Pham, "Modeling the Dependent Competing Risks With Multiple Degradation Processes and Random Shock Using Time-Varying Copulas," IEEE Trans. Reliability, vol. 6, pp. 3-22, 202. [27] V. K. Kaishev, D. S. Dimitrova, and S. Haberman, "Modelling the oint distribution of competing risks survival times using copula functions," Insurance: Mathematics and Economics, vol. 4, pp , [28] J. K. Sari, M. J. Newby, A. C. Brombacher, and L. C. Tang, "Bivariate constant stress degradation model: LED lighting system reliability estimation with two-stage modelling," Quality and Reliability Engineering International, vol. 25, pp , [29] R. B. Nelsen, An introduction to copulas, 2nd ed. New York: Springer, [30] I. Zezula, "On multivariate Gaussian copulas," Journal of Statistical Planning and Inference, vol. 39, pp , [3] W. Q. Meeker and L. A. Escobar, Statistical methods for reliability data. New York: Wiley, 998. [32] P. X. K. Song, "Multivariate dispersion models generated from Gaussian copula," Scandinavian Journal of Statistics, vol. 27, pp , Jun [33] G. Casella and R. L. Berger, Statistical Inference, 2nd ed. Pacific Grove, Calif.: Duxbury Press, 200. [34] Å. Svensson, "Asymptotic Estimation in Counting Processes with Parametric Intensities Based on One Realization," Scandinavian Journal of Statistics, vol. 7, pp ,

35 [35] R. Rebonato and P. Jäckel, "The most general methodology for creating a valid correlation matrix for risk management and option pricing purposes," The Journal of Risk, vol. 2, pp. 7-28, 999. QingyuYang received the B.S., and M.S. degrees in Automatic Control and Intelligent System from the University of Science and Technology of China in 2000, and 2003 respectively; and the M.S. degree in Statistics, and the Ph.D. degree in Industrial Engineering, from the University of Iowa in 2007, and 2008 respectively. Currently, he is an Assistant Professor in the Department of Industrial and System Engineering at Wayne State University. His research interests include statistical data analysis, reliability and quality, and complex system modeling. He is a member of INFORMS and IIE. Nailong Zhang is a Ph.D. candidate in the Department of Industrial and Systems Engineering at Wayne State University. He received a B.S. degree in Mechanical Engineering from Harbin Institute of Technology, Harbin, China, in His research interests include reliability & maintenance engineering and optimization, especially reliability modeling and analysis of complex systems under competing risks. Yili Hong received a BS in statistics (2004) from University of Science and Technology of China; and MS, and PhD degrees in statistics (2005, 2009) from Iowa State University. He is currently an Assistant Professor in the Department of Statistics at Virginia Tech. His research interests include reliability data analysis, and engineering statistics. 35

Bivariate Degradation Modeling Based on Gamma Process

Bivariate Degradation Modeling Based on Gamma Process Bivariate Degradation Modeling Based on Gamma Process Jinglun Zhou Zhengqiang Pan Member IAENG and Quan Sun Abstract Many highly reliable products have two or more performance characteristics (PCs). The

More information

Simultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions

Simultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions Statistics Preprints Statistics 10-2014 Simultaneous Prediction Intervals for the (Log)- Location-Scale Family of Distributions Yimeng Xie Virginia Tech Yili Hong Virginia Tech Luis A. Escobar Louisiana

More information

STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL. A Thesis. Presented to the. Faculty of. San Diego State University

STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL. A Thesis. Presented to the. Faculty of. San Diego State University STATISTICAL INFERENCE IN ACCELERATED LIFE TESTING WITH GEOMETRIC PROCESS MODEL A Thesis Presented to the Faculty of San Diego State University In Partial Fulfillment of the Requirements for the Degree

More information

Multivariate Distribution Models

Multivariate Distribution Models Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is

More information

A Tool for Evaluating Time-Varying-Stress Accelerated Life Test Plans with Log-Location- Scale Distributions

A Tool for Evaluating Time-Varying-Stress Accelerated Life Test Plans with Log-Location- Scale Distributions Statistics Preprints Statistics 6-2010 A Tool for Evaluating Time-Varying-Stress Accelerated Life Test Plans with Log-Location- Scale Distributions Yili Hong Virginia Tech Haiming Ma Iowa State University,

More information

Financial Econometrics and Volatility Models Copulas

Financial Econometrics and Volatility Models Copulas Financial Econometrics and Volatility Models Copulas Eric Zivot Updated: May 10, 2010 Reading MFTS, chapter 19 FMUND, chapters 6 and 7 Introduction Capturing co-movement between financial asset returns

More information

Statistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions

Statistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS040) p.4828 Statistical Inference on Constant Stress Accelerated Life Tests Under Generalized Gamma Lifetime Distributions

More information

Time-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 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 information

The Relationship Between Confidence Intervals for Failure Probabilities and Life Time Quantiles

The Relationship Between Confidence Intervals for Failure Probabilities and Life Time Quantiles Statistics Preprints Statistics 2008 The Relationship Between Confidence Intervals for Failure Probabilities and Life Time Quantiles Yili Hong Iowa State University, yili_hong@hotmail.com William Q. Meeker

More information

Estimation of Operational Risk Capital Charge under Parameter Uncertainty

Estimation of Operational Risk Capital Charge under Parameter Uncertainty Estimation of Operational Risk Capital Charge under Parameter Uncertainty Pavel V. Shevchenko Principal Research Scientist, CSIRO Mathematical and Information Sciences, Sydney, Locked Bag 17, North Ryde,

More information

Package CopulaRegression

Package CopulaRegression Type Package Package CopulaRegression Title Bivariate Copula Based Regression Models Version 0.1-5 Depends R (>= 2.11.0), MASS, VineCopula Date 2014-09-04 Author, Daniel Silvestrini February 19, 2015 Maintainer

More information

1 Appendix A: Matrix Algebra

1 Appendix A: Matrix Algebra Appendix A: Matrix Algebra. Definitions Matrix A =[ ]=[A] Symmetric matrix: = for all and Diagonal matrix: 6=0if = but =0if 6= Scalar matrix: the diagonal matrix of = Identity matrix: the scalar matrix

More information

Bayesian Life Test Planning for the Log-Location- Scale Family of Distributions

Bayesian Life Test Planning for the Log-Location- Scale Family of Distributions Statistics Preprints Statistics 3-14 Bayesian Life Test Planning for the Log-Location- Scale Family of Distributions Yili Hong Virginia Tech Caleb King Virginia Tech Yao Zhang Pfizer Global Research and

More information

Statistical Inference and Methods

Statistical Inference and Methods Department of Mathematics Imperial College London d.stephens@imperial.ac.uk http://stats.ma.ic.ac.uk/ das01/ 31st January 2006 Part VI Session 6: Filtering and Time to Event Data Session 6: Filtering and

More information

Point and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples

Point and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples 90 IEEE TRANSACTIONS ON RELIABILITY, VOL. 52, NO. 1, MARCH 2003 Point and Interval Estimation for Gaussian Distribution, Based on Progressively Type-II Censored Samples N. Balakrishnan, N. Kannan, C. T.

More information

INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION

INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION Pak. J. Statist. 2017 Vol. 33(1), 37-61 INVERTED KUMARASWAMY DISTRIBUTION: PROPERTIES AND ESTIMATION A. M. Abd AL-Fattah, A.A. EL-Helbawy G.R. AL-Dayian Statistics Department, Faculty of Commerce, AL-Azhar

More information

Frailty Models and Copulas: Similarities and Differences

Frailty Models and Copulas: Similarities and Differences Frailty Models and Copulas: Similarities and Differences KLARA GOETHALS, PAUL JANSSEN & LUC DUCHATEAU Department of Physiology and Biometrics, Ghent University, Belgium; Center for Statistics, Hasselt

More information

Probability and Stochastic Processes

Probability and Stochastic Processes Probability and Stochastic Processes A Friendly Introduction Electrical and Computer Engineers Third Edition Roy D. Yates Rutgers, The State University of New Jersey David J. Goodman New York University

More information

Constant Stress Partially Accelerated Life Test Design for Inverted Weibull Distribution with Type-I Censoring

Constant Stress Partially Accelerated Life Test Design for Inverted Weibull Distribution with Type-I Censoring Algorithms Research 013, (): 43-49 DOI: 10.593/j.algorithms.01300.0 Constant Stress Partially Accelerated Life Test Design for Mustafa Kamal *, Shazia Zarrin, Arif-Ul-Islam Department of Statistics & Operations

More information

A new condition based maintenance model with random improvements on the system after maintenance actions: Optimizing by monte carlo simulation

A new condition based maintenance model with random improvements on the system after maintenance actions: Optimizing by monte carlo simulation ISSN 1 746-7233, England, UK World Journal of Modelling and Simulation Vol. 4 (2008) No. 3, pp. 230-236 A new condition based maintenance model with random improvements on the system after maintenance

More information

Distribution Fitting (Censored Data)

Distribution Fitting (Censored Data) Distribution Fitting (Censored Data) Summary... 1 Data Input... 2 Analysis Summary... 3 Analysis Options... 4 Goodness-of-Fit Tests... 6 Frequency Histogram... 8 Comparison of Alternative Distributions...

More information

Failure modeling and maintenance optimization for a railway line

Failure modeling and maintenance optimization for a railway line Failure modeling and maintenance optimization for a railway line Per Hokstad *, Helge Langseth SINTEF Technology and Society, Department of Safety and Reliability, N-7465 Trondheim, Norway Bo H. Lindqvist,

More information

Strategic Allocation of Test Units in an Accelerated Degradation Test Plan

Strategic Allocation of Test Units in an Accelerated Degradation Test Plan Strategic Allocation of Test Units in an Accelerated Degradation Test Plan Zhi-Sheng YE, Qingpei HU, and Dan YU Department of Industrial and Systems Engineering, National University of Singapore Academy

More information

Multistate Modeling and Applications

Multistate Modeling and Applications Multistate Modeling and Applications Yang Yang Department of Statistics University of Michigan, Ann Arbor IBM Research Graduate Student Workshop: Statistics for a Smarter Planet Yang Yang (UM, Ann Arbor)

More information

Some New Aspects of Dose-Response Models with Applications to Multistage Models Having Parameters on the Boundary

Some New Aspects of Dose-Response Models with Applications to Multistage Models Having Parameters on the Boundary Some New Aspects of Dose-Response Models with Applications to Multistage Models Having Parameters on the Boundary Bimal Sinha Department of Mathematics & Statistics University of Maryland, Baltimore County,

More information

Integrating Quality and Inspection for the Optimal Lot-sizing Problem with Rework, Minimal Repair and Inspection Time

Integrating Quality and Inspection for the Optimal Lot-sizing Problem with Rework, Minimal Repair and Inspection Time Proceedings of the 2011 International Conference on Industrial Engineering and Operations Management Kuala Lumpur, Malaysia, January 22 24, 2011 Integrating Quality and Inspection for the Optimal Lot-sizing

More information

Design of Optimal Bayesian Reliability Test Plans for a Series System

Design of Optimal Bayesian Reliability Test Plans for a Series System Volume 109 No 9 2016, 125 133 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://wwwijpameu ijpameu Design of Optimal Bayesian Reliability Test Plans for a Series System P

More information

Dependence. MFM Practitioner Module: Risk & Asset Allocation. John Dodson. September 11, Dependence. John Dodson. Outline.

Dependence. MFM Practitioner Module: Risk & Asset Allocation. John Dodson. September 11, Dependence. John Dodson. Outline. MFM Practitioner Module: Risk & Asset Allocation September 11, 2013 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y

More information

A novel repair model for imperfect maintenance

A 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 information

The Instability of Correlations: Measurement and the Implications for Market Risk

The Instability of Correlations: Measurement and the Implications for Market Risk The Instability of Correlations: Measurement and the Implications for Market Risk Prof. Massimo Guidolin 20254 Advanced Quantitative Methods for Asset Pricing and Structuring Winter/Spring 2018 Threshold

More information

Bayesian Life Test Planning for the Weibull Distribution with Given Shape Parameter

Bayesian Life Test Planning for the Weibull Distribution with Given Shape Parameter Statistics Preprints Statistics 10-8-2002 Bayesian Life Test Planning for the Weibull Distribution with Given Shape Parameter Yao Zhang Iowa State University William Q. Meeker Iowa State University, wqmeeker@iastate.edu

More information

Chapter 15. System Reliability Concepts and Methods. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University

Chapter 15. System Reliability Concepts and Methods. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Chapter 15 System Reliability Concepts and Methods William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on

More information

Marginal Specifications and a Gaussian Copula Estimation

Marginal Specifications and a Gaussian Copula Estimation Marginal Specifications and a Gaussian Copula Estimation Kazim Azam Abstract Multivariate analysis involving random variables of different type like count, continuous or mixture of both is frequently required

More information

Correlation: Copulas and Conditioning

Correlation: Copulas and Conditioning Correlation: Copulas and Conditioning This note reviews two methods of simulating correlated variates: copula methods and conditional distributions, and the relationships between them. Particular emphasis

More information

Dependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.

Dependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline. Practitioner Course: Portfolio Optimization September 10, 2008 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y ) (x,

More information

Modelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich

Modelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich Modelling Dependence with Copulas and Applications to Risk Management Filip Lindskog, RiskLab, ETH Zürich 02-07-2000 Home page: http://www.math.ethz.ch/ lindskog E-mail: lindskog@math.ethz.ch RiskLab:

More information

When is a copula constant? A test for changing relationships

When is a copula constant? A test for changing relationships When is a copula constant? A test for changing relationships Fabio Busetti and Andrew Harvey Bank of Italy and University of Cambridge November 2007 usetti and Harvey (Bank of Italy and University of Cambridge)

More information

A Dynamic Programming Approach for Sequential Preventive Maintenance Policies with Two Failure Modes

A Dynamic Programming Approach for Sequential Preventive Maintenance Policies with Two Failure Modes Chapter 1 A Dynamic Programming Approach for Sequential Preventive Maintenance Policies with Two Failure Modes Hiroyuki Okamura 1, Tadashi Dohi 1 and Shunji Osaki 2 1 Department of Information Engineering,

More information

Optimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests

Optimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests International Journal of Performability Engineering, Vol., No., January 24, pp.3-4. RAMS Consultants Printed in India Optimum Test Plan for 3-Step, Step-Stress Accelerated Life Tests N. CHANDRA *, MASHROOR

More information

Introduction to Maximum Likelihood Estimation

Introduction to Maximum Likelihood Estimation Introduction to Maximum Likelihood Estimation Eric Zivot July 26, 2012 The Likelihood Function Let 1 be an iid sample with pdf ( ; ) where is a ( 1) vector of parameters that characterize ( ; ) Example:

More information

Load-strength Dynamic Interaction Principle and Failure Rate Model

Load-strength Dynamic Interaction Principle and Failure Rate Model International Journal of Performability Engineering Vol. 6, No. 3, May 21, pp. 25-214. RAMS Consultants Printed in India Load-strength Dynamic Interaction Principle and Failure Rate Model LIYANG XIE and

More information

Comparison study between MCMC-based and weight-based Bayesian methods for identification of joint distribution

Comparison study between MCMC-based and weight-based Bayesian methods for identification of joint distribution Struct Multidisc Optim (2010) 42:823 833 DOI 10.1007/s00158-010-0539-1 RESEARCH PAPER Comparison study between MCMC-based and weight-based Bayesian methods for identification of joint distribution Yoojeong

More information

How to select a good vine

How to select a good vine Universitetet i Oslo ingrihaf@math.uio.no International FocuStat Workshop on Focused Information Criteria and Related Themes, May 9-11, 2016 Copulae Regular vines Model selection and reduction Limitations

More information

Bayesian Methods for Accelerated Destructive Degradation Test Planning

Bayesian Methods for Accelerated Destructive Degradation Test Planning Statistics Preprints Statistics 11-2010 Bayesian Methods for Accelerated Destructive Degradation Test Planning Ying Shi Iowa State University William Q. Meeker Iowa State University, wqmeeker@iastate.edu

More information

Estimation Under Multivariate Inverse Weibull Distribution

Estimation Under Multivariate Inverse Weibull Distribution Global Journal of Pure and Applied Mathematics. ISSN 097-768 Volume, Number 8 (07), pp. 4-4 Research India Publications http://www.ripublication.com Estimation Under Multivariate Inverse Weibull Distribution

More information

Step-Stress Models and Associated Inference

Step-Stress Models and Associated Inference Department of Mathematics & Statistics Indian Institute of Technology Kanpur August 19, 2014 Outline Accelerated Life Test 1 Accelerated Life Test 2 3 4 5 6 7 Outline Accelerated Life Test 1 Accelerated

More information

Reliability Analysis of Tampered Failure Rate Load-Sharing k-out-of-n:g Systems

Reliability Analysis of Tampered Failure Rate Load-Sharing k-out-of-n:g Systems Reliability Analysis of Tampered Failure Rate Load-Sharing k-out-of-n:g Systems Suprasad V. Amari Relex Software Corporation 540 Pellis Road Greensburg, PA 15601 USA Krishna B. Misra RAMS Consultants 71

More information

MISCELLANEOUS TOPICS RELATED TO LIKELIHOOD. Copyright c 2012 (Iowa State University) Statistics / 30

MISCELLANEOUS TOPICS RELATED TO LIKELIHOOD. Copyright c 2012 (Iowa State University) Statistics / 30 MISCELLANEOUS TOPICS RELATED TO LIKELIHOOD Copyright c 2012 (Iowa State University) Statistics 511 1 / 30 INFORMATION CRITERIA Akaike s Information criterion is given by AIC = 2l(ˆθ) + 2k, where l(ˆθ)

More information

Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories

Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories Entropy 2012, 14, 1784-1812; doi:10.3390/e14091784 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories Lan Zhang

More information

Bayesian Methods for Estimating the Reliability of Complex Systems Using Heterogeneous Multilevel Information

Bayesian Methods for Estimating the Reliability of Complex Systems Using Heterogeneous Multilevel Information Statistics Preprints Statistics 8-2010 Bayesian Methods for Estimating the Reliability of Complex Systems Using Heterogeneous Multilevel Information Jiqiang Guo Iowa State University, jqguo@iastate.edu

More information

Overview of Extreme Value Theory. Dr. Sawsan Hilal space

Overview of Extreme Value Theory. Dr. Sawsan Hilal space Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate

More information

A CONDITION-BASED MAINTENANCE MODEL FOR AVAILABILITY OPTIMIZATION FOR STOCHASTIC DEGRADING SYSTEMS

A CONDITION-BASED MAINTENANCE MODEL FOR AVAILABILITY OPTIMIZATION FOR STOCHASTIC DEGRADING SYSTEMS A CONDITION-BASED MAINTENANCE MODEL FOR AVAILABILITY OPTIMIZATION FOR STOCHASTIC DEGRADING SYSTEMS Abdelhakim Khatab, Daoud Ait-Kadi, Nidhal Rezg To cite this version: Abdelhakim Khatab, Daoud Ait-Kadi,

More information

Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach

Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score

More information

Simulating Exchangeable Multivariate Archimedean Copulas and its Applications. Authors: Florence Wu Emiliano A. Valdez Michael Sherris

Simulating Exchangeable Multivariate Archimedean Copulas and its Applications. Authors: Florence Wu Emiliano A. Valdez Michael Sherris Simulating Exchangeable Multivariate Archimedean Copulas and its Applications Authors: Florence Wu Emiliano A. Valdez Michael Sherris Literatures Frees and Valdez (1999) Understanding Relationships Using

More information

Evaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions

Evaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions Evaluating the value of structural heath monitoring with longitudinal performance indicators and hazard functions using Bayesian dynamic predictions C. Xing, R. Caspeele, L. Taerwe Ghent University, Department

More information

Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011

Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011 Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011 Outline Ordinary Least Squares (OLS) Regression Generalized Linear Models

More information

Structural Reliability

Structural Reliability Structural Reliability Thuong Van DANG May 28, 2018 1 / 41 2 / 41 Introduction to Structural Reliability Concept of Limit State and Reliability Review of Probability Theory First Order Second Moment Method

More information

Multivariate Non-Normally Distributed Random Variables

Multivariate Non-Normally Distributed Random Variables Multivariate Non-Normally Distributed Random Variables An Introduction to the Copula Approach Workgroup seminar on climate dynamics Meteorological Institute at the University of Bonn 18 January 2008, Bonn

More information

Multivariate Assays With Values Below the Lower Limit of Quantitation: Parametric Estimation By Imputation and Maximum Likelihood

Multivariate Assays With Values Below the Lower Limit of Quantitation: Parametric Estimation By Imputation and Maximum Likelihood Multivariate Assays With Values Below the Lower Limit of Quantitation: Parametric Estimation By Imputation and Maximum Likelihood Robert E. Johnson and Heather J. Hoffman 2* Department of Biostatistics,

More information

COMPOSITE RELIABILITY MODELS FOR SYSTEMS WITH TWO DISTINCT KINDS OF STOCHASTIC DEPENDENCES BETWEEN THEIR COMPONENTS LIFE TIMES

COMPOSITE RELIABILITY MODELS FOR SYSTEMS WITH TWO DISTINCT KINDS OF STOCHASTIC DEPENDENCES BETWEEN THEIR COMPONENTS LIFE TIMES COMPOSITE RELIABILITY MODELS FOR SYSTEMS WITH TWO DISTINCT KINDS OF STOCHASTIC DEPENDENCES BETWEEN THEIR COMPONENTS LIFE TIMES Jerzy Filus Department of Mathematics and Computer Science, Oakton Community

More information

Estimation of Quantiles

Estimation of Quantiles 9 Estimation of Quantiles The notion of quantiles was introduced in Section 3.2: recall that a quantile x α for an r.v. X is a constant such that P(X x α )=1 α. (9.1) In this chapter we examine quantiles

More information

Chapter 17. Failure-Time Regression Analysis. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University

Chapter 17. Failure-Time Regression Analysis. William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Chapter 17 Failure-Time Regression Analysis William Q. Meeker and Luis A. Escobar Iowa State University and Louisiana State University Copyright 1998-2008 W. Q. Meeker and L. A. Escobar. Based on the authors

More information

Condition based maintenance optimization using neural network based health condition prediction

Condition based maintenance optimization using neural network based health condition prediction Condition based maintenance optimization using neural network based health condition prediction Bairong Wu a,b, Zhigang Tian a,, Mingyuan Chen b a Concordia Institute for Information Systems Engineering,

More information

Correlation & Dependency Structures

Correlation & Dependency Structures Correlation & Dependency Structures GIRO - October 1999 Andrzej Czernuszewicz Dimitris Papachristou Why are we interested in correlation/dependency? Risk management Portfolio management Reinsurance purchase

More information

The comparative studies on reliability for Rayleigh models

The comparative studies on reliability for Rayleigh models Journal of the Korean Data & Information Science Society 018, 9, 533 545 http://dx.doi.org/10.7465/jkdi.018.9..533 한국데이터정보과학회지 The comparative studies on reliability for Rayleigh models Ji Eun Oh 1 Joong

More information

An Analysis of Record Statistics based on an Exponentiated Gumbel Model

An Analysis of Record Statistics based on an Exponentiated Gumbel Model Communications for Statistical Applications and Methods 2013, Vol. 20, No. 5, 405 416 DOI: http://dx.doi.org/10.5351/csam.2013.20.5.405 An Analysis of Record Statistics based on an Exponentiated Gumbel

More information

A Simulation Study on Confidence Interval Procedures of Some Mean Cumulative Function Estimators

A 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 information

Quality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2

Quality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2 Quality Control Using Inferential Statistics In Weibull Based Reliability Analyses S. F. Duffy 1 and A. Parikh 2 1 Cleveland State University 2 N & R Engineering www.inl.gov ASTM Symposium on Graphite

More information

Chapter 2 Inference on Mean Residual Life-Overview

Chapter 2 Inference on Mean Residual Life-Overview Chapter 2 Inference on Mean Residual Life-Overview Statistical inference based on the remaining lifetimes would be intuitively more appealing than the popular hazard function defined as the risk of immediate

More information

Modelling Dependent Credit Risks

Modelling Dependent Credit Risks Modelling Dependent Credit Risks Filip Lindskog, RiskLab, ETH Zürich 30 November 2000 Home page:http://www.math.ethz.ch/ lindskog E-mail:lindskog@math.ethz.ch RiskLab:http://www.risklab.ch Modelling Dependent

More information

Modeling 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 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 information

A Goodness-of-fit Test for Copulas

A Goodness-of-fit Test for Copulas A Goodness-of-fit Test for Copulas Artem Prokhorov August 2008 Abstract A new goodness-of-fit test for copulas is proposed. It is based on restrictions on certain elements of the information matrix and

More information

A measure of radial asymmetry for bivariate copulas based on Sobolev norm

A measure of radial asymmetry for bivariate copulas based on Sobolev norm A measure of radial asymmetry for bivariate copulas based on Sobolev norm Ahmad Alikhani-Vafa Ali Dolati Abstract The modified Sobolev norm is used to construct an index for measuring the degree of radial

More information

Fall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.

Fall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A. 1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n

More information

Lifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution

Lifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution Lifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution Daniel Alai Zinoviy Landsman Centre of Excellence in Population Ageing Research (CEPAR) School of Mathematics, Statistics

More information

Fleet Maintenance Simulation With Insufficient Data

Fleet Maintenance Simulation With Insufficient Data Fleet Maintenance Simulation With Insufficient Data Zissimos P. Mourelatos Mechanical Engineering Department Oakland University mourelat@oakland.edu Ground Robotics Reliability Center (GRRC) Seminar 17

More information

Introducing the Normal Distribution

Introducing the Normal Distribution Department of Mathematics Ma 3/103 KC Border Introduction to Probability and Statistics Winter 2017 Lecture 10: Introducing the Normal Distribution Relevant textbook passages: Pitman [5]: Sections 1.2,

More information

Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution

Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution Journal of Probability and Statistical Science 14(), 11-4, Aug 016 Some Theoretical Properties and Parameter Estimation for the Two-Sided Length Biased Inverse Gaussian Distribution Teerawat Simmachan

More information

Solutions of the Financial Risk Management Examination

Solutions of the Financial Risk Management Examination Solutions of the Financial Risk Management Examination Thierry Roncalli January 9 th 03 Remark The first five questions are corrected in TR-GDR and in the document of exercise solutions, which is available

More information

Reliability Engineering I

Reliability 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

Introduction to Reliability Theory (part 2)

Introduction 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 information

Web-based Supplementary Material for A Two-Part Joint. Model for the Analysis of Survival and Longitudinal Binary. Data with excess Zeros

Web-based Supplementary Material for A Two-Part Joint. Model for the Analysis of Survival and Longitudinal Binary. Data with excess Zeros Web-based Supplementary Material for A Two-Part Joint Model for the Analysis of Survival and Longitudinal Binary Data with excess Zeros Dimitris Rizopoulos, 1 Geert Verbeke, 1 Emmanuel Lesaffre 1 and Yves

More information

e 4β e 4β + e β ˆβ =0.765

e 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 information

Calibration Estimation of Semiparametric Copula Models with Data Missing at Random

Calibration Estimation of Semiparametric Copula Models with Data Missing at Random Calibration Estimation of Semiparametric Copula Models with Data Missing at Random Shigeyuki Hamori 1 Kaiji Motegi 1 Zheng Zhang 2 1 Kobe University 2 Renmin University of China Econometrics Workshop UNC

More information

Identification of marginal and joint CDFs using Bayesian method for RBDO

Identification of marginal and joint CDFs using Bayesian method for RBDO Struct Multidisc Optim (2010) 40:35 51 DOI 10.1007/s00158-009-0385-1 RESEARCH PAPER Identification of marginal and joint CDFs using Bayesian method for RBDO Yoojeong Noh K. K. Choi Ikjin Lee Received:

More information

Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems

Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems Analysis of the AIC Statistic for Optimal Detection of Small Changes in Dynamic Systems Jeremy S. Conner and Dale E. Seborg Department of Chemical Engineering University of California, Santa Barbara, CA

More information

Copulas. MOU Lili. December, 2014

Copulas. MOU Lili. December, 2014 Copulas MOU Lili December, 2014 Outline Preliminary Introduction Formal Definition Copula Functions Estimating the Parameters Example Conclusion and Discussion Preliminary MOU Lili SEKE Team 3/30 Probability

More information

Warranty Prediction Based on Auxiliary Use-rate Information

Warranty Prediction Based on Auxiliary Use-rate Information Statistics Preprints Statistics 2009 Warranty Prediction Based on Auxiliary Use-rate Information Yili Hong Virginia Polytechnic Institute and State University William Q. Meeker Iowa State University, wqmeeker@iastate.edu

More information

Notes largely based on Statistical Methods for Reliability Data by W.Q. Meeker and L. A. Escobar, Wiley, 1998 and on their class notes.

Notes 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 information

Semi-parametric predictive inference for bivariate data using copulas

Semi-parametric predictive inference for bivariate data using copulas Semi-parametric predictive inference for bivariate data using copulas Tahani Coolen-Maturi a, Frank P.A. Coolen b,, Noryanti Muhammad b a Durham University Business School, Durham University, Durham, DH1

More information

Bayesian Regression Linear and Logistic Regression

Bayesian Regression Linear and Logistic Regression When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we

More information

A copula goodness-of-t approach. conditional probability integral transform. Daniel Berg 1 Henrik Bakken 2

A copula goodness-of-t approach. conditional probability integral transform. Daniel Berg 1 Henrik Bakken 2 based on the conditional probability integral transform Daniel Berg 1 Henrik Bakken 2 1 Norwegian Computing Center (NR) & University of Oslo (UiO) 2 Norwegian University of Science and Technology (NTNU)

More information

Appendix A: The time series behavior of employment growth

Appendix A: The time series behavior of employment growth Unpublished appendices from The Relationship between Firm Size and Firm Growth in the U.S. Manufacturing Sector Bronwyn H. Hall Journal of Industrial Economics 35 (June 987): 583-606. Appendix A: The time

More information

Statistical Inference Using Progressively Type-II Censored Data with Random Scheme

Statistical Inference Using Progressively Type-II Censored Data with Random Scheme International Mathematical Forum, 3, 28, no. 35, 1713-1725 Statistical Inference Using Progressively Type-II Censored Data with Random Scheme Ammar M. Sarhan 1 and A. Abuammoh Department of Statistics

More information

Estimation for inverse Gaussian Distribution Under First-failure Progressive Hybird Censored Samples

Estimation for inverse Gaussian Distribution Under First-failure Progressive Hybird Censored Samples Filomat 31:18 (217), 5743 5752 https://doi.org/1.2298/fil1718743j Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Estimation for

More information

LOGISTIC REGRESSION Joseph M. Hilbe

LOGISTIC REGRESSION Joseph M. Hilbe LOGISTIC REGRESSION Joseph M. Hilbe Arizona State University Logistic regression is the most common method used to model binary response data. When the response is binary, it typically takes the form of

More information

COM336: Neural Computing

COM336: Neural Computing COM336: Neural Computing http://www.dcs.shef.ac.uk/ sjr/com336/ Lecture 2: Density Estimation Steve Renals Department of Computer Science University of Sheffield Sheffield S1 4DP UK email: s.renals@dcs.shef.ac.uk

More information

A Regression Model for the Copula Graphic Estimator

A Regression Model for the Copula Graphic Estimator Discussion Papers in Economics Discussion Paper No. 11/04 A Regression Model for the Copula Graphic Estimator S.M.S. Lo and R.A. Wilke April 2011 2011 DP 11/04 A Regression Model for the Copula Graphic

More information

Using copulas to model time dependence in stochastic frontier models

Using copulas to model time dependence in stochastic frontier models Using copulas to model time dependence in stochastic frontier models Christine Amsler Michigan State University Artem Prokhorov Concordia University November 2008 Peter Schmidt Michigan State University

More information

Songklanakarin Journal of Science and Technology SJST R1 Sukparungsee

Songklanakarin Journal of Science and Technology SJST R1 Sukparungsee Songklanakarin Journal of Science and Technology SJST-0-0.R Sukparungsee Bivariate copulas on the exponentially weighted moving average control chart Journal: Songklanakarin Journal of Science and Technology

More information