Predicting cancer involvement of genes from heterogeneous data

Background: Systematic approaches for identifying proteins involved in different types of cancer are needed. Experimental techniques such as microarrays are being used to characterize cancer, but validating their results can be a laborious task. Computational approaches are used to prioritize between genes putatively involved in cancer, usually based on further analyzing experimental data. Results: We implemented a systematic method using the PIANA software that predicts cancer involvement of genes by integrating heterogeneous datasets. Specifically, we produced lists of genes likely to be involved in cancer by relying on: (i) protein-protein interactions; (ii) differential expression data; and (iii) structural and functional properties of cancer genes. The integrative approach that combines multiple sources of data obtained positive predictive values ranging from 23% (on a list of 811 genes) to 73% (on a list of 22 genes), outperforming the use of any of the data sources alone. We analyze a list of 20 cancer gene predictions, finding that most of them have been recently linked to cancer in literature. Conclusion: Our approach to identifying and prioritizing candidate cancer genes can be used to produce lists of genes likely to be involved in cancer. Our results suggest that differential expression studies yielding high numbers of candidate cancer genes can be filtered using protein interaction networks.

Updated survival data of the phase iii clinical trial of novottf-100a versus best standard chemotherapy for recurrent glioblastoma

NovoTTF-100A (TTF) is a portable device delivering low-intensity, intermediate-frequency, alternating electric fields using noninvasive, disposable scalp electrodes. TTF interferes with tumor cell division, and it has been approved by the US Food and Drug Administration (FDA) for the treatment of recurrent glioblastoma (rGBM) based on data from a phase III trial. This presentation describes the updated survival data 2 years after completing recruitment. Adults with rGBM (KPS ≥ 70) were randomized (stratified by surgery and center) to either continuous TTF (20-24 h/day, 7 days/week) or efficacious chemotherapy based on best physician choice (BPC). The primary endpoint was overall survival (OS), and secondary endpoints were PFS6, 1-year survival, and QOL. Patients were randomized (28 US and European centers) to either TTF alone (n ¼ 120) or BPC (n ¼ 117). Patient characteristics were balanced, median age was 54 years (range, 23-80 years), and median KPS was 80 (range, 50-100). One quarter of the patients had debulking surgery, and over half of the patients were at their second or later recurrence. OS in the intent-to-treat (ITT) population was equivalent in TTF versus BPC patients (median OS, 6.6vs. 6.0 months; n ¼ 237; p ¼ 0.26; HR ¼ 0.86). With a median follow-up of 33.6 months, long-term survival in the TTF group was higher than that in the BPC group at 2, 3, and 4 years of follow-up (9.3% vs. 6.6%; 8.4% vs. 1.4%; 8.4% vs. 0.0%, respectively). Analysis of patients who received at least one treatment course demonstrated a survival benefit for TTF patients compared to BPC patients (median OS, 7.8 vs. 6.0 months; n ¼ 93 vs. n ¼ 117; p ¼ 0.012; HR ¼ 0.69). In this group, 1-year survival was 28% vs. 20%, and PFS6 was 26.2% vs. 15.2% (p ¼ 0.034). TTF, a noninvasive, novel cancer treatment modality shows significant therapeutic efficacy with promising long-term survival results. The impact of TTF was more pronounced when comparing only patients who received the minimal treatment course. A large-scale phase III trial in newly diagnosed GBM is ongoing.

Molecular phylogeny and evolution of Sorex shrews (Soricidae: insectivora) inferred from mitochondrial DNA sequence data.

Shrews of the genus Sorex are characterized by a Holarctic distribution, and relationships among extant taxa have never been fully resolved. Phylogenies have been proposed based on morphological, karyological, and biochemical comparisons, but these analyses often produced controversial and contradictory results. Phylogenetic analyses of partial mitochondrial cytochrome b gene sequences (1011 bp) were used to examine the relationships among 27 Sorex species. The molecular data suggest that Sorex comprises two major monophyletic lineages, one restricted mostly to the New World and one with a primarily Palearctic distribution. Furthermore, several sister-species relationships are revealed by the analysis. Based on the split between the Soricinae and Crocidurinae subfamilies, we used a 95% confidence interval for both the calibration of a molecular clock and the subsequent calculation of major diversification events within the genus Sorex. Our analysis does not support an unambiguous acceleration of the molecular clock in shrews, the estimated rate being similar to other estimates of mammalian mitochondrial clocks. In addition, the data presented here indicate that estimates from the fossil record greatly underestimate divergence dates among Sorex taxa.

@neurIST infrastructure for advanced disease management through integration of heterogeneous data, computing, and complex processing services

The increasing volume of data describing humandisease processes and the growing complexity of understanding, managing, and sharing such data presents a huge challenge for clinicians and medical researchers. This paper presents the@neurIST system, which provides an infrastructure for biomedical research while aiding clinical care, by bringing together heterogeneous data and complex processing and computing services. Although @neurIST targets the investigation and treatment of cerebral aneurysms, the system’s architecture is generic enough that it could be adapted to the treatment of other diseases.Innovations in @neurIST include confining the patient data pertaining to aneurysms inside a single environment that offers cliniciansthe tools to analyze and interpret patient data and make use of knowledge-based guidance in planning their treatment. Medicalresearchers gain access to a critical mass of aneurysm related data due to the system’s ability to federate distributed informationsources. A semantically mediated grid infrastructure ensures that both clinicians and researchers are able to seamlessly access andwork on data that is distributed across multiple sites in a secure way in addition to providing computing resources on demand forperforming computationally intensive simulations for treatment planning and research.

What is the value of joint processing of pilots and data in block-fading channels?

The spectral efficiency achievable with joint processing of pilot and data symbol observations is compared with that achievable through the conventional (separate) approach of first estimating the channel on the basis of the pilot symbols alone, and subsequently detecting the datasymbols. Studied on the basis of a mutual information lower bound, joint processing is found to provide a non-negligible advantage relative to separate processing, particularly for fast fading. It is shown that, regardless of the fading rate, only a very small number of pilot symbols (at most one per transmit antenna and per channel coherence interval) shouldbe transmitted if joint processing is allowed.

What is "clinical data"? Why and how can they be collected during field surveys on medicinal plants?

ETHNOPHARMACOLOGICAL RELEVANCE: "Reverse pharmacology", also called "bedside-to-bench" or "field to pharmacy" approach, is a research process starting with documentation of clinical outcome as observed by patients with different therapeutic regimens. The treatment most significantly associated with cure is selected for future studies: first, clinical safety and efficacy; then in vivo and vitro studies. Some clinical data, i.e. details on patient status and progress, can be collected during ethnobotanical surveys; they will help clinical researchers and, once effectiveness and safety are established, will also help users of traditional medicine make safer and more effective choices. To gather clinical data successfully, ethnopharmacologists need to be backed by an appropriate team of specialists in medicine and epidemiology. Ethnopharmacologists can also gather important data on traditional medicine safety. MATERIALS AND METHODS: The first step is to create a consensus on the meaning of "clinical data", their interest and importance. An understanding of why "a cure is not a proof of effectiveness" is a starting point to avoid faulty interpretation of the clinical observations. RESULTS: Experience showed that, with the "bedside-to-bench" approach, a treatment derived from traditional recipe can be scientifically validated (in terms of safety and effectiveness) with a cost of less than a million euros, thus providing an end-product that is affordable, available and sustainable. CONCLUSIONS: With rigorous clinical study results, medicinal plant users gain the possibility to refine heath strategies. The field surveyor may gain a better relationship with the population, once she/he is seen as bringing information useful for the quality of care in the community.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Field Experiments of Current Concrete Pavement Surface Characteristics Practices: Iowa Data Collection and Analysis, December 2005

One of the most important issues in portland cement concrete pavement research today is surface characteristics. The issue is one of balancing surface texture construction with the need for durability, skid resistance, and noise reduction. The National Concrete Pavement Technology Center at Iowa State University, in conjunction with the Federal Highway Administration, American Concrete Pavement Association, International Grinding and Grooving Association, Iowa Highway Research Board, and other states, have entered into a three-part National Surface Characteristics Program to resolve the balancing problem. As a portion of Part 2, this report documents the construction of 18 separate pavement surfaces for use in the first level of testing for the national project. It identifies the testing to be done and the limitations observed in the construction process. The results of the actual tests will be included in the subsequent national study reports.

Estimating the basic reproductive number from viral sequence data.

Epidemiological processes leave a fingerprint in the pattern of genetic structure of virus populations. Here, we provide a new method to infer epidemiological parameters directly from viral sequence data. The method is based on phylogenetic analysis using a birth-death model (BDM) rather than the commonly used coalescent as the model for the epidemiological transmission of the pathogen. Using the BDM has the advantage that transmission and death rates are estimated independently and therefore enables for the first time the estimation of the basic reproductive number of the pathogen using only sequence data, without further assumptions like the average duration of infection. We apply the method to genetic data of the HIV-1 epidemic in Switzerland.

The Exotic Palaeo-Thethys Terrane in Central Iran: New geological data from Anarak, Jandaq and Posht-e-Badam areas

Résumé Le « terrane » d'Anarak-Jandak occupe une position géologique clé au nord-ouest du Microcontinent Centre-East Iranien (CE1M), connecté avec le Bloc du Grand Kavir et la ceinture métamorphique de Sanandaj-Sirjan. Nous discutons ici l'origine de ces différentes unités, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, pour conclure finalement de leur affinité paléotéthysienne. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur-Dévonien inférieur, pour se terminer au Trias par la collision des blocs Cimmériens dérivé du Gondwana avec le Bloc du Turan d'affinité asiatique (événement Eocimmérien). La plus importante unité métamorphique affleurant au sud-ouest de la région de Jandak-Anarak-Kaboudan est une épaisse séquence silicoclastique à grains fins contenant des blocs ophiolitiques (marginal-sea-type), et des associations basalte-gabbro à signatures géochimiques de type supra-subduction. Dans la région de Nakhlak, nous avons daté ces gabbros par la méthode U-Pb à 387f0.11 Ma ; les roches métamorphiques pélitiques ont donné des âges de refroidissement Ar-Ar pour la muscovite de 320 à 333 Ma. Ce complexe d'accrétion "varisque" a été métamorphisé dans le faciès schiste vert-amphibolite au cours de l'accrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma par la méthode U/Pb), qui affleure aujourd'hui à l'extrémité nord-ouest du terrane d'Anarak-Jandak . La subduction vers le nord de l'océan Paléotéthys depuis le Paléazoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de guyots (Anarak, Kaboudan, et Meraji Seamounts) et de hauts sous-marins, entrés en collision oblique avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries {âges Ar-Ar de 280 à 230 Ma). De plus, le magmatisme bimodal de Chah Gorbeh est caractérisé d'une part par des roches de type trondjémite-gabbros (262 Ma), d'autre part par des laves en coussin de type basaltes alcalins-rhyolites; ces roches magmatiques ont recoupé l'ophiolite d'Anarak lors de la mise en place de cette dernière dans la fosse interne de subduction. Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, i1 a été accrété le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge probable Triasique. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé dans les dépôts infra-arc Dévonien supérieur-Carbonifère de Godar-e-Siah, ainsi que dans la succession d'avant-arc de Nakhlak. Pendant l'intervalle Paléozoïque supérieur-Trias, la région de Jandak a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, elle-même comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites d'arc à collisionnel datés à 215±15 Ma. Dans la région de Yazd, témoin de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur; il en a été de même pour les dépôts de plate-forme Paléozoïque supérieur. L'érosion, qui dans ce dernier cas a atteint le Permien, pourrait être liée au bombement flexural de la marge passive. La collision finale n'a pas induit de déformations trop importantes, et se caractérise par la mise en place de nappes sur la marge passive. Cet événement est scellé par des dépôts molassique du Lias. D'un point de vue régional, la zone s'étendant actuellement de la Mer Noire au Pamir a été soumise à six épisodes d'extension-compression du Jurassique inférieur (début du l'ouverture en position arrière-arc de la Néotéthys) à l'Eocène moyen. Par exemple, le terrane d'AnarakJandak, probablement situé entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de sa patrie d'origine au début du Crétacé supérieur. Des preuves de cet événement se retrouvent dans les séries de plate-forme de Khur (préservation de séries syn-rift puis de marge passive). Les ophiolites de Nain et de Sabzevar sont de plus interprétée comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation par la plaque indienne de l'Eurasie a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent Iranien et de la formation du CEIM. Cette rotation est responsable du transport du terrane d'Anarak-Jandak vers sa position actuelle en Iran Central, et de la dislocation de Terranes de moindre importance, comme le bloc de Posht-e Badam. Depuis le Miocène supérieur, et à la suite de la collision entre l'Arabie et l'Iran, le ternane d'Anarak-Jandak a subi des déformations liées à l'activité d'une zone de cisaillement dextre parallèle à la suture du Zagros, à l'arrière de l'arc magmatique d'Uromieh-Dokhtar. Résumé large public Le Microcontinent Centre-Est Iranien occupe une position géologique clé au centre de l'Iran. Les différentes unités qui le composent, reliées jusqu'à présent à des épisodes orogéniques d'âge Précambrien à Paléozoïque inférieur, sont maintenant rajeunies et liés à la fermeture de l'océean Paléotéthys. Leur histoire commence par un épisode de rifting d'âge Ordovicien supérieur à Dévonien inférieur, pour se terminer au Trias par la collision des- blocs Cimmériens, dérivés du Gondwana, avec le Bloc du Turan d'affinité asiatique. Dans la marge active asiatique de la Paléotéthys, nous avons daté les restes d'un océan marginal à 387±0.11 Ma. Ce complexe d'accrétion a été métamorphisé au cours de la réaccrétion de la ceinture granitique d'Airekan, d'âge Cambrien inférieur (549±15 Ma), qui affleure aujourd'hui à l'extrémité nord-ouest du « terrane » d'Anarak-Jandak correspondant à la plus grande partie de la région étudiée. Le matériel issu de l'arc magmatique de la Paléotéthys est très bien préservé et daté du Dévonien supérieur-Carbonifère. Pendant l'intervalle Paléozoïque supérieur-Trias, la région a été soumise à un régime extensif de type bassin d'arrière-arc, dont un témoin pourrait être la ceinture ophiolitique d'Arusan, comparable aux écailles ophiolitiques d'Aghdarband au nord-est de l'Iran. Cet ensemble métamorphique est recoupé par des granites datés à 215±15 Ma. La subduction vers le nord de l'océan Paléotéthys depuis le Paléozoïque supérieur jusqu'au Trias, a permis l'accumulation de grandes quantités de matériel océanique dans la zone de subduction. Par exemple, une succession de volcans sous-marins, entrés en collision avec le prisme d'accrétion, est à l'origine d'un léger métamorphisme de type HP qui affecte ces séries (280 à 230 Ma). Quant au prisme d'accrétion de Doshakh, d'âge essentiellement Permien supérieur, il a été mis en place le long de la marge continentale et métamorphisé dans le faciès schiste vert. La fermeture de la Paléotéthys s'enregistre finalement par la sédimentation dans le bassin d'avant pays du flysch de Bayazeh, d'âge Triasique. Dans la région de Yazd, on trouve les témoins de la marge passive Cimmérienne, la sédimentation syn-rift Silurienne à Dévonienne inférieure a été interrompue pendant l'intervalle Trias moyen-Trias supérieur, marqué par la flexuration de la marge passive lorsqu'elle rentra en collision avec la marge active asiatique. Cet événement est scellé par des dépôts molassique à charbon du Lias. Le «terrane» d'Anarak-Jandak, probablement situé à l'origine entre le Kopeh Dagh et la plate-forme nord Afghane, s'est complètement détaché de cette région au début du Crétacé supérieur lors de l'ouverture d'un bassin d'arrière-arc, engendré, cette fois, par la subduction de l'océan Néotéthys situé au sud des blocs cimmériens. Des preuves de cet événement se retrouvent dans les séries syn-rift, puis de marge passive de Khour. Les ophiolites de Nain et de Sabzevar sont interprétées comme un témoin de l'existence de ce bassin d'arrière-arc. Dans l'intervalle Eocène-Oligocène, l'indentation de l'Eurasie par la plaque indienne a été contemporaine de la rotation horaire de fragments de l'ancien microcontinent centre-Iranien. Cette rotation de près de 90° est responsable du transport du « terrane » d'Anarak-Jandak vers sa position actuelle. Abstract The Anarak-Jandaq terrane occupies a strategic geological situation at the north-western part of the Central-East Iranian Microcontinent (CEIM) and in connection with the Great Kavir Block and Sanandaj-Sirjan metamorphic belt. Our recent findings redefine the origin of these mentioned areas so far attributed to the Precambrian-Early Palaeozoic orogenic episodes, to be now directly related to the tectonic evolution of the Palaeo-Tethys Ocean, commenced by Late Ordovician-Early Devonian rifting events and terminated in the Triassic by the Eocimmerian tectonic event due to the collision of the Cimmerian blocks with the Asiatic Turan block. The most distributed metamorphic unit that is exposed from the south-west of Jandaq to the Anarak and Kaboudan areas is a thick and fine grain siliciclastic sequence accompanied by marginal-sea-basin ophiolitic blocks including basalt-gabbro association with supra-subduction-geochemical signature. These gabbros in the Nakhlak area were dated by U/Pb method at 387.6 ± 0.11 Ma and the metamorphic pelitic rocks yielded a range of 320 to 333 Ma muscovite-cooling ages based on 40Ar/39 Ar method. This "Variscan" accretionary complex was metamorphosed in greenschist-amphibolite facies during accretion to the Lower Cambrian Airekan granitic belt (549 ± 15 Ma by U/Pb method) that crops out at the northwestern edge of the Anarak-Jandaq terrane. Continued northward subduction of the Palaeo-Tethys Ocean during the entire Late Palaeozoic-Middle Triassic brought huge amount of oceanic material to the subduction zone. One chain of Carboniferous-Triassic oceanic rises and seamounts (the Anarak, Kaboudan, and Meraji Seamounts) obliquely collided with the accretionary wedge and created a mild HP metamorphic event (280-230 Ma based on 40Ar/39Ar results). Bimodal magmatism of the Chah Gorbeh area is characterized by a 262 Ma trondjemite-gabbro as well as pillow alkalibasalts-rhyolites which intruded the Anarak ophiolite when it was being emplaced within the inner-wall trench. The mainly Late Permian-Triassic Doshakh wedge was accreted along the continent and metamorphosed under lower greenschist facies and the probable Triassic Bayazeh flysch filled the foreland basin during the final closure. The Palaeo-Tethys magmatic arc products have been well preserved in the Late Devonian-Carboniferous Godar-e-Siah intra-arc deposits and the Triassic Nakhlak fore-arc succession. During the Late Palaeozoic-Triassic times, the Jandaq area has been affected by back-arc extension and probably the Arusan ophiolitic belt is the remnant of this narrow basin comparable to the Aqdarband ophiolitic remnant in north-east Iran. This metamorphic belt was intruded by 215 ± 15 Ma arc to collisional granites. In the passive margin of the Cimmerian block, on the Yazd region, the Silurian-Early Devonian syn-rift succession as well as the nearly continuous Upper Palaeozoic platform-type deposition was interrupted during the Middle to Late Triassic time, local erosion down to Devonian levels may be related to flexural bulge erosion. The collision event was not so strong to generate intensive deformation but was accompanied by some nappe thrusting onto the passive margin. It is finally unconformably covered by Liassic continental molassic deposits. Related to the onset of Neo-Tethyan back-arc opening in Early Jurassic to Mid-Eocene times, six periods of extensional-compressional events have differently influenced an elongated area, extending from the West Black Sea to Pamir. The Anarak-Jandaq terrane which was situated somewhere in this affected area, probably between the Kopeh Dagh and North Afghan platform, was completely detached from its source at the beginning of the Late Cretaceous

Focused time-lapse inversion of radio and audio magnetotelluric data

Geoelectrical techniques are widely used to monitor groundwater processes, while surprisingly few studies have considered audio (AMT) and radio (RMT) magnetotellurics for such purposes. In this numerical investigation, we analyze to what extent inversion results based on AMT and RMT monitoring data can be improved by (1) time-lapse difference inversion; (2) incorporation of statistical information about the expected model update (i.e., the model regularization is based on a geostatistical model); (3) using alternative model norms to quantify temporal changes (i.e., approximations of l(1) and Cauchy norms using iteratively reweighted least-squares), (4) constraining model updates to predefined ranges (i.e., using Lagrange Multipliers to only allow either increases or decreases of electrical resistivity with respect to background conditions). To do so, we consider a simple illustrative model and a more realistic test case related to seawater intrusion. The results are encouraging and show significant improvements when using time-lapse difference inversion with non l(2) model norms. Artifacts that may arise when imposing compactness of regions with temporal changes can be suppressed through inequality constraints to yield models without oscillations outside the true region of temporal changes. Based on these results, we recommend approximate l(1)-norm solutions as they can resolve both sharp and smooth interfaces within the same model. (C) 2012 Elsevier B.V. All rights reserved.

Measures of fit in multiple correspondence analysis of crisp and fuzzy coded data

When continuous data are coded to categorical variables, two types of coding are possible: crisp coding in the form of indicator, or dummy, variables with values either 0 or 1; or fuzzy coding where each observation is transformed to a set of "degrees of membership" between 0 and 1, using co-called membership functions. It is well known that the correspondence analysis of crisp coded data, namely multiple correspondence analysis, yields principal inertias (eigenvalues) that considerably underestimate the quality of the solution in a low-dimensional space. Since the crisp data only code the categories to which each individual case belongs, an alternative measure of fit is simply to count how well these categories are predicted by the solution. Another approach is to consider multiple correspondence analysis equivalently as the analysis of the Burt matrix (i.e., the matrix of all two-way cross-tabulations of the categorical variables), and then perform a joint correspondence analysis to fit just the off-diagonal tables of the Burt matrix - the measure of fit is then computed as the quality of explaining these tables only. The correspondence analysis of fuzzy coded data, called "fuzzy multiple correspondence analysis", suffers from the same problem, albeit attenuated. Again, one can count how many correct predictions are made of the categories which have highest degree of membership. But here one can also defuzzify the results of the analysis to obtain estimated values of the original data, and then calculate a measure of fit in the familiar percentage form, thanks to the resultant orthogonal decomposition of variance. Furthermore, if one thinks of fuzzy multiple correspondence analysis as explaining the two-way associations between variables, a fuzzy Burt matrix can be computed and the same strategy as in the crisp case can be applied to analyse the off-diagonal part of this matrix. In this paper these alternative measures of fit are defined and applied to a data set of continuous meteorological variables, which are coded crisply and fuzzily into three categories. Measuring the fit is further discussed when the data set consists of a mixture of discrete and continuous variables.