gems: An R Package for Simulating from Disease Progression Models

Mathematical models of disease progression predict disease outcomes and are useful epidemiological tools for planners and evaluators of health interventions. The R package gems is a tool that simulates disease progression in patients and predicts the effect of different interventions on patient outcome. Disease progression is represented by a series of events (e.g., diagnosis, treatment and death), displayed in a directed acyclic graph. The vertices correspond to disease states and the directed edges represent events. The package gems allows simulations based on a generalized multistate model that can be described by a directed acyclic graph with continuous transition-specific hazard functions. The user can specify an arbitrary hazard function and its parameters. The model includes parameter uncertainty, does not need to be a Markov model, and may take the history of previous events into account. Applications are not limited to the medical field and extend to other areas where multistate simulation is of interest. We provide a technical explanation of the multistate models used by gems, explain the functions of gems and their arguments, and show a sample application.

Fluid intake and all-cause mortality, cardiovascular mortality, and kidney function: a population-based longitudinal cohort study

BACKGROUND Drinking eight glasses of fluid or water each day is widely believed to improve health, but evidence is sparse and conflicting. We aimed to investigate the association between fluid consumption and long-term mortality and kidney function. METHODS We conducted a longitudinal analysis within a prospective, population-based cohort study of 3858 men and women aged 49 years or older residing in Australia. Daily fluid intake from food and beverages not including water was measured using a food frequency questionnaire. We did multivariable adjusted Cox proportional hazard models for all-cause and cardiovascular mortality and a boot-strapping procedure for estimated glomerular filtration rate (eGFR). RESULTS Upper and lower quartiles of daily fluid intake corresponded to >3 L and <2 L, respectively. During a median follow-up of 13.1 years (total 43 093 years at risk), 1127 deaths (26.1 per 1000 years at risk) including 580 cardiovascular deaths (13.5 per 1000 years at risk) occurred. Daily fluid intake (per 250 mL increase) was not associated with all-cause [adjusted hazard ratio (HR) 0.99 (95% CI 0.98-1.01)] or cardiovascular mortality [HR 0.98 (95% CI 0.95-1.01)]. Overall, eGFR reduced by 2.2 mL/min per 1.73 m(2) (SD 10.9) in the 1207 (31%) participants who had repeat creatinine measurements and this was not associated with fluid intake [adjusted regression coefficient 0.06 mL/min/1.73 m(2) per 250 mL increase (95% CI -0.03 to 0.14)]. CONCLUSIONS Fluid intake from food and beverages excluding water is not associated with improved kidney function or reduced mortality.

Prognostic relevance of coronary collateral function: confounded or causal relationship?

OBJECTIVE To expand the limited information on the prognostic impact of quantitatively obtained collateral function in patients with coronary artery disease (CAD) and to estimate causality of such a relation. DESIGN Prospective cohort study with long-term observation of clinical outcome. SETTING University Hospital. PATIENTS One thousand one hundred and eighty-one patients with chronic stable CAD undergoing 1771 quantitative, coronary pressure-derived collateral flow index measurements, as obtained during a 1-min coronary balloon occlusion (CFI is the ratio between mean distal coronary occlusive pressure and mean aortic pressure both subtracted by central venous pressure). Subgroup of 152 patients included in randomised trials on the longitudinal effect of different arteriogenic protocols on CFI. INTERVENTIONS Collection of long-term follow-up information on clinical outcome. MAIN OUTCOME MEASURES All-cause mortality and major adverse cardiac events. RESULTS Cumulative 15-year survival rate was 48% in patients with CFI<0.25 and 65% in the group with CFI≥0.25 (p=0.0057). Cumulative 10-year survival rate was 75% in patients without arteriogenic therapy and 88% (p=0.0482) in the group with arteriogenic therapy and showing a significant increase in CFI at follow-up. By proportional hazard analysis, the following variables predicted increased all-cause mortality: age, low CFI, left ventricular end-diastolic pressure and number of vessels with CAD. CONCLUSIONS A well-functioning coronary collateral circulation independently predicts lowered mortality in patients with chronic CAD. This relation appears to be causal, because augmented collateral function by arteriogenic therapy is associated with prolonged survival.

EGFR inhibitors erlotinib and lapatinib ameliorate epidermal blistering in pemphigus vulgaris in a non-linear, V-shaped relationship.

Novel insights into intra-cellular signalling involved in pemphigus vulgaris (PV), an autoimmune blistering disease of skin and mucous membranes, are now revealing new therapeutic approaches such as the chemical inhibition of PV-associated signals in conjunction with standard immunosuppressive therapy. However, extensive inhibition of signalling molecules that are required for normal tissue function and integrity may hamper this approach. Using a neonatal PV mouse model, we demonstrate that epidermal blistering can be prevented in a dose-dependent manner by clinically approved EGFR inhibitors erlotinib and lapatinib, but only up to approximately 50% of normal EGFR activity. At lower EGFR activity, blisters again aggravated and were highly exacerbated in mice with a conditional deletion of EGFR. Statistical analysis of the relation between EGFR activity and the extent of skin blistering revealed the best fit with a non-linear, V-shaped curve with a median break point at 52% EGFR activity (P = 0.0005). Moreover, lapatinib (a dual EGFR/ErbB2 inhibitor) but not erlotinib significantly reduced blistering in the oral cavity, suggesting that signalling mechanisms differ between PV predilection sites. Our results demonstrate that future clinical trials evaluating EGFR/ErbB2 inhibitors in PV patients must select treatment doses that retain a specific level of signal molecule activity. These findings may also be of relevance for cancer patients treated with EGFR inhibitors, for whom skin lesions due to extensive EGFR inhibition represent a major threat.

Thyroid function within the normal range and risk of coronary heart disease: an individual participant data analysis of 14 cohorts

IMPORTANCE Some experts suggest that serum thyrotropin levels in the upper part of the current reference range should be considered abnormal, an approach that would reclassify many individuals as having mild hypothyroidism. Health hazards associated with such thyrotropin levels are poorly documented, but conflicting evidence suggests that thyrotropin levels in the upper part of the reference range may be associated with an increased risk of coronary heart disease (CHD). OBJECTIVE To assess the association between differences in thyroid function within the reference range and CHD risk. DESIGN, SETTING, AND PARTICIPANTS Individual participant data analysis of 14 cohorts with baseline examinations between July 1972 and April 2002 and with median follow-up ranging from 3.3 to 20.0 years. Participants included 55,412 individuals with serum thyrotropin levels of 0.45 to 4.49 mIU/L and no previously known thyroid or cardiovascular disease at baseline. EXPOSURES Thyroid function as expressed by serum thyrotropin levels at baseline. MAIN OUTCOMES AND MEASURES Hazard ratios (HRs) of CHD mortality and CHD events according to thyrotropin levels after adjustment for age, sex, and smoking status. RESULTS Among 55,412 individuals, 1813 people (3.3%) died of CHD during 643,183 person-years of follow-up. In 10 cohorts with information on both nonfatal and fatal CHD events, 4666 of 48,875 individuals (9.5%) experienced a first-time CHD event during 533,408 person-years of follow-up. For each 1-mIU/L higher thyrotropin level, the HR was 0.97 (95% CI, 0.90-1.04) for CHD mortality and 1.00 (95% CI, 0.97-1.03) for a first-time CHD event. Similarly, in analyses by categories of thyrotropin, the HRs of CHD mortality (0.94 [95% CI, 0.74-1.20]) and CHD events (0.97 [95% CI, 0.83-1.13]) were similar among participants with the highest (3.50-4.49 mIU/L) compared with the lowest (0.45-1.49 mIU/L) thyrotropin levels. Subgroup analyses by sex and age group yielded similar results. CONCLUSIONS AND RELEVANCE Thyrotropin levels within the reference range are not associated with risk of CHD events or CHD mortality. This finding suggests that differences in thyroid function within the population reference range do not influence the risk of CHD. Increased CHD risk does not appear to be a reason for lowering the upper thyrotropin reference limit.

Normal Thyroid Function and the Risk of Atrial Fibrillation: the Rotterdam Study

CONTEXT Hyperthyroidism is an established risk factor for atrial fibrillation (AF), but information concerning the association with variations within the normal range of thyroid function and subgroups at risk is lacking. OBJECTIVE This study aimed to investigate the association between normal thyroid function and AF prospectively and explore potential differential risk patterns. DESIGN, SETTING, AND PARTICIPANTS From the Rotterdam Study we included 9166 participants ≥ 45 y with TSH and/or free T4 (FT4) measurements and AF assessment (1997-2012 median followup, 6.8 y), with 399 prevalent and 403 incident AF cases. MAIN OUTCOME MEASURES Outcome measures were 3-fold: 1) hazard ratios (HRs) for the risk of incident AF by Cox proportional-hazards models, 2) 10-year absolute risks taking competing risk of death into account, and 3) discrimination ability of adding FT4 to the CHARGE-AF simple model, an established prediction model for AF. RESULTS Higher FT4 levels were associated with higher risks of AF (HR 1.63, 95% confidence interval, 1.19-2.22), when comparing those in the highest quartile to those in lowest quartile. Absolute 10-year risks increased with higher FT4 in participants ≤ 65 y from 1-9% and from 6-12% in subjects ≥ 65 y. Discrimination of the prediction model improved when adding FT4 to the simple model (c-statistic, 0.722 vs 0.729; P = .039). TSH levels were not associated with AF. CONCLUSIONS There is an increased risk of AF with higher FT4 levels within the normal range, especially in younger subjects. Adding FT4 to the simple model slightly improved discrimination of risk prediction.