Spatial weights : constructing weight-compatible exchange matrices from proximity matrices

Exchange matrices represent spatial weights as symmetric probability distributions on pairs of regions, whose margins yield regional weights, generally well-specified and known in most contexts. This contribution proposes a mechanism for constructing exchange matrices, derived from quite general symmetric proximity matrices, in such a way that the margin of the exchange matrix coincides with the regional weights. Exchange matrices generate in turn diffusive squared Euclidean dissimilarities, measuring spatial remoteness between pairs of regions. Unweighted and weighted spatial frameworks are reviewed and compared, regarding in particular their impact on permutation and normal tests of spatial autocorrelation. Applications include tests of spatial autocorrelation with diagonal weights, factorial visualization of the network of regions, multivariate generalizations of Moran's I, as well as "landscape clustering", aimed at creating regional aggregates both spatially contiguous and endowed with similar features.

Adaptive statistical iterative reconstruction reduces patient radiation dose in neuroradiology CT studies.

INTRODUCTION: Adaptive statistical iterative reconstruction (ASIR) can decrease image noise, thereby generating CT images of comparable diagnostic quality with less radiation. The purpose of this study is to quantify the effect of systematic use of ASIR versus filtered back projection (FBP) for neuroradiology CT protocols on patients' radiation dose and image quality. METHODS: We evaluated the effect of ASIR on six types of neuroradiologic CT studies: adult and pediatric unenhanced head CT, adult cervical spine CT, adult cervical and intracranial CT angiography, adult soft tissue neck CT with contrast, and adult lumbar spine CT. For each type of CT study, two groups of 100 consecutive studies were retrospectively reviewed: 100 studies performed with FBP and 100 studies performed with ASIR/FBP blending factor of 40 %/60 % with appropriate noise indices. The weighted volume CT dose index (CTDIvol), dose-length product (DLP) and noise were recorded. Each study was also reviewed for image quality by two reviewers. Continuous and categorical variables were compared by t test and free permutation test, respectively. RESULTS: For adult unenhanced brain CT, CT cervical myelography, cervical and intracranial CT angiography and lumbar spine CT both CTDIvol and DLP were lowered by up to 10.9 % (p < 0.001), 17.9 % (p = 0.005), 20.9 % (p < 0.001), and 21.7 % (p = 0.001), respectively, by using ASIR compared with FBP alone. Image quality and noise were similar for both FBP and ASIR. CONCLUSION: We recommend routine use of iterative reconstruction for neuroradiology CT examinations because this approach affords a significant dose reduction while preserving image quality.

Local optima networks of hard combinatorial landscapes

Combinatorial optimization involves finding an optimal solution in a finite set of options; many everyday life problems are of this kind. However, the number of options grows exponentially with the size of the problem, such that an exhaustive search for the best solution is practically infeasible beyond a certain problem size. When efficient algorithms are not available, a practical approach to obtain an approximate solution to the problem at hand, is to start with an educated guess and gradually refine it until we have a good-enough solution. Roughly speaking, this is how local search heuristics work. These stochastic algorithms navigate the problem search space by iteratively turning the current solution into new candidate solutions, guiding the search towards better solutions. The search performance, therefore, depends on structural aspects of the search space, which in turn depend on the move operator being used to modify solutions. A common way to characterize the search space of a problem is through the study of its fitness landscape, a mathematical object comprising the space of all possible solutions, their value with respect to the optimization objective, and a relationship of neighborhood defined by the move operator. The landscape metaphor is used to explain the search dynamics as a sort of potential function. The concept is indeed similar to that of potential energy surfaces in physical chemistry. Borrowing ideas from that field, we propose to extend to combinatorial landscapes the notion of the inherent network formed by energy minima in energy landscapes. In our case, energy minima are the local optima of the combinatorial problem, and we explore several definitions for the network edges. At first, we perform an exhaustive sampling of local optima basins of attraction, and define weighted transitions between basins by accounting for all the possible ways of crossing the basins frontier via one random move. Then, we reduce the computational burden by only counting the chances of escaping a given basin via random kick moves that start at the local optimum. Finally, we approximate network edges from the search trajectory of simple search heuristics, mining the frequency and inter-arrival time with which the heuristic visits local optima. Through these methodologies, we build a weighted directed graph that provides a synthetic view of the whole landscape, and that we can characterize using the tools of complex networks science. We argue that the network characterization can advance our understanding of the structural and dynamical properties of hard combinatorial landscapes. We apply our approach to prototypical problems such as the Quadratic Assignment Problem, the NK model of rugged landscapes, and the Permutation Flow-shop Scheduling Problem. We show that some network metrics can differentiate problem classes, correlate with problem non-linearity, and predict problem hardness as measured from the performances of trajectory-based local search heuristics.

Divergent and narrower climatic niches characterize polyploid species of European primroses in Primula sect. Aleuritia

Aim It is hypothesized that the ecological niches of polyploids should be both distinct and broader than those of diploids - characteristics that might have allowed the successful colonization of open habitats by polyploids during the Pleistocene glacial cycles. Here, we test these hypotheses by quantifying and comparing the ecological niches and niche breadths of a group of European primroses. Location Europe. Methods We gathered georeferenced data of four related species in Primula sect. Aleuritia at different ploidy levels (diploid, tetraploid, hexaploid and octoploid) and used seven bioclimatic variables to quantify niche overlap between species by applying a series of univariate and multivariate analyses combined with modelling techniques. We also employed permutation-based tests to evaluate niche similarity between the four species. Niche breadth for each species was evaluated both in the multivariate environmental space and in geographical space. Results The four species differed significantly from each other in mono-dimensional comparisons of climatological variables and occupied distinct habitats in the multi-dimensional environmental space. The majority of the permutation-based tests either indicated that the four species differed significantly in their habitat preferences and ecological niches or did not support significant niche similarity. Furthermore, our results revealed narrower niche breadths and geographical ranges in species of P. sect. Aleuritia at higher ploidy levels. Main conclusions The detected ecological differentiation between the four species of P. sect. Aleuritia at different ploidy levels is consistent with the hypothesis that polyploids occupy distinct ecological niches that differ from those of their diploid relative. Contrary to expectations, we find that polyploid species of P. sect. Aleuritia occupy narrower environmental and geographical spaces than their diploid relative. These results on the ecological niches of closely related polyploid and diploid species highlight factors that potentially contribute to the evolution and distribution of polyploid species.

A multi-contrast MRI study of microstructural brain damage in patients with mild cognitive impairment.

OBJECTIVES: The aim of this study was to investigate pathological mechanisms underlying brain tissue alterations in mild cognitive impairment (MCI) using multi-contrast 3 T magnetic resonance imaging (MRI). METHODS: Forty-two MCI patients and 77 healthy controls (HC) underwent T1/T2* relaxometry as well as Magnetization Transfer (MT) MRI. Between-groups comparisons in MRI metrics were performed using permutation-based tests. Using MRI data, a generalized linear model (GLM) was computed to predict clinical performance and a support-vector machine (SVM) classification was used to classify MCI and HC subjects. RESULTS: Multi-parametric MRI data showed microstructural brain alterations in MCI patients vs HC that might be interpreted as: (i) a broad loss of myelin/cellular proteins and tissue microstructure in the hippocampus (p ≤ 0.01) and global white matter (p < 0.05); and (ii) iron accumulation in the pallidus nucleus (p ≤ 0.05). MRI metrics accurately predicted memory and executive performances in patients (p ≤ 0.005). SVM classification reached an accuracy of 75% to separate MCI and HC, and performed best using both volumes and T1/T2*/MT metrics. CONCLUSION: Multi-contrast MRI appears to be a promising approach to infer pathophysiological mechanisms leading to brain tissue alterations in MCI. Likewise, parametric MRI data provide powerful correlates of cognitive deficits and improve automatic disease classification based on morphometric features.

Weakened functional connectivity in patients with psychogenic non-epileptic seizures (PNES) converges on basal ganglia.

BACKGROUND: Psychogenic non-epileptic seizures (PNES) are involuntary paroxysmal events that are unaccompanied by epileptiform EEG discharges. We hypothesised that PNES are a disorder of distributed brain networks resulting from their functional disconnection.The disconnection may underlie a dissociation mechanism that weakens the influence of unconsciously presented traumatising information but exerts maladaptive effects leading to episodic failures of behavioural control manifested by psychogenic 'seizures'. METHODS: To test this hypothesis, we compared functional connectivity (FC) derived from resting state high-density EEGs of 18 patients with PNES and 18 age-matched and gender-matched controls. To this end, the EEGs were transformed into source space using the local autoregressive average inverse solution. FC was estimated with a multivariate measure of lagged synchronisation in the θ, α and β frequency bands for 66 brain sites clustered into 18 regions. A multiple comparison permutation test was applied to deduce significant between-group differences in inter-regional and intraregional FC. RESULTS: The significant effect of PNES-a decrease in lagged FC between the basal ganglia and limbic, prefrontal, temporal, parietal and occipital regions-was found in the α band. CONCLUSION: We believe that this finding reveals a possible neurobiological substrate of PNES, which explains both attenuation of the effect of potentially disturbing mental representations and the occurrence of PNES episodes. By improving understanding of the aetiology of this condition, our results suggest a potential refinement of diagnostic criteria and management principles.