The McDonald extended distribution: properties and applications

We study a five-parameter lifetime distribution called the McDonald extended exponential model to generalize the exponential, generalized exponential, Kumaraswamy exponential and beta exponential distributions, among others. We obtain explicit expressions for the moments and incomplete moments, quantile and generating functions, mean deviations, Bonferroni and Lorenz curves and Gini concentration index. The method of maximum likelihood and a Bayesian procedure are adopted for estimating the model parameters. The applicability of the new model is illustrated by means of a real data set.

A SIMULATION STUDY OF THE SIZE AND POWER OF SOME TWO-SAMPLE TESTS FOR UNEQUALLY CENSORED SURVIVAL DATA (RANK TESTS, LOGRANK, PROPORTIONAL HAZARDS, WILCOXON TESTS, EXPONENTIAL DISTRIBUTION)

Sizes and power of selected two-sample tests of the equality of survival distributions are compared by simulation for small samples from unequally, randomly-censored exponential distributions. The tests investigated include parametric tests (F, Score, Likelihood, Asymptotic), logrank tests (Mantel, Peto-Peto), and Wilcoxon-Type tests (Gehan, Prentice). Equal sized samples, n = 18, 16, 32 with 1000 (size) and 500 (power) simulation trials, are compared for 16 combinations of the censoring proportions 0%, 20%, 40%, and 60%. For n = 8 and 16, the Asymptotic, Peto-Peto, and Wilcoxon tests perform at nominal 5% size expectations, but the F, Score and Mantel tests exceeded 5% size confidence limits for 1/3 of the censoring combinations. For n = 32, all tests showed proper size, with the Peto-Peto test most conservative in the presence of unequal censoring. Powers of all tests are compared for exponential hazard ratios of 1.4 and 2.0. There is little difference in power characteristics of the tests within the classes of tests considered. The Mantel test showed 90% to 95% power efficiency relative to parametric tests. Wilcoxon-type tests have the lowest relative power but are robust to differential censoring patterns. A modified Peto-Peto test shows power comparable to the Mantel test. For n = 32, a specific Weibull-exponential comparison of crossing survival curves suggests that the relative powers of logrank and Wilcoxon-type tests are dependent on the scale parameter of the Weibull distribution. Wilcoxon-type tests appear more powerful than logrank tests in the case of late-crossing and less powerful for early-crossing survival curves. Guidelines for the appropriate selection of two-sample tests are given. ^

Bayesian and Frequentist Two-Sample Predictions of the Inverse Weibull Model Based on Generalized Order Statistics

2000 Mathematics Subject Classification: 62E16,62F15, 62H12, 62M20.

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.

Fatigue life and failure modes of crowns systems with a modified framework design

Objectives: To evaluate the effect of framework design on the fatigue life and failure modes of metal ceramic (MC, Ni-Cr alloy core, VMK 95 porcelain veneer), glass-infiltrated alumina (ICA, In-Ceram Alumina/VM7), and veneered yttria-stabilized tetragonal zirconia polycrystals (Y-TZP, IPSe.max ZirCAD/IPS e.max,) crowns. Methods: Sixty composite resin tooth replicas of a prepared maxillary first molar were produced to receive crowns systems of a standard (MCs, ICAs, and Y-TZPs, n = 10 each) or a modified framework design (MCm, ICAm, and Y-TZPm, n = 10 each). Fatigue loading was delivered with a spherical steel indenter (3.18 mm radius) on the center of the occlusal surface using r-ratio fatigue (30-300 N) until completion of 10(6) cycles or failure. Fatigue was interrupted every 125,000 cycles for damage evaluation. Weibull distribution fits and contour plots were used for examining differences between groups. Failure mode was evaluated by light polarized and SEM microscopy. Results: Weibull analysis showed the highest fatigue life for MC crowns regardless of framework design. No significant difference (confidence bound overlaps) was observed between ICA and Y-TZP with or without framework design modification. Y-TZPm crowns presented fatigue life in the range of MC crowns. No porcelain veneer fracture was observed in the MC groups, whereas ICAs presented bulk fracture and ICAm failed mainly through the veneer. Y-TZP crowns failed through chipping within the veneer, without core fractures. Conclusions: Framework design modification did not improve the fatigue life of the crown systems investigated. Y-TZPm crowns showed comparable fatigue life to MC groups. Failure mode varied according to crown system. (C) 2010 Elsevier Ltd. All rights reserved.

Organização e gestão da manutencão de ferramentas numa empresa de produção de peças metálicas

Dissertação de mestrado Engenharia e Gestão da Qualidade

Determining the best population-level alcohol consumption model and its impact on estimates of alcohol-attributable harms.

BACKGROUND: The goals of our study are to determine the most appropriate model for alcohol consumption as an exposure for burden of disease, to analyze the effect of the chosen alcohol consumption distribution on the estimation of the alcohol Population- Attributable Fractions (PAFs), and to characterize the chosen alcohol consumption distribution by exploring if there is a global relationship within the distribution. METHODS: To identify the best model, the Log-Normal, Gamma, and Weibull prevalence distributions were examined using data from 41 surveys from Gender, Alcohol and Culture: An International Study (GENACIS) and from the European Comparative Alcohol Study. To assess the effect of these distributions on the estimated alcohol PAFs, we calculated the alcohol PAF for diabetes, breast cancer, and pancreatitis using the three above-named distributions and using the more traditional approach based on categories. The relationship between the mean and the standard deviation from the Gamma distribution was estimated using data from 851 datasets for 66 countries from GENACIS and from the STEPwise approach to Surveillance from the World Health Organization. RESULTS: The Log-Normal distribution provided a poor fit for the survey data, with Gamma and Weibull distributions providing better fits. Additionally, our analyses showed that there were no marked differences for the alcohol PAF estimates based on the Gamma or Weibull distributions compared to PAFs based on categorical alcohol consumption estimates. The standard deviation of the alcohol distribution was highly dependent on the mean, with a unit increase in alcohol consumption associated with a unit increase in the mean of 1.258 (95% CI: 1.223 to 1.293) (R2 = 0.9207) for women and 1.171 (95% CI: 1.144 to 1.197) (R2 = 0. 9474) for men. CONCLUSIONS: Although the Gamma distribution and the Weibull distribution provided similar results, the Gamma distribution is recommended to model alcohol consumption from population surveys due to its fit, flexibility, and the ease with which it can be modified. The results showed that a large degree of variance of the standard deviation of the alcohol consumption Gamma distribution was explained by the mean alcohol consumption, allowing for alcohol consumption to be modeled through a Gamma distribution using only average consumption.

Caracterización selvícola de las masas monoespecíficas de pino laricio en Cataluña

En el presente trabajo se ha abordado, a partir de la información suministrada por el Segundo Inventario Forestal Nacional, la caracterización de las masas monoespecíficas de pino laricio en Cataluña (NE de España), formalizándola en una tipología de enfoque silvogenético. Basándose en los resultados obtenidos del análisis factorial realizado sobre variables dasométricas, los categoremas de la tipología se articularon sobre aspectos concernientes a la distribución diamétrica y espesura del rodal, identificándose nueve tipos. El carácter irregularizado, presente en cinco de ellos aunque con rasgos diferentes, se examinó a través de la distribución de Weibull truncada. La tipología, que se maneja computando tan sólo el nivel de regenerado y los diámetros de los pies inventariables, contempla los distintos estadios evolutivos y permite diagnosticar situaciones comprometidas en cuanto a la persistencia de la masa, bien causadas por deficientes estados de espesura o por falta de regenerado.

Size Effect and Notch Size Effect in Metal Fatigue

It is a well known phenomenon that the constant amplitude fatigue limit of a large component is lower than the fatigue limit of a small specimen made of the same material. In notched components the opposite occurs: the fatigue limit defined as the maximum stress at the notch is higher than that achieved with smooth specimens. These two effects have been taken into account in most design handbooks with the help of experimental formulas or design curves. The basic idea of this study is that the size effect can mainly be explained by the statistical size effect. A component subjected to an alternating load can be assumed to form a sample of initiated cracks at the end of the crack initiation phase. The size of the sample depends on the size of the specimen in question. The main objective of this study is to develop a statistical model for the estimation of this kind of size effect. It was shown that the size of a sample of initiated cracks shall be based on the stressed surface area of the specimen. In case of varying stress distribution, an effective stress area must be calculated. It is based on the decreasing probability of equally sized initiated cracks at lower stress level. If the distribution function of the parent population of cracks is known, the distribution of the maximum crack size in a sample can be defined. This makes it possible to calculate an estimate of the largest expected crack in any sample size. The estimate of the fatigue limit can now be calculated with the help of the linear elastic fracture mechanics. In notched components another source of size effect has to be taken into account. If we think about two specimens which have similar shape, but the size is different, it can be seen that the stress gradient in the smaller specimen is steeper. If there is an initiated crack in both of them, the stress intensity factor at the crack in the larger specimen is higher. The second goal of this thesis is to create a calculation method for this factor which is called the geometric size effect. The proposed method for the calculation of the geometric size effect is also based on the use of the linear elastic fracture mechanics. It is possible to calculate an accurate value of the stress intensity factor in a non linear stress field using weight functions. The calculated stress intensity factor values at the initiated crack can be compared to the corresponding stress intensity factor due to constant stress. The notch size effect is calculated as the ratio of these stress intensity factors. The presented methods were tested against experimental results taken from three German doctoral works. Two candidates for the parent population of initiated cracks were found: the Weibull distribution and the log normal distribution. Both of them can be used successfully for the prediction of the statistical size effect for smooth specimens. In case of notched components the geometric size effect due to the stress gradient shall be combined with the statistical size effect. The proposed method gives good results as long as the notch in question is blunt enough. For very sharp notches, stress concentration factor about 5 or higher, the method does not give sufficient results. It was shown that the plastic portion of the strain becomes quite high at the root of this kind of notches. The use of the linear elastic fracture mechanics becomes therefore questionable.

Thinning effect on plant growth of pruned eucalypt clone

A pruned stand of eucalypt clone underwent five thinning treatments with the removal of different proportion of the planted trees, at different ages: a) 0% - unthinned, b) 35% at 55 months, c) 35% at 81 months, d) 70% at 81 months, removing sprouts in the thinned plant stumps and, e) 70% at 81 months, without coppice sprouts removal. By the age of 141 months, the Weibull distribution showed higher number of trees in the smallest diameter classes for the unthinned treatment. The 70% thinning, with thinned coppice sprouts removal, presented higher number of individuals in the largest diameter classes. Height and yield were the smallest with the removal of 70% of the trees at 81 months, maintaining coppice sprouts. The afterthinning periodic annual increment was greater by thinning 35% of the trees at 55 months resulting in greater number of trees in the largest diameter classes as compared to the other treatments. Yield was higher for the unthinned treatment. The results of this study indicated that thinning 70% of the trees at the age of 81 months, with coppice sprout removal, could be recommended to obtain trees of larger diameter for multiproduct.

Assessment of lettuce quality during storage at low relative humidity using Global Stability Index methodology

During postharvest, lettuce is usually exposed to adverse conditions (e.g. low relative humidity) that reduce the vegetable quality. In order to evaluate its shelf life, a great number of quality attributes must be analyzed, which requires careful experimental design, and it is time consuming. In this study, the modified Global Stability Index method was applied to estimate the quality of butter lettuce at low relative humidity during storage discriminating three lettuce zones (internal, middle, and external). The results indicated that the most relevant attributes were: the external zone - relative water content, water content , ascorbic acid, and total mesophilic counts; middle zone - relative water content, water content, total chlorophyll, and ascorbic acid; internal zone - relative water content, bound water, water content, and total mesophilic counts. A mathematical model that takes into account the Global Stability Index and overall visual quality for each lettuce zone was proposed. Moreover, the Weibull distribution was applied to estimate the maximum vegetable storage time which was 5, 4, and 3 days for the internal, middle, and external zone, respectively. When analyzing the effect of storage time for each lettuce zone, all the indices evaluated in the external zone of lettuce presented significant differences (p < 0.05). For both, internal and middle zones, the attributes presented significant differences (p < 0.05), except for water content and total chlorophyll.

Mathematical modeling of microwave dried celery leaves and determination of the effective moisture diffusivities and activation energy

Celery (Apium graveolens L. var. secalinum Alef) leaves with 50±0.07 g weight and 91.75±0.15% humidity (~11.21 db) were dried using 8 different microwave power densities ranging between 1.8-20 W g-1, until the humidity fell down to 8.95±0.23% (~0.1 db). Microwave drying processes were completed between 5.5 and 77 min depending on the microwave power densities. In this study, measured values were compared with predicted values obtained from twenty thin layer drying theoretical, semi-empirical and empirical equations with a new thin layer drying equation. Within applied microwave power density; models whose coefficient and correlation (R²) values are highest were chosen as the best models. Weibull distribution model gave the most suitable predictions at all power density. At increasing microwave power densities, the effective moisture diffusivity values ranged from 1.595 10-10 to 6.377 10-12 m2 s-1. The activation energy was calculated using an exponential expression based on Arrhenius equation. The linear relationship between the drying rate constant and effective moisture diffusivity gave the best fit.

Évaluation de la modélisation et des prévisions de la vitesse du vent menant à l'estimation de la production d'énergie annuelle d'une turbine éolienne

Suite à un stage avec la compagnie Hatch, nous possédons des jeux de données composés de séries chronologiques de vitesses de vent mesurées à divers sites dans le monde, sur plusieurs années. Les ingénieurs éoliens de la compagnie Hatch utilisent ces jeux de données conjointement aux banques de données d’Environnement Canada pour évaluer le potentiel éolien afin de savoir s’il vaut la peine d’installer des éoliennes à ces endroits. Depuis quelques années, des compagnies offrent des simulations méso-échelle de vitesses de vent, basées sur divers indices environnementaux de l’endroit à évaluer. Les ingénieurs éoliens veulent savoir s’il vaut la peine de payer pour ces données simulées, donc si celles-ci peuvent être utiles lors de l’estimation de la production d’énergie éolienne et si elles pourraient être utilisées lors de la prévision de la vitesse du vent long terme. De plus, comme l’on possède des données mesurées de vitesses de vent, l’on en profitera pour tester à partir de diverses méthodes statistiques différentes étapes de l’estimation de la production d’énergie. L’on verra les méthodes d’extrapolation de la vitesse du vent à la hauteur d’une turbine éolienne et l’on évaluera ces méthodes à l’aide de l’erreur quadratique moyenne. Aussi, on étudiera la modélisation de la vitesse du vent par la distributionWeibull et la variation de la distribution de la vitesse dans le temps. Finalement, l’on verra à partir de la validation croisée et du bootstrap si l’utilisation de données méso-échelle est préférable à celle de données des stations de référence, en plus de tester un modèle où les deux types de données sont utilisées pour prédire la vitesse du vent. Nous testerons la méthodologie globale présentement utilisée par les ingénieurs éoliens pour l’estimation de la production d’énergie d’un point de vue statistique, puis tenterons de proposer des changements à cette méthodologie, qui pourraient améliorer l’estimation de la production d’énergie annuelle.

Characterization of Probability Distributions Using the Residual Entropy Function

The present study on the characterization of probability distributions using the residual entropy function. The concept of entropy is extensively used in literature as a quantitative measure of uncertainty associated with a random phenomenon. The commonly used life time models in reliability Theory are exponential distribution, Pareto distribution, Beta distribution, Weibull distribution and gamma distribution. Several characterization theorems are obtained for the above models using reliability concepts such as failure rate, mean residual life function, vitality function, variance residual life function etc. Most of the works on characterization of distributions in the reliability context centers around the failure rate or the residual life function. The important aspect of interest in the study of entropy is that of locating distributions for which the shannon’s entropy is maximum subject to certain restrictions on the underlying random variable. The geometric vitality function and examine its properties. It is established that the geometric vitality function determines the distribution uniquely. The problem of averaging the residual entropy function is examined, and also the truncated form version of entropies of higher order are defined. In this study it is established that the residual entropy function determines the distribution uniquely and that the constancy of the same is characteristics to the geometric distribution

Exercise sheet 5

Exercises and solutions in LaTex