A frailty modeling approach for parental effects in animal breeding

Survival models involving frailties are commonly applied in studies where correlated event time data arise due to natural or artificial clustering. In this paper we present an application of such models in the animal breeding field. Specifically, a mixed survival model with a multivariate correlated frailty term is proposed for the analysis of data from over 3611 Brazilian Nellore cattle. The primary aim is to evaluate parental genetic effects on the trait length in days that their progeny need to gain a commercially specified standard weight gain. This trait is not measured directly but can be estimated from growth data. Results point to the importance of genetic effects and suggest that these models constitute a valuable data analysis tool for beef cattle breeding.

A log-extended Weibull regression model

A bathtub-shaped failure rate function is very useful in survival analysis and reliability studies. The well-known lifetime distributions do not have this property. For the first time, we propose a location-scale regression model based on the logarithm of an extended Weibull distribution which has the ability to deal with bathtub-shaped failure rate functions. We use the method of maximum likelihood to estimate the model parameters and some inferential procedures are presented. We reanalyze a real data set under the new model and the log-modified Weibull regression model. We perform a model check based on martingale-type residuals and generated envelopes and the statistics AIC and BIC to select appropriate models. (C) 2009 Elsevier B.V. All rights reserved.

Formulating Mixed Models for Experiments, Including Longitudinal Experiments

Mixed models have become important in analyzing the results of experiments, particularly those that require more complicated models (e.g., those that involve longitudinal data). This article describes a method for deriving the terms in a mixed model. Our approach extends an earlier method by Brien and Bailey to explicitly identify terms for which autocorrelation and smooth trend arising from longitudinal observations need to be incorporated in the model. At the same time we retain the principle that the model used should include, at least, all the terms that are justified by the randomization. This is done by dividing the factors into sets, called tiers, based on the randomization and determining the crossing and nesting relationships between factors. The method is applied to formulate mixed models for a wide range of examples. We also describe the mixed model analysis of data from a three-phase experiment to investigate the effect of time of refinement on Eucalyptus pulp from four different sources. Cubic smoothing splines are used to describe differences in the trend over time and unstructured covariance matrices between times are found to be necessary.

A generalized modified Weibull distribution for lifetime modeling

A four parameter generalization of the Weibull distribution capable of modeling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone as well as non-monotone failure rates, which are quite common in lifetime problems and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the Weibull, extreme value, exponentiated Weibull, generalized Rayleigh and modified Weibull distributions, among others. We derive two infinite sum representations for its moments. The density of the order statistics is obtained. The method of maximum likelihood is used for estimating the model parameters. Also, the observed information matrix is obtained. Two applications are presented to illustrate the proposed distribution. (C) 2008 Elsevier B.V. All rights reserved.

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.

On estimation and influence diagnostics for zero-inflated negative binomial regression models

The zero-inflated negative binomial model is used to account for overdispersion detected in data that are initially analyzed under the zero-Inflated Poisson model A frequentist analysis a jackknife estimator and a non-parametric bootstrap for parameter estimation of zero-inflated negative binomial regression models are considered In addition an EM-type algorithm is developed for performing maximum likelihood estimation Then the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and some ways to perform global influence analysis are derived In order to study departures from the error assumption as well as the presence of outliers residual analysis based on the standardized Pearson residuals is discussed The relevance of the approach is illustrated with a real data set where It is shown that zero-inflated negative binomial regression models seems to fit the data better than the Poisson counterpart (C) 2010 Elsevier B V All rights reserved

Interval breast cancers in an Australian mammographic screening program

Objective: To determine the incidence of interval cancers which occurred in the first 12 months after mammographic screening at a mammographic screening service. Design: Retrospective analysis of data obtained by crossmatching the screening Service and the New South Wales Central Cancer Registry databases. Setting: The Central & Eastern Sydney Service of BreastScreen NSW. Participants: Women aged 40-69 years at first screen, who attended for their first or second screen between 1 March 1988 and 31 December 1992. Main outcome measures: Interval-cancer rates per 10 000 screens and as a proportion of the underlying incidence of breast cancer (as estimated by the underlying rate in the total NSW population). Results: The 12-month interval-cancer incidence per 10 000 screens was 4.17 for the 40-49 years age group (95% confidence interval [CI], 1.35-9.73) and 4.64 for the 50-69 years age group (95% CI, 2.47-7.94). Proportional incidence rates were 30.1% for the 40-49 years age group (95% CI, 9.8-70.3) and 22% for the 50-69 years age group (95% CI, 11.7-37.7). There was no significant difference between the proportional incidence rate for the 50-69 years age group for the Central & Eastern Sydney Service and those of major successful overseas screening trials. Conclusion: Screening quality was acceptable and should result in a significant mortality reduction in the screened population. Given the small number of cancers involved, comparison of interval-cancer statistics of mammographic screening programs with trials requires age-specific or age-adjusted data, and consideration of confidence intervals of both program and trial data.

Cohort trends in youth suicide in Australia 1964-1997

Objective: This paper examines trends in the rate of suicide among young Australians aged 15-24 years from 1964 to 1997 and presents an age-period-cohort analysis of these trends. Method: Study design consisted of an age-period-cohort analysis of suicide mortality in Australian youth aged between 15 and 24 for the years 1964-1997 inclusive. Data sources were Australian Bureau of Statistics data on: numbers of deaths due to suicide by gender and age at death; and population at risk in each of eight birth cohorts (1940-1944, 1945-1949, 1950-1954, 1955-1959, 1960-1964, 1965-1969, 1970-1974, and 1975-1979). Main outcome measures were population rates of deaths among males and females in each birth cohort attributed to suicide in each year 1964-1997. Results: The rate of suicide deaths among Australian males aged 15-24 years increased from 8.7 per 100 000 in 1964 to 30.9 per 100 000 in 1997, with the rate among females changing little over the period, from 5.2 per 100 000 in 1964 to 7.1 per 100 000 in 1997. While the rate of deaths attributed to suicide increased over the birth cohorts, analyses revealed that these increases were largely due to period effects, with suicide twice as likely among those aged 15-24 years in 1985-1997 than between 1964 and 1969. Conclusions: The rate of youth suicide in Australia has increased since 1964, particularly among males. This increase can largely be attributed to period effects rather than to a cohort effect and has been paralleled by an increased rate of youth suicides internationally and by an increase in other psychosocial problems including psychiatric illness, criminal offending and substance use disorders.

Effects of phonophoresis with Arnica montana onto acute inflammatory process in rat skeletal muscles: An experimental study

This study aimed at verifying the effects of phonophoresis associated with Arnica montana on the acute phase of an inflammatory muscle lesion. Forty Wistar male rats (300 +/- 50 g), of which the Tibialis Anterior muscle was surgically lesioned, were divided into four groups (n = 10 each): control group received no treatment; the ultrasound group (US) was treated in pulsed mode with 1-MHz frequency, 0.5 W/cm(2) intensity (spatial and temporal average - SATA), duty cycle of 1: 2 (2 ms on, 4 ms off, 50%), time of application 3 min per session, one session per day, for 3 days; the phonophoresis or ultrasound plus arnica (US+A) group was treated with arnica with the same US parameters plus arnica gel; and the arnica group (A) was submitted to massage with arnica gel, also for 3 min, once a day, for 3 days. Treatment started 24 h after the surgical lesion. On the 4th day after lesion creation, animals were sacrificed and sections of the lesioned, inflamed muscle were removed for quantitative (mononuclear and polymorphonuclear cell count) and qualitative histological analysis. Collected data from the 4 groups were statistically analyzed and the significance level set at p < 0.05. Results show higher mononuclear cell density in all three treated groups with no significant difference between them, but values were significantly different (p < 0.0001) when compared to control group`s. As to polymorphonuclear cell density, significant differences were found between control group (p = 0.0134) and US, US+A and A groups; the arnica group presented lesser density of polymorphonuclear cells when compared (p = 0.0134) to the other groups. No significant difference was found between US and US+A groups. While the massage with arnica gel proved to be an effective anti-inflammatory on acute muscle lesion in topic use, these results point to ineffectiveness of Arnica montana phonophoresis, US having seemingly checked or minimized its anti-inflammatory effect. (C) 2008 Elsevier B. V. All rights reserved.

Recent trends in pay equity: Beyond the aggregate statistics

The earnings gap between men and women has remained comparatively stable at an aggregate level over the 1990s in Australia. From one perspective, this is a reminder of the considerable difficulty of addressing wage differentials once the most overt forms of wage discrimination have been removed, and of the limited impact of most policy initiatives. From another, it may be seen as evidence that dire predictions about the effects of decentralisation on the earnings gap have failed to materialise. In this paper, I use Australian Bureau of Statistics data to show that a number of different trends are evident underneath the relatively static picture shown by the aggregate statistics, particularly as wage dispersion has increased. The data suggest not only that the prospects for pay equity are far from benign, but also that in the current labour market the issue of gender pay inequality cannot be effectively addressed separately from wage inequality more generally.

Effect of changing data collection parameters on statistical motor unit number estimates

The effect of number of samples and selection of data for analysis on the calculation of surface motor unit potential (SMUP) size in the statistical method of motor unit number estimates (MUNE) was determined in 10 normal subjects and 10 with amyotrophic lateral sclerosis (ALS). We recorded 500 sequential compound muscle action potentials (CMAPs) at three different stable stimulus intensities (10–50% of maximal CMAP). Estimated mean SMUP sizes were calculated using Poisson statistical assumptions from the variance of 500 sequential CMAP obtained at each stimulus intensity. The results with the 500 data points were compared with smaller subsets from the same data set. The results using a range of 50–80% of the 500 data points were compared with the full 500. The effect of restricting analysis to data between 5–20% of the CMAP and to standard deviation limits was also assessed. No differences in mean SMUP size were found with stimulus intensity or use of different ranges of data. Consistency was improved with a greater sample number. Data within 5% of CMAP size gave both increased consistency and reduced mean SMUP size in many subjects, but excluded valid responses present at that stimulus intensity. These changes were more prominent in ALS patients in whom the presence of isolated SMUP responses was a striking difference from normal subjects. Noise, spurious data, and large SMUP limited the Poisson assumptions. When these factors are considered, consistent statistical MUNE can be calculated from a continuous sequence of data points. A 2 to 2.5 SD or 10% window are reasonable methods of limiting data for analysis. Muscle Nerve 27: 320–331, 2003

Guimarães residents’ perceptions towards tourism impacts: a cluster analysis

The present study was designed to assess and segment local residents with respect to their perceived impacts of Guimarães tourism development. The residents of this municipality (located in the northern part of Portugal) are quite strong in their support to tourism. However, they do not keep a homogeneous perception of tourism impacts. A clusters analysis using data from a survey of 400 Guimarães residents’ has revealed the existence of three clusters, according the different degrees of perceived tourism impacts: the Skeptics - moderate in relation to the benefits (averages range from 2.89-3.74) and the ones more concerned with its costs (averages range from 2.86-3.74); the Moderately optimistic - very optimistic about the benefits of tourism (averages range from 3.74-4.51) and conscious of the costs (averages range from 2.71-3.49); the Enthusiasts - very optimistic about tourism benefits (averages range from 2.92-4.52) and little worried about its costs (averages range from 1.78-3.26). Following the data from the survey, the findings are discussed and a few conclusions are extracted.

Sistema de apoio à decisão espacial para análise do impacte ambiental da dispersão de poluentes atmosféricos

Neste artigo é apresentado um Sistema de Apoio à Decisão Espacial (SADE) onde os decisores podem facilmente definir diferentes tipos de problemas espaciais recorrendo a diferentes categorias de objetos, pré-definidas ou a definir, associando- lhes características com ou sem dependência espacial, e indicando formas de interferência (impactos) entre essas caracte- rísticas/propriedades. A análise espacial para determinação ou avaliação de configurações alternativas para a localização de diferentes tipos de ocorrências espaciais será feita através da utilização interativa do SADE de acordo com conjuntos de regras intrínsecas aos vários elementos gráficos (objetos, categorias, características, impactos) utilizados na definição dos problemas. O teste à generalidade representativa e analítica do SADE proposto é efectuado recorrendo a um problema concreto, suficientemente específico e complexo, relativo à aplicação de modelos gaussianos para o estudo da dispersão atmosférica de eventuais poluentes resultantes do tratamento de resíduos sólidos. A região em estudo está limitada, neste exemplo, ao município de Coimbra, Portugal. Para este município estão acessíveis, e são utilizados, os dados demográficos ao nível da secção de voto (censos oficiais) e, como tal, é possível a realização de um estudo realístico do impacto com populações humanas.

Analisador de vibrações de dois canais baseado em computador

O presente trabalho teve como principal objectivo o desenvolvimento de um analisador de vibrações de dois canais baseado em computador, para a realização de diagnóstico no âmbito do controlo de condição de máquinas. Foi desenvolvida uma aplicação num computador comum, no software LabVIEW, que através de transdutores de aceleração do tipo MEMS conectados via USB, faz a recolha de dados de vibração e procede ao seu processamento e apresentação ao utilizador. As ferramentas utilizadas para o processamento de dados são ferramentas comuns encontradas em vários analisadores de vibrações disponíveis no mercado. Estas podem ser: gráficos de espectro de frequência, sinal no tempo, cascata ou valores de nível global de vibração, entre outras. Apesar do analisador desenvolvido não apresentar inovação nas ferramentas de análise adoptadas, este pretende ser distinguido pelo baixo custo, simplicidade e carácter didáctico. Este trabalho vem evidenciar as vantagens, desvantagens e potencialidades de um analisador desta natureza. São tiradas algumas conclusões quanto à sua capacidade de diagnóstico de avarias, capacidades como ferramenta didáctica, sensores utilizados e linguagem de programação escolhida. Como conclusões principais, o trabalho revela que os sensores escolhidos não são os indicados para efectuar o diagnóstico de avarias em ambiente industrial, contudo são ideais para tornar este analisador numa boa ferramenta didáctica e de treino.

Análise de fiabilidade de sistemas e componentes de helicópteros

Este trabalho ocorre face à necessidade da empresa Helisuporte ter uma perspectiva a nível de fiabilidade das suas aeronaves. Para isso, foram traçados como objectivos de estudo a criação de uma base de dados de anomalias; identificação de sistemas e componentes problemáticos; caracterização dos mesmos, avaliar a condição de falha e, com isto, apresentar soluções de controlo de anomalias. Assim, foi desenvolvida uma metodologia que proporciona tratamento de dados com recurso a uma análise não-paramétrica, tendo sido escolhida a estatística de amostra. Esta irá permitir a identificação dos sistemas problemáticos e seus componentes anómalos. Efectuado o tratamento de dados, passamos para a caracterização fiabilística desses componentes, assumindo o tempo de operação e a vida útil específica de cada um. Esta foi possível recorrendo ao cálculo do nível de fiabilidade, MTBF, MTBUR e taxa de avarias. De modo a identificar as diferentes anomalias e caracterizar o “know-how” da equipa de manutenção, implementou-se a análise de condição de falha, mais propriamente a análise dos modos e efeitos de falha. Tendo isso em atenção, foi construído um encadeamento lógico simples, claro e eficaz, face a uma frota complexa. Implementada essa metodologia e analisados os resultados podemos afirmar que os objectivos foram alcançados, concluindo-se que os valores de fiabilidade que caracterizam alguns dos componentes das aeronaves pertencentes à frota em estudo não correspondem ao esperado e idealizado como referência de desempenho dos mesmos. Assim, foram sugeridas alterações no manual de manutenção de forma a melhorar estes índices. Com isto conseguiu-se desenvolver, o que se poderá chamar de, “fiabilidade na óptica do utilizador”.