Post-operative stress hyperglycemia is a predictor of mortality in liver transplantation.

BACKGROUND: A significant association is known between increased glycaemic variability and mortality in critical patients. To ascertain whether glycaemic profiles during the first week after liver transplantation might be associated with long-term mortality in these patients, by analysing whether diabetic status modified this relationship. METHOD: Observational long-term survival study includes 642 subjects undergoing liver transplantation from July 1994 to July 2011. Glucose profiles, units of insulin and all variables with influence on mortality are analysed using joint modelling techniques. RESULTS: Patients registered a survival rate of 85% at 1 year and 65% at 10 years, without differences in mortality between patients with and without diabetes. In glucose profiles, however, differences were observed between patients with and without diabetes: patients with diabetes registered lower baseline glucose values, which gradually rose until reaching a peak on days 2-3 and then subsequently declined, diabetic subjects started from higher values which gradually decreased across the first week. Patients with diabetes showed an association between mortality and age, Model for End-Stage Liver Disease score (MELD) score and hepatitis C virus; among non-diabetic patients, mortality was associated with age, body mass index, malignant aetiology, red blood cell requirements and parenteral nutrition. Glucose profiles were observed to be statistically associated with mortality among patients without diabetes (P = 0.022) but not among patients who presented with diabetes prior to transplantation (P = 0.689). CONCLUSIONS: Glucose profiles during the first week after liver transplantation are different in patients with and without diabetes. While glucose profiles are associated with long-term mortality in patients without diabetes, after adjusting for potential confounding variables such as age, cause of transplantation, MELD, nutrition, immunosuppressive drugs, and units of insulin administered, this does not occur among patients with diabetes.

Bootstrap-based procedures for inference in nonparametric receiver-operating characteristic curve regression analysis.

Prior to using a diagnostic test in a routine clinical setting, the rigorous evaluation of its diagnostic accuracy is essential. The receiver-operating characteristic curve is the measure of accuracy most widely used for continuous diagnostic tests. However, the possible impact of extra information about the patient (or even the environment) on diagnostic accuracy also needs to be assessed. In this paper, we focus on an estimator for the covariate-specific receiver-operating characteristic curve based on direct regression modelling and nonparametric smoothing techniques. This approach defines the class of generalised additive models for the receiver-operating characteristic curve. The main aim of the paper is to offer new inferential procedures for testing the effect of covariates on the conditional receiver-operating characteristic curve within the above-mentioned class. Specifically, two different bootstrap-based tests are suggested to check (a) the possible effect of continuous covariates on the receiver-operating characteristic curve and (b) the presence of factor-by-curve interaction terms. The validity of the proposed bootstrap-based procedures is supported by simulations. To facilitate the application of these new procedures in practice, an R-package, known as npROCRegression, is provided and briefly described. Finally, data derived from a computer-aided diagnostic system for the automatic detection of tumour masses in breast cancer is analysed.

Studying the relationship between a woman's reproductive lifespan and age at menarche using a Bayesian multivariate structured additive distributional regression model.

Studies addressing breast cancer risk factors have been looking at trends relative to age at menarche and menopause. These studies point to a downward trend of age at menarche and an upward trend for age at menopause, meaning an increase of a woman's reproductive lifespan cycle. In addition to studying the effect of the year of birth on the expectation of age at menarche and a woman's reproductive lifespan, it is important to understand how a woman's cohort affects the correlation between these two variables. Since the behavior of age at menarche and menopause may vary with the geographic location of a woman's residence, the spatial effect of the municipality where a woman resides needs to be considered. Thus, a Bayesian multivariate structured additive distributional regression model is proposed in order to analyze how a woman's municipality and year of birth affects a woman's age of menarche, her lifespan cycle, and the correlation of the two. The data consists of 212,517 postmenopausal women, born between 1920 and 1965, who attended the breast cancer screening program in the central region of Portugal.

Glycemic Variability and Its Association With Demographics and Lifestyles in a General Adult Population.

OBJECTIVE: The objective was to investigate glycemic variability indices in relation to demographic factors and common environmental lifestyles in a general adult population. METHODS: The A Estrada Glycation and Inflammation Study is a cross-sectional study covering 1516 participants selected by sampling of the population aged 18 years and over. A subsample of 622 individuals participated in the Glycation project, which included continuous glucose monitoring procedures. Five glycemic variability indices were analyzed, that is, SD, MAGE, MAG, CONGA1, and MODD. RESULTS: Participants had a mean age of 48 years, 62% were females, and 12% had been previously diagnosed with diabetes. In the population without diabetes, index distributions were not normal but skewed to the right. Distributional regression models that adjusted for age, gender, BMI, alcohol intake, smoking status, and physical activity confirmed that all indices were positively and independently associated with fasting glucose levels and negatively with heavy drinking. SD, MAGE, and CONGA1 were positively associated with aging, and MAG was negatively associated with BMI. None of the GVI studied were influenced by physical activity. Age-group-specific reference values are given for the indices. CONCLUSIONS: This study yielded age-specific reference values for glucose variability indices in a general adult population. Significant increases were observed with aging. Heavy drinking of more than 140 g/week was associated with significant decreases in variability indices. No differences were found between males and females. These normative ranges provide a guide for clinical care, and may offer an alternative treatment target among persons with diabetes.

Confidence intervals for the symmetry point: an optimal cutpoint in continuous diagnostic tests.

Continuous diagnostic tests are often used for discriminating between healthy and diseased populations. For this reason, it is useful to select an appropriate discrimination threshold. There are several optimality criteria: the North-West corner, the Youden index, the concordance probability and the symmetry point, among others. In this paper, we focus on the symmetry point that maximizes simultaneously the two types of correct classifications. We construct confidence intervals for this optimal cutpoint and its associated specificity and sensitivity indexes using two approaches: one based on the generalized pivotal quantity and the other on empirical likelihood. We perform a simulation study to check the practical behaviour of both methods and illustrate their use by means of three real biomedical datasets on melanoma, prostate cancer and coronary artery disease.

Structured additive regression modeling of age of menarche and menopause in a breast cancer screening program.

Breast cancer risk is believed to be associated with several reproductive factors, such as early menarche and late menopause. This study is based on the registries of the first time a woman enters the screening program, and presents a spatio-temporal analysis of the variables age of menarche and age of menopause along with other reproductive and socioeconomic factors. The database was provided by the Portuguese Cancer League (LPCC), a private nonprofit organization dealing with multiple issues related to oncology of which the Breast Cancer Screening Program is one of its main activities. The registry consists of 259,652 records of women who entered the screening program for the first time between 1990 and 2007 (45-69-year age group). Structured Additive Regression (STAR) models were used to explore spatial and temporal correlations with a wide range of covariates. These models are flexible enough to deal with a variety of complex datasets, allowing us to reveal possible relationships among the variables considered in this study. The analysis shows that early menarche occurs in younger women and in municipalities located in the interior of central Portugal. Women living in inland municipalities register later ages for menopause, and those born in central Portugal after 1933 show a decreasing trend in the age of menopause. Younger ages of menarche and late menopause are observed in municipalities with a higher purchasing power index. The analysis performed in this study portrays the time evolution of the age of menarche and age of menopause and their spatial characterization, adding to the identification of factors that could be of the utmost importance in future breast cancer incidence research.

smoothHR: an R package for pointwise nonparametric estimation of hazard ratio curves of continuous predictors.

The Cox proportional hazards regression model has become the traditional choice for modeling survival data in medical studies. To introduce flexibility into the Cox model, several smoothing methods may be applied, and approaches based on splines are the most frequently considered in this context. To better understand the effects that each continuous covariate has on the outcome, results can be expressed in terms of splines-based hazard ratio (HR) curves, taking a specific covariate value as reference. Despite the potential advantages of using spline smoothing methods in survival analysis, there is currently no analytical method in the R software to choose the optimal degrees of freedom in multivariable Cox models (with two or more nonlinear covariate effects). This paper describes an R package, called smoothHR, that allows the computation of pointwise estimates of the HRs--and their corresponding confidence limits--of continuous predictors introduced nonlinearly. In addition the package provides functions for choosing automatically the degrees of freedom in multivariable Cox models. The package is available from the R homepage. We illustrate the use of the key functions of the smoothHR package using data from a study on breast cancer and data on acute coronary syndrome, from Galicia, Spain.

Selecting variables in non-parametric regression models for binary response. An application to the computerized detection of breast cancer.

In many biomedical applications, interest lies in being able to distinguish between two possible states of a given response variable, depending on the values of certain continuous predictors. If the number of predictors, p, is high, or if there is redundancy among them, it then becomes important to decide on the selection of the best subset of predictors that will be able to obtain the models with greatest discrimination capacity. With this aim in mind, logistic generalized additive models were considered and receiver operating characteristic (ROC) curves were applied in order to determine and compare the discriminatory capacity of such models. This study sought to develop bootstrap-based tests that allow for the following to be ascertained: (a) the optimal number q < or = p of predictors; and (b) the model or models including q predictors, which display the largest AUC (area under the ROC curve). A simulation study was conducted to verify the behaviour of these tests. Finally, the proposed method was applied to a computer-aided diagnostic system dedicated to early detection of breast cancer.

Multi-state models for the analysis of time-to-event data.

The experience of a patient in a survival study may be modelled as a process with two states and one possible transition from an "alive" state to a "dead" state. In some studies, however, the "alive" state may be partitioned into two or more intermediate (transient) states, each of which corresponding to a particular stage of the illness. In such studies, multi-state models can be used to model the movement of patients among the various states. In these models issues, of interest include the estimation of progression rates, assessing the effects of individual risk factors, survival rates or prognostic forecasting. In this article, we review modelling approaches for multi-state models, and we focus on the estimation of quantities such as the transition probabilities and survival probabilities. Differences between these approaches are discussed, focussing on possible advantages and disadvantages for each method. We also review the existing software currently available to fit the various models and present new software developed in the form of an R library to analyse such models. Different approaches and software are illustrated using data from the Stanford heart transplant study and data from a study on breast cancer conducted in Galicia, Spain.

Categorical variables, interactions and generalized additive models. Applications in computer-aided diagnosis systems.

Recently, the generalized additive models (GAMs) have been presented as a novel statistical approach to distinguish lesion/non-lesion in computer-aided diagnosis (CAD) systems. In this paper, we present an extension of the GAM that allows for the introduction of factors and their interactions with continuous variables, for reducing false positives in a CAD system for detecting clustered microcalcifications in digital mammograms. The results obtained have shown an increase in the sensitivity from 83.12% to 85.71%, while the false positive rate was drastically reduced from 1.46 to 0.74 false detections per image.