Influence diagnostics for polyhazard models in the presence of covariates

In this paper, we present various diagnostic methods for polyhazard models. Polyhazard models are a flexible family for fitting lifetime data. Their main advantage over the single hazard models, such as the Weibull and the log-logistic models, is to include a large amount of nonmonotone hazard shapes, as bathtub and multimodal curves. Some influence methods, such as the local influence and total local influence of an individual are derived, analyzed and discussed. A discussion of the computation of the likelihood displacement as well as the normal curvature in the local influence method are presented. Finally, an example with real data is given for illustration.

The exponentiated generalized gamma distribution with application to lifetime data

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.

Comparison of HRT-3 glaucoma probability score and subjective stereophotograph assessment for prediction of progression in glaucoma

PURPOSE. To assess whether baseline Glaucoma Probability Score (GPS; HRT-3; Heidelberg Engineering, Dossenheim, Germany) results are predictive of progression in patients with suspected glaucoma. The GPS is a new feature of the confocal scanning laser ophthalmoscope that generates an operator-independent, three-dimensional model of the optic nerve head and gives a score for the probability that this model is consistent with glaucomatous damage. METHODS. The study included 223 patients with suspected glaucoma during an average follow-up of 63.3 months. Included subjects had a suspect optic disc appearance and/or elevated intraocular pressure, but normal visual fields. Conversion was defined as development of either repeatable abnormal visual fields or glaucomatous deterioration in the appearance of the optic disc during the study period. The association between baseline GPS and conversion was investigated by Cox regression models. RESULTS. Fifty-four (24.2%) eyes converted. In multivariate models, both higher values of GPS global and subjective stereophotograph assessment ( larger cup-disc ratio and glaucomatous grading) were predictive of conversion: adjusted hazard ratios (95% CI): 1.31 (1.15 - 1.50) per 0.1 higher global GPS, 1.34 (1.12 - 1.62) per 0.1 higher CDR, and 2.34 (1.22 - 4.47) for abnormal grading, respectively. No significant differences ( P > 0.05 for all comparisons) were found between the c-index values ( equivalent to area under ROC curve) for the multivariate models (0.732, 0.705, and 0.699, respectively). CONCLUSIONS. GPS values were predictive of conversion in our population of patients with suspected glaucoma. Further, they performed as well as subjective assessment of the optic disc. These results suggest that GPS could potentially replace stereophotograph as a tool for estimating the likelihood of conversion to glaucoma.

A critical view on biaxial and short-beam uniaxial flexural strength tests applied to resin composites using Weibull, fractographic and finite element analyses

Objective. To evaluate the biaxial and short-beam uniaxial strength tests applied to resin composites based upon their Weibull parameters, fractographic features and stress distribution. Methods. Disk- (15 mm x 1 mm) and beam-shaped specimens (10 mm x 2 mm x 1 mm) of three commercial composites (Concept/Vigodent, CA; Heliomolar/Ivoclar-Vivadent, HE; Z250/3M ESPE, FZ) were prepared. After 48h dry storage at 37 degrees C, disks and beams were submitted to piston-on-three-balls (BI) and three-point bending (UNI) tests, respectively. Data were analyzed by Weibull statistics. Fractured surfaces were observed under stereomicroscope and scanning electron microscope. Maximum principal stress (sigma(1)) distribution was determined by finite element analysis (FEA). Maximum sigma(1-BI) and sigma(1-UNI) were compared to FZ strengths calculated by applying the average failure loads to the analytical equations (sigma(a-BI) and sigma(a-UNI)). Results. For BI, characteristic strengths were: 169.9a (FZ), 122.4b (CA) and 104.8c (HE), and for UNI were: 160.3a (FZ), 98.2b (CA) and 91.6b (HE). Weibull moduli ( m) were similar within the same test. CA and HE presented statistically higher m for BI. Surface pores ( BI) and edge flaws ( UNI) were the most frequent fracture origins. sigma(1-BI) was 14% lower than sigma(a-BI.) sigma(1-UNI) was 43% higher than sigma(a-UNI). Significance. Compared to the short-beam uniaxial test, the biaxial test detected more differences among composites and displayed less data scattering for two of the tested materials. Also, biaxial strength was closer to the material`s strength estimated by FEA. (C) 2009 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Effect of ion exchange on strength and slow crack growth of a dental porcelain

Objectives. To determine the effect of ion exchange on slow crack growth (SCG) parameters (n, stress corrosion susceptibility coefficient, and sigma(f0), scaling parameter) and Weibull parameters (m, Weibull modulus, and sigma(0), characteristic strength) of a dental porcelain. Methods. 160 porcelain discs were fabricated according to manufacturer`s instructions, polished through 1 mu m and divided into two groups: GC (control) and GI (submitted to an ion exchange procedure using a KNO(3) paste at 470 degrees C for 15 min). SCG parameters were determined by biaxial flexural strength test in artificial saliva at 37 degrees C using five constant stress rates (n =10). 20 specimens of each group were tested at 1 MPa/s to determine Weibull parameters. The SPT diagram was constructed using the least-squares fit of the strength data versus probability of failure. Results. Mean values of m and sigma(0) (95% confidence interval), n and sigma(f0) (standard deviation) were, respectively: 13.8 (10.1-18.8) and 60.4 (58.5 - 62.2), 24.1 (2.5) and 58.1 (0.01) for GC and 7.4 (5.3 -10.0) and 136.8 (129.1-144.7), 36.7 (7.3) and 127.9 (0.01) for GI. Fracture stresses (MPa) calculated using the SPT diagram for lifetimes of 1 day, 1 year and 10 years (at a 5% failure probability) were, respectively, 31.8, 24.9 and 22.7 for GC and 71.2, 60.6 and 56.9 for GI. Significance. For the porcelain tested, the ion exchange process improved strength and resistance to SCG, however, the material`s reliability decreased. The predicted fracture stress at 5% failure probability for a lifetime of 10 years was also higher for the ion treated group. (C) 009 Academy of Dental Materials. Published by Elsevier Ltd. All rights reserved.

Failure modes of Y-TZP crowns at different cusp inclines

Objectives To compare the reliability of the disto-facial (DF) and mesio-lingual (ML) cusps of an anatomically correct zirconia (Y-TZP) crown system The research hypotheses tested were (1) fatigue reliability and failure mode are similar for the ML and DF cusps, (2) failure mode of one cusp does not affect the failure of the other Methods The average dimensions of a mandibular first molar crown were imported into CAD software, a tooth preparation was modelled by 1 5 mm marginal high reduction of proximal walls and occlusal surface by 2 0 mm The CAD-based tooth preparation was milled and used as a die to fabricate crowns (n = 14) with porcelain veneer on a 0 5 mm Y-TZP core. Crowns were cemented on composite reproductions of the tooth preparation The crowns were step-stress mouth motion fatigued with sliding (0 7 mm) a tungsten-carbide indenter of 6 25 mm diameter down on the inclines of either the DF or ML cusps Use level probability Weibull curve with use stress of 200 N and the reliability for completion of a mission of 50,000 cycles at 200 N load were calculated Results Reliability for a 200 N at 50,000 cycles mission was not different between tested cusps SEM imaging showed large cohesive failures within the veneer for the ML and smaller for the DF Fractures originated from the contact area regardless of the cusp loaded Conclusion No significant difference on fatigue reliability was observed between the DF compared to the ML cusp Fracture of one cusp did not affect the other (c) 2010 Elsevier Ltd All rights reserved

Fatigue and damage accumulation of veneer porcelain pressed on Y-TZP

Objectives: This study compared the reliability and fracture patterns of zirconia cores veneered with pressable porcelain submitted to either axial or off-axis sliding contact fatigue. Methods: Forty-two Y-TZP plates (12 mm x 12 mm x 0.5 mm) veneered with pressable porcelain (12 mm x 12 mm x 1.2 mm) and adhesively luted to water aged composite resin blocks (12 mm x 12 mm x 4 mm) were stored in water at least 7 days prior to testing. Profiles for step-stress fatigue (ratio 3:2:1) were determined from single load to fracture tests (n = 3). Fatigue loading was delivered on specimen either on axial (n = 18) or off-axis 30 degrees angulation (n = 18) to simulate posterior tooth cusp inclination creating a 0.7 mm slide. Single load and fatigue tests utilized a 6.25 mm diameter WC indenter. Specimens were inspected by means of polarized-light microscope and SEM. Use level probability Weibull curves were plotted with 2-sided 90% confidence bounds (CB) and reliability for missions of 50,000 cycles at 200 N (90% CB) were calculated. Results: The calculated Weibull Beta was 3.34 and 2.47 for axial and off-axis groups, respectively, indicating that fatigue accelerated failure in both loading modes. The reliability data for a mission of 50,000 cycles at 200 N load with 90% CB indicates no difference between loading groups. Deep penetrating cone cracks reaching the core-veneer interface were observed in both groups. Partial cones due to the sliding component were observed along with the cone cracking for the off-axis group. No Y-TZP core fractures were observed. Conclusions: Reliability was not significantly different between axial and off-axis mouth-motion fatigued pressed over Y-TZP cores, but incorporation of sliding resulted in more aggressive damage on the veneer. (C) 2009 Elsevier Ltd. All rights reserved.

Reliability of Metalloceramic and Zirconia-based Ceramic Crowns

Despite the increasing utilization of all-ceramic crown systems, their mechanical performance relative to that of metal ceramic restorations (MCR) has yet to be determined. This investigation tested the hypothesis that MCR present higher reliability over two Y-TZP all-ceramic crown systems under mouth-motion fatigue conditions. A CAD-based tooth preparation with the average dimensions of a mandibular first molar was used as a master die to fabricate all restorations. One 0.5-mm Pd-Ag and two Y-TZP system cores were veneered with 1.5 mm porcelain. Crowns were cemented onto aged (60 days in water) composite (Z100, 3M/ESPE) reproductions of the die. Mouth-motion fatigue was performed, and use level probability Weibull curves were determined. Failure modes of all systems included chipping or fracture of the porcelain veneer initiating at the indentation site. Fatigue was an acceleration factor for all-ceramic systems, but not for the MCR system. The latter presented significantly higher reliability under mouth-motion cyclic mechanical testing.

Using non-homogeneous Poisson models with multiple change-points to estimate the number of ozone exceedances in Mexico City

In this paper, we consider some non-homogeneous Poisson models to estimate the probability that an air quality standard is exceeded a given number of times in a time interval of interest. We assume that the number of exceedances occurs according to a non-homogeneous Poisson process (NHPP). This Poisson process has rate function lambda(t), t >= 0, which depends on some parameters that must be estimated. We take into account two cases of rate functions: the Weibull and the Goel-Okumoto. We consider models with and without change-points. When the presence of change-points is assumed, we may have the presence of either one, two or three change-points, depending of the data set. The parameters of the rate functions are estimated using a Gibbs sampling algorithm. Results are applied to ozone data provided by the Mexico City monitoring network. In a first instance, we assume that there are no change-points present. Depending on the adjustment of the model, we assume the presence of either one, two or three change-points. Copyright (C) 2009 John Wiley & Sons, Ltd.

Infant mortality model for lifetime data

In this paper we introduce a parametric model for handling lifetime data where an early lifetime can be related to the infant-mortality failure or to the wear processes but we do not know which risk is responsible for the failure. The maximum likelihood approach and the sampling-based approach are used to get the inferences of interest. Some special cases of the proposed model are studied via Monte Carlo methods for size and power of hypothesis tests. To illustrate the proposed methodology, we introduce an example consisting of a real data set.

A Bayesian Analysis in the Presence of Covariates for Multivariate Survival Data: An example of Application

In this paper, we introduce a Bayesian analysis for survival multivariate data in the presence of a covariate vector and censored observations. Different ""frailties"" or latent variables are considered to capture the correlation among the survival times for the same individual. We assume Weibull or generalized Gamma distributions considering right censored lifetime data. We develop the Bayesian analysis using Markov Chain Monte Carlo (MCMC) methods.

Estimating the number of ozone peaks in Mexico City using a non-homogeneous Poisson model

In this paper, we consider the problem of estimating the number of times an air quality standard is exceeded in a given period of time. A non-homogeneous Poisson model is proposed to analyse this issue. The rate at which the Poisson events occur is given by a rate function lambda(t), t >= 0. This rate function also depends on some parameters that need to be estimated. Two forms of lambda(t), t >= 0 are considered. One of them is of the Weibull form and the other is of the exponentiated-Weibull form. The parameters estimation is made using a Bayesian formulation based on the Gibbs sampling algorithm. The assignation of the prior distributions for the parameters is made in two stages. In the first stage, non-informative prior distributions are considered. Using the information provided by the first stage, more informative prior distributions are used in the second one. The theoretical development is applied to data provided by the monitoring network of Mexico City. The rate function that best fit the data varies according to the region of the city and/or threshold that is considered. In some cases the best fit is the Weibull form and in other cases the best option is the exponentiated-Weibull. Copyright (C) 2007 John Wiley & Sons, Ltd.

The complementary exponential power lifetime model

In this paper we propose a new lifetime distribution which can handle bathtub-shaped unimodal increasing and decreasing hazard rate functions The model has three parameters and generalizes the exponential power distribution proposed by Smith and Bain (1975) with the inclusion of an additional shape parameter The maximum likelihood estimation procedure is discussed A small-scale simulation study examines the performance of the likelihood ratio statistics under small and moderate sized samples Three real datasets Illustrate the methodology (C) 2010 Elsevier B V All rights reserved

A new family of generalized distributions

Kumaraswamy [Generalized probability density-function for double-bounded random-processes, J. Hydrol. 462 (1980), pp. 79-88] introduced a distribution for double-bounded random processes with hydrological applications. For the first time, based on this distribution, we describe a new family of generalized distributions (denoted with the prefix `Kw`) to extend the normal, Weibull, gamma, Gumbel, inverse Gaussian distributions, among several well-known distributions. Some special distributions in the new family such as the Kw-normal, Kw-Weibull, Kw-gamma, Kw-Gumbel and Kw-inverse Gaussian distribution are discussed. We express the ordinary moments of any Kw generalized distribution as linear functions of probability weighted moments (PWMs) of the parent distribution. We also obtain the ordinary moments of order statistics as functions of PWMs of the baseline distribution. We use the method of maximum likelihood to fit the distributions in the new class and illustrate the potentiality of the new model with an application to real data.

A flexible model for survival data with a cure rate: a Bayesian approach

In this paper we deal with a Bayesian analysis for right-censored survival data suitable for populations with a cure rate. We consider a cure rate model based on the negative binomial distribution, encompassing as a special case the promotion time cure model. Bayesian analysis is based on Markov chain Monte Carlo (MCMC) methods. We also present some discussion on model selection and an illustration with a real dataset.