Semiparametric Normal Transformation Models for Spatially Correlated Survival Data

There is an emerging interest in modeling spatially correlated survival data in biomedical and epidemiological studies. In this paper, we propose a new class of semiparametric normal transformation models for right censored spatially correlated survival data. This class of models assumes that survival outcomes marginally follow a Cox proportional hazard model with unspecified baseline hazard, and their joint distribution is obtained by transforming survival outcomes to normal random variables, whose joint distribution is assumed to be multivariate normal with a spatial correlation structure. A key feature of the class of semiparametric normal transformation models is that it provides a rich class of spatial survival models where regression coefficients have population average interpretation and the spatial dependence of survival times is conveniently modeled using the transformed variables by flexible normal random fields. We study the relationship of the spatial correlation structure of the transformed normal variables and the dependence measures of the original survival times. Direct nonparametric maximum likelihood estimation in such models is practically prohibited due to the high dimensional intractable integration of the likelihood function and the infinite dimensional nuisance baseline hazard parameter. We hence develop a class of spatial semiparametric estimating equations, which conveniently estimate the population-level regression coefficients and the dependence parameters simultaneously. We study the asymptotic properties of the proposed estimators, and show that they are consistent and asymptotically normal. The proposed method is illustrated with an analysis of data from the East Boston Ashma Study and its performance is evaluated using simulations.

Checking Assumptions in Latent Class Regression Models via a Markov Chain Monte Carlo Estimation Approach: An Application to Depression and Socio-Economic Status

Latent class regression models are useful tools for assessing associations between covariates and latent variables. However, evaluation of key model assumptions cannot be performed using methods from standard regression models due to the unobserved nature of latent outcome variables. This paper presents graphical diagnostic tools to evaluate whether or not latent class regression models adhere to standard assumptions of the model: conditional independence and non-differential measurement. An integral part of these methods is the use of a Markov Chain Monte Carlo estimation procedure. Unlike standard maximum likelihood implementations for latent class regression model estimation, the MCMC approach allows us to calculate posterior distributions and point estimates of any functions of parameters. It is this convenience that allows us to provide the diagnostic methods that we introduce. As a motivating example we present an analysis focusing on the association between depression and socioeconomic status, using data from the Epidemiologic Catchment Area study. We consider a latent class regression analysis investigating the association between depression and socioeconomic status measures, where the latent variable depression is regressed on education and income indicators, in addition to age, gender, and marital status variables. While the fitted latent class regression model yields interesting results, the model parameters are found to be invalid due to the violation of model assumptions. The violation of these assumptions is clearly identified by the presented diagnostic plots. These methods can be applied to standard latent class and latent class regression models, and the general principle can be extended to evaluate model assumptions in other types of models.

Body Mass Index, Cigarette Smoking, and Alcohol Consumption and Cancers of the Oral Cavity, Pharynx, and Larynx: Modeling Odds Ratios in Pooled Case-Control Data

study-specific results, their findings should be interpreted with caution

Evolutionary dynamics on rugged fitness landscapes: Exact dynamics and information theoretical aspects

The parallel mutation-selection evolutionary dynamics, in which mutation and replication are independent events, is solved exactly in the case that the Malthusian fitnesses associated to the genomes are described by the random energy model (REM) and by a ferromagnetic version of the REM. The solution method uses the mapping of the evolutionary dynamics into a quantum Ising chain in a transverse field and the Suzuki-Trotter formalism to calculate the transition probabilities between configurations at different times. We find that in the case of the REM landscape the dynamics can exhibit three distinct regimes: pure diffusion or stasis for short times, depending on the fitness of the initial configuration, and a spin-glass regime for large times. The dynamic transition between these dynamical regimes is marked by discontinuities in the mean-fitness as well as in the overlap with the initial reference sequence. The relaxation to equilibrium is described by an inverse time decay. In the ferromagnetic REM, we find in addition to these three regimes, a ferromagnetic regime where the overlap and the mean-fitness are frozen. In this case, the system relaxes to equilibrium in a finite time. The relevance of our results to information processing aspects of evolution is discussed.

A theoretical model of the evolution of virulence in sexually transmitted HIV/AIDS

INTRODUCTION: The evolution of virulence in host-parasite relationships has been the subject of several publications. In the case of HIV virulence, some authors suggest that the evolution of HIV virulence correlates with the rate of acquisition of new sexual partners. In contrast some other authors argue that the level of HIV virulence is independent of the sexual activity of the host population. METHODS: Provide a mathematical model for the study of the potential influence of human sexual behaviour on the evolution of virulence of HIV is provided. RESULTS: The results indicated that, when the probability of acquisition of infection is a function both of the sexual activity and of the virulence level of HIV strains, the evolution of HIV virulence correlates positively with the rate of acquisition of new sexual partners. CONCLUSION: It is concluded that in the case of a host population with a low (high) rate of exchange of sexual partners the evolution of HIV virulence is such that the less (more) virulent strain prevails.

Economic burden of expected epidemiological changes in diseases related to tobacco, Mexico

OBJECTIVE: To determine health care costs and economic burden of epidemiological changes in diseases related to tobacco consumption. METHODS: A time-series analysis in Mexico (1994-2005) was carried out on seven health interventions: chronic obstructive pulmonary diseases, lung cancer with and without surgical intervention, asthma in smokers and non-smokers, full treatment course with nicotine gum, and full treatment course with nicotine patch. According with Box-Jenkins methodology, probabilistic models were developed to forecast the expected changes in the epidemiologic profile and the expected changes in health care services required for selected interventions. Health care costs were estimated following the instrumentation methods and validated with consensus technique. RESULTS: A comparison of the economic impact in 2006 vs. 2008 showed 20-90% increase in expected cases depending on the disease (p<0.05), and 25-93% increase in financial requirements (p<0.01). The study data suggest that changes in the demand for health services for patients with respiratory diseases related to tobacco consumption will continue showing an increasing trend. CONCLUSIONS: In economic terms, the growing number of cases expected during the study period indicates a process of internal competition and adds an element of intrinsic competition in the management of preventive and curative interventions. The study results support the assumption that if preventive programs remain unchanged, the increasing demands for curative health care may cause great financial and management challenges to the health care system of middle-income countries like Mexico.

A mathematical model to describe fat oxidation kinetics during graded exercise.

PURPOSE: The purpose of this study was to develop a mathematical model (sine model, SIN) to describe fat oxidation kinetics as a function of the relative exercise intensity [% of maximal oxygen uptake (%VO2max)] during graded exercise and to determine the exercise intensity (Fatmax) that elicits maximal fat oxidation (MFO) and the intensity at which the fat oxidation becomes negligible (Fatmin). This model included three independent variables (dilatation, symmetry, and translation) that incorporated primary expected modulations of the curve because of training level or body composition. METHODS: Thirty-two healthy volunteers (17 women and 15 men) performed a graded exercise test on a cycle ergometer, with 3-min stages and 20-W increments. Substrate oxidation rates were determined using indirect calorimetry. SIN was compared with measured values (MV) and with other methods currently used [i.e., the RER method (MRER) and third polynomial curves (P3)]. RESULTS: There was no significant difference in the fitting accuracy between SIN and P3 (P = 0.157), whereas MRER was less precise than SIN (P < 0.001). Fatmax (44 +/- 10% VO2max) and MFO (0.37 +/- 0.16 g x min(-1)) determined using SIN were significantly correlated with MV, P3, and MRER (P < 0.001). The variable of dilatation was correlated with Fatmax, Fatmin, and MFO (r = 0.79, r = 0.67, and r = 0.60, respectively, P < 0.001). CONCLUSIONS: The SIN model presents the same precision as other methods currently used in the determination of Fatmax and MFO but in addition allows calculation of Fatmin. Moreover, the three independent variables are directly related to the main expected modulations of the fat oxidation curve. SIN, therefore, seems to be an appropriate tool in analyzing fat oxidation kinetics obtained during graded exercise.

CD8 kinetically promotes ligand binding to the T-cell antigen receptor.

The mechanism of CD8 cooperation with the TCR in antigen recognition was studied on live T cells. Fluorescence correlation measurements yielded evidence of the presence of two TCR and CD8 subpopulations with different lateral diffusion rate constants. Independently, evidence for two subpopulations was derived from the experimentally observed two distinct association phases of cognate peptide bound to class I MHC (pMHC) tetramers and the T cells. The fast phase rate constant ((1.7 +/- 0.2) x 10(5) M(-1) s(-1)) was independent of examined cell type or MHC-bound peptides' structure. Its value was much faster than that of the association of soluble pMHC and TCR ((7.0 +/- 0.3) x 10(3) M(-1) s(-1)), and close to that of the association of soluble pMHC with CD8 ((1-2) x 10(5) M(-1) s(-1)). The fast binding phase disappeared when CD8-pMHC interaction was blocked by a CD8-specific mAb. The latter rate constant was slowed down approximately 10-fold after cells treatment with methyl-beta-cyclodextrin. These results suggest that the most efficient pMHC-cell association route corresponds to a fast tetramer binding to a colocalized CD8-TCR subpopulation, which apparently resides within membrane rafts: the reaction starts by pMHC association with the CD8. This markedly faster step significantly increases the probability of pMHC-TCR encounters and thereby promotes pMHC association with CD8-proximal TCR. The slow binding phase is assigned to pMHC association with a noncolocalized CD8-TCR subpopulation. Taken together with results of cytotoxicity assays, our data suggest that the colocalized, raft-associated CD8-TCR subpopulation is the one capable of inducing T-cell activation.

Predicting hematologic toxicity in patients undergoing radioimmunotherapy with 90Y-ibritumomab tiuxetan or 131I-tositumomab.

This study aimed at identifying clinical factors for predicting hematologic toxicity after radioimmunotherapy with (90)Y-ibritumomab tiuxetan or (131)I-tositumomab in clinical practice. Hematologic data were available from 14 non-Hodgkin lymphoma patients treated with (90)Y-ibritumomab tiuxetan and 18 who received (131)I-tositumomab. The percentage baseline at nadir and 4 wk post nadir and the time to nadir were selected as the toxicity indicators for both platelets and neutrophils. Multiple linear regression analysis was performed to identify significant predictors (P < 0.05) of each indicator. For both platelets and neutrophils, pooled and separate analyses of (90)Y-ibritumomab tiuxetan and (131)I-tositumomab data yielded the time elapsed since the last chemotherapy as the only significant predictor of the percentage baseline at nadir. The extent of bone marrow involvement was not a significant factor in this study, possibly because of the short time elapsed since the last chemotherapy of the 7 patients with bone marrow involvement. Because both treatments were designed to deliver a comparable bone marrow dose, this factor also was not significant. None of the 14 factors considered was predictive of the time to nadir. The R(2) value for the model predicting percentage baseline at nadir was 0.60 for platelets and 0.40 for neutrophils. This model predicted the platelet and neutrophil toxicity grade to within ±1 for 28 and 30 of the 32 patients, respectively. For the 7 patients predicted with grade I thrombocytopenia, 6 of whom had actual grade I-II, dosing might be increased to improve treatment efficacy. The elapsed time since the last chemotherapy can be used to predict hematologic toxicity and customize the current dosing method in radioimmunotherapy.

Microsatellites can be misleading: an empirical and simulation study.

It has been long recognized that highly polymorphic genetic markers can lead to underestimation of divergence between populations when migration is low. Microsatellite loci, which are characterized by extremely high mutation rates, are particularly likely to be affected. Here, we report genetic differentiation estimates in a contact zone between two chromosome races of the common shrew (Sorex araneus), based on 10 autosomal microsatellites, a newly developed Y-chromosome microsatellite, and mitochondrial DNA. These results are compared to previous data on proteins and karyotypes. Estimates of genetic differentiation based on F- and R-statistics are much lower for autosomal microsatellites than for all other genetic markers. We show by simulations that this discrepancy stems mainly from the high mutation rate of microsatellite markers for F-statistics and from deviations from a single-step mutation model for R-statistics. The sex-linked genetic markers show that all gene exchange between races is mediated by females. The absence of male-mediated gene flow most likely results from male hybrid sterility.

Plasma aldosterone is independently associated with the metabolic syndrome.

The aim of this study was to analyze the associations of plasma aldosterone and plasma renin activity with the metabolic syndrome and each of its components. We analyzed data from a family based study in the Seychelles made up of 356 participants (160 men and 196 women) from 69 families of African descent. In multivariable models, plasma aldosterone was associated positively (P < 0.05) with blood pressure in older individuals (interaction with age, P < 0.05) and with waist circumference in men (interaction with sex, P < 0.05) and negatively with high-density lipoprotein cholesterol, in particular in individuals with elevated urinary potassium excretion (interaction with urinary potassium, P < 0.05); plasma renin activity was significantly associated with triglycerides and fasting blood glucose. Plasma aldosterone, but not plasma renin activity, was associated with the metabolic syndrome per se, independently of the association with its separate components. The observation that plasma renin activity was associated with some components of the metabolic syndrome, whereas plasma aldosterone was associated with other components of the metabolic syndrome, suggests different underlying mechanisms. These findings reinforce previous observations suggesting that aldosterone is associated with several cardiovascular risk factors and also suggest that aldosterone might contribute to the increased cardiovascular disease risk in individuals of African descent with the metabolic syndrome.