Wavelet analysis residual kriging vs. neural network residual kriging

This paper deals with the problem of spatial data mapping. A new method based on wavelet interpolation and geostatistical prediction (kriging) is proposed. The method - wavelet analysis residual kriging (WARK) - is developed in order to assess the problems rising for highly variable data in presence of spatial trends. In these cases stationary prediction models have very limited application. Wavelet analysis is used to model large-scale structures and kriging of the remaining residuals focuses on small-scale peculiarities. WARK is able to model spatial pattern which features multiscale structure. In the present work WARK is applied to the rainfall data and the results of validation are compared with the ones obtained from neural network residual kriging (NNRK). NNRK is also a residual-based method, which uses artificial neural network to model large-scale non-linear trends. The comparison of the results demonstrates the high quality performance of WARK in predicting hot spots, reproducing global statistical characteristics of the distribution and spatial correlation structure.

Detection of Optimal Models in Parameter Space with Support Vector Machines

The paper proposes an approach aimed at detecting optimal model parameter combinations to achieve the most representative description of uncertainty in the model performance. A classification problem is posed to find the regions of good fitting models according to the values of a cost function. Support Vector Machine (SVM) classification in the parameter space is applied to decide if a forward model simulation is to be computed for a particular generated model. SVM is particularly designed to tackle classification problems in high-dimensional space in a non-parametric and non-linear way. SVM decision boundaries determine the regions that are subject to the largest uncertainty in the cost function classification, and, therefore, provide guidelines for further iterative exploration of the model space. The proposed approach is illustrated by a synthetic example of fluid flow through porous media, which features highly variable response due to the parameter values' combination.

BME-based uncertainty assessment of the Chernobyl fallout

The vast territories that have been radioactively contaminated during the 1986 Chernobyl accident provide a substantial data set of radioactive monitoring data, which can be used for the verification and testing of the different spatial estimation (prediction) methods involved in risk assessment studies. Using the Chernobyl data set for such a purpose is motivated by its heterogeneous spatial structure (the data are characterized by large-scale correlations, short-scale variability, spotty features, etc.). The present work is concerned with the application of the Bayesian Maximum Entropy (BME) method to estimate the extent and the magnitude of the radioactive soil contamination by 137Cs due to the Chernobyl fallout. The powerful BME method allows rigorous incorporation of a wide variety of knowledge bases into the spatial estimation procedure leading to informative contamination maps. Exact measurements (?hard? data) are combined with secondary information on local uncertainties (treated as ?soft? data) to generate science-based uncertainty assessment of soil contamination estimates at unsampled locations. BME describes uncertainty in terms of the posterior probability distributions generated across space, whereas no assumption about the underlying distribution is made and non-linear estimators are automatically incorporated. Traditional estimation variances based on the assumption of an underlying Gaussian distribution (analogous, e.g., to the kriging variance) can be derived as a special case of the BME uncertainty analysis. The BME estimates obtained using hard and soft data are compared with the BME estimates obtained using only hard data. The comparison involves both the accuracy of the estimation maps using the exact data and the assessment of the associated uncertainty using repeated measurements. Furthermore, a comparison of the spatial estimation accuracy obtained by the two methods was carried out using a validation data set of hard data. Finally, a separate uncertainty analysis was conducted that evaluated the ability of the posterior probabilities to reproduce the distribution of the raw repeated measurements available in certain populated sites. The analysis provides an illustration of the improvement in mapping accuracy obtained by adding soft data to the existing hard data and, in general, demonstrates that the BME method performs well both in terms of estimation accuracy as well as in terms estimation error assessment, which are both useful features for the Chernobyl fallout study.

Cytogenetic response to imatinib in chronic myeloid leukaemia patients: concentration dependency and associated influences under field-conditions

Introduction: Imatinib, a first-line drug for chronic myeloid leukaemia (CML), has been increasingly proposed for therapeutic drug monitoring (TDM), as trough concentrations >=1000 ng/ml (Cmin) have been associated with improved molecular and complete cytogenetic response (CCyR). The pharmacological monitoring project of EUTOS (European Treatment and Outcome Study) was launched to validate retrospectively the correlation between Cmin and response in a large population of patients followed by central TDM in Bordeaux.¦Methods: 1898 CML patients with first TDM 0-9 years after imatinib initiation, providing cytogenetic data along with demographic and comedication (37%) information, were included. Individual Cmin, estimated by non-linear regression (NONMEM), was adjusted to initial standard dose (400 mg/day) and stratified at 1000 ng/ml. Kaplan-Meier estimates of overall cumulative CCyR rates (stratified by sex, age, comedication and Cmin) were compared using asymptotic logrank k-sample test for interval-censored data. Differences in Cmin were assessed by Wilcoxon test.¦Results: There were no significant differences in overall cumulative CCyR rates between Cmin strata, sex and comedication with P-glycoprotein inhibitors/inducers or CYP3A4 inhibitors (p >0.05). Lower rates were observed in 113 young patients <30 years (p = 0.037; 1-year rates: 43% vs 60% in older patients), as well as in 29 patients with CYP3A4 inducers (p = 0.001, 1-year rates: 40% vs 66% without). Higher rates were observed in 108 patients on organic-cation-transporter-1 (hOCT-1) inhibitors (p = 0.034, 1-year rates: 83% vs 56% without). Considering 1-year CCyR rates, a trend towards better response for Cmin above 1000 ng/ml was observed: 64% (95%CI: 60-69%) vs 59% (95%CI: 56-61%). Median Cmin (400 mg/day) was significantly reduced in male patients (732 vs 899ng/ml, p <0.001), young patients <30 years (734 vs 802 ng/ml, p = 0.037) and under CYP3A4 inducers (758 vs 859 ng/ml, p = 0.022). Under hOCT-1 inhibitors, Cmin was increased (939 vs 827 ng/ml, p = 0.038).¦Conclusion: Based on observational TDM data, the impact of imatinib Cmin >1000 ng/ml on CCyR was not salient. Young CML patients (<30 years) and patients taking CYP3A4 inducers probably need close monitoring and possibly higher imatinib doses, due to lower Cmin along with lower CCyR rates. Patients taking hOCT-1 inhibitors seem in contrast to have improved CCyR response rates. The precise role for imatinib TDM remains to be established prospectively.

My older sibling was drunk - younger siblings' drunkenness in relation to parental monitoring and the parent-adolescent relationship.

This study explored the links between having older siblings who get drunk, satisfaction with the parent-adolescent relationship, parental monitoring, and adolescents' risky drinking. Regression models were conducted based on a national representative sample of 3725 8th to 10th graders in Switzerland (mean age 15.0, SD = .93) who indicated having older siblings. Results showed that both parental factors and older siblings' drinking behaviour shape younger siblings' frequency of risky drinking. Parental monitoring showed a linear dose-response relationship, and siblings' influence had an additive effect. There was a non-linear interaction effect between parent-adolescent relationship and older sibling's drunkenness. The findings suggest that, apart from avoiding an increasingly unsatisfactory relationship with their children, parental monitoring appears to be important in preventing risky drinking by their younger children, even if the older sibling drinks in such a way. However, a satisfying relationship with parents does not seem to be sufficient to counterbalance older siblings' influence.

How quickly should we titrate antihypertensive medication? Systematic review modelling blood pressure response from trial data.

Context There are no evidence syntheses available to guide clinicians on when to titrate antihypertensive medication after initiation. Objective To model the blood pressure (BP) response after initiating antihypertensive medication. Data sources electronic databases including Medline, Embase, Cochrane Register and reference lists up to December 2009. Study selection Trials that initiated antihypertensive medication as single therapy in hypertensive patients who were either drug naive or had a placebo washout from previous drugs. Data extraction Office BP measurements at a minimum of two weekly intervals for a minimum of 4 weeks. An asymptotic approach model of BP response was assumed and non-linear mixed effects modelling used to calculate model parameters. Results and conclusions Eighteen trials that recruited 4168 patients met inclusion criteria. The time to reach 50% of the maximum estimated BP lowering effect was 1 week (systolic 0.91 weeks, 95% CI 0.74 to 1.10; diastolic 0.95, 0.75 to 1.15). Models incorporating drug class as a source of variability did not improve fit of the data. Incorporating the presence of a titration schedule improved model fit for both systolic and diastolic pressure. Titration increased both the predicted maximum effect and the time taken to reach 50% of the maximum (systolic 1.2 vs 0.7 weeks; diastolic 1.4 vs 0.7 weeks). Conclusions Estimates of the maximum efficacy of antihypertensive agents can be made early after starting therapy. This knowledge will guide clinicians in deciding when a newly started antihypertensive agent is likely to be effective or not at controlling BP.

Statistical learning theory for geospatial data. Case study: Aral sea

In recent years there has been an explosive growth in the development of adaptive and data driven methods. One of the efficient and data-driven approaches is based on statistical learning theory (Vapnik 1998). The theory is based on Structural Risk Minimisation (SRM) principle and has a solid statistical background. When applying SRM we are trying not only to reduce training error ? to fit the available data with a model, but also to reduce the complexity of the model and to reduce generalisation error. Many nonlinear learning procedures recently developed in neural networks and statistics can be understood and interpreted in terms of the structural risk minimisation inductive principle. A recent methodology based on SRM is called Support Vector Machines (SVM). At present SLT is still under intensive development and SVM find new areas of application (www.kernel-machines.org). SVM develop robust and non linear data models with excellent generalisation abilities that is very important both for monitoring and forecasting. SVM are extremely good when input space is high dimensional and training data set i not big enough to develop corresponding nonlinear model. Moreover, SVM use only support vectors to derive decision boundaries. It opens a way to sampling optimization, estimation of noise in data, quantification of data redundancy etc. Presentation of SVM for spatially distributed data is given in (Kanevski and Maignan 2004).

Deletion of delta-opioid receptor in mice alters skin differentiation and delays wound healing.

In addition to their well-known antinociceptive action, opioids can modulate non-neuronal functions, such as immune activity and physiology of different cell types. Several findings suggest that the delta-opioid receptor (DOR) and its endogenous ligands (enkephalins) are important players in cell differentiation and proliferation. Here we show the expression of DOR in mouse skin and human skin cultured fibroblasts and keratinocytes using RT-PCR. In DOR knock-out (KO) mice, a phenotype of thinner epidermis and higher expression of cell differentiation marker cytokeratin 10 (CK 10) were observed compared with wild type (WT). Using a burn wound model, significant wound healing delay (about 2 days) and severe epidermal hypertrophy were shown at the wound margin of DOR KO mice. This wound healing delay was further investigated by immunohistochemistry using markers for proliferation, differentiation, re-epithelialization, and dermal repair (CK 6, CK 10, and collagen IV). The levels of all these markers were increased in wounds of KO mice compared with WT. During the wound healing, the epidermal thickness in KO mice augments faster and exceeds that of the WT by day 3. These results suggest an essential role of DOR in skin differentiation, proliferation, and migration, factors that are important for wound healing.

The Smiling Project: development of a new program for gait and balance training

Introduction: The SMILING project, a multicentric project fundedby the European Union, aims to develop a new gait and balance trainingprogram to prevent falls in older persons. The program includes the"SMILING shoe", an innovative device that generates mechanical perturbationwhile walking by changing the soles' inclination. Induced perturbationschallenge subjects' balance and force them to react to avoidfalls. By training specifically the complex motor reactions used to maintainbalance when walking on irregular ground, the program will improvesubjects' ability to react in situation of unsteadiness and reduce theirrisk of falling. Methods: The program will be evaluated in a multicentric,cross-over randomized controlled trial. Overall, 112 subjects (aged≥65 years, ≥1 falls, POMA score 22-26/28) will be enrolled. Subjectswill be randomised in 2 groups: group A begin the training with active"SMILING shoes", group B with inactive dummy shoes. After 4 weeksof training, group A and B will exchange the shoes. Supervised trainingsessions (30 minutes twice a week for 8 weeks) include walkingtasks of progressive difficulties.To avoid a learning effect, "SMILINGshoes" perturbations will be generated in a non-linear and chaotic way.Gait performance, fear of falling, and acceptability of the program willbe assessed. Conclusion: The SMILING program is an innovative interventionfor falls prevention in older persons based on gait and balancetraining using chaotic perturbations. Because of the easy use of the"SMILING shoes", this program could be used in various settings, suchas geriatric clinics or at home.

Saturable metabolism of continuous high-dose ifosfamide with mesna and GM-CSF: a pharmacokinetic study in advanced sarcoma patients. Swiss Group for Clinical Cancer Research (SAKK).

BACKGROUND: The aim of this study was to assess the pharmacology, toxicity and activity of high-dose ifosfamide mesna +/- GM-CSF administered by a five-day continuous infusion at a total ifosfamide dose of 12-18 g/m2 in adult patients with advanced sarcomas. PATIENTS AND METHODS: Between January 1991 and October 1992 32 patients with advanced or metastatic sarcoma were entered the study. Twenty-seven patients were pretreated including twenty-three with prior ifosfamide at less than 8 g/m2 total dose/cycle. In 25 patients (27 cycles) extensive pharmacokinetic analyses were performed. RESULTS: The area under the plasma concentration-time curve (AUC) for ifosfamide increased linearly with dose while the AUC's of the metabolites measured in plasma by thin-layer chromatography did not increase with dose, particularly that of the active metabolite isophosphoramide mustard. Furthermore the AUC of the inactive carboxymetabolite did not increase with dose. Interpatient variability of pharmacokinetic parameters was high. Dose-limiting toxicity was myelosuppression at 18 g/m2 total dose with grade 4 neutropenia in five of six patients and grade 4 thrombocytopenia in four of six patients. Therefore the maximum tolerated dose was considered to be 18 g/m2 total dose. There was one CR and eleven PR in twenty-nine evaluable patients (overall response rate 41%). CONCLUSION: Both the activation and inactivation pathways of ifosfamide are non-linear and saturable at high-doses although the pharmacokinetics of the parent drug itself are dose linear. Ifosfamide doses greater than 14-16 g/m2 per cycle appear to result in a relative decrease of the active metabolite isophosphoramide mustard. These data suggest a dose-dependent saturation or even inhibition of ifosfamide metabolism by increasing high dose ifosfamide and suggest the need for further metabolic studies.

Dynamic changes in brain functional connectivity during concurrent dual-task performance.

This study investigated the spatial, spectral, temporal and functional proprieties of functional brain connections involved in the concurrent execution of unrelated visual perception and working memory tasks. Electroencephalography data was analysed using a novel data-driven approach assessing source coherence at the whole-brain level. Three connections in the beta-band (18-24 Hz) and one in the gamma-band (30-40 Hz) were modulated by dual-task performance. Beta-coherence increased within two dorsofrontal-occipital connections in dual-task conditions compared to the single-task condition, with the highest coherence seen during low working memory load trials. In contrast, beta-coherence in a prefrontal-occipital functional connection and gamma-coherence in an inferior frontal-occipitoparietal connection was not affected by the addition of the second task and only showed elevated coherence under high working memory load. Analysis of coherence as a function of time suggested that the dorsofrontal-occipital beta-connections were relevant to working memory maintenance, while the prefrontal-occipital beta-connection and the inferior frontal-occipitoparietal gamma-connection were involved in top-down control of concurrent visual processing. The fact that increased coherence in the gamma-connection, from low to high working memory load, was negatively correlated with faster reaction time on the perception task supports this interpretation. Together, these results demonstrate that dual-task demands trigger non-linear changes in functional interactions between frontal-executive and occipitoparietal-perceptual cortices.

Breakpoints in ventilation, cerebral and muscle oxygenation, and muscle activity during an incremental cycling exercise.

The aim of this study was to locate the breakpoints of cerebral and muscle oxygenation and muscle electrical activity during a ramp exercise in reference to the first and second ventilatory thresholds. Twenty-five cyclists completed a maximal ramp test on an electromagnetically braked cycle-ergometer with a rate of increment of 25 W/min. Expired gazes (breath-by-breath), prefrontal cortex and vastus lateralis (VL) oxygenation [Near-infrared spectroscopy (NIRS)] together with electromyographic (EMG) Root Mean Square (RMS) activity for the VL, rectus femoris (RF), and biceps femoris (BF) muscles were continuously assessed. There was a non-linear increase in both cerebral deoxyhemoglobin (at 56 ± 13% of the exercise) and oxyhemoglobin (56 ± 8% of exercise) concomitantly to the first ventilatory threshold (57 ± 6% of exercise, p > 0.86, Cohen's d < 0.1). Cerebral deoxyhemoglobin further increased (87 ± 10% of exercise) while oxyhemoglobin reached a plateau/decreased (86 ± 8% of exercise) after the second ventilatory threshold (81 ± 6% of exercise, p < 0.05, d > 0.8). We identified one threshold only for muscle parameters with a non-linear decrease in muscle oxyhemoglobin (78 ± 9% of exercise), attenuation in muscle deoxyhemoglobin (80 ± 8% of exercise), and increase in EMG activity of VL (89 ± 5% of exercise), RF (82 ± 14% of exercise), and BF (85 ± 9% of exercise). The thresholds in BF and VL EMG activity occurred after the second ventilatory threshold (p < 0.05, d > 0.6). Our results suggest that the metabolic and ventilatory events characterizing this latter cardiopulmonary threshold may affect both cerebral and muscle oxygenation levels, and in turn, muscle recruitment responses.

Accelerated Microstructure Imaging via Convex Optimization (AMICO) from diffusion MRI data.

Microstructure imaging from diffusion magnetic resonance (MR) data represents an invaluable tool to study non-invasively the morphology of tissues and to provide a biological insight into their microstructural organization. In recent years, a variety of biophysical models have been proposed to associate particular patterns observed in the measured signal with specific microstructural properties of the neuronal tissue, such as axon diameter and fiber density. Despite very appealing results showing that the estimated microstructure indices agree very well with histological examinations, existing techniques require computationally very expensive non-linear procedures to fit the models to the data which, in practice, demand the use of powerful computer clusters for large-scale applications. In this work, we present a general framework for Accelerated Microstructure Imaging via Convex Optimization (AMICO) and show how to re-formulate this class of techniques as convenient linear systems which, then, can be efficiently solved using very fast algorithms. We demonstrate this linearization of the fitting problem for two specific models, i.e. ActiveAx and NODDI, providing a very attractive alternative for parameter estimation in those techniques; however, the AMICO framework is general and flexible enough to work also for the wider space of microstructure imaging methods. Results demonstrate that AMICO represents an effective means to accelerate the fit of existing techniques drastically (up to four orders of magnitude faster) while preserving accuracy and precision in the estimated model parameters (correlation above 0.9). We believe that the availability of such ultrafast algorithms will help to accelerate the spread of microstructure imaging to larger cohorts of patients and to study a wider spectrum of neurological disorders.

Somewhere between illusion of control and powerlessness : trying to situate the pathological gambler's locus of control

It has been suggested that pathological gamblers develop illusory perceptions of control regarding the outcome of the games and should express higher Internal and Chance locus of control. A sample of 48 outpatients diagnosed with pathological gambling disorder who participated in this ex post facto study, completed the Internality, Powerful Others, and Chance scale, the South Oaks Gambling Screen questionnaire, and the Beck Depression Inventory. Results for the locus of control measure were compared with a reference group. Pathological gamblers scored higher than the reference group on the Chance locus of control, which increased with the severity of cases. Moreover, Internal locus of control did show a curvilinear relationship with the severity of cases. Pathological gamblers have specific locus of control scores that vary in function of the severity, in a linear fashion or a non-linear fashion according to the scale. This effect might be caused by competition between "illusion of control" and the tendency to attribute adverse consequence of gambling to external causes.

Estimating the time since discharge of spent cartridges: a logical approach fro interpreting the evidence

Estimating the time since discharge of a spent cartridge or a firearm can be useful in criminal situa-tions involving firearms. The analysis of volatile gunshot residue remaining after shooting using solid-phase microextraction (SPME) followed by gas chromatography (GC) was proposed to meet this objective. However, current interpretative models suffer from several conceptual drawbacks which render them inadequate to assess the evidential value of a given measurement. This paper aims to fill this gap by proposing a logical approach based on the assessment of likelihood ratios. A probabilistic model was thus developed and applied to a hypothetical scenario where alternative hy-potheses about the discharge time of a spent cartridge found on a crime scene were forwarded. In order to estimate the parameters required to implement this solution, a non-linear regression model was proposed and applied to real published data. The proposed approach proved to be a valuable method for interpreting aging-related data.