An approximation algorithm for #k-SAT

"Vegeu el resum a l'inici del document del fitxer adjunt"

Imaging EEG synchronization in schizophrenia patients with S-estimator

Recently a new measure of the cooperative behavior of simultaneous time series was introduced (Carmeli et al. NeuroImage 2005). This measure called S-estimator is defined from the embedding dimension in a state space. S-estimator quantifies the amount of synchronization within a data set by comparing the actual dimensionality of the set with the expected full dimensionality of the asynchronous set. It has the advantage of being a multivariate measure over traditionally used in systems neuroscience bivariate measures of synchronization. Multivariate measures of synchronization are of particular interest for applications in the field of modern multichannel EEG research, since they easily allow mapping of local and/or regional synchronization and are compatible with other imaging techniques. We applied Sestimator to the analysis of EEG synchronization in schizophrenia patients vs. matched controls. The whole-head mapping with S-estimator revealed a specific pattern of local synchronization in schizophrenia patients. The differences in the landscape of synchronization included decreased local synchronization in the territories over occipital and midline areas and increased synchronization over temporal areas. In frontal areas, the S-estimator revealed a tendency for an asymmetry: decreased S-values over the left hemisphere were adjacent to increased values over the right hemisphere. Separate calculations showed reproducibility of this pattern across the main EEG frequency bands. The maintenance of the same synchronization landscape across EEG frequencies probably implies the structural changes in the cortical circuitry of schizophrenia patients. These changes are regionally specific and suggest that schizophrenia is a misconnectivity rather than hypo- or hyper-connectivity disorder.

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.

Abnormal EEG Synchronization Correlates with Demyelination in Alzheimer's Disease

Introduction : The pathological processes caused by Alzheimer's disease (AD) supposedly disrupt communication between and within the distributed cortical networks due to the dysfunction/loss of synapses and myelination breakdown. Indeed, recently (Knyazeva et al. 2008), we have revealed the whole-head topography of EEG synchronization specific to AD. Here we analyze whether and how these abnormalities of synchronization are related to the demyelination of cortico-cortical fibers. Methods : Fifteen newly diagnosed AD patients (CDR 0.5-1) and 15 controls matched for age, participated in the study. Their multichannel (128) EEGs were recorded during 3-5 min at rest. They were submitted to the multivariate phase synchronization (MPS) analysis for mapping regional synchronization. To obtain individual whole-head maps, the MPS was computed for each sensor considering its 2nd nearest topographical neighbors. Separate calculations were performed for the delta, theta, alpha-1/−2, and beta-1/−2 EEG bands. The same subjects were scanned on a 3 Tesla Philips scanner. The protocol included a high-resolution T1-weighted sequence and a Magnetization Transfer Imaging (MTI) acquisition. For each subject, we defined a 3mm thick layer of white matter exactly below the cortical gray matter. The magnetization transfer ratio (MTR) - an estimator of myelination - was calculated for this layer in 39 Brodmann-defined ROIs per hemisphere. To assess the between-group differences, we used a permutation version of Hotelling's T2 test or two-sample T-test (Pcorrected <0.05). For correlation analysis, Spearman Rank Correlation was calculated. Results : In AD patients, we have found an abnormal landscape of synchronization characterized by a decrease in MPS over the fronto-temporal region of the left hemisphere and an increase over the temporo-parieto-occipital regions bilaterally. Also, we have shown a widespread decrease in regional MTR in the AD patients for all the areas excluding motor, premotor, and primary sensory ones. Assuming that AD-related changes in synchronization are associated with demyelination, we hypothesized a correlation between the regional MTR values and MPS values in the hypo- and hyper-synchronized clusters. We found that MPS in the left fronto-temporal hypo-synchronized cluster directly correlates with myelination in BA42-46 of the left hemisphere: the lower the myelination in individual patients, the lower the EEG synchronization. By contrast, in the posterior hyper-synchronized cluster, MPS inversely correlated with myelination, i.e., the lower the myelination, the higher the synchronization. This posterior hyper-synchronization, more characteristic for early-onset AD, probably, results from the initial effect of the disease on cortical inhibition, reducing cortical capacity for decoupling irrelevant connections. Remarkably, it showed different topography of correlations in early- vs. late-onset patients. In the early-onset patients, hyper-synchronization was mainly related to demyelination in posterior BAs, the effect being significant in all the EEG frequency bands. In the late-onset patients, widely distributed correlations were significant for the EEG delta band, suggesting an interaction between the cerebral manifestations of AD and the age of its onset, i.e., topographically selective impairment of cortical inhibition in early-onset AD vs. its wide-spread weakening in old age. Conclusions : Overall, our results document that the degradation of white matter is a significant factor of AD pathogenesis leading to functional dysconnection, the latter being reflected in EEG synchronization abnormalities.

Neonatal and adult ICU ventilators to provide ventilation in neonates, infants, and children: a bench model study.

BACKGROUND: Using a bench test model, we investigated the hypothesis that neonatal and/or adult ventilators equipped with neonatal/pediatric modes currently do not reliably administer pressure support (PS) in neonatal or pediatric patient groups in either the absence or presence of air leaks. METHODS: PS was evaluated in 4 neonatal and 6 adult ventilators using a bench model to evaluate triggering, pressurization, and cycling in both the absence and presence of leaks. Delivered tidal volumes were also assessed. Three patients were simulated: a preterm infant (resistance 100 cm H2O/L/s, compliance 2 mL/cm H2O, inspiratory time of the patient [TI] 400 ms, inspiratory effort 1 and 2 cm H2O), a full-term infant (resistance 50 cm H2O/L/s, compliance 5 mL/cm H2O, TI 500 ms, inspiratory effort 2 and 4 cm H2O), and a child (resistance 30 cm H2O/L/s, compliance 10 mL/cm H2O, TI 600 ms, inspiratory effort 5 and 10 cm H2O). Two PS levels were tested (10 and 15 cm H2O) with and without leaks and with and without the leak compensation algorithm activated. RESULTS: Without leaks, only 2 neonatal ventilators and one adult ventilator had trigger delays under a given predefined acceptable limit (1/8 TI). Pressurization showed high variability between ventilators. Most ventilators showed TI in excess high enough to seriously impair patient-ventilator synchronization (> 50% of the TI of the subject). In some ventilators, leaks led to autotriggering and impairment of ventilation performance, but the influence of leaks was generally lower in neonatal ventilators. When a noninvasive ventilation algorithm was available, this was partially corrected. In general, tidal volume was calculated too low by the ventilators in the presence of leaks; the noninvasive ventilation algorithm was able to correct this difference in only 2 adult ventilators. CONCLUSIONS: No ventilator performed equally well under all tested conditions for all explored parameters. However, neonatal ventilators tended to perform better in the presence of leaks. These findings emphasize the need to improve algorithms for assisted ventilation modes to better deal with situations of high airway resistance, low pulmonary compliance, and the presence of leaks.

An iterative Multiscale Finite-Volume algorithm converging to the exact solution

The multiscale finite volume (MsFV) method has been developed to efficiently solve large heterogeneous problems (elliptic or parabolic); it is usually employed for pressure equations and delivers conservative flux fields to be used in transport problems. The method essentially relies on the hypothesis that the (fine-scale) problem can be reasonably described by a set of local solutions coupled by a conservative global (coarse-scale) problem. In most cases, the boundary conditions assigned for the local problems are satisfactory and the approximate conservative fluxes provided by the method are accurate. In numerically challenging cases, however, a more accurate localization is required to obtain a good approximation of the fine-scale solution. In this paper we develop a procedure to iteratively improve the boundary conditions of the local problems. The algorithm relies on the data structure of the MsFV method and employs a Krylov-subspace projection method to obtain an unconditionally stable scheme and accelerate convergence. Two variants are considered: in the first, only the MsFV operator is used; in the second, the MsFV operator is combined in a two-step method with an operator derived from the problem solved to construct the conservative flux field. The resulting iterative MsFV algorithms allow arbitrary reduction of the solution error without compromising the construction of a conservative flux field, which is guaranteed at any iteration. Since it converges to the exact solution, the method can be regarded as a linear solver. In this context, the schemes proposed here can be viewed as preconditioned versions of the Generalized Minimal Residual method (GMRES), with a very peculiar characteristic that the residual on the coarse grid is zero at any iteration (thus conservative fluxes can be obtained).

Proposal of a parallel architecture for a motion detection algorithm

This paper proposes a parallel architecture for estimation of the motion of an underwater robot. It is well known that image processing requires a huge amount of computation, mainly at low-level processing where the algorithms are dealing with a great number of data. In a motion estimation algorithm, correspondences between two images have to be solved at the low level. In the underwater imaging, normalised correlation can be a solution in the presence of non-uniform illumination. Due to its regular processing scheme, parallel implementation of the correspondence problem can be an adequate approach to reduce the computation time. Taking into consideration the complexity of the normalised correlation criteria, a new approach using parallel organisation of every processor from the architecture is proposed

Speculative parallelization of multipath radiosity algorithm

In computer graphics, global illumination algorithms take into account not only the light that comes directly from the sources, but also the light interreflections. This kind of algorithms produce very realistic images, but at a high computational cost, especially when dealing with complex environments. Parallel computation has been successfully applied to such algorithms in order to make it possible to compute highly-realistic images in a reasonable time. We introduce here a speculation-based parallel solution for a global illumination algorithm in the context of radiosity, in which we have taken advantage of the hierarchical nature of such an algorithm

A Monte Carlo-Based Fiber Tracking Algorithm using Diffusion Tensor MRI

Diffusion tensor magnetic resonance imaging, which measures directional information of water diffusion in the brain, has emerged as a powerful tool for human brain studies. In this paper, we introduce a new Monte Carlo-based fiber tracking approach to estimate brain connectivity. One of the main characteristics of this approach is that all parameters of the algorithm are automatically determined at each point using the entropy of the eigenvalues of the diffusion tensor. Experimental results show the good performance of the proposed approach

Constraint algorithm for k-presymplectic Hamiltonian systems. Application to singular field theories

The k-symplectic formulation of field theories is especially simple, since only tangent and cotangent bundles are needed in its description. Its defining elements show a close relationship with those in the symplectic formulation of mechanics. It will be shown that this relationship also stands in the presymplectic case. In a natural way,one can mimick the presymplectic constraint algorithm to obtain a constraint algorithmthat can be applied to k-presymplectic field theory, and more particularly to the Lagrangian and Hamiltonian formulations offield theories defined by a singular Lagrangian, as well as to the unified Lagrangian-Hamiltonian formalism (Skinner--Rusk formalism) for k-presymplectic field theory. Two examples of application of the algorithm are also analyzed.

Risk assessment of recurrence in sporadic retinoblastoma using a molecular-based algorithm

Le rétinoblastome (Rb) est une tumeur provenant des cellules rétiniennes progénitrices des photorécepteurs. C'est la tumeur pédiatrique maligne la plus fréquente avec une incidence par naissance évaluée entre 1/15Ό00 et 1/20Ό00. Les enfants atteints de Rb sont diagnostiqué dans leur grande majorité avant l'âge de 4 ans, soit le temps nécessaire à la différentiation et à la maturation des photorécepteurs et donc à la disparition de la cellule d'origine du Rb. La survie du patient, la sauvegarde oculaire et le pronostic visuel restent excellents pour autant que le traitement ne soit pas différé. Dans sa variante non héréditaire (60%) le Rb est toujours unilatéral et sporadique. Le Rb héréditaire de transmission dominante autosomique (40%), se décline sous toutes les formes, familiale (10%) ou sporadique (30%), que l'atteinte soit unilatérale ou bilatérale. La majorité des mutations causales sont uniques et distribuées de façon aléatoire sur la totalité du gène RB1 sans région prédisposante. La détection de ces mutations est couteuse et chronophage, tout en présentant un taux de détection relativement bas; surtout dans les cas de Rb sporadiques unilatéraux. Dans le but d'identifier les patients présentant un risque réel de développer un Rb, et de réduire le nombre d'examens sous narcose requis pour le dépistage de la maladie chez les sujets à risque, nous avons développé une stratégie sensible, rapide, efficace et peu couteuse basée sur une analyse de l'haplotype intragénique. Cet algorithme prend en compte a) la perte d'hétérozygotie intratumorale du gène RB1, b) l'origine paternelle préférentielle des nouvelles mutations germinales et c) un risque a priori dérivé des données empiriques de Vogel. Pendant la période allant de janvier 1994 à décembre 2006, nous avons comparé l'apparition de nouveau Rb parmi la fratrie et la descendance de patient atteints au nombre de nouveaux cas attendus calculé par notre algorithme. 134 familles ont été étudiées. L'analyse moléculaire a été effectuée chez 570 personnes dont 99 patients âgés de moins de 4 ans et donc à risque de développer un Rb. Parmi cette cohorte, nous avons observé l'apparition d'un cas de Rb, alors que les risques cumulés a posteriori calculé par notre algorithme prédisait l'apparition de 1.77 nouveau cas. Dans cette étude, nous avons pu valider notre algorithme prédisant la récurrence de Rb chez les parents de 1er degré de patients atteints. Cet outil devrait grandement faciliter le conseil génétique ainsi que le suivi des patients à risque de développer un Rb, surtout dans les cas ou le séquençage direct du gène RB1 n'est pas disponible ou est resté non informatif. - Purpose: Most RBI mutations are unique and distributed throughout the RBI gene. Their detection can be time-consuming and the yield especially low in cases of conservatively-treated sporadic unilateral retinoblas-toma (Rb) patients. In order to identify patients with true risk of developing Rb, and to reduce the number of unnecessary examinations under anesthesia in all other cases, we developed a universal sensitive, efficient and cost-effective strategy based on intragenic haplotype analysis. Methods: This algorithm allows the calculation of the a posteriori risk of developing Rb and takes into account (a) RBI loss of heterozygosity in tumors, (b) preferential paternal origin of new germline mutations, (c) a priori risk derived from empirical data by Vogel, and (d) disease penetrance of 90% in most cases. We report the occurrence of Rb in first degree relatives of patients with sporadic Rb who visited the Jules Gonin Eye Hospital, Lausanne, Switzerland, from January 1994 to December 2006 compared to expected new cases of Rb using our algorithm. Results: A total of 134 families with sporadic Rb were enrolled; testing was performed in 570 individuals and 99 patients younger than 4 years old were identified. We observed one new case of Rb. Using our algorithm, the cumulated total a posteriori risk of recurrence was 1.77. Conclusions: This is the first time that linkage analysis has been validated to monitor the risk of recurrence in sporadic Rb. This should be a useful tool in genetic counseling, especially when direct RBI screening for mutations leaves a negative result or is unavailable.

Glutathione Precursor, N-Acetyl-Cysteine, Modulates EEG Synchronization in Schizophrenia Patients

Background: Glutathione (GSH) dysregulation at the gene, protein and functional levels observed in schizophrenia patients, and schizophrenia-like anomalies in GSH deficit experimental models, suggest that genetic glutathione synthesis impairments represent one major risk factor for the disease (Do et al., 2009). In a randomized, double blind, placebo controlled, add-on clinical trial of 140 patients, the GSH precursor N-Acetyl-Cysteine (NAC, 2g/day, 6 months) significantly improved the negative symptoms and reduced sideeffects due to antipsychotics (Berk et al., 2008). In a subset of patients (n=7), NAC (2g/day, 2 months, cross-over design) also improved auditory evoked potentials, the NMDA-dependent mismatch negativity (Lavoie et al, 2008). Methods: To determine whether increased GSH levels would modulate the topography of functional brain connectivity, we applied a multivariate phase synchronization (MPS) estimator (Knyazeva et al, 2008) to dense-array EEGs recorded during rest with eyes closed at the protocol onset, the point of crossover, and at its end. Results: The whole-head imaging revealed a specific synchronization landscape in NAC compared to placebo condition. In particular, NAC increased MPS over frontal and left temporal regions in a frequency-specific manner. The topography and direction of MPS changes were similar and robust in all 7 patients. Moreover, these changes correlated with the changes in the Liddle's score of disorganization, thus linking EEG synchronization to the improvement of the clinical picture. Conclusions: The data suggest an important pathway towards new therapeutic strategies that target GSH dysregulation in schizophrenia. They also show the utility of MPS mapping as a marker of treatment efficacy.

A Label space reduction algorithm for P2MP LSPs using asymmetric tunnels

The aim of traffic engineering is to optimise network resource utilization. Although several works on minimizing network resource utilization have been published, few works have focused on LSR label space. This paper proposes an algorithm that uses MPLS label stack features in order to reduce the number of labels used in LSPs forwarding. Some tunnelling methods and their MPLS implementation drawbacks are also discussed. The algorithm described sets up the NHLFE tables in each LSR, creating asymmetric tunnels when possible. Experimental results show that the algorithm achieves a large reduction factor in the label space. The work presented here applies for both types of connections: P2MP and P2P