Robust Nonlinear Mapping of Soil Contamination Using Support Regression

Spatial data analysis mapping and visualization is of great importance in various fields: environment, pollution, natural hazards and risks, epidemiology, spatial econometrics, etc. A basic task of spatial mapping is to make predictions based on some empirical data (measurements). A number of state-of-the-art methods can be used for the task: deterministic interpolations, methods of geostatistics: the family of kriging estimators (Deutsch and Journel, 1997), machine learning algorithms such as artificial neural networks (ANN) of different architectures, hybrid ANN-geostatistics models (Kanevski and Maignan, 2004; Kanevski et al., 1996), etc. All the methods mentioned above can be used for solving the problem of spatial data mapping. Environmental empirical data are always contaminated/corrupted by noise, and often with noise of unknown nature. That's one of the reasons why deterministic models can be inconsistent, since they treat the measurements as values of some unknown function that should be interpolated. Kriging estimators treat the measurements as the realization of some spatial randomn process. To obtain the estimation with kriging one has to model the spatial structure of the data: spatial correlation function or (semi-)variogram. This task can be complicated if there is not sufficient number of measurements and variogram is sensitive to outliers and extremes. ANN is a powerful tool, but it also suffers from the number of reasons. of a special type ? multiplayer perceptrons ? are often used as a detrending tool in hybrid (ANN+geostatistics) models (Kanevski and Maignank, 2004). Therefore, development and adaptation of the method that would be nonlinear and robust to noise in measurements, would deal with the small empirical datasets and which has solid mathematical background is of great importance. The present paper deals with such model, based on Statistical Learning Theory (SLT) - Support Vector Regression. SLT is a general mathematical framework devoted to the problem of estimation of the dependencies from empirical data (Hastie et al, 2004; Vapnik, 1998). SLT models for classification - Support Vector Machines - have shown good results on different machine learning tasks. The results of SVM classification of spatial data are also promising (Kanevski et al, 2002). The properties of SVM for regression - Support Vector Regression (SVR) are less studied. First results of the application of SVR for spatial mapping of physical quantities were obtained by the authorsin for mapping of medium porosity (Kanevski et al, 1999), and for mapping of radioactively contaminated territories (Kanevski and Canu, 2000). The present paper is devoted to further understanding of the properties of SVR model for spatial data analysis and mapping. Detailed description of the SVR theory can be found in (Cristianini and Shawe-Taylor, 2000; Smola, 1996) and basic equations for the nonlinear modeling are given in section 2. Section 3 discusses the application of SVR for spatial data mapping on the real case study - soil pollution by Cs137 radionuclide. Section 4 discusses the properties of the modelapplied to noised data or data with outliers.

Micro-ordinateurs de poche : du gadget à l'irremplaçable [Handheld computers: the gadget that became an invaluable tool]

Les ordinateurs de poche ont le potentiel de devenir un outil de travail indispensable pour le médecin praticien. De part une facilité d'utilisation toujours plus grande, un format ultra-portable, une grande capacité de stockage de l'information et de synchronisation à de multiples bases de données, ils permettent l'accès en tout temps et littéralement au lit du malade à une gamme complète de ressources médicales. Cet article se veut une lllustration des différentes possibilités d'utilisation de ces nouveaux outils, au travers de situations cliniques banales rencontrées par tout praticien.

MRS glucose mapping and PET joining forces: re-evaluation of the lumped constant in the rat brain under isoflurane anaesthesia.

Although numerous positron emission tomography (PET) studies with (18) F-fluoro-deoxyglucose (FDG) have reported quantitative results on cerebral glucose kinetics and consumption, there is a large variation between the absolute values found in the literature. One of the underlying causes is the inconsistent use of the lumped constants (LCs), the derivation of which is often based on multiple assumptions that render absolute numbers imprecise and errors hard to quantify. We combined a kinetic FDG-PET study with magnetic resonance spectroscopic imaging (MRSI) of glucose dynamics in Sprague-Dawley rats to obtain a more comprehensive view of brain glucose kinetics and determine a reliable value for the LC under isoflurane anaesthesia. Maps of Tmax /CMRglc derived from MRSI data and Tmax determined from PET kinetic modelling allowed to obtain an LC-independent CMRglc . The LC was estimated to range from 0.33 ± 0.07 in retrosplenial cortex to 0.44 ± 0.05 in hippocampus, yielding CMRglc between 62 ± 14 and 54 ± 11 μmol/min/100 g, respectively. These newly determined LCs for four distinct areas in the rat brain under isoflurane anaesthesia provide means of comparing the growing amount of FDG-PET data available from translational studies.

Interest rates mapping

The present study deals with the analysis and mapping of Swiss franc interest rates. Interest rates depend on time and maturity, defining term structure of the interest rate curves (IRC). In the present study IRC are considered in a two-dimensional feature space - time and maturity. Exploratory data analysis includes a variety of tools widely used in econophysics and geostatistics. Geostatistical models and machine learning algorithms (multilayer perceptron and Support Vector Machines) were applied to produce interest rate maps. IR maps can be used for the visualisation and pattern perception purposes, to develop and to explore economical hypotheses, to produce dynamic asset-liability simulations and for financial risk assessments. The feasibility of an application of interest rates mapping approach for the IRC forecasting is considered as well. (C) 2008 Elsevier B.V. All rights reserved.

Detection and mapping of a lethal locus in a eucalyptus hybrid population

The objective of this work was to verify the existence of a lethal locus in a eucalyptus hybrid population, and to quantify the segregation distortion in the linkage group 3 of the Eucalyptus genome. A E. grandis x E. urophylla hybrid population, which segregates for rust resistance, was genotyped with 19 microsatellite markers belonging to linkage group 3 of the Eucalyptus genome. To quantify the segregation distortion, maximum likelihood (ML) models, specific to outbreeding populations, were used. These models consider the observed marker genotypes and the lethal locus viability as parameters. The ML solutions were obtained using the expectation‑maximization algorithm. A lethal locus in the linkage group 3 was verified and mapped, with high confidence, between the microssatellites EMBRA 189 e EMBRA 122. This lethal locus causes an intense gametic selection from the male side. Its map position is 25 cM from the locus which controls the rust resistance in this population.

Detection and chemical profiling of medicine counterfeits by Raman spectroscopy and chemometrics

Raman spectroscopy combined with chemometrics has recently become a widespread technique for the analysis of pharmaceutical solid forms. The application presented in this paper is the investigation of counterfeit medicines. This increasingly serious issue involves networks that are an integral part of industrialized organized crime. Efficient analytical tools are consequently required to fight against it. Quick and reliable authentication means are needed to allow the deployment of measures from the company and the authorities. For this purpose a method in two steps has been implemented here. The first step enables the identification of pharmaceutical tablets and capsules and the detection of their counterfeits. A nonlinear classification method, the Support Vector Machines (SVM), is computed together with a correlation with the database and the detection of Active Pharmaceutical Ingredient (API) peaks in the suspect product. If a counterfeit is detected, the second step allows its chemical profiling among former counterfeits in a forensic intelligence perspective. For this second step a classification based on Principal Component Analysis (PCA) and correlation distance measurements is applied to the Raman spectra of the counterfeits.

The role of interval nodes in sentinel lymph node mapping and dissection for melanoma patients.

In sentinel node (SN) biopsy, an interval SN is defined as a lymph node or group of lymph nodes located between the primary melanoma and an anatomically well-defined lymph node group directly draining the skin. As shown in previous reports, these interval SNs seem to be at the same metastatic risk as are SNs in the usual, classic areas. This study aimed to review the incidence, lymphatic anatomy, and metastatic risk of interval SNs. METHODS: SN biopsy was performed at a tertiary center by a single surgical team on a cohort of 402 consecutive patients with primary melanoma. The triple technique of localization was used-that is, lymphoscintigraphy, blue dye, and gamma-probe. Otolaryngologic melanoma and mucosal melanoma were excluded from this analysis. SNs were examined by serial sectioning and immunohistochemistry. All patients with metastatic SNs were recommended to undergo a radical selective lymph node dissection. RESULTS: The primary locations of the melanomas included the trunk (188), an upper limb (67), or a lower limb (147). Overall, 97 (24.1%) of the 402 SNs were metastatic. Interval SNs were observed in 18 patients, in all but 2 of whom classic SNs were also found. The location of the primary was truncal in 11 (61%) of the 18, upper limb in 5, and lower limb in 2. One patient with a dorsal melanoma had drainage exclusively in a cervicoscapular area that was shown on removal to contain not lymph node tissue but only a blue lymph channel without tumor cells. Apart from the interval SN, 13 patients had 1 classic SN area and 3 patients 2 classic SN areas. Of the 18 patients, 2 had at least 1 metastatic interval SN and 2 had a classic SN that was metastatic; overall, 4 (22.2%) of 18 patients were node-positive. CONCLUSION: We found that 2 of 18 interval SNs were metastatic: This study showed that preoperative lymphoscintigraphy must review all known lymphatic areas in order to exclude an interval SN.

Bone mineral density (BMD), micro-architecture estimation (TBS) and vertebral fracture assessment (VFA) extracted from a single DXA device in combination with clinical risk factors improve significantly the identification of women at high risk of fracture.

Osteoporosis (OP) is a systemic skeletal disease characterized by a low bone mineral density (BMD) and a micro-architectural (MA) deterioration. Clinical risk factors (CRF) are often used as a MA approximation. MA is yet evaluable in daily practice by the trabecular bone score (TBS) measure. TBS is very simple to obtain, by reanalyzing a lumbar DXA-scan. TBS has proven to have diagnosis and prognosis values, partially independent of CRF and BMD. The aim of the OsteoLaus cohort is to combine in daily practice the CRF and the information given by DXA (BMD, TBS and vertebral fracture assessment (VFA)) to better identify women at high fracture risk. The OsteoLaus cohort (1400 women 50 to 80 years living in Lausanne, Switzerland) started in 2010. This study is derived from the cohort COLAUS who started in Lausanne in 2003. The main goal of COLAUS is to obtain information on the epidemiology and genetic determinants of cardiovascular risk in 6700 men and women. CRF for OP, bone ultrasound of the heel, lumbar spine and hip BMD, VFA by DXA and MA evaluation by TBS are recorded in OsteoLaus. Preliminary results are reported. We included 631 women: mean age 67.4 ± 6.7 years, BMI 26.1 ± 4.6, mean lumbar spine BMD 0.943 ± 0.168 (T-score − 1.4 SD), and TBS 1.271 ± 0.103. As expected, correlation between BMD and site matched TBS is low (r2 = 0.16). Prevalence of VFx grade 2/3, major OP Fx and all OP Fx is 8.4%, 17.0% and 26.0% respectively. Age- and BMI-adjusted ORs (per SD decrease) are 1.8 (1.2-2.5), 1.6 (1.2-2.1), and 1.3 (1.1-1.6) for BMD for the different categories of fractures and 2.0 (1.4-3.0), 1.9 (1.4-2.5), and 1.4 (1.1-1.7) for TBS respectively. Only 32 to 37% of women with OP Fx have a BMD < − 2.5 SD or a TBS < 1.200. If we combine a BMD < − 2.5 SD or a TBS < 1.200, 54 to 60% of women with an osteoporotic Fx are identified. As in the already published studies, these preliminary results confirm the partial independence between BMD and TBS. More importantly, a combination of TBS subsequent to BMD increases significantly the identification of women with prevalent OP Fx which would have been misclassified by BMD alone. For the first time we are able to have complementary information about fracture (VFA), density (BMD), micro- and macro architecture (TBS and HAS) from a simple, low ionizing radiation and cheap device: DXA. Such complementary information is very useful for the patient in the daily practice and moreover will likely have an impact on cost effectiveness analysis.

Evaluating spine micro-architectural texture (via TBS) discriminates major osteoporotic fractures from controls both as well as and independent of site matched BMD: the Eastern European TBS study.

The aim of the study was to assess the clinical performance of the model combining areal bone mineral density (aBMD) at spine and microarchitecural texture (TBS) for the detection of the osteoporotic fracture. The Eastern European Study is a multicenter study (Serbia, Bulgaria, Romania and Ukraine) evaluating the role of TBS in routine clinical practice as a complement to aBMD. All scans were acquired on Hologic Discovery and GE Prodigy densitometers in a routine clinical manner. The additional clinical values of aBMD and TBS were analyzed using a two steps classification tree approach (aBMD followed by TBS tertiles) for all type of osteoporotic fracture (All-OP Fx). Sensitivity, specificity and accuracy of fracture detection as well as the Net Reclassification Index (NRI) were calculated. This study involves 1031 women subjects aged 45 and older recruited in east European countries. Clinical centers were cross-calibrated in terms of BMD and TBS. As expected, areal BMD (aBMD) at spine and TBS were only moderately correlated (r (2) = 0.19). Prevalence rate for All-OP Fx was 26 %. Subjects with fracture have significant lower TBS and aBMD than subjects without fracture (p < 0.01). TBS remains associated with the fracture even after adjustment for age and aBMD with an OR of 1.27 [1.07-1.51]. When using aBMD T-score of -2.5 and the lowest TBS tertile thresholds, both BMD and TBS were similar in terms of sensitivity (35 vs. 39 %), specificity (78 vs. 80 %) and accuracy (64 vs. 66 %). aBMD and TBS combination, induced a significant improvement in sensitivity (+28 %) and accuracy (+17 %) compared to aBMD alone whereas a moderate improvement was observed in terms of specificity (+9 %). The overall combination gain was 36 % as expressed using the NRI. aBMD and TBS combination decrease significantly the number of subjects needed to diagnose from 7 for aBMD alone to 2. In a multi-centre Eastern European cohort, we have shown that the use of TBS in addition to the aBMD permit to reclassified correctly more than one-third of the overall subjects. Furthermore, the number of subjects needed to diagnose fell to 2 subjects. Economical studies have to be performed to evaluate the gain induced by the use of TBS for the healthcare system.

Quantitative multi-parameter mapping of R1, PD(*), MT, and R2(*) at 3T: a multi-center validation.

Multi-center studies using magnetic resonance imaging facilitate studying small effect sizes, global population variance and rare diseases. The reliability and sensitivity of these multi-center studies crucially depend on the comparability of the data generated at different sites and time points. The level of inter-site comparability is still controversial for conventional anatomical T1-weighted MRI data. Quantitative multi-parameter mapping (MPM) was designed to provide MR parameter measures that are comparable across sites and time points, i.e., 1 mm high-resolution maps of the longitudinal relaxation rate (R1 = 1/T1), effective proton density (PD(*)), magnetization transfer saturation (MT) and effective transverse relaxation rate (R2(*) = 1/T2(*)). MPM was validated at 3T for use in multi-center studies by scanning five volunteers at three different sites. We determined the inter-site bias, inter-site and intra-site coefficient of variation (CoV) for typical morphometric measures [i.e., gray matter (GM) probability maps used in voxel-based morphometry] and the four quantitative parameters. The inter-site bias and CoV were smaller than 3.1 and 8%, respectively, except for the inter-site CoV of R2(*) (<20%). The GM probability maps based on the MT parameter maps had a 14% higher inter-site reproducibility than maps based on conventional T1-weighted images. The low inter-site bias and variance in the parameters and derived GM probability maps confirm the high comparability of the quantitative maps across sites and time points. The reliability, short acquisition time, high resolution and the detailed insights into the brain microstructure provided by MPM makes it an efficient tool for multi-center imaging studies.

ADC mapping of chronic cerebral hypoperfusion induced by carotid artery stenosis.

BACKGROUND: Carotid artery stenosis is associated with the occurrence of acute and chronic ischemic lesions that increase with age in the elderly population. Diffusion Imaging and ADC mapping may be an appropriate method to investigate patients with chronic hypoperfusion consecutive to carotid stenosis. This non-invasive technique allows to investigate brain integrity and structure, in particular hypoperfusion induced by carotid stenosis diseases. The aim of this study was to evaluate the impact of a carotid stenosis on the parenchyma using ADC mapping. METHODS: Fifty-nine patients with symptomatic (33) and asymptomatic (26) carotid stenosis were recruited from our multidisciplinary consultation. Both groups demonstrated a similar degree of stenosis. All patients underwent MRI of the brain including diffusion-weighted MR imaging with ADC mapping. Regions of interest were defined in the anterior and posterior paraventricular regions both ipsilateral and contralateral to the stenosis (anterior circulation). The same analysis was performed for the thalamic and occipital regions (posterior circulation). RESULTS: ADC values of the affected vascular territory were significantly higher on the side of the stenosis in the periventricular anterior (P<0.001) and posterior (P<0.01) area. There was no difference between ipsilateral and contralateral ADC values in the thalamic and occipital regions. CONCLUSIONS: We have shown that carotid stenosis is associated with significantly higher ADC values in the anterior circulation, probably reflecting an impact of chronic hypoperfusion on the brain parenchyma in symptomatic and asymptomatic patients. This is consistent with previous data in the literature.