Latent residual analysis in binary regression with skewed link

Model diagnostics is an integral part of model determination and an important part of the model diagnostics is residual analysis. We adapt and implement residuals considered in the literature for the probit, logistic and skew-probit links under binary regression. New latent residuals for the skew-probit link are proposed here. We have detected the presence of outliers using the residuals proposed here for different models in a simulated dataset and a real medical dataset.

Retrospective Clinical Study of 988 Dual Acid-Etched Implants Placed in Grafted and Native Bone for Single-Tooth Replacement

Purpose: To evaluate the influence of sex, implant characteristics, and bone grafting on the survival rate of dual acid-etched (DAE) implants. Materials and Methods: Patients treated with internal-hex DAE implants for single-tooth replacement in a military dental clinic between January 2005 and December 2010 were included in this study. Clinical data related to implant characteristics, implant location, presence of grafted bone, and implant failures were collected. The primary outcome was implant loss. The survival rate was analyzed using the Kaplan-Meier method. Cox regression modeling was used to determine which factors would predict implant failure. Results: DAE implants were evaluated in a total of 988 patients (80.3% men). Twenty-four (2.4%) implants failed, most were cylindric (54.2%) with regular platforms (70.8%) and were 10 mm long (58.3%). The failure rate was 2.4% for the anterior maxilla, 3.3% for the posterior maxilla, 1.6% for the anterior mandible, and 2.0% for posterior mandible. The cumulative survival rate was 97.6%. The failure rate was 8.8% in implants placed after sinus augmentation, 7.3% in bone block-grafted areas, and 1.6% in native bone. Based on multivariable analysis (Cox regression), sinus augmentation and bone block grafting had a statistically significant effect on implant failure; the hazard ratios were 5.5 and 4.6, respectively. Conclusion: The results revealed that DAE implants had high survival rates, and no influence of sex, location, shape, diameter, or length on failure rates could be observed. However, a significant association was observed between failure and presence of bone graft in the implant area. Int J Oral Maxillofac Implants 2012;27:1243-1248

A log-Birnbaum-Saunders regression model with asymmetric errors

This paper introduces a skewed log-Birnbaum-Saunders regression model based on the skewed sinh-normal distribution proposed by Leiva et al. [A skewed sinh-normal distribution and its properties and application to air pollution, Comm. Statist. Theory Methods 39 (2010), pp. 426-443]. Some influence methods, such as the local influence and generalized leverage, are presented. Additionally, we derived the normal curvatures of local influence under some perturbation schemes. An empirical application to a real data set is presented in order to illustrate the usefulness of the proposed model.

Covariates of high-risk human papillomavirus (HPV) infections are distinct for incident CIN1, CIN2 and CIN3 as disclosed by competing-risks regression models

Background: In addition to the oncogenic human papillomavirus (HPV), several cofactors are needed in cervical carcinogenesis, but whether the HPV covariates associated with incident i) CIN1 are different from those of incident ii) CIN2 and iii) CIN3 needs further assessment. Objectives: To gain further insights into the true biological differences between CIN1, CIN2 and CIN3, we assessed HPV covariates associated with incident CIN1, CIN2, and CIN3. Study Design and Methods: HPV covariates associated with progression to CIN1, CIN2 and CIN3 were analysed in the combined cohort of the NIS (n = 3,187) and LAMS study (n = 12,114), using competing-risks regression models (in panel data) for baseline HR-HPV-positive women (n = 1,105), who represent a sub-cohort of all 1,865 women prospectively followed-up in these two studies. Results: Altogether, 90 (4.8%), 39 (2.1%) and 14 (1.4%) cases progressed to CIN1, CIN2, and CIN3, respectively. Among these baseline HR-HPV-positive women, the risk profiles of incident GIN I, CIN2 and CIN3 were unique in that completely different HPV covariates were associated with progression to CIN1, CIN2 and CIN3, irrespective which categories (non-progression, CIN1, CIN2, CIN3 or all) were used as competing-risks events in univariate and multivariate models. Conclusions: These data confirm our previous analysis based on multinomial regression models implicating that distinct covariates of HR-HPV are associated with progression to CIN1, CIN2 and CIN3. This emphasises true biological differences between the three grades of GIN, which revisits the concept of combining CIN2 with CIN3 or with CIN1 in histological classification or used as a common end-point, e.g., in HPV vaccine trials.

The log-exponentiated generalized gamma regression model for censored data

For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution [Cordeiro et al. The exponentiated generalized gamma distribution with application to lifetime data, J. Statist. Comput. Simul. 81 (2011), pp. 827-842.] that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub- and unimodal-shaped hazard functions. A second advantage is that it includes classical distributions reported in the lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful by means of two applications to real data.

AGE INFLUENCE ON THE HEART RATE BEHAVIOR ON THE REST-EXERCICIO TRANSITION: AN ANALYSIS BY DELTAS AND LINEAR REGRESSION

Background: Changes in heart rate during rest-exercise transition can be characterized by the application of mathematical calculations, such as deltas 0-10 and 0-30 seconds to infer on the parasympathetic nervous system and linear regression and delta applied to data range from 60 to 240 seconds to infer on the sympathetic nervous system. The objective of this study was to test the hypothesis that young and middle-aged subjects have different heart rate responses in exercise of moderate and intense intensity, with different mathematical calculations. Methods: Seven middle-aged men and ten young men apparently healthy were subject to constant load tests (intense and moderate) in cycle ergometer. The heart rate data were submitted to analysis of deltas (0-10, 0-30 and 60-240 seconds) and simple linear regression (60-240 seconds). The parameters obtained from simple linear regression analysis were: intercept and slope angle. We used the Shapiro-Wilk test to check the distribution of data and the "t" test for unpaired comparisons between groups. The level of statistical significance was 5%. Results: The value of the intercept and delta 0-10 seconds was lower in middle age in two loads tested and the inclination angle was lower in moderate exercise in middle age. Conclusion: The young subjects present greater magnitude of vagal withdrawal in the initial stage of the HR response during constant load exercise and higher speed of adjustment of sympathetic response in moderate exercise.

Survival Analysis of Patients with Heart Failure: Implications of Time-Varying Regression Effects in Modeling Mortality

Background: Several models have been designed to predict survival of patients with heart failure. These, while available and widely used for both stratifying and deciding upon different treatment options on the individual level, have several limitations. Specifically, some clinical variables that may influence prognosis may have an influence that change over time. Statistical models that include such characteristic may help in evaluating prognosis. The aim of the present study was to analyze and quantify the impact of modeling heart failure survival allowing for covariates with time-varying effects known to be independent predictors of overall mortality in this clinical setting. Methodology: Survival data from an inception cohort of five hundred patients diagnosed with heart failure functional class III and IV between 2002 and 2004 and followed-up to 2006 were analyzed by using the proportional hazards Cox model and variations of the Cox's model and also of the Aalen's additive model. Principal Findings: One-hundred and eighty eight (188) patients died during follow-up. For patients under study, age, serum sodium, hemoglobin, serum creatinine, and left ventricular ejection fraction were significantly associated with mortality. Evidence of time-varying effect was suggested for the last three. Both high hemoglobin and high LV ejection fraction were associated with a reduced risk of dying with a stronger initial effect. High creatinine, associated with an increased risk of dying, also presented an initial stronger effect. The impact of age and sodium were constant over time. Conclusions: The current study points to the importance of evaluating covariates with time-varying effects in heart failure models. The analysis performed suggests that variations of Cox and Aalen models constitute a valuable tool for identifying these variables. The implementation of covariates with time-varying effects into heart failure prognostication models may reduce bias and increase the specificity of such models.

The negative binomial-beta Weibull regression model to predict the cure of prostate cancer

In this article, for the first time, we propose the negative binomial-beta Weibull (BW) regression model for studying the recurrence of prostate cancer and to predict the cure fraction for patients with clinically localized prostate cancer treated by open radical prostatectomy. The cure model considers that a fraction of the survivors are cured of the disease. The survival function for the population of patients can be modeled by a cure parametric model using the BW distribution. We derive an explicit expansion for the moments of the recurrence time distribution for the uncured individuals. The proposed distribution can be used to model survival data when the hazard rate function is increasing, decreasing, unimodal and bathtub shaped. Another advantage is that the proposed model includes as special sub-models some of the well-known cure rate models discussed in the literature. We derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes. We analyze a real data set for localized prostate cancer patients after open radical prostatectomy.

A log-linear regression model for the beta-Birnbaum-Saunders distribution with censored data

The beta-Birnbaum-Saunders (Cordeiro and Lemonte, 2011) and Birnbaum-Saunders (Birnbaum and Saunders, 1969a) distributions have been used quite effectively to model failure times for materials subject to fatigue and lifetime data. We define the log-beta-Birnbaum-Saunders distribution by the logarithm of the beta-Birnbaum-Saunders distribution. Explicit expressions for its generating function and moments are derived. We propose a new log-beta-Birnbaum-Saunders regression model that can be applied to censored data and be used more effectively in survival analysis. We obtain the maximum likelihood estimates of the model parameters for censored data and investigate influence diagnostics. The new location-scale regression model is modified for the possibility that long-term survivors may be presented in the data. Its usefulness is illustrated by means of two real data sets. (C) 2011 Elsevier B.V. All rights reserved.

Preditores da realização de consultas odontológicas de rotina e por problema em pré-escolares

OBJETIVO: Estimar a prevalência do uso de serviços odontológicos por pré- escolares e fatores associados. MÉTODOS: Estudo transversal com 1.129 crianças de cinco anos de idade da Coorte de Nascimentos de Pelotas 2004, RS, de setembro de 2009 a janeiro de 2010. Registrou-se o uso de serviço odontológico pelo menos uma vez na vida e o motivo para a primeira consulta odontológica da criança. As categorias do desfecho foram: ter feito a primeira consulta por rotina, para resolver um problema ou nunca ter ido ao dentista. Os exames bucais e as entrevistas foram realizados nos domicílios. Aspectos socioeconômicos e variáveis independentes ligadas à mãe e à criança foram analisados por meio de regressão logística multinomial. RESULTADOS: A prevalência de uso por qualquer motivo foi 37,0%. Os principais preditores para consulta de rotina foram nível econômico mais elevado, mãe com maior escolaridade e ter recebido orientação sobre prevenção. Principais preditores para consulta por problema foram ter sentido dor nos últimos seis meses, mãe com maior escolaridade e ter recebido orientação sobre prevenção. Cerca de 45,0% das mães receberam orientação de como prevenir cárie, principalmente fornecida por dentistas. Filhos de mães com história de maior aderência a programas de saúde tiveram maior probabilidade de ter feito uma consulta odontológica de rotina. CONCLUSÕES: A taxa de utilização dos serviços odontológicos por pré- escolares foi inferior às de consultas médicas (puericultura). Além da renda e da escolaridade, comportamentos maternos têm papel importante no uso por rotina. Relato de dor nos últimos seis meses e número elevado de dentes afetados por cárie, independentemente dos demais fatores, estiveram associados ao uso para resolver problema. É necessária a integração de ações de saúde bucal nos programas materno-infantis.

Competing regression models for longitudinal data

The choice of an appropriate family of linear models for the analysis of longitudinal data is often a matter of concern for practitioners. To attenuate such difficulties, we discuss some issues that emerge when analyzing this type of data via a practical example involving pretestposttest longitudinal data. In particular, we consider log-normal linear mixed models (LNLMM), generalized linear mixed models (GLMM), and models based on generalized estimating equations (GEE). We show how some special features of the data, like a nonconstant coefficient of variation, may be handled in the three approaches and evaluate their performance with respect to the magnitude of standard errors of interpretable and comparable parameters. We also show how different diagnostic tools may be employed to identify outliers and comment on available software. We conclude by noting that the results are similar, but that GEE-based models may be preferable when the goal is to compare the marginal expected responses.

Risk Factors and Survival Impact of Primary Graft Dysfunction After Lung Transplantation in a Single Institution

Background. Lung transplantation has become a standard procedure for some end-stage lung diseases, but primary graft dysfunction (PGD) is an inherent problem that impacts early and late outcomes. The aim of this study was to define the incidence, risk factors, and impact of mechanical ventilation time on mortality rates among a retrospective cohort of lung transplantations performed in a single institution. Methods. We performed a retrospective study of 118 lung transplantations performed between January 2003 and July 2010. The most severe form of PGD (grade III) as defined at 48 and 72 hours was examined for risk factors by multivariable logistic regression models using donor, recipient, and transplant variables. Results. The overall incidence of PGD at 48 hours was 19.8%, and 15.4% at 72 hours. According multivariate analysis, risk factors associated with PGD were donor smoking history for 48 hours (adjusted odds ratio [OR], 4.83; 95% confidence interval [CI], 1.236-18.896; P = .022) and older donors for 72 hours (adjusted OR, 1.046; 95% CI, 0.997-1.098; P = .022). The operative mortality was 52.9% among patients with PGD versus 20.3% at 48 hours (P = .012). At 72 hours, the mortality rate was 58.3% versus 21.2% (P = .013). The 90-days mortality was also higher among patients with PGD. The mechanical ventilation time was longer in patients with PGD III at 48 hours namely, a mean time of 72 versus 24 hours (P = .001). When PGD was defined at 72 hours, the mean ventilation time was even longer, namely 151 versus 24 hours (P < .001). The mean overall survival for patients who developed PGD at 48 hours was 490.9 versus 1665.5 days for subjects without PGD (P = .001). Considering PGD only at 72 hours, the mean survival was 177.7 days for the PGD group and 1628.9 days for the other patients (P < .001). Conclusion. PGD showed an important impacts on operative and 90-day mortality rates, mechanical ventilation time, and overall survival among lung transplant patients. PGD at 72 hours was a better predictor of lung transplant outcomes than at 48 hours. The use of donors with a smoking history or of advanced age were risk factors for the development of PGD.

Accuracy of positron emission tomography/computed tomography and clinical assessment in the detection of complete rectal tumor regression after neoadjuvant chemoradiation Long-Term Results of a Prospective Trial (National Clinical Trial 00254683)

BACKGROUND: Neoadjuvant chemoradiation (CRT) therapy may result in significant tumor regression in patients with rectal cancer. Patients who develop complete tumor regression have been managed by treatment strategies that are alternatives to standard total mesorectal excision. Therefore, assessment of tumor response with positron emission tomography/computed tomography (PET/CT) after neoadjuvant treatment may offer relevant information for the selection of patients to receive alternative treatment strategies. METHODS: Patients with clinical T2 (cT2) through cT4NxM0 rectal adenocarcinoma were included prospectively. Neoadjuvant therapy consisted of 54 grays of radiation and 5-fluorouracil-based chemotherapy. Baseline PET/CT studies were obtained before CRT followed by PET/CT studies at 6 weeks and 12 weeks after the completion of CRT. Clinical assessment was performed at 12 weeks after CRT completion. PET/CT results were compared with clinical and pathologic data. RESULTS: In total, 99 patients were included in the study. Twenty-three patients were complete responders (16 had a complete clinical response, and 7 had a complete pathologic response). The PET/CT response evaluation at 12 weeks indicated that 18 patients had a complete response, and 81 patients had an incomplete response. There were 5 false-negative and 10 false-positive PET/CT results. PET/CT for the detection of residual cancer had 93% sensitivity, 53% specificity, a 73% negative predictive value, an 87% positive predictive value, and 85% accuracy. Clinical assessment alone resulted in an accuracy of 91%. PET/CT information may have detected misdiagnoses made by clinical assessment alone, improving overall accuracy to 96%. CONCLUSIONS: Assessment of tumor response at 12 weeks after CRT completion with PET/CT imaging may provide a useful additional tool with good overall accuracy for the selection of patients who may avoid unnecessary radical resection after achieving a complete clinical response. Cancer 2012;35013511. (C) 2011 American Cancer Society.

Are mental health problems and depression associated with bruxism in children?

Renner AC, da Silva AAM, Rodriguez JDM, Simoes VMF, Barbieri MA, Bettiol H, Thomaz EBAF, Saraiva MC. Are mental health problems and depression associated with bruxism in children? Community Dent Oral Epidemiol 2011. (C) 2011 John Wiley & Sons A/S Abstract Objectives: Previous studies have found an association between bruxism and emotional and behavioral problems in children, but reported data are inconsistent. The objective of this study was to estimate the prevalence of bruxism, and of its components clenching and grinding, and its associations with mental problems and depression. Methods: Data from two Brazilian birth cohorts were analyzed: one from 869 children in Ribeirao Preto RP (Sao Paulo), a more developed city, and the other from 805 children in Sao Luis SL (Maranhao). Current bruxism evaluated by means of a questionnaire applied to the parents/persons responsible for the children was defined when the habit of tooth clenching during daytime and/or tooth grinding at night still persisted until the time of the assessment. Additionally, the lifetime prevalence of clenching during daytime only and grinding at night only was also evaluated. Mental health problems were investigated using the Strength and Difficulties Questionnaire (SDQ) and depression using the Childrens Depression Inventory (CDI). Analyses were carried out for each city: with the SDQ subscales (emotional symptoms, conduct problems, peer problems, attention/hyperactivity disorder), with the total score (sum of the subscales), and with the CDI. These analyses were performed considering different response variables: bruxism, clenching only, and grinding only. The risks were estimated using a Poisson regression model. Statistical inferences were based on 95% confidence intervals (95% CI). Results: There was a high prevalence of current bruxism: 28.7% in RP and 30.0% in SL. The prevalence of clenching was 20.3% in RP and 18.8% in SL, and grinding was found in 35.7% of the children in RP and 39.1% in SL. Multivariable analysis showed a significant association of bruxism with emotional symptoms and total SDQ score in both cities. When analyzed separately, teeth clenching was associated with emotional symptoms, peer problems, and total SDQ score; grinding was significantly associated with emotional symptoms and total SDQ score in RP and SL. Female sex appeared as a protective factor for bruxism, and for clenching and grinding in RP. Furthermore, maternal employment outside the home and white skin color of children were associated with increased prevalence of teeth clenching in SL. Conclusions: Mental health problems were associated with bruxism, with teeth clenching only and grinding at night only. No association was detected between depression and bruxism, neither clenching nor grinding. But it is necessary to be cautious regarding the inferences from some of our results.

A general class of zero-or-one inflated beta regression models

This paper proposes a general class of regression models for continuous proportions when the data contain zeros or ones. The proposed class of models assumes that the response variable has a mixed continuous-discrete distribution with probability mass at zero or one. The beta distribution is used to describe the continuous component of the model, since its density has a wide range of different shapes depending on the values of the two parameters that index the distribution. We use a suitable parameterization of the beta law in terms of its mean and a precision parameter. The parameters of the mixture distribution are modeled as functions of regression parameters. We provide inference, diagnostic, and model selection tools for this class of models. A practical application that employs real data is presented. (C) 2011 Elsevier B.V. All rights reserved.