From EEG dependency multichannel matching pursuit to sparse topographic EEG decomposition

In this work, we present a multichannel EEG decomposition model based on an adaptive topographic time-frequency approximation technique. It is an extension of the Matching Pursuit algorithm and called dependency multichannel matching pursuit (DMMP). It takes the physiologically explainable and statistically observable topographic dependencies between the channels into account, namely the spatial smoothness of neighboring electrodes that is implied by the electric leadfield. DMMP decomposes a multichannel signal as a weighted sum of atoms from a given dictionary where the single channels are represented from exactly the same subset of a complete dictionary. The decomposition is illustrated on topographical EEG data during different physiological conditions using a complete Gabor dictionary. Further the extension of the single-channel time-frequency distribution to a multichannel time-frequency distribution is given. This can be used for the visualization of the decomposition structure of multichannel EEG. A clustering procedure applied to the topographies, the vectors of the corresponding contribution of an atom to the signal in each channel produced by DMMP, leads to an extremely sparse topographic decomposition of the EEG.

How does spatial extent of fMRI datasets affect independent component analysis decomposition?

Spatial independent component analysis (sICA) of functional magnetic resonance imaging (fMRI) time series can generate meaningful activation maps and associated descriptive signals, which are useful to evaluate datasets of the entire brain or selected portions of it. Besides computational implications, variations in the input dataset combined with the multivariate nature of ICA may lead to different spatial or temporal readouts of brain activation phenomena. By reducing and increasing a volume of interest (VOI), we applied sICA to different datasets from real activation experiments with multislice acquisition and single or multiple sensory-motor task-induced blood oxygenation level-dependent (BOLD) signal sources with different spatial and temporal structure. Using receiver operating characteristics (ROC) methodology for accuracy evaluation and multiple regression analysis as benchmark, we compared sICA decompositions of reduced and increased VOI fMRI time-series containing auditory, motor and hemifield visual activation occurring separately or simultaneously in time. Both approaches yielded valid results; however, the results of the increased VOI approach were spatially more accurate compared to the results of the decreased VOI approach. This is consistent with the capability of sICA to take advantage of extended samples of statistical observations and suggests that sICA is more powerful with extended rather than reduced VOI datasets to delineate brain activity.

DECOMPOSITION OF REGRESSION ESTIMATORS TO EXPLORE THE INFLUENCE OF "UNMEASURED" TIME-VARYING CONFOUNDERS

In environmental epidemiology, exposure X and health outcome Y vary in space and time. We present a method to diagnose the possible influence of unmeasured confounders U on the estimated effect of X on Y and to propose several approaches to robust estimation. The idea is to use space and time as proxy measures for the unmeasured factors U. We start with the time series case where X and Y are continuous variables at equally-spaced times and assume a linear model. We define matching estimator b(u)s that correspond to pairs of observations with specific lag u. Controlling for a smooth function of time, St, using a kernel estimator is roughly equivalent to estimating the association with a linear combination of the b(u)s with weights that involve two components: the assumptions about the smoothness of St and the normalized variogram of the X process. When an unmeasured confounder U exists, but the model otherwise correctly controls for measured confounders, the excess variation in b(u)s is evidence of confounding by U. We use the plot of b(u)s versus lag u, lagged-estimator-plot (LEP), to diagnose the influence of U on the effect of X on Y. We use appropriate linear combination of b(u)s or extrapolate to b(0) to obtain novel estimators that are more robust to the influence of smooth U. The methods are extended to time series log-linear models and to spatial analyses. The LEP plot gives us a direct view of the magnitude of the estimators for each lag u and provides evidence when models did not adequately describe the data.

A guide to sparkology: the taxonomy of elementary cellular Ca2+ signaling events

Since the discovery of the Ca(2+) spark as an elementary event of cellular Ca(2+) signaling almost 15 years ago, the family of newly described Ca(2+) signal entities has been ever growing. While scientists working in Ca(2+) signaling may have maintained an overview over the specifics of this nomenclature, those outside the field often make the complaint that they feel hopelessly lost. With the present review we collect and summarize systematic information on the many Ca(2+) signaling events described in a variety of tissues and cells, and we emphasize why and how each of them has its own importance. Most of these signals are taking place in the cytosol of the respective cells, but several events have been recorded from intracellular organelles as well, where they may serve their own specific functions. Finally, we also try to convey an integrated view as to why cellular microdomain signaling is of fundamental biological importance.

Enumerating constrained elementary flux vectors of metabolic networks

The concept of elementary vector is generalised to the case where the steady-state space of the metabolic network is not a flux cone but is a general polyhedron due to further inhomogeneous constraints on the flows through some of the reactions. On one hand, this allows to selectively enumerate elementary modes which satisfy certain optimality criteria and this can yield a large computational gain compared with full enumeration. On the other hand, in contrast to the single optimum found by executing a linear program, this enables a comprehensive description of the set of alternate optima often encountered in flux balance analysis. The concepts are illustrated on a metabolic network model of human cardiac mitochondria.

Observation and identification of metabolites emerging during postmortem decomposition of brain tissue by means of in situ 1H-magnetic resonance spectroscopy

Postmortem decomposition of brain tissue was investigated by (1)H-magnetic resonance spectroscopy (MRS) in a sheep head model and selected human cases. Aiming at the eventual estimation of postmortem intervals in forensic medicine, this study focuses on the characterization and identification of newly observed metabolites. In situ single-voxel (1)H-MRS at 1.5 T was complemented by multidimensional homo- and heteronuclear high-resolution NMR spectroscopy of an extract of sheep brain tissue. The inclusion of spectra of model solutions in the program LC Model confirmed the assignments in situ. The first postmortem phase was characterized mainly by changes in the concentrations of metabolites usually observed in vivo and by the appearance of previously reported decay products. About 3 days postmortem, new metabolites, including free trimethylammonium, propionate, butyrate, and iso-butyrate, started to appear in situ. Since the observed metabolites and the time course is comparable in sheep and human brain tissue, the model system seems to be appropriate.

Precise pose recovery of distal locking holes from single calibrated fluoroscopic image via a novel variable decomposition approach

This paper presents a novel variable decomposition approach for pose recovery of the distal locking holes using single calibrated fluoroscopic image. The problem is formulated as a model-based optimal fitting process, where the control variables are decomposed into two sets: (a) the angle between the nail axis and its projection on the imaging plane, and (b) the translation and rotation of the geometrical model of the distal locking hole around the nail axis. By using an iterative algorithm to find the optimal values of the latter set of variables for any given value of the former variable, we reduce the multiple-dimensional model-based optimal fitting problem to a one-dimensional search along a finite interval. We report the results of our in vitro experiments, which demonstrate that the accuracy of our approach is adequate for successful distal locking of intramedullary nails.

Effects of elevated atmospheric CO2 and O3 on wood density, antomical proerties and decomposition of Northern Hardwoods

Anthropogenic activities continue to drive atmospheric CO2 and O3 concentrations to levels higher than during the pre-industrial era. Accumulating evidence indicates that both elevated CO2 and elevated O3 could modify the quantity and biochemistry of woody plant biomass. Anatomical properties of woody plants are largely influenced by the activity of the cambium and the growth characteristics of wood cells, which are in turn influenced by a range of environmental factors. Hence, alterations in the concentrations of atmospheric CO2 and / or O3 could also impact wood anatomical properties. Many fungi derive their metabolic resources for growth from plant litter, including woody tissue, and therefore modifications in the quantity, biochemistry and anatomical properties of woody plants in response to elevated CO2 and / or O3 could impact the community of wood-decaying fungi and rates of wood decomposition. Consequently carbon and nutrient cycling and productivity of terrestrial ecosystem could also be impacted. Alterations in wood structure and biochemistry of woody plants could also impact wood density and subsequently impact wood quality. This dissertation examined the long term effects of elevated CO2 and / or O3 on wood anatomical properties, wood density, wood-decaying fungi and wood decomposition of northern hardwood tree species at the Aspen Free-Air CO2 and O3 Enrichment (Aspen FACE) project, near Rhinelander, WI, USA. Anatomical properties of wood varied significantly with species and aspen genotypes and radial position within the stem. Elevated CO2 did not have significant effects on wood anatomical properties in trembling aspen, paper birch or sugar maple, except for marginally increasing (P < 0.1) the number of vessels per square millimeter. Elevated O3 marginally or significantly altered vessel lumen diameter, cell wall area and vessel lumen area proportions depending on species and radial position. In line with the modifications in the anatomical properties, elevated CO2 and O3, alone, significantly modified wood density but effects were species and / or genotype specific. However, the effects of elevated CO2 and O3, alone, on wood anatomical properties and density were ameliorated when in combination. Wood species had a much greater impact on the wood-decaying fungal community and initial wood decomposition rate than did growth or decomposition of wood in elevated CO2 and / or O3. Polyporales, Agaricales, and Russulales were the dominant orders of fungi isolated. Based on the current results, future higher levels of CO2 and O3 may have moderate effects on wood quality of northern hardwoods, but for utilization purposes these may not be considered significant. However, wood-decaying fungal community composition and decomposition of northern hardwoods may be altered via shifts in species and / or genotype composition under future higher levels of CO2 and O3.

HIGH PERFORMANCE, LOW COST SUBSPACE DECOMPOSITION AND POLYNOMIAL ROOTING FOR REAL TIME DIRECTION OF ARRIVAL ESTIMATION: ANALYSIS AND IMPLEMENTATION

This thesis develops high performance real-time signal processing modules for direction of arrival (DOA) estimation for localization systems. It proposes highly parallel algorithms for performing subspace decomposition and polynomial rooting, which are otherwise traditionally implemented using sequential algorithms. The proposed algorithms address the emerging need for real-time localization for a wide range of applications. As the antenna array size increases, the complexity of signal processing algorithms increases, making it increasingly difficult to satisfy the real-time constraints. This thesis addresses real-time implementation by proposing parallel algorithms, that maintain considerable improvement over traditional algorithms, especially for systems with larger number of antenna array elements. Singular value decomposition (SVD) and polynomial rooting are two computationally complex steps and act as the bottleneck to achieving real-time performance. The proposed algorithms are suitable for implementation on field programmable gated arrays (FPGAs), single instruction multiple data (SIMD) hardware or application specific integrated chips (ASICs), which offer large number of processing elements that can be exploited for parallel processing. The designs proposed in this thesis are modular, easily expandable and easy to implement. Firstly, this thesis proposes a fast converging SVD algorithm. The proposed method reduces the number of iterations it takes to converge to correct singular values, thus achieving closer to real-time performance. A general algorithm and a modular system design are provided making it easy for designers to replicate and extend the design to larger matrix sizes. Moreover, the method is highly parallel, which can be exploited in various hardware platforms mentioned earlier. A fixed point implementation of proposed SVD algorithm is presented. The FPGA design is pipelined to the maximum extent to increase the maximum achievable frequency of operation. The system was developed with the objective of achieving high throughput. Various modern cores available in FPGAs were used to maximize the performance and details of these modules are presented in detail. Finally, a parallel polynomial rooting technique based on Newton’s method applicable exclusively to root-MUSIC polynomials is proposed. Unique characteristics of root-MUSIC polynomial’s complex dynamics were exploited to derive this polynomial rooting method. The technique exhibits parallelism and converges to the desired root within fixed number of iterations, making this suitable for polynomial rooting of large degree polynomials. We believe this is the first time that complex dynamics of root-MUSIC polynomial were analyzed to propose an algorithm. In all, the thesis addresses two major bottlenecks in a direction of arrival estimation system, by providing simple, high throughput, parallel algorithms.