Missing data mechanisms and their implications on the analysis of categorical data

We review some issues related to the implications of different missing data mechanisms on statistical inference for contingency tables and consider simulation studies to compare the results obtained under such models to those where the units with missing data are disregarded. We confirm that although, in general, analyses under the correct missing at random and missing completely at random models are more efficient even for small sample sizes, there are exceptions where they may not improve the results obtained by ignoring the partially classified data. We show that under the missing not at random (MNAR) model, estimates on the boundary of the parameter space as well as lack of identifiability of the parameters of saturated models may be associated with undesirable asymptotic properties of maximum likelihood estimators and likelihood ratio tests; even in standard cases the bias of the estimators may be low only for very large samples. We also show that the probability of a boundary solution obtained under the correct MNAR model may be large even for large samples and that, consequently, we may not always conclude that a MNAR model is misspecified because the estimate is on the boundary of the parameter space.

Improved testing inference in mixed linear models

Mixed linear models are commonly used in repeated measures studies. They account for the dependence amongst observations obtained from the same experimental unit. Often, the number of observations is small, and it is thus important to use inference strategies that incorporate small sample corrections. In this paper, we develop modified versions of the likelihood ratio test for fixed effects inference in mixed linear models. In particular, we derive a Bartlett correction to such a test, and also to a test obtained from a modified profile likelihood function. Our results generalize those in [Zucker, D.M., Lieberman, O., Manor, O., 2000. Improved small sample inference in the mixed linear model: Bartlett correction and adjusted likelihood. Journal of the Royal Statistical Society B, 62,827-838] by allowing the parameter of interest to be vector-valued. Additionally, our Bartlett corrections allow for random effects nonlinear covariance matrix structure. We report simulation results which show that the proposed tests display superior finite sample behavior relative to the standard likelihood ratio test. An application is also presented and discussed. (C) 2008 Elsevier B.V. All rights reserved.

A new class of survival regression models with heavy-tailed errors: robustness and diagnostics

Birnbaum-Saunders models have largely been applied in material fatigue studies and reliability analyses to relate the total time until failure with some type of cumulative damage. In many problems related to the medical field, such as chronic cardiac diseases and different types of cancer, a cumulative damage caused by several risk factors might cause some degradation that leads to a fatigue process. In these cases, BS models can be suitable for describing the propagation lifetime. However, since the cumulative damage is assumed to be normally distributed in the BS distribution, the parameter estimates from this model can be sensitive to outlying observations. In order to attenuate this influence, we present in this paper BS models, in which a Student-t distribution is assumed to explain the cumulative damage. In particular, we show that the maximum likelihood estimates of the Student-t log-BS models attribute smaller weights to outlying observations, which produce robust parameter estimates. Also, some inferential results are presented. In addition, based on local influence and deviance component and martingale-type residuals, a diagnostics analysis is derived. Finally, a motivating example from the medical field is analyzed using log-BS regression models. Since the parameter estimates appear to be very sensitive to outlying and influential observations, the Student-t log-BS regression model should attenuate such influences. The model checking methodologies developed in this paper are used to compare the fitted models.

Phase transitions in biodegradable films based on blends of gelatin and poly (vinyl alcohol)

The aim of this work was to study the effect of the hydrolysis degree (HD) and the concentration (C PVA) of two types of poly (vinyl alcohol) (PVA) and the effect of the type and the concentration of plasticizers on the phase properties of biodegradable films based on blends of gelatin and PVA, using a response-surface methodology. The films were made by casting and the studied properties were their glass (Tg) and melting (Tm) transition temperatures, which were determined by diferential scanning calorimetry (DSC). For the data obtained on the first scan, the fitting of the linear model was statistically significant and predictive only for the second melting temperature. In this case, the most important effect on the second Tm of the first scan was due to the HD of the PVA. In relation to the second scan, the linear model could be fit to Tg data with only two statistically significant parameters. Both the PVA and plasticizer concentrations had an important effect on Tg. Concerning the second Tm of the second scan, the linear model was fit to data with two statistically significant parameters, namely the HD and the plasticizer concentration. But, the most important effect was provoked by the HD of the PVA.

Tendência temporal e desigualdades na mortalidade por diarreias em menores de 5 anos

OBJETIVO: Analisar a tendência da mortalidade por diarreia entre menores de 5 anos, no município de Osasco (SP), entre 1980 e 2000. MÉTODOS: Trata-se de estudo observacional com dois delineamentos. Um descritivo, que toma o indivíduo como unidade do estudo, e outro ecológico, analisando agregado populacional que incluiu análise de séries temporais. A fonte de dados foi o sistema de informação de mortalidade do Estado de São Paulo e censos de 1980, 1991 e 2000. Descreveu-se a variação sazonal e para a análise de tendência aplicaram-se modelos log lineares de regressão polinomiais, utilizando-se variáveis sociodemográficas da criança e da mãe. Foram analisadas a evolução de indicadores sociodemográficos do município de 1980 a 2000, as taxas médias de mortalidade por diarreia nos menores de 5 anos e seus diferenciais por distrito nos anos 90. RESULTADOS: Dos 1.360 óbitos, 94,3 e 75,3% atingiram, respectivamente, menores de 1 ano e de 6 meses. O declínio da mortalidade foi de 98,3%, com deslocamento da sazonalidade do verão para o outono. A mediana da idade elevou-se de 2 meses nos primeiros períodos para 3 meses no último. O resíduo de óbitos manteve-se entre filhos de mães de 20 a 29 anos e escolaridade < 8 anos. O risco relativo entre o distrito mais atingido e a taxa média do município diminuiu de 3,4 para 1,3 do primeiro para o segundo quinquênio dos anos 90. CONCLUSÃO: Nossos resultados apontam uma elevação da idade mais vulnerável e a provável mudança do agente mais frequentemente associado ao óbito por diarreia.

Prevalência de hipertensão arterial em militares jovens e fatores associados

OBJETIVO: Estimar a prevalência de hipertensão arterial entre militares jovens e fatores associados. MÉTODOS: Estudo transversal realizado com amostra de 380 militares do sexo masculino de 19 e 35 anos de idade em uma unidade da Força Aérea Brasileira em São Paulo, SP, entre 2000 e 2001. Os pontos de corte para hipertensão foram: >140mmHg para pressão sistólica e > 90mmHg para pressão diastólica. As variáveis estudadas incluíram fatores de risco e de proteção para hipertensão, como características comportamentais e nutricionais. Para análise das associações, utilizou-se regressão linear generalizada múltipla, com família binomial e ligação logarítmica, obtendo-se razões de prevalências com intervalo de 90% de confiança e seleção hierarquizada das variáveis. RESULTADOS: A prevalência de hipertensão arterial foi de 22% (IC 90%: 21;29). No modelo final da regressão múltipla verificou-se prevalência de hipertensão 68% maior entre os ex-fumantes em relação aos não fumantes (IC 90%: 1,13;2,50). Entre os indivíduos com sobrepeso (índice de massa corporal - IMC de 25 a 29kg/m2) e com obesidade (IMC>29kg/m2) as prevalências foram, respectivamente, 75% (IC 90%: 1,23;2,50) e 178% (IC 90%: 1,82;4,25) maiores do que entre os eutróficos. Entre os que praticavam atividade física regular, comparado aos que não praticavam, a prevalência foi 52% menor (IC 90%: 0,30;0,90). CONCLUSÕES: Ser ex-fumante e ter sobrepeso ou obesidade foram situações de risco para hipertensão, enquanto que a prática regular de atividade física foi fator de proteção em militares jovens.

Bovine gene polymorphisms related to fat deposition and meat tenderness

Leptin, thyroglobulin and diacylglycerol O-acyltransferase play important roles in fat metabolism. Fat deposition has an influence on meat quality and consumers' choice. The aim of this study was to determine allele and genotype frequencies of polymorphisms of the bovine genes, which encode leptin (LEP), thyroglobulin (TG) and diacylglycerol O-acyltransferase (DGAT1). A further objective was to establish the effects of these polymorphisms on meat characteristics. We genotyped 147 animals belonging to the Nelore (Bos indicus), Canchim (5/8 Bos taurus + 3/8 Bos indicus), Rubia Gallega X Nelore (1/2 Bos taurus + 1/2 Bos indicus), Brangus Three-way cross (9/16 Bos taurus + 7/16 Bos indicus) and Braunvieh Three-way cross (3/4 Bos taurus + 1/4 Bos indicus) breeds. Backfat thickness, total lipids, marbling score, ribeye area and shear force were fitted, using the General Linear Model (GLM) procedure of the SAS software. The least square means of genotypes and genetic groups were compared using Tukey's test. Allele frequencies vary among the genetic groups, depending on Bos indicus versus Bos taurus influence. The LEP polymorphism segregates in pure Bos indicus Nelore animals, which is a new finding. The T allele of TG is fixed in Nelore, and DGAT1 segregates in all groups, but the frequency of allele A is lower in Nelore animals. The results showed no association between the genotypes and traits studied, but a genetic group effect on these traits was found. So, the genetic background remains relevant for fat deposition and meat tenderness, but the gene markers developed for Bos taurus may be insufficient for Bos indicus.

Tendência temporal e desigualdades na mortalidade por diarreias em

Objetivo: Analisar a tendência da mortalidade por diarreia entre menores de 5 anos, no município de Osasco (SP), entre 1980 e 2000. Métodos: Trata-se de estudo observacional com dois delineamentos. Um descritivo, que toma o indivíduo como unidade do estudo, e outro ecológico, analisando agregado populacional que incluiu análise de séries temporais. A fonte de dados foi o sistema de informação de mortalidade do Estado de São Paulo e censos de 1980, 1991 e 2000. Descreveu-se a variação sazonal e para a análise de tendência aplicaram-se modelos log lineares de regressão polinomiais, utilizando-se variáveis sociodemográficas da criança e da mãe. Foram analisadas a evolução de indicadores sociodemográficos do município de 1980 a 2000, as taxas médias de mortalidade por diarreia nos menores de 5 anos e seus diferenciais por distrito nos anos 90. Resultados: Dos 1.360 óbitos, 94,3 e 75,3% atingiram, respectivamente, menores de 1 ano e de 6 meses. O declínio da mortalidade foi de 98,3%, com deslocamento da sazonalidade do verão para o outono. A mediana da idade elevou-se de 2 meses nos primeiros períodos para 3 meses no último. O resíduo de óbitos manteve-se entre filhos de mães de 20 a 29 anos e escolaridade < 8 anos. O risco relativo entre o distrito mais atingido e a taxa média do município diminuiu de 3,4 para 1,3 do primeiro para o segundo quinquênio dos anos 90. Conclusão: Nossos resultados apontam uma elevação da idade mais vulnerável e a provável mudança do agente mais frequentemente associado ao óbito por diarreia

Prevalência de sobrepeso e obesidade em crianças de dois a seis anos matriculadas em escolas particulares no município de São Paulo

OBJETIVO: verificar a prevalência de sobrepeso e obesidade segundo sexo e idade em crianças de 2 a 6 anos de idade, alunos de escolas particulares no município de São Paulo. MÉTODO: foram realizadas medidas de peso e de altura para verificação do estado nutricional de oitocentos e seis crianças de ambos os sexos. Para a classificação do estado nutricional das crianças foram utilizadas as curvas de percentis do Índice de Massa Corporal (IMC = Peso (kg) / Altura² (cm)) para idade, conforme padrão de referência do Multicentre Growth Study, recomendado pela Organização Mundial de Saúde que classifica como sobrepeso valores de percentis > 85 e < 97 e para a obesidade valores < 97. Para análise da relação entre sexo, idade da criança e estado nutricional utilizou-se modelo linear generalizado de regressão múltipla (glm) com ligação logarítmica e família binomial, que permite, diretamente, a estimação das razões de prevalências. A prevalência de sobrepeso+obesidade foi 37,2 por cento para o sexo masculino e 33,4 por cento para o sexo feminino. A razão de prevalência (RP) mostrou que não existe diferença significativa entre obesidade e sobrepeso+obesidade para sexo e idade. CONCLUSÃO: observaram-se prevalências de sobrepeso e de obesidade superiores às prevalências médias da população brasileira. Os resultados encontrados neste estudo reforçam a preocupação com a obesidade infantil que aparentemente vem crescendo, em idades mais precoces como dos pré-escolares

Opportunity for catch-up HPV vaccination in young women after first delivery

Background Early age at first delivery has been identified as a risk factor for high-risk HPV-type infection and cervical cancer development. Methods A cross-sectional study was carried out in a large public maternity hospital in Sao Paulo, Brazil. During June 2006 to February 2007, 301 women aged 15-24 years who gave birth to their first child were recruited between 43 and 60 days after delivery. Detection of HPV DNA in cervical specimens was performed using a standardised PCR protocol with PGMY09/11 primers. The association of selected factors with HPV infection was assessed by using a Generalised Linear Model. Results HPV DNA was detected in 58.5% (95% CI 52.7% to 64.0%) of the enrolled young women. The most common types of HPV found were: HPV16, HPV51, HPV52, HPV58 and HPV71. The overall prevalence of HPV types targeted by the HPV prophylactic vaccines was: HPV 16-12.0%, HPV 18-2.3% and HPV 6 and 11 4.3%. In the multivariate analysis, only age (inversely, p for trend=0.02) and smoking habits were independently associated with HPV infection. Conclusions The findings show that these young primiparous women had high cervical HPV prevalence, suggesting that this is a high-risk group for cervical cancer development. Nevertheless, 17.3% were positive for any of the four HPV types included in HPV vaccines (HPV6, 11, 16 or 18), with 13.3% positive for HPV 16 or 18 and only 1.0% having both vaccine related-oncogenic HPV types. Thus, young primiparous women could benefit from catch-up HPV vaccination programmes.

Coactivation of the shoulder and arm muscles during closed kinetic chain exercises on an unstable surface

Introduction: The purpose of this study was to compare the electromyography index of muscle coactivation of the following muscle pairs: posterior deltoid and pectoralis major (PD/PM); triceps brachii and biceps brachii (TB/BB); and serratus anterior and upper trapezius (SA/UT) during three different closed kinetic chain exercises (wall-press, bench-press and push-up) on an unstable surface at the maximal load. Methods: A total of 20 healthy sedentary men participated in the study. Integral linear values were obtained from three sustained contractions of six seconds each for the three proposed exercises. Mean coactivation index values were compared using the mixed-effects linear model, with a five percent significance level. Results: Electromyography indexes of muscle coactivation showed significant differences for the PD/PM and TB/BB muscle pairs. No differences were found between exercises for the SA/UT muscle pair. Conclusion: Our results seem to differ from those of previous studies, which reported that the similarity in exercises performed is responsible for the comparable muscle activation levels.

Evaluating Patellar Kinematics Through Magnetic Resonance Imaging During Open- and Closed-Kinetic-Chain Exercises

Purpose: To evaluate patellar kinematics of volunteers Without knee pain at rest and during isometric contraction in open- and closed-kinetic-chain exercises. Methods: Twenty individuals took part in this study. All were submitted to magnetic resonance imaging (MRI) during rest and voluntary isometric contraction (VIC) in the open anti closed kinetic chain at 15 degrees, 30 degrees, and 45 degrees of knee flexion. Through MRI and using medical e-film software, the following measurements were evaluated: sulcus angle, patellar-tilt angle, and bisect offset. The mixed-effects linear model was used for comparison between knee positions, between rest and isometric contractions, and between (he exercises. Results: Data analysis revealed that the sulcus angle decreased as knee flexion increased and revealed increases with isometric contractions in both the open and closed kinetic chain for all knee-flexion angles. The patellar-tilt angle decreased with isometric contractions in both the open and closed kinetic chain for every knee position. However, in the closed kinetic chain, patellar tilt increased significantly with the knee flexed at 15 degrees. The bisect offset increased with the knee flexed at 15 degrees during isometric contractions and decreased as knee flexion increased during both exercises. Conclusion: VIC in the last degrees of knee extension may compromise patellar dynamics. On the other hand, it is possible to favor patellar stability by performing muscle contractions with the knee flexed at 30 degrees and 45 degrees in either the open or closed kinetic chain.

Viability of fibroblasts cultured under nutritional stress irradiated with red laser, infrared laser, and red light-emitting diode

Phototherapy is noninvasive, painless and has no known side effect. However, for its incorporation into clinical practice, more well-designed studies are necessary to define optimal parameters for its application. The viability of fibroblasts cultured under nutritional stress irradiated with either a red laser, an infrared laser, or a red light-emitting diode (LED) was analyzed. Irradiation parameters were: red laser (660 nm, 40 mW, 1 W/cm(2)), infrared laser (780 nm, 40 mW, 1 W/cm(2)), and red LED (637 +/- 15 nm, 40 mW, 1 W/cm(2)). All applications were punctual and performed with a spot with 0.4 mm(2) of diameter for 4 or 8 s. The Kruskal-Wallis test and analysis of variance of the general linear model (p <= 0.05) were used for statistical analysis. After 72 h, phototherapy with low-intensity laser and LED showed no toxicity at the cellular level. It even stimulated methylthiazol tetrazolium assay (MTT) conversion and neutral red uptake of fibroblasts cultured under nutritional stress, especially in the group irradiated with infrared laser (p = 0.004 for MTT conversion and p < 0.001 for neutral red uptake). Considering the parameters and protocol of phototherapy used, it can be concluded that phototherapy stimulated the viability of fibroblasts cultured under nutritional deficit resembling those found in traumatized tissue in which cell viability is reduced. (C) 2011 Society of Photo-Optical Instrumentation Engineers (SPIE). [DOI: 10.1117/1.3602850]

Prevalence, magnitude, and methods of rapid weight loss among judo competitors

ARTIOLI, G. G., B. GUALANO, E. FRANCHINI, F. B. SCAGLIUSI, M. TAKESIAN, M. FUCHS, and A. H. LANCHA. Prevalence, Magnitude, and Methods of Rapid Weight Loss among Judo Competitors. Med. Sci. Sports Exerc., Vol. 42, No. 3, pp. 436-442, 2010. Purpose: To identify the prevalence, magnitude, and methods of rapid weight loss among judo competitors. Methods: Athletes (607 males and 215 females; age = 19.3 +/- 5.3 yr, weight = 70 +/- 7.5 kg, height = 170.6 +/- 9.8 cm) completed a previously validated questionnaire developed to evaluate rapid weight loss in judo athletes, which provides a score. The higher the score obtained, the more aggressive the weight loss behaviors. Data were analyzed using descriptive statistics and frequency analyses. Mean scores obtained in the questionnaire were used to compare specific groups of athletes using, when appropriate, Mann-Whitney U-test or general linear model one-way ANOVA followed by Tamhane post hoc test. Results: Eighty-six percent of athletes reported that have already lost weight to compete. When heavyweights are excluded, this percentage rises to 89%. Most athletes reported reductions of up to 5% of body weight (mean +/- SD: 2.5 +/- 2.3%). The most weight ever lost was 2%-5%, whereas a great part of athletes reported reductions of 5%-10% (mean +/- SD: 6 +/- 4%). The number of reductions underwent in a season was 3 +/- 5. The reductions usually occurred within 7 +/- 7 d. Athletes began cutting weight at 12.6 +/- 6.1 yr. No significant differences were found in the score obtained by male versus female athletes as well as by athletes from different weight classes. Elite athletes scored significantly higher in the questionnaire than nonelite. Athletes who began cutting weight earlier also scored higher than those who began later. Conclusions: Rapid weight loss is highly prevalent in judo competitors. The level of aggressiveness in weight management behaviors seems to not be influenced by the gender or by the weight class, but it seems to be influenced by competitive level and by the age at which athletes began cutting weight.

The influence of environmental variables on capybara (Hydrochoerus hydrochaeris: Rodentia, Hydrochoeridae) detectability in anthropogenic environments of southeastern Brazil

Capybaras were monitored weekly from 1998 to 2006 by counting individuals in three anthropogenic environments (mixed agricultural fields, forest and open areas) of southeastern Brazil in order to examine the possible influence of environmental variables (temperature, humidity, wind speed, precipitation and global radiation) on the detectability of this species. There was consistent seasonality in the number of capybaras in the study area, with a specific seasonal pattern in each area. Log-linear models were fitted to the sample counts of adult capybaras separately for each sampled area, with an allowance for monthly effects, time trends and the effects of environmental variables. Log-linear models containing effects for the months of the year and a quartic time trend were highly significant. The effects of environmental variables on sample counts were different in each type of environment. As environmental variables affect capybara detectability, they should be considered in future species survey/monitoring programs.