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.

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.

Using Dynamic Programming Decomposition for Revenue Management with Opaque Products

Opaque products enable service providers to hide specific characteristics of their service fulfillment from the customer until after purchase. Prominent examples include internet-based service providers selling airline tickets without defining details, such as departure time or operating airline, until the booking has been made. Owing to the resulting flexibility in resource utilization, the traditional revenue management process needs to be modified. In this paper, we extend dynamic programming decomposition techniques widely used for traditional revenue management to develop an intuitive capacity control approach that allows for the incorporation of opaque products. In a simulation study, we show that the developed approach significantly outperforms other well-known capacity control approaches adapted to the opaque product setting. Based on the approach, we also provide computational examples of how the share of opaque products as well as the degree of opacity can influence the results.

The role of peat decomposition in patterned mires: a case study from the central Swiss Alps

A number of hydrological, botanical, macro- and micro-climatological processes are involved in the formation of patterned peatlands. La Grande Tsa at 2336 m a.s.l. is probably the highest bog in the central Swiss Alps and is unique in its pattern. In two of five pools there is in the contact zone between the basal peat and the overlying gyttja an unconformity in the depth-age models based on radiocarbon dates. Palynostratigraphies of cores from a ridge and a pool confirm the occurrence of an unconformity in the contact zone. We conclude that deepening of the pools results from decomposition of peat. The fact that the dated unconformities in the two pools and the unconformity in the ridge-core all fall within the Bronze Age suggest they were caused by events external to the bog. We hypothesize that early transhumance resulted in anthropogenic lowering of the timberline, which resulted in a reduction in the leaf-area index and evapotranspiration, and in higher water levels and thus pool formation.