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.

In conspectu prope totius urbis: an application of different visual methods at the ager Tarraconensis landscape

We present in this paper the results of the application of several visual methods on a group of locations, dated between VI and I centuries BC, of the ager Tarraconensis (Tarragona, Spain) a Hinterland of the roman colony of Tarraco. The difficulty in interpreting the diverse results in a combined way has been resolved by means of the use of statistical methods, such as Principal Components Analysis (PCA) and K-means clustering analysis. These methods have allowed us to carry out site classifications in function of the landscape's visual structure that contains them and of the visual relationships that could be given among them.

Normative Topographic ERP Analyses of Speed of Speech Processing and Grammar Before and After Grammatical Treatment.

Typically developing (TD) preschoolers and age-matched preschoolers with specific language impairment (SLI) received event-related potentials (ERPs) to four monosyllabic speech sounds prior to treatment and, in the SLI group, after 6 months of grammatical treatment. Before treatment, the TD group processed speech sounds faster than the SLI group. The SLI group increased the speed of their speech processing after treatment. Posttreatment speed of speech processing predicted later impairment in comprehending phrase elaboration in the SLI group. During the treatment phase, change in speed of speech processing predicted growth rate of grammar in the SLI group.

Sugarcane response to nitrogen rates, measured by a canopy reflectance sensor

Abstract:The objective of this work was to evaluate whether a canopy sensor is capable of estimating sugarcane response to N, as well as to propose strategies for handling the data generated by this device during the decision-making process for crop N fertilization. Four N rate-response experiments were carried out, with N rates varying from 0 to 240 kg ha-1. Two evaluations with the canopy sensor were performed when the plants reached average stalk height of 0.3 and 0.5 m. Only two experiments showed stalk yield response to N rates. The canopy sensor was able to identify the crop response to different N rates and the relationship of the nutrient with sugarcane yield. The response index values obtained from the canopy sensor readings were useful in assessing sugarcane response to the applied N rate. Canopy reflectance sensors can help to identify areas responsive to N fertilization and, therefore, improve sugarcane fertilizer management.

Detection of Body Dump Sites and Clandestine Burials: a GIS-Based Landscape Approach

Forensic archaeology has been of inestimable help in the location and excavation of clandestine burials (Hunter 1996a:16-17, 1996b; Killam 2004; Levine et al. 1984). Landscape archaeology techniques have been adapted and widely applied for such purposes. In this article, an adaptation of one of the most common applications of archaeological GIS, that is predictive site location, is applied to the detection of body dump sites and clandestine burial sites, bringing together in the process the fields of landscape archaeology, forensic sciences, and GIS.

Linking landscape- and population ecology for large population management modelling: the case of Ibex (capra ibex) in Switzerland

La biologie de la conservation est communément associée à la protection de petites populations menacées d?extinction. Pourtant, il peut également être nécessaire de soumettre à gestion des populations surabondantes ou susceptibles d?une trop grande expansion, dans le but de prévenir les effets néfastes de la surpopulation. Du fait des différences tant quantitatives que qualitatives entre protection des petites populations et contrôle des grandes, il est nécessaire de disposer de modèles et de méthodes distinctes. L?objectif de ce travail a été de développer des modèles prédictifs de la dynamique des grandes populations, ainsi que des logiciels permettant de calculer les paramètres de ces modèles et de tester des scénarios de gestion. Le cas du Bouquetin des Alpes (Capra ibex ibex) - en forte expansion en Suisse depuis sa réintroduction au début du XXème siècle - servit d?exemple. Cette tâche fut accomplie en trois étapes : En premier lieu, un modèle de dynamique locale, spécifique au Bouquetin, fut développé : le modèle sous-jacent - structuré en classes d?âge et de sexe - est basé sur une matrice de Leslie à laquelle ont été ajoutées la densité-dépendance, la stochasticité environnementale et la chasse de régulation. Ce modèle fut implémenté dans un logiciel d?aide à la gestion - nommé SIM-Ibex - permettant la maintenance de données de recensements, l?estimation automatisée des paramètres, ainsi que l?ajustement et la simulation de stratégies de régulation. Mais la dynamique d?une population est influencée non seulement par des facteurs démographiques, mais aussi par la dispersion et la colonisation de nouveaux espaces. Il est donc nécessaire de pouvoir modéliser tant la qualité de l?habitat que les obstacles à la dispersion. Une collection de logiciels - nommée Biomapper - fut donc développée. Son module central est basé sur l?Analyse Factorielle de la Niche Ecologique (ENFA) dont le principe est de calculer des facteurs de marginalité et de spécialisation de la niche écologique à partir de prédicteurs environnementaux et de données d?observation de l?espèce. Tous les modules de Biomapper sont liés aux Systèmes d?Information Géographiques (SIG) ; ils couvrent toutes les opérations d?importation des données, préparation des prédicteurs, ENFA et calcul de la carte de qualité d?habitat, validation et traitement des résultats ; un module permet également de cartographier les barrières et les corridors de dispersion. Le domaine d?application de l?ENFA fut exploré par le biais d?une distribution d?espèce virtuelle. La comparaison à une méthode couramment utilisée pour construire des cartes de qualité d?habitat, le Modèle Linéaire Généralisé (GLM), montra qu?elle était particulièrement adaptée pour les espèces cryptiques ou en cours d?expansion. Les informations sur la démographie et le paysage furent finalement fusionnées en un modèle global. Une approche basée sur un automate cellulaire fut choisie, tant pour satisfaire aux contraintes du réalisme de la modélisation du paysage qu?à celles imposées par les grandes populations : la zone d?étude est modélisée par un pavage de cellules hexagonales, chacune caractérisée par des propriétés - une capacité de soutien et six taux d?imperméabilité quantifiant les échanges entre cellules adjacentes - et une variable, la densité de la population. Cette dernière varie en fonction de la reproduction et de la survie locale, ainsi que de la dispersion, sous l?influence de la densité-dépendance et de la stochasticité. Un logiciel - nommé HexaSpace - fut développé pour accomplir deux fonctions : 1° Calibrer l?automate sur la base de modèles de dynamique (par ex. calculés par SIM-Ibex) et d?une carte de qualité d?habitat (par ex. calculée par Biomapper). 2° Faire tourner des simulations. Il permet d?étudier l?expansion d?une espèce envahisseuse dans un paysage complexe composé de zones de qualité diverses et comportant des obstacles à la dispersion. Ce modèle fut appliqué à l?histoire de la réintroduction du Bouquetin dans les Alpes bernoises (Suisse). SIM-Ibex est actuellement utilisé par les gestionnaires de la faune et par les inspecteurs du gouvernement pour préparer et contrôler les plans de tir. Biomapper a été appliqué à plusieurs espèces (tant végétales qu?animales) à travers le Monde. De même, même si HexaSpace fut initialement conçu pour des espèces animales terrestres, il pourrait aisément être étndu à la propagation de plantes ou à la dispersion d?animaux volants. Ces logiciels étant conçus pour, à partir de données brutes, construire un modèle réaliste complexe, et du fait qu?ils sont dotés d?une interface d?utilisation intuitive, ils sont susceptibles de nombreuses applications en biologie de la conservation. En outre, ces approches peuvent également s?appliquer à des questions théoriques dans les domaines de l?écologie des populations et du paysage.<br/><br/>Conservation biology is commonly associated to small and endangered population protection. Nevertheless, large or potentially large populations may also need human management to prevent negative effects of overpopulation. As there are both qualitative and quantitative differences between small population protection and large population controlling, distinct methods and models are needed. The aim of this work was to develop theoretical models to predict large population dynamics, as well as computer tools to assess the parameters of these models and to test management scenarios. The alpine Ibex (Capra ibex ibex) - which experienced a spectacular increase since its reintroduction in Switzerland at the beginning of the 20th century - was used as paradigm species. This task was achieved in three steps: A local population dynamics model was first developed specifically for Ibex: the underlying age- and sex-structured model is based on a Leslie matrix approach with addition of density-dependence, environmental stochasticity and culling. This model was implemented into a management-support software - named SIM-Ibex - allowing census data maintenance, parameter automated assessment and culling strategies tuning and simulating. However population dynamics is driven not only by demographic factors, but also by dispersal and colonisation of new areas. Habitat suitability and obstacles modelling had therefore to be addressed. Thus, a software package - named Biomapper - was developed. Its central module is based on the Ecological Niche Factor Analysis (ENFA) whose principle is to compute niche marginality and specialisation factors from a set of environmental predictors and species presence data. All Biomapper modules are linked to Geographic Information Systems (GIS); they cover all operations of data importation, predictor preparation, ENFA and habitat suitability map computation, results validation and further processing; a module also allows mapping of dispersal barriers and corridors. ENFA application domain was then explored by means of a simulated species distribution. It was compared to a common habitat suitability assessing method, the Generalised Linear Model (GLM), and was proven better suited for spreading or cryptic species. Demography and landscape informations were finally merged into a global model. To cope with landscape realism and technical constraints of large population modelling, a cellular automaton approach was chosen: the study area is modelled by a lattice of hexagonal cells, each one characterised by a few fixed properties - a carrying capacity and six impermeability rates quantifying exchanges between adjacent cells - and one variable, population density. The later varies according to local reproduction/survival and dispersal dynamics, modified by density-dependence and stochasticity. A software - named HexaSpace - was developed, which achieves two functions: 1° Calibrating the automaton on the base of local population dynamics models (e.g., computed by SIM-Ibex) and a habitat suitability map (e.g. computed by Biomapper). 2° Running simulations. It allows studying the spreading of an invading species across a complex landscape made of variously suitable areas and dispersal barriers. This model was applied to the history of Ibex reintroduction in Bernese Alps (Switzerland). SIM-Ibex is now used by governmental wildlife managers to prepare and verify culling plans. Biomapper has been applied to several species (both plants and animals) all around the World. In the same way, whilst HexaSpace was originally designed for terrestrial animal species, it could be easily extended to model plant propagation or flying animals dispersal. As these softwares were designed to proceed from low-level data to build a complex realistic model and as they benefit from an intuitive user-interface, they may have many conservation applications. Moreover, theoretical questions in the fields of population and landscape ecology might also be addressed by these approaches.

Modelling bird species richness with neural networks for forest landscape management in NE Spain

Para preservar la biodiversidad de los ecosistemas forestales de la Europa mediterránea en escenarios actuales y futuros de cambio global mediante una gestión forestal sostenible es necesario determinar cómo influye el medio ambiente y las propias características de los bosques sobre la biodiversidad que éstos albergan. Con este propósito, se analizó la influencia de diferentes factores ambientales y de estructura y composición del bosque sobre la riqueza de aves forestales a escala 1 × 1 km en Cataluña (NE de España). Se construyeron modelos univariantes y multivariantes de redes neuronales para respectivamente explorar la respuesta individual a las variables y obtener un modelo parsimonioso (ecológicamente interpretable) y preciso. La superficie de bosque (con una fracción de cabida cubierta superior a 5%), la fracción de cabida cubierta media, la temperatura anual y la precipitación estival medias fueron los mejores predictores de la riqueza de aves forestales. La red neuronal multivariante obtenida tuvo una buena capacidad de generalización salvo en las localidades con una mayor riqueza. Además, los bosques con diferentes grados de apertura del dosel arbóreo, más maduros y más diversos en cuanto a su composición de especies arbóreas se asociaron de forma positiva con una mayor riqueza de aves forestales. Finalmente, se proporcionan directrices de gestión para la planificación forestal que permitan promover la diversidad ornítica en esta región de la Europa mediterránea.

Food traditions and landscape histories of the Indian Ocean World : theoretical and methodological reflections

Environmental histories of plant exchanges have largely centred on their eco- nomic importance in international trade and on their ecological and social impacts in the places where they were introduced. Yet few studies have at- tempted to examine how plants brought from elsewhere become incorporated over time into the regional cultures of material life and agricultural landscapes. This essay considers the theoretical and methodological problems in inves- tigating the environmental history, diversity and distribution of food plants transferred across the Indian Ocean over several millennia. It brings together concepts of creolisation, syncretism, and hybridity to outline a framework for understanding how biotic exchanges and diffusions have been translated into regional landscape histories through food traditions, ritual practices and articu- lation of cultural identity. We use the banana plant - which underwent early domestication across New Guinea, South-east Asia and peninsular India and reached East Africa roughly two thousand years ago - as an example for il- lustrating the diverse patterns of incorporation into the cultural symbolism, material life and regional landscapes of the Indian Ocean World. We show that this cultural evolutionary approach allows new historical insights to emerge and enriches ongoing debates regarding the antiquity of the plant's diffusion from South-east Asia to Africa.

Brain energy metabolism measured by (13) C magnetic resonance spectroscopy in vivo upon infusion of [3-(13) C]lactate.

The brain uses lactate produced by glycolysis as an energy source. How lactate originated from the blood stream is used to fuel brain metabolism is not clear. The current study measures brain metabolic fluxes and estimates the amount of pyruvate that becomes labeled in glial and neuronal compartments upon infusion of [3-(13) C]lactate. For that, labeling incorporation into carbons of glutamate and glutamine was measured by (13) C magnetic resonance spectroscopy at 14.1 T and analyzed with a two-compartment model of brain metabolism to estimate rates of mitochondrial oxidation, glial pyruvate carboxylation, and the glutamate-glutamine cycle as well as pyruvate fractional enrichments. Extracerebral lactate at supraphysiological levels contributes at least two-fold more to replenish the neuronal than the glial pyruvate pools. The rates of mitochondrial oxidation in neurons and glia, pyruvate carboxylase, and glutamate-glutamine cycles were similar to those estimated by administration of (13) C-enriched glucose, the main fuel of brain energy metabolism. These results are in agreement with primary utilization of exogenous lactate in neurons rather than astrocytes. © 2014 Wiley Periodicals, Inc.

Aortic blood pressure measured via EIT: investigation of different measurement settings.

Electrical impedance tomography (EIT) allows the measurement of intra-thoracic impedance changes related to cardiovascular activity. As a safe and low-cost imaging modality, EIT is an appealing candidate for non-invasive and continuous haemodynamic monitoring. EIT has recently been shown to allow the assessment of aortic blood pressure via the estimation of the aortic pulse arrival time (PAT). However, finding the aortic signal within EIT image sequences is a challenging task: the signal has a small amplitude and is difficult to locate due to the small size of the aorta and the inherent low spatial resolution of EIT. In order to most reliably detect the aortic signal, our objective was to understand the effect of EIT measurement settings (electrode belt placement, reconstruction algorithm). This paper investigates the influence of three transversal belt placements and two commonly-used difference reconstruction algorithms (Gauss-Newton and GREIT) on the measurement of aortic signals in view of aortic blood pressure estimation via EIT. A magnetic resonance imaging based three-dimensional finite element model of the haemodynamic bio-impedance properties of the human thorax was created. Two simulation experiments were performed with the aim to (1) evaluate the timing error in aortic PAT estimation and (2) quantify the strength of the aortic signal in each pixel of the EIT image sequences. Both experiments reveal better performance for images reconstructed with Gauss-Newton (with a noise figure of 0.5 or above) and a belt placement at the height of the heart or higher. According to the noise-free scenarios simulated, the uncertainty in the analysis of the aortic EIT signal is expected to induce blood pressure errors of at least ± 1.4 mmHg.

Land use changes, landscape ecology and their socioeconomic driving forces in the Spanish Mediterranean coast (the Maresme County, 1850-2005)

[cat] Utilitzem un conjunt de mètriques del paisatge per estudiar l'evolució a llarg termini seguida en una típica zona costanera del Mediterrani des de 1850 fins a 2005, que mostren una greu deterioració del medi ambient entre 1950 i 2005. Les principals forces motores d'aquesta degradació del paisatge han estat el creixement urbà experimentat a les antigues zones agrícoles situades a les planes litorals, juntament amb l'abandonament i la reforestació dels vessants dels pujols interceptats per àrees residencials de baixa densitat, carreteres i altres infraestructures lineals. Duem a terme una anàlisi estadística de redundància (RDA) amb la finalitat d'identificar els que considerem com alguns agents rectors socioeconòmics i polítics d'última instància d'aquests impactes ambientals. Els resultats confirmen les nostres hipòtesis interpretatives, que són que: 1) els canvis en les cobertes i usos del sòl determinen canvis en les propietats dels paisatge, tant estructurals com funcionals; 2) aquests canvis no es produeixen per atzar, sinó que estan relacionats amb factors geogràfics i forces socioeconòmiques i polítiques.

Land use changes, landscape ecology and their socioeconomic driving forces in the Spanish Mediterranean coast (the Maresme County, 1850-2005)

[cat] Utilitzem un conjunt de mètriques del paisatge per estudiar l'evolució a llarg termini seguida en una típica zona costanera del Mediterrani des de 1850 fins a 2005, que mostren una greu deterioració del medi ambient entre 1950 i 2005. Les principals forces motores d'aquesta degradació del paisatge han estat el creixement urbà experimentat a les antigues zones agrícoles situades a les planes litorals, juntament amb l'abandonament i la reforestació dels vessants dels pujols interceptats per àrees residencials de baixa densitat, carreteres i altres infraestructures lineals. Duem a terme una anàlisi estadística de redundància (RDA) amb la finalitat d'identificar els que considerem com alguns agents rectors socioeconòmics i polítics d'última instància d'aquests impactes ambientals. Els resultats confirmen les nostres hipòtesis interpretatives, que són que: 1) els canvis en les cobertes i usos del sòl determinen canvis en les propietats dels paisatge, tant estructurals com funcionals; 2) aquests canvis no es produeixen per atzar, sinó que estan relacionats amb factors geogràfics i forces socioeconòmiques i polítiques.

The Ovarian Cancer Chemokine Landscape Is Conducive to Homing of Vaccine-Primed and CD3/CD28-Costimulated T Cells Prepared for Adoptive Therapy.

PURPOSE: Chemokines are implicated in T-cell trafficking. We mapped the chemokine landscape in advanced stage ovarian cancer and characterized the expression of cognate receptors in autologous dendritic cell (DC)-vaccine primed T cells in the context of cell-based immunotherapy. EXPERIMENTAL DESIGN: The expression of all known human chemokines in patients with primary ovarian cancer was analyzed on two independent microarray datasets and validated on tissue microarray. Peripheral blood T cells from five HLA-A2 patients with recurrent ovarian cancer, who previously received autologous tumor DC vaccine, underwent CD3/CD28 costimulation and expansion ex vivo. Tumor-specific T cells were identified by HER2/neu pentamer staining and were evaluated for the expression and functionality of chemokine receptors important for homing to ovarian cancer. RESULTS: The chemokine landscape of ovarian cancer is heterogeneous with high expression of known lymphocyte-recruiting chemokines (CCL2, CCL4, and CCL5) in tumors with intraepithelial T cells, whereas CXCL10, CXCL12, and CXCL16 are expressed quasi-universally, including in tumors lacking tumor-infiltrating T cells. DC-vaccine primed T cells were found to express the cognate receptors for the above chemokines. Ex vivo CD3/CD28 costimulation and expansion of vaccine-primed Tcells upregulated CXCR3 and CXCR4, and enhanced their migration toward universally expressed chemokines in ovarian cancer. CONCLUSIONS: DC-primed tumor-specific T cells are armed with the appropriate receptors to migrate toward universal ovarian cancer chemokines, and these receptors are further upregulated by ex vivo CD3/CD28 costimulation, which render T cells more fit for migrating toward these chemokines. Clin Cancer Res; 21(12); 2840-50. ©2015 AACR.