Therapeutic drug monitoring of imatinib: Bayesian and alternative methods to predict trough levels

Background: The imatinib trough plasma concentration (C(min)) correlates with clinical response in cancer patients. Therapeutic drug monitoring (TDM) of plasma C(min) is therefore suggested. In practice, however, blood sampling for TDM is often not performed at trough. The corresponding measurement is thus only remotely informative about C(min) exposure. Objectives: The objectives of this study were to improve the interpretation of randomly measured concentrations by using a Bayesian approach for the prediction of C(min), incorporating correlation between pharmacokinetic parameters, and to compare the predictive performance of this method with alternative approaches, by comparing predictions with actual measured trough levels, and with predictions obtained by a reference method, respectively. Methods: A Bayesian maximum a posteriori (MAP) estimation method accounting for correlation (MAP-ρ) between pharmacokinetic parameters was developed on the basis of a population pharmacokinetic model, which was validated on external data. Thirty-one paired random and trough levels, observed in gastrointestinal stromal tumour patients, were then used for the evaluation of the Bayesian MAP-ρ method: individual C(min) predictions, derived from single random observations, were compared with actual measured trough levels for assessment of predictive performance (accuracy and precision). The method was also compared with alternative approaches: classical Bayesian MAP estimation assuming uncorrelated pharmacokinetic parameters, linear extrapolation along the typical elimination constant of imatinib, and non-linear mixed-effects modelling (NONMEM) first-order conditional estimation (FOCE) with interaction. Predictions of all methods were finally compared with 'best-possible' predictions obtained by a reference method (NONMEM FOCE, using both random and trough observations for individual C(min) prediction). Results: The developed Bayesian MAP-ρ method accounting for correlation between pharmacokinetic parameters allowed non-biased prediction of imatinib C(min) with a precision of ±30.7%. This predictive performance was similar for the alternative methods that were applied. The range of relative prediction errors was, however, smallest for the Bayesian MAP-ρ method and largest for the linear extrapolation method. When compared with the reference method, predictive performance was comparable for all methods. The time interval between random and trough sampling did not influence the precision of Bayesian MAP-ρ predictions. Conclusion: Clinical interpretation of randomly measured imatinib plasma concentrations can be assisted by Bayesian TDM. Classical Bayesian MAP estimation can be applied even without consideration of the correlation between pharmacokinetic parameters. Individual C(min) predictions are expected to vary less through Bayesian TDM than linear extrapolation. Bayesian TDM could be developed in the future for other targeted anticancer drugs and for the prediction of other pharmacokinetic parameters that have been correlated with clinical outcomes.

Binary classification of dyslipidemia from the waist-to-hip ratio and body mass index: a comparison of linear, logistic, and CART models.

BACKGROUND: We sought to improve upon previously published statistical modeling strategies for binary classification of dyslipidemia for general population screening purposes based on the waist-to-hip circumference ratio and body mass index anthropometric measurements. METHODS: Study subjects were participants in WHO-MONICA population-based surveys conducted in two Swiss regions. Outcome variables were based on the total serum cholesterol to high density lipoprotein cholesterol ratio. The other potential predictor variables were gender, age, current cigarette smoking, and hypertension. The models investigated were: (i) linear regression; (ii) logistic classification; (iii) regression trees; (iv) classification trees (iii and iv are collectively known as "CART"). Binary classification performance of the region-specific models was externally validated by classifying the subjects from the other region. RESULTS: Waist-to-hip circumference ratio and body mass index remained modest predictors of dyslipidemia. Correct classification rates for all models were 60-80%, with marked gender differences. Gender-specific models provided only small gains in classification. The external validations provided assurance about the stability of the models. CONCLUSIONS: There were no striking differences between either the algebraic (i, ii) vs. non-algebraic (iii, iv), or the regression (i, iii) vs. classification (ii, iv) modeling approaches. Anticipated advantages of the CART vs. simple additive linear and logistic models were less than expected in this particular application with a relatively small set of predictor variables. CART models may be more useful when considering main effects and interactions between larger sets of predictor variables.

Synchronization of EEG: bivariate and multivariate measures.

Synchronization behavior of electroencephalographic (EEG) signals is important for decoding information processing in the human brain. Modern multichannel EEG allows a transition from traditional measurements of synchronization in pairs of EEG signals to whole-brain synchronization maps. The latter can be based on bivariate measures (BM) via averaging over pair-wise values or, alternatively, on multivariate measures (MM), which directly ascribe a single value to the synchronization in a group. In order to compare BM versus MM, we applied nine different estimators to simulated multivariate time series with known parameters and to real EEGs.We found widespread correlations between BM and MM, which were almost frequency-independent for all the measures except coherence. The analysis of the behavior of synchronization measures in simulated settings with variable coupling strength, connection probability, and parameter mismatch showed that some of them, including S-estimator, S-Renyi, omega, and coherence, aremore sensitive to linear interdependences,while others, like mutual information and phase locking value, are more responsive to nonlinear effects. Onemust consider these properties together with the fact thatMM are computationally less expensive and, therefore, more efficient for the large-scale data sets than BM while choosing a synchronization measure for EEG analysis.

Association between pre-diagnostic circulating vitamin D concentration and risk of colorectal cancer in European populations:a nested case-control study

Objective. To examine the association between pre-diagnostic circulating vitamin D concentration, dietary intake of vitamin D and calcium, and the risk of colorectal cancer in European populations. Design Nested case-control study. Setting. The study was conducted within the EPIC study, a cohort of more than 520 000 participants from 10 western European countries. Participants: 1248 cases of incident colorectal cancer, which developed after enrolment into the cohort, were matched to 1248 controls. Main outcome measures. Circulating vitamin D concentration (25-hydroxy-vitamin-D, 25-(OH)D) was measured by enzyme immunoassay. Dietary and lifestyle data were obtained from questionnaires. Incidence rate ratios and 95% confidence intervals for the risk of colorectal cancer by 25-(OH)D concentration and levels of dietary calcium and vitamin D intake were estimated from multivariate conditional logistic regression models, with adjustment for potential dietary and other confounders. Results. 25-(OH)D concentration showed a strong inverse linear dose-response association with risk of colorectal cancer (P for trend <0.001). Compared with a pre-defined mid-level concentration of 25-(OH)D (50.0-75.0 nmol/l), lower levels were associated with higher colorectal cancer risk (<25.0 nmol/l: incidence rate ratio 1.32 (95% confidence interval 0.87 to 2.01); 25.0-49.9 nmol/l: 1.28 (1.05 to 1.56), and higher concentrations associated with lower risk (75.0-99.9 nmol/l: 0.88 (0.68 to 1.13); ≥100.0 nmol/l: 0.77 (0.56 to 1.06)). In analyses by quintile of 25-(OH)D concentration, patients in the highest quintile had a 40% lower risk of colorectal cancer than did those in the lowest quintile (P<0.001). Subgroup analyses showed a strong association for colon but not rectal cancer (P for heterogeneity=0.048). Greater dietary intake of calcium was associated with a lower colorectal cancer risk. Dietary vitamin D was not associated with disease risk. Findings did not vary by sex and were not altered by corrections for season or month of blood donation. Conclusions The results of this large observational study indicate a strong inverse association between levels of pre-diagnostic 25-(OH)D concentration and risk of colorectal cancer in western European populations. Further randomised trials are needed to assess whether increases in circulating 25-(OH)D concentration can effectively decrease the risk of colorectal cancer.

Factors associated with improvement in disability-adjusted life years in patients with HIV/AIDS

Background The epidemic of HIV/AIDS and treatments that have emerged to alleviate, have brought about a shift in the burden of disease from death to quality of life/disability. The aim was to determine which factors are associated with improvements in the level of health of male and female patients with HIV/AIDS in Andalusia, in terms of disability-adjusted life years. Methods Descriptive study based on a sample group of 8800 people on the Andalusian AIDS register between 1983 and 2004. Dependent variables: Life lost due to premature mortality (YLL), years lost due to disability (YLD) and disability-adjusted life years (DALY). Independent variables: vital state, sex, age at the time of diagnosis, age at the time of death, transmission category, province of residence, AIDS-indicator disease and the period of diagnosis. A bivariate analysis was carried out to find out if the health level variables changed in accordance with the independent variables. Using the independent variables which had a statistically significant link with the level of health variables, a multivariate linear regression model, disaggregated by gender, was constructed. Results Amongst the women, we found a model which explained the level of health of 64.9%: a link was found between a higher level of health (lower DALYs) and not intravenous drug use, the province of residence, being diagnosed during the HAART era and older age at the time of diagnosis. Amongst the men, we found a model which explained the level of health of 64.4%: a link was found between a higher level of health (lower DALYs) and intravenous drug use, the province of residence, being diagnosed during the HAART era and older age at the time of diagnosis. Conclusion A higher level of health (lower DALY) amongst both men and women was found to be linked to not be intravenous drug user, the province of residence, being diagnosed during the HAART era and older age at the time of diagnosis.

Macronutrient composition of the diet and prospective weight change in participants of the EPIC-PANACEA study.

BACKGROUND The effect of the macronutrient composition of the usual diet on long term weight maintenance remains controversial. METHODS 373,803 subjects aged 25-70 years were recruited in 10 European countries (1992-2000) in the PANACEA project of the EPIC cohort. Diet was assessed at baseline using country-specific validated questionnaires and weight and height were measured at baseline and self-reported at follow-up in most centers. The association between weight change after 5 years of follow-up and the iso-energetic replacement of 5% of energy from one macronutrient by 5% of energy from another macronutrient was assessed using multivariate linear mixed-models. The risk of becoming overweight or obese after 5 years was investigated using multivariate Poisson regressions stratified according to initial Body Mass Index. RESULTS A higher proportion of energy from fat at the expense of carbohydrates was not significantly associated with weight change after 5 years. However, a higher proportion of energy from protein at the expense of fat was positively associated with weight gain. A higher proportion of energy from protein at the expense of carbohydrates was also positively associated with weight gain, especially when carbohydrates were rich in fibre. The association between percentage of energy from protein and weight change was slightly stronger in overweight participants, former smokers, participants ≥60 years old, participants underreporting their energy intake and participants with a prudent dietary pattern. Compared to diets with no more than 14% of energy from protein, diets with more than 22% of energy from protein were associated with a 23-24% higher risk of becoming overweight or obese in normal weight and overweight subjects at baseline. CONCLUSION Our results show that participants consuming an amount of protein above the protein intake recommended by the American Diabetes Association may experience a higher risk of becoming overweight or obese during adult life.

Mirror image atrial dilatation in adult patients with atrial fibrillation and congenital heart disease.

BACKGROUND: Atrial fibrillation (AF) is largely regarded to be initiated from left atrial (LA) dilatation, with subsequent dilatation of the right atrium (RA) in those who progress to chronic AF. We hypothesized that in adult patients with right-sided congenital heart disease (CHD) and AF, RA dilatation will predominate with subsequent dilatation of the left atrium, as a mirror image. METHODS: Adult patients with diagnosis of right-sided, ASD or left-sided CHD who had undergone an echocardiographic study and electrocardiographic recording in 2007 were included. RA and LA area were measured from the apical view. AF was diagnosed from a 12-lead electrocardiogram or Holter recording. A multivariate logistic regression model was used to identify predictors of AF and linear regression models were performed to measure relationship between RA and LA area and AF. RESULTS: A total of 291 patients were included in the study. Multivariate analysis showed that age (p=0.0001), RA (p=0.025) and LA area (p=0.0016) were significantly related to AF. In patients with pure left-sided pathologies, there was progressive and predominant LA dilatation that paralleled the development of AF from none to paroxysmal to chronic AF. In patients with pure right-sided pathologies, there was a mirror image of progressive and predominant RA dilatation with the development of AF. CONCLUSION: We observed a mirror image atrial dilatation in patients with right sided disease and AF. This may provide novel mechanistic insight as to the origin of AF in these patients and deserves further studying in the form of targeted electrophysiological studies.

TRAF1/C5 but not PTPRC variants are potential predictors of rheumatoid arthritis response to anti-tumor necrosis factor therapy.

BACKGROUND The aim of our work was to replicate, in a Southern European population, the association reported in Northern populations between PTPRC locus and response to anti-tumor necrosis factor (anti-TNF) treatment in rheumatoid arthritis (RA). We also looked at associations between five RA risk alleles and treatment response. METHODS We evaluated associations between anti-TNF treatment responses assessed by DAS28 change and by EULAR response at six months in 383 Portuguese patients. Univariate and multivariate linear and logistic regression analyses were performed. In a second step to confirm our findings, we pooled our population with 265 Spanish patients. RESULTS No association was found between PTPRC rs10919563 allele and anti-TNF treatment response, neither in Portuguese modeling for several clinical variables nor in the overall population combining Portuguese and Spanish patients. The minor allele for RA susceptibility, rs3761847 SNP in TRAF1/C5 region, was associated with a poor response in linear and logistic univariate and multivariate regression analyses. No association was observed with the other allellic variants. Results were confirmed in the pooled analysis. CONCLUSION This study did not replicate the association between PTPRC and the response to anti-TNF treatment in our Southern European population. We found that TRAF1/C5 risk RA variants potentially influence anti-TNF treatment response.

Genotype 3 is associated with accelerated fibrosis progression in chronic hepatitis C.

BACKGROUND/AIMS: While several risk factors for the histological progression of chronic hepatitis C have been identified, the contribution of HCV genotypes to liver fibrosis evolution remains controversial. The aim of this study was to assess independent predictors for fibrosis progression. METHODS: We identified 1189 patients from the Swiss Hepatitis C Cohort database with at least one biopsy prior to antiviral treatment and assessable date of infection. Stage-constant fibrosis progression rate was assessed using the ratio of fibrosis Metavir score to duration of infection. Stage-specific fibrosis progression rates were obtained using a Markov model. Risk factors were assessed by univariate and multivariate regression models. RESULTS: Independent risk factors for accelerated stage-constant fibrosis progression (>0.083 fibrosis units/year) included male sex (OR=1.60, [95% CI 1.21-2.12], P<0.001), age at infection (OR=1.08, [1.06-1.09], P<0.001), histological activity (OR=2.03, [1.54-2.68], P<0.001) and genotype 3 (OR=1.89, [1.37-2.61], P<0.001). Slower progression rates were observed in patients infected by blood transfusion (P=0.02) and invasive procedures or needle stick (P=0.03), compared to those infected by intravenous drug use. Maximum likelihood estimates (95% CI) of stage-specific progression rates (fibrosis units/year) for genotype 3 versus the other genotypes were: F0-->F1: 0.126 (0.106-0.145) versus 0.091 (0.083-0.100), F1-->F2: 0.099 (0.080-0.117) versus 0.065 (0.058-0.073), F2-->F3: 0.077 (0.058-0.096) versus 0.068 (0.057-0.080) and F3-->F4: 0.171 (0.106-0.236) versus 0.112 (0.083-0.142, overall P<0.001). CONCLUSIONS: This study shows a significant association of genotype 3 with accelerated fibrosis using both stage-constant and stage-specific estimates of fibrosis progression rates. This observation may have important consequences for the management of patients infected with this genotype.

Variational calculation of vibrational linear and nonlinear optical properties

A variational approach for reliably calculating vibrational linear and nonlinear optical properties of molecules with large electrical and/or mechanical anharmonicity is introduced. This approach utilizes a self-consistent solution of the vibrational Schrödinger equation for the complete field-dependent potential-energy surface and, then, adds higher-level vibrational correlation corrections as desired. An initial application is made to static properties for three molecules of widely varying anharmonicity using the lowest-level vibrational correlation treatment (i.e., vibrational Møller-Plesset perturbation theory). Our results indicate when the conventional Bishop-Kirtman perturbation method can be expected to break down and when high-level vibrational correlation methods are likely to be required. Future improvements and extensions are discussed

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.

Generalized linear and generalized additive models in studies of species distributions: setting the scene

An important statistical development of the last 30 years has been the advance in regression analysis provided by generalized linear models (GLMs) and generalized additive models (GAMs). Here we introduce a series of papers prepared within the framework of an international workshop entitled: Advances in GLMs/GAMs modeling: from species distribution to environmental management, held in Riederalp, Switzerland, 6-11 August 2001.We first discuss some general uses of statistical models in ecology, as well as provide a short review of several key examples of the use of GLMs and GAMs in ecological modeling efforts. We next present an overview of GLMs and GAMs, and discuss some of their related statistics used for predictor selection, model diagnostics, and evaluation. Included is a discussion of several new approaches applicable to GLMs and GAMs, such as ridge regression, an alternative to stepwise selection of predictors, and methods for the identification of interactions by a combined use of regression trees and several other approaches. We close with an overview of the papers and how we feel they advance our understanding of their application to ecological modeling.

Equivalence of piecewise-linear approximation and Lagrangian relaxation for network revenue management

The network revenue management (RM) problem arises in airline, hotel, media,and other industries where the sale products use multiple resources. It can be formulatedas a stochastic dynamic program but the dynamic program is computationallyintractable because of an exponentially large state space, and a number of heuristicshave been proposed to approximate it. Notable amongst these -both for their revenueperformance, as well as their theoretically sound basis- are approximate dynamic programmingmethods that approximate the value function by basis functions (both affinefunctions as well as piecewise-linear functions have been proposed for network RM)and decomposition methods that relax the constraints of the dynamic program to solvesimpler dynamic programs (such as the Lagrangian relaxation methods). In this paperwe show that these two seemingly distinct approaches coincide for the network RMdynamic program, i.e., the piecewise-linear approximation method and the Lagrangianrelaxation method are one and the same.

Asymptotic robust inferences in the analysis of mean and covariance structures

Structural equation models are widely used in economic, socialand behavioral studies to analyze linear interrelationships amongvariables, some of which may be unobservable or subject to measurementerror. Alternative estimation methods that exploit different distributionalassumptions are now available. The present paper deals with issues ofasymptotic statistical inferences, such as the evaluation of standarderrors of estimates and chi--square goodness--of--fit statistics,in the general context of mean and covariance structures. The emphasisis on drawing correct statistical inferences regardless of thedistribution of the data and the method of estimation employed. A(distribution--free) consistent estimate of $\Gamma$, the matrix ofasymptotic variances of the vector of sample second--order moments,will be used to compute robust standard errors and a robust chi--squaregoodness--of--fit squares. Simple modifications of the usual estimateof $\Gamma$ will also permit correct inferences in the case of multi--stage complex samples. We will also discuss the conditions under which,regardless of the distribution of the data, one can rely on the usual(non--robust) inferential statistics. Finally, a multivariate regressionmodel with errors--in--variables will be used to illustrate, by meansof simulated data, various theoretical aspects of the paper.

Step count is associated with lower nighttime systolic blood pressure and increased dipping.

BACKGROUND: Higher nighttime blood pressure (BP) and the loss of nocturnal dipping of BP are associated with an increased risk for cardiovascular events. However, the determinants of the loss of nocturnal BP dipping are only beginning to be understood. We investigated whether different indicators of physical activity were associated with the loss of nocturnal dipping of BP. METHODS: We conducted a cross-sectional study of 103 patients referred for 24-hour ambulatory monitoring of BP. We measured these patients' step count (SC), active energy expenditure (AEE), and total energy expenditure simultaneously, using actigraphs. RESULTS: In our study population of 103 patients, most of whom were hypertensive, SC and AEE were associated with nighttime systolic BP in univariate (SC, r = -0.28, P < 0.01; AEE, r = -0.20, P = 0.046) and multivariate linear regression analyses (SC, coefficient beta = -5.37, P < 0.001; AEE, coefficient beta = -0.24, P < 0.01). Step count was associated with both systolic (r = 0.23, P = 0.018) and diastolic (r = 0.20, P = 0.045) BP dipping. Nighttime systolic BP decreased progressively across the categories of sedentary, moderately active, and active participants (125mm Hg, 116mm Hg, 112mm Hg, respectively; P = 0.002). The degree of BP dipping of BP increased progressively across the same three categories of activity (respectively 8.9%, 14.6%, and 18.6%, P = 0.002, for systolic BP and respectively 12.8%, 18.1%, and 22.2%, P = 0.006, for diastolic BP). CONCLUSIONS: Step count is continuously associated with nighttime systolic BP and with the degree of BP dipping independently of 24-hour mean BP. The combined use of an actigraph for measuring indicators of physical activity and a device for 24-hour measurement of ambulatory BP may help identify patients at increased risk for cardiovascular events in whom increased physical activity toward higher target levels may be recommended.