Spatial-temporal analysis of mortality using splines

A method for constructing contour maps in disease mapping using smoothing splines is presented. Smoothing errors are discussed and error maps suggested. The approach is illustrated for mortality data on cardiovascular diseases in northwestern Germany (Lower Saxony and Bremen) 1970-1979.     

Spatio-temporal models with errors in covariates: mapping Ohio lung cancer mortality

In estimating spatial disease patterns, as well as in related assessments of environmental equity, regional morbidity and mortality rate maps are widely used. Hierarchical Bayes methods are increasingly popular tools for creating such maps, since they permit smoothing of the fitted rates toward spatially local mean values, with more unreliable estimates (those arising in low-population regions) receiving more smoothing. In this paper we blend methods for spatial-temporal mapping with those for handling errors in covariates in a single hierarchical model framework. Estimated posterior distributions for the resulting highly-parameterized models are obtained via Markov chain Monte Carlo (MCMC) methods, which also play a key role in our approach to model evaluation and selection. We apply our approach to a data set of county-specific lung cancer rates in the state of Ohio during the period 1968-1988. Our model uses age-adjusted death rates, and incorporates recent information regarding smoking prevalence, population density, and the socio-economic status of the counties. This information is critical to understanding the role played by a certain depleted uranium fuel processing facility on the elevated lung cancer rates in the counties that neighbour it.     

The role of visual experience in the development of columns in cat visual cortex

The role of experience in the development of the cerebral cortex has long been controversial. Patterned visual experience in the cat begins when the eyes open about a week after birth. Cortical maps for orientation and ocular dominance in the primary visual cortex of cats were found to be present by 2 weeks. Early pattern vision appeared unimportant because these cortical maps developed identically until nearly 3 weeks of age, whether or not the eyes were open. The na?ve maps were powerfully dominated by the contralateral eye, and experience was needed for responses to the other eye to become strong, a process unlikely to be strictly Hebbian. With continued visual deprivation, responses to both eyes deteriorated, with a time course parallel to the well-known critical period of cortical plasticity. The basic structure of cortical maps is therefore innate, but experience is essential for specific features of these maps, as well as for maintaining the responsiveness and selectivity of cortical neurons.     

Development of spatiotemporal receptive fields of simple cells: II. Simulation and analysis

In part I of this article a correlation based model for the developmental process of spatiotemporal receptive fields has been introduced. In this model the development is described as an activity-dependent competition between four types of input from the lateral geniculate nucleus onto a cortical cell, viz. non-lagged ON and OFF and lagged ON and OFF inputs. In the present paper simulation results and a first analysis are presented for this model. We study the developmental process both before and after eye-opening and compare the results with experimental data from reverse correlation measurements. The outcome of the developmental process is determined mainly by the spatial and the temporal correlations between the different inputs. In particular, if the mean correlation between non-lagged and lagged inputs is weak, receptive fields with a widely varying degree of direction selectivity emerge. However, spatiotemporal receptive fields may show rotation of their preferred orientation as a function of response delay. Even if the mean correlation between two types of temporal input is not weak, direction-selective receptive fields may emerge because of an intracortical interaction between different cortical maps. In an environment of moving lines or gratings, direction-selective receptive fields develop only if the distribution of the directions of motion presented during development shows some anisotropy. In this case, a continuous map of preferred direction is also shown to develop.     

A comprehensive appraisal of 241Am in soils around Rocky Flats, Colorado

Soils east of Rocky Flats Plant (RFETS) near Golden, Colorado, were contaminated with actinides because of accidental release of oils laden with plutonium isotopes. Consequently, these soils were contaminated by 241Am due to radioactive decay of 241Pu (t1/2 = 14.4 y). A spatial analysis of 241Am activity in soils east of RFETS was conducted to elucidate the magnitude and the mode of 241Am dispersion in the soil environment. 241Am activity of 178 soil samples ranged from 0.037 Bq kg-1 to 10,004 Bq kg-1 with a mean of 214 Bq kg-1, median of 7.28 Bq kg-1, standard deviation of 942 Bq kg-1, and a coefficient of variation of 4.3. Spatial analysis of 241Am in soils around RFETS was conducted using indicator kriging, which is a nonparametric technique especially suitable to model a conditional cumulative distribution function (ccdf) of highly skewed environmental data such as 241Am. The ccdf was used to generate an E-type (mean of the conditional cdf) surface. The resulted surfaces were consistent with the hypothesis that the westerly winds were the dominant mechanism of americium dispersal. The spatial distribution and dispersal mechanisms of 241Am were similar to those of 239+240Pu. The ccdf was also used to construct probability of exceedence maps of 241Am in soils. For the purpose of this report two threshold values for the probability maps were selected: (1) the mean measured background activity of 241Am (0.4 Bq kg-1), and (2) the programmatic preliminary remediation goal for residential occupancy scenario (87.7 Bq kg-1). The probability-of-exceedance maps provide estimates of spatial uncertainty associated with each threshold. The E-type maps in conjunction with the probability-of-exceedance maps provide a robust framework for future cleanup options and land use decisions.     

The retinotopy of visual spatial attention

We used high-field (3T) functional magnetic resonance imaging (fMRI) to label cortical activity due to visual spatial attention, relative to flattened cortical maps of the retinotopy and visual areas from the same human subjects. In the main task, the visual stimulus remained constant, but covert visual spatial attention was varied in both location and load. In each of the extrastriate retinotopic areas, we found MR increases at the representations of the attended target. Similar but smaller increases were found in V1. Decreased MR levels were found in the same cortical locations when attention was directed at retinotopically different locations. In and surrounding area MT+, MR increases were lateralized but not otherwise retinotopic. At the representation of eccentricities central to that of the attended targets, prominent MR decreases occurred during spatial attention.     

Pre-processing variance reducing techniques in multispectral positron emission tomography

Stochastic fluctuations and systematic errors severely restrict the potential of multispectral acquisition to improve scatter correction by energy-dependent processing in high-resolution positron emission tomography (PET). To overcome this limitation, three pre-processing approaches which reduce stochastic fluctuations and systematic errors without degrading spatial resolution were investigated: statistical variance was reduced by smoothing acquired data in energy space, systematic errors due to nonuniform detector efficiency were minimized by normalizing the data in the spatial domain and the overall variance was further reduced by selecting an optimal pre-processing sequence. Selection of the best protocol to reduce stochastic fluctuations entailed comparisons between four smoothing algorithms (prior constrained (PC) smoothing, weighted smoothing (WS), ideal low-pass filtering (ILF) and mean median (MM) smoothing) and permutations of three pre-processing procedures (smoothing, normalization and subtraction of random events). Results demonstrated that spectral smoothing by WS, ILF and MM efficiently reduces the statistical variance in both the energy and spatial domains without observable spatial resolution loss. The ILF algorithm was found to be the most convenient in terms of simplicity and efficiency. Regardless of the position of subtraction of randoms in the sequence, reduction of the systematic errors by normalization followed by spectral smoothing to suppress statistical noise produced the best results. However, subtraction of random events first in the sequence reduces computation load by half since the need to pre-process this distribution before subtraction is removed. In summary, normalizing data in the spatial domain and smoothing data in energy space are essential steps required to reduce systematic errors and statistical variance independently without degrading spatial resolution of multispectral PET data.     

Comparison of Ordinary Kriging, Regression Kriging, and Cokriging Techniques to Estimate Soil Salinity Using LANDSAT Images

The objectives of this study are: (1) to evaluate the LANDSAT best band combinations to estimate soil salinity with different crop types; (2) to compare ordinary kriging, regression kriging, and cokriging techniques to generate accurate soil salinity maps when applied to LANDSAT images; and (3) to compare the performance of different crop types: alfalfa, cantaloupe, corn, and wheat as indicators of soil salinity. This study was conducted in an area in the southern part of the Arkansas River Basin in Colorado. Six LANDSAT images acquired during the years: 2000, 2001, 2003, 2004, 2005, and 2006 in conjunction with field data were used to estimate soil salinity in the study area. The optimal subsets of band combinations from the LANDSAT images that correlate best with the soil salinity data sets were selected. Ordinary kriging, regression kriging, and cokriging were applied to 2,914 soil salinity data points collected in alfalfa, cantaloupe, corn, and wheat fields in conjunction with the selected LANDSAT image band combination subsets. Ordinary least-squares (OLSs) were used to regress the correlated band combinations to generate a soil salinity surface. The residuals of the OLS multiple regression model were kriged and combined with the soil salinity surface generated using the OLS multiple regression model to produce the final soil salinity surface of the regression kriging model. The same LANDSAT band combinations used with the regression kriging technique were used as secondary data variables with the cokriging technique, while the soil salinity data was used as a primary variable. The results show that the best band combinations for estimating soil salinity with different crops are as follows: alfalfa [red, near infrared, and normalized difference vegetation index (NDVI)]; cantaloupe (blue and green); corn (near, thermal, and NDVI); and wheat (blue and thermal). The performance of the different geostatistical models used in this study is: (1) ordinary kriging; (2) regression kriging; and (3) cokriging. Estimation of soil salinity works best for corn, then wheat, cantaloupe, and alfalfa.     

Metrics for cortical map organization and lateralization

Cerebral lateralization refers to the poorly understood fact that some functions are better controlled by one side of the brain than the other (e.g. handedness, language). Of particular concern here are the asymmetries apparent in cortical topographic maps that can be demonstrated electrophysiologically in mirror-image locations of the cerebral cortex. In spite of great interest in issues surrounding cerebral lateralization, methods for measuring the degree of organization and asymmetry in cortical maps are currently quite limited. In this paper, several measures are developed and used to assess the degree of organization, lateralization, and mirror symmetry in topographic map formation. These measures correct for large constant displacements as well as curving of maps. The behavior of the measures is tested on several topographic maps obtained by self-organization of an initially random artificial neural network model of a bihemispheric brain, and the results are compared with subjective assessments made by humans.     

Genetic features of the spatial organization of the human cerebral cortex

The within-pair similarity of the topographical maps of the spatial synchronization of the cerebral cortical potentials was studied in 11 pairs of monozygotic twins and 20 genetically unrelated subjects in the resting state and during four types of intellectual activity. It was demonstrated that the degree of similarity of the topographical maps is higher in the resting state in the monozygotic twins than in the genetically unrelated subjects. No frequency band of the EEG which exerts a special influence on the formation of the high degree of similarity of the spatial synchronization of the potentials was found. The maximal similarity in the topographical maps of the monozygotic twins is observed in the anterior regions of the left hemisphere and in the temporoparietal regions of the right. Differences between groups are recorded in both regions; however they are especially distinct in the right hemisphere. The data obtained make it possible to draw an inference regarding the influence of a genetic factor not only on individual components of the EEG (which had been previously demonstrated by a number of authors), but also on the systemic organization of the cortical processes. The right-hemispheric activity apparently relates to the most controllable genetic processes. The results of the analysis of the similarity of topographical maps during various types of activity confirm this: the coefficient of similarity reaches the level of significance in the majority of pairs of monozygotic twins only in a "right-hemispheric" test (simultaneous analysis of nonverbal material).     

A computer simulation study of cortical imaging from scalp potentials

In this paper, computer simulation studies were conducted to test the feasibility of imaging brain electrical activity from the scalp electroencephalograms. The inhomogeneous three-concentric-sphere head model was used to represent the head volume conductor. Closed spherical dipole layers, consisting of several thousand uniformly distributed dipoles, were used to reconstruct the cortical potential maps corresponding to neuronal sources located inside the brain. Simulation results indicate that the present procedure can image both cortical and deep sources, and for the cortical sources, a spatial resolution as high as 1.2 cm can be achieved.     

Effects of diffusion anisotropy on lesion delineation in a rat model of cerebral ischemia

The effects of white and gray matter diffusion anisotropy on ischemic lesion delineation have been studied in the rat model of middle cerebral artery occlusion. Apparent diffusion coefficient (ADC) maps obtained by conventional pulsed gradient spin echo diffusion-weighted imaging (PGSE-DWI) were compared with maps of the trace of the diffusion tensor in both normal and occluded animals. Diffusion tensor trace maps were derived from the average of the ADC maps from three separate experiments with diffusion weighting along three orthogonal axes, and also from a single-scan method. A marked degree of diffusion anisotropy was observed in both cortical gray matter and white matter from ADC maps of the control animals. In the occluded animals, the systematic effects of anisotropy on ADC and lesion area influenced the delineation of the ischemic territory in the PGSE-DWI ADC maps. However, the two trace methods eliminated these effects and gave consistent ischemic lesion depiction, despite the use of differing diffusion times in the two measurements.     

Geographic Variability of Rainfall Erosivity Estimation and Impact on Construction Site Erosion Control Design

The Revised Universal Soil Loss Equation (RUSLE) is used often by erosion control planners to estimate the soil loss from urban construction sites when sizing sediment ponds and determining the soil loss under vegetative mats. This project used the existing, complete-year rainfall record for 27 sites in the state of Pennsylvania to compare the USDA isoerodent maps to the annual rainfall erosivity, R, values calculated using the USEPA equations for the National Resource Conservation Service Type II rainfall. The USDA and USEPA maps showed a general trend of increasing median annual R from west-to-east and north-to-south. A trend analysis relating the median R values calculated during this project to geographic location had similar, large-scale geographic trends as the USDA and USEPA maps. However, the R values more closely followed a combination of the annual rainfall pattern and topography (the Appalachian mountains bisect the state). Two case studies of the impacts of these calculations were developed to show the impact of using different values of R on the design of sediment ponds and predicting vegetation establishment. The results of these scenarios indicate that the source of data to predict R can affect the frequency and cost of sediment pond maintenance and may under-predict the protection level required of a vegetated erosion control mat.     

Generalization of prism adaptation.

Prism exposure produces 2 kinds of adaptive response. Recalibration is ordinary strategic remapping of spatially coded movement commands to rapidly reduce performance error. Realignment is the extraordinary process of transforming spatial maps to bring the origins of coordinate systems into correspondence. Realignment occurs when spatial discordance signals noncorrespondence between spatial maps. In Experiment 1, generalization of recalibration aftereffects from prism exposure to postexposure depended upon the similarity of target pointing limb postures. Realignment aftereffects generalized to the spatial maps involved in exposure. In Experiment 2, the 2 kinds of aftereffects were measured for 3 test positions, one of which was the exposure training position. Recalibration aftereffects generalized nonlinearly, while realignment aftereffects generalized linearly, replicating Bedford (1989, 1993a) using a more familiar prism adaptation paradigm. Recalibration and realignment require methods for distinguishing their relative contribution to prism adaptation. (PsycINFO Database Record (c) 2011 APA, all rights reserved)     

Anisotropy in the Spatial Distribution of Roller-Measured Soil Stiffness

The geostatistical analysis of roller-measured soil properties (from continuous compaction control and intelligent compaction) is required for advanced quality control/quality assurance of earthwork and asphalt compaction. This paper explores the existence of anisotropy in the spatial distribution of roller-measured soil stiffness and the effect of anisotropy on kriging. Field testing was conducted to collect roller measurement value (MV) data over typical roadway embankment evaluation areas and on a large square area to enable a robust investigation of anisotropy. The semivariogram analysis of the field data clearly indicates that range anisotropy exists. The spatial distributions of roller MV data are different in the longitudinal x-direction than in the transverse y-direction. Magnitudes of range anisotropy (x-range/y-range) varied from 2.4 to over 5. The observed range anisotropy is not due to the roller measurement system; rather, it is likely due to the directional nature of earthwork construction activities and to alignment geometry. The influence of anisotropy on kriging was found to be significant when considering the use of kriged data in earthwork specifications. The error introduced by not accounting for anisotropy in kriging varied from 5 to 17% when considering pass to pass or layer to layer map analysis. The anisotropy in the spatial distribution of roller MV data should be factored into kriging and other geostatistical analyses. For typical earthwork area geometries, the roller mapping procedure requires a slight modification to determine the y-range and anisotropy ratio.     

Sperm subpopulations in boar (Sus scrofa) and gazelle (Gazella dama mhorr) semen as revealed by pattern analysis of computer-assisted motility assessments

PURPOSE: Combining computational blood flow modeling with three-dimensional medical imaging provides a new approach for studying links between hemodynamic factors and arterial disease. Although this provides patient-specific hemodynamic information, it is subject to several potential errors. This study quantifies some of these errors and identifies optimal reconstruction methodologies. METHODS: A carotid artery bifurcation phantom of known geometry was imaged using a commercial magnetic resonance (MR) imager. Three-dimensional models were reconstructed from the images using several reconstruction techniques, and steady and unsteady blood flow simulations were performed. The carotid bifurcation from a healthy, human volunteer was then imaged in vivo, and geometric models were reconstructed. RESULTS: Reconstructed models of the phantom showed good agreement with the gold standard geometry, with a mean error of approximately 15% between the computed wall shear stress fields. Reconstructed models of the in vivo carotid bifurcation were unacceptably noisy, unless lumenal profile smoothing and approximating surface splines were used. CONCLUSIONS: All reconstruction methods gave acceptable results for the phantom model, but in vivo models appear to require smoothing. If proper attention is paid to smoothing and geometric fidelity issues, models reconstructed from MR images appear to be suitable for use in computational studies of in vivo hemodynamics.     

A voxel-based method for the statistical analysis of gray and white matter density applied to schizophrenia

We describe a novel technique for characterizing regional cerebral gray and white matter differences in structural magnetic resonance images by the application of methods derived from functional imaging. The technique involves automatic scalp-editing of images followed by segmentation, smoothing, and spatial normalization to a symmetrical template brain in stereotactic Talairach space. The basic idea is (i) to convert structural magnetic resonance image data into spatially normalized images of gray (or white) matter density, effected by segmenting the images and smoothing, and then (ii) to use Statistical Parametric Mapping to make inferences about the relationship between gray (or white) matter density and symptoms (or other pathophysiological measures) in a regionally specific fashion. Because the whole brain sum of gray (or white) matter indices is treated as a confound, the analysis reduces to a characterization of relative gray (or white) matter density on a voxel by voxel basis. We suggest that this is a powerful approach to voxel-based statistical anatomy. Using the technique, we constructed maps of the regional cerebral gray and white matter density correlates of syndrome scores (distinct psychotic symptoms) in a group of 15 schizophrenic patients. There was a negative correlation between the score for the reality distortion syndrome and regional gray matter density in the left superior temporal lobe (P = 0.01) and regional white matter density in the corpus callosum (P < 0.001). These abnormalities may be associated with functional changes predisposing to auditory hallucinations and delusions. This method permits the detection of structural differences within the entire brain (as opposed to selected regions of interest) and may be of value in the investigation of structural gray and white matter abnormalities in a variety of brain diseases.     

A method for assessing the accuracy of intersubject registration of the human brain using anatomic landmarks

Several groups have developed methods for registering an individual's 3D MRI by deforming a standard template. This achievement leads to many possibilities for segmentation and morphology that will impact nuclear medical research in areas such as activation and receptor studies. Accordingly, there is a need for methods that can assess the accuracy of intersubject registration. We have developed a method based on a set of 128 anatomic landmarks per hemisphere, both cortical and subcortical, that allows assessment of both global and local transformation accuracy. We applied our method to compare the accuracy of two standard methods of intersubject registration, AIR 3.0 with fifth-order polynomial warping and the Talairach stereotaxic transformation (Talairach and Tournoux, 1988). SPGR MRI's (256 x 256 x 160) of six normal subjects (age 18-24 years) were derformed to match a standard template volume. To assess registration accuracy the landmarks were located on both the template volume and the transformed volumes by an experienced neuroanatomist. The resulting list of coordinates was analyzed graphically and by ANOVA to compare the accuracy of the two methods and the results of the manual analysis. ANOVA performed over all 128 landmarks showed that the Woods method was more accurate than Talairach (left hemisphere F = 2.8, P < 0.001 and right hemisphere F =2.4, P < 0.006). The Woods method provided a better brain surface transformation than did Talairach (F = 18.0, P < 0.0001), but as expected there was a smaller difference for subcortical structures and both had an accuracy <1 mm for the majority of subcortical landmarks. Overall, both the Woods and Talairach method located about 70% of landmarks with an error of 3 mm or less. More striking differences were noted for landmark accuracy 相似文献   

Topographical and temporal specificity of human intraoperative optical intrinsic signals

The goal of this study was to determine the topographical and temporal specificity of neuronal and vascular responses using an intraoperative optical technique (iOIS). The face, thumb, index, and middle fingers were stimulated individually to obtain separate maps of cortical activation. Peak optical responses provided unique, non-overlapping cortical brain maps. Non-peak signals were more dispersed and produced overlapping responses from different digits. Peak iOIS responses colocalized with electrocortical stimulation mapping and evoked potentials. Temporally, we observed statistically significant specificity corresponding to sequential cortical activation during early optical signals (500-1750 ms), but later perfusion responses were non-specific. To our knowledge, this is the first report of either topographical specificity in overlapping spatial patterns, and/or temporal specificity in early perfusion profiles. These results therefore may have significant implications for other perfusion dependent functional imaging techniques.