Computation of multiple correspondence analysis, with code in R

The generalization of simple correspondence analysis, for two categorical variables, to multiple correspondence analysis where they may be three or more variables, is not straighforward, both from a mathematical and computational point of view. In this paper we detail the exact computational steps involved in performing a multiple correspondence analysis, including the special aspects of adjusting the principal inertias to correct the percentages of inertia, supplementary points and subset analysis. Furthermore, we give the algorithm for joint correspondence analysis where the cross-tabulations of all unique pairs of variables are analysed jointly. The code in the R language for every step of the computations is given, as well as the results of each computation.

Inundações urbanas: desafios ao ordenamento do território. O caso da Cidade da Praia (Cabo Verde)

Pela sua posição geográfica e particulares condições climáticas, devidas à sua inserção na faixa saheliana com caraterísticas de marcada aridez, Cabo Verde é um arquipélago com condições naturais adversas, pautado principalmente pela seca prolongada. Muitas vezes esta conjuntura é interpolada por curtos períodos de fortes chuvadas que podem originar cheias e inundações nos principais centros urbanos. Os eventos ocorridos revelam consequências graves, desde prejuízos na agricultura, perda de animais, destruição de infra-estruturas, perda de bens materiais e, mesmo, vítimas humanas mortais. Este estudo tem como objetivos principais: (i) perceber os problemas e desafios que se colocam à cidade da Praia (maior centro urbano do país, com forte crescimento e expansão urbana) perante situações de inundação; (ii) contribuir para o maior conhecimento das causas e consequências dessas inundações; (iii) definir quais as áreas de maior suscetibilidade às cheias e quais as que possuem um maior risco potencial. Optou-se por uma metodologia integrada, através do levantamento bibliográfico, cartográfico, numérico e percetivo (com base em entrevistas e inquéritos). Para o estudo das bacias hidrográficas foram calculados índices morfométricos, definidas classes de permeabilidade do substrato geológico e aplicado o método multicritério de Reis (2011) para a definição das áreas suscetíveis às cheias. Analisaram-se as precipitações máximas diárias anuais e respetivos períodos de retorno, com a aplicação do método de Gumbel. A análise de notícias de jornais, referentes ao período compreendido entre 1980 e 2011, foi fundamental para o conhecimento da distribuição espácio-temporal dos eventos perigosos de inundação em Cabo Verde e na cidade da Praia. Os resultados obtidos revelam um significativo grau de suscetibilidade às cheias na cidade da Praia. As áreas de maior risco potencial às inundações encontram-se no setor central da cidade, resultante da conjugação da convergência do escoamento das três ribeiras principais, da elevada densidade populacional e de construção desordenada nos leitos de cheia e nas áreas deprimidas, onde se acumulam as águas.

Inundações Urbanas - Desafios ao Ordenamento do Território - O

Pela sua posição geográfica e particulares condições climáticas, devidas à sua inserção na faixa saheliana com caraterísticas de marcada aridez, Cabo Verde é um arquipélago com condições naturais adversas, pautado principalmente pela seca prolongada. Muitas vezes esta conjuntura é interpolada por curtos períodos de fortes chuvadas que podem originar cheias e inundações nos principais centros urbanos. Os eventos ocorridos revelam consequências graves, desde prejuízos na agricultura, perda de animais, destruição de infra-estruturas, perda de bens materiais e, mesmo, vítimas humanas mortais. Este estudo tem como objetivos principais: (i) perceber os problemas e desafios que se colocam à cidade da Praia (maior centro urbano do país, com forte crescimento e expansão urbana) perante situações de inundação; (ii) contribuir para o maior conhecimento das causas e consequências dessas inundações; (iii) definir quais as áreas de maior suscetibilidade às cheias e quais as que possuem um maior risco potencial. Optou-se por uma metodologia integrada, através do levantamento bibliográfico, cartográfico, numérico e percetivo (com base em entrevistas e inquéritos). Para o estudo das bacias hidrográficas foram calculados índices morfométricos, definidas classes de permeabilidade do substrato geológico e aplicado o método multicritério de Reis (2011) para a definição das áreas suscetíveis às cheias. Analisaram-se as precipitações máximas diárias anuais e respetivos períodos de retorno, com a aplicação do método de Gumbel. A análise de notícias de jornais, referentes ao período compreendido entre 1980 e 2011, foi fundamental para o conhecimento da distribuição espácio-temporal dos eventos perigosos de inundação em Cabo Verde e na cidade da Praia. Os resultados obtidos revelam um significativo grau de suscetibilidade às cheias na cidade da Praia. As áreas de maior risco potencial às inundações encontram-se no setor central da cidade, resultante da conjugação da convergência do escoamento das três ribeiras principais, da elevada densidade populacional e de construção desordenada nos leitos de cheia e nas áreas deprimidas, onde se acumulam as águas.

Surviving in a competitive spatial market: The threshold capture model

Most facility location decision models ignore the fact that for a facility to survive it needs a minimum demand level to cover costs. In this paper we present a decision model for a firm thatwishes to enter a spatial market where there are several competitors already located. This market is such that for each outlet there is a demand threshold level that has to be achievedin order to survive. The firm wishes to know where to locate itsoutlets so as to maximize its market share taking into account the threshold level. It may happen that due to this new entrance, some competitors will not be able to meet the threshold and therefore will disappear. A formulation is presented together with a heuristic solution method and computational experience.

Positional information by differential endocytosis splits auxin response to drive Arabidopsis root meristem growth.

In the Arabidopsis root meristem, polar auxin transport creates a transcriptional auxin response gradient that peaks at the stem cell niche and gradually decreases as stem cell daughters divide and differentiate [1-3]. The amplitude and extent of this gradient are essential for both stem cell maintenance and root meristem growth [4, 5]. To investigate why expression of some auxin-responsive genes, such as the essential root meristem growth regulator BREVIS RADIX (BRX) [6], deviates from this gradient, we combined experimental and computational approaches. We created cellular-level root meristem models that accurately reproduce distribution of nuclear auxin activity and allow dynamic modeling of regulatory processes to guide experimentation. Expression profiles deviating from the auxin gradient could only be modeled after intersection of auxin activity with the observed differential endocytosis pattern and positive autoregulatory feedback through plasma-membrane-to-nucleus transfer of BRX. Because BRX is required for expression of certain auxin response factor targets, our data suggest a cell-type-specific endocytosis-dependent input into transcriptional auxin perception. This input sustains expression of a subset of auxin-responsive genes across the root meristem's division and transition zones and is essential for meristem growth. Thus, the endocytosis pattern provides specific positional information to modulate auxin response.

Development of LRFD Procedures for Bridge Pile Foundations in Iowa, September 2011

In response to the mandate on Load and Resistance Factor Design (LRFD) implementations by the Federal Highway Administration (FHWA) on all new bridge projects initiated after October 1, 2007, the Iowa Highway Research Board (IHRB) sponsored these research projects to develop regional LRFD recommendations. The LRFD development was performed using the Iowa Department of Transportation (DOT) Pile Load Test database (PILOT). To increase the data points for LRFD development, develop LRFD recommendations for dynamic methods, and validate the results of LRFD calibration, 10 full-scale field tests on the most commonly used steel H-piles (e.g., HP 10 x 42) were conducted throughout Iowa. Detailed in situ soil investigations were carried out, push-in pressure cells were installed, and laboratory soil tests were performed. Pile responses during driving, at the end of driving (EOD), and at re-strikes were monitored using the Pile Driving Analyzer (PDA), following with the CAse Pile Wave Analysis Program (CAPWAP) analysis. The hammer blow counts were recorded for Wave Equation Analysis Program (WEAP) and dynamic formulas. Static load tests (SLTs) were performed and the pile capacities were determined based on the Davisson’s criteria. The extensive experimental research studies generated important data for analytical and computational investigations. The SLT measured load displacements were compared with the simulated results obtained using a model of the TZPILE program and using the modified borehole shear test method. Two analytical pile setup quantification methods, in terms of soil properties, were developed and validated. A new calibration procedure was developed to incorporate pile setup into LRFD.

From subdiffusion to superdiffusion of particles on solid surfaces

We present a numerical and partially analytical study of classical particles obeying a Langevin equation that describes diffusion on a surface modeled by a two-dimensional potential. The potential may be either periodic or random. Depending on the potential and the damping, we observe superdiffusion, large-step diffusion, diffusion, and subdiffusion. Superdiffusive behavior is associated with low damping and is in most cases transient, albeit often long. Subdiffusive behavior is associated with highly damped particles in random potentials. In some cases subdiffusive behavior persists over our entire simulation and may be characterized as metastable. In any case, we stress that this rich variety of behaviors emerges naturally from an ordinary Langevin equation for a system described by ordinary canonical Maxwell-Boltzmann statistics.

Dynamics of the two-dimensional gonihedric spin model

In this paper, we study dynamical aspects of the two-dimensional (2D) gonihedric spin model using both numerical and analytical methods. This spin model has vanishing microscopic surface tension and it actually describes an ensemble of loops living on a 2D surface. The self-avoidance of loops is parametrized by a parameter ¿. The ¿=0 model can be mapped to one of the six-vertex models discussed by Baxter, and it does not have critical behavior. We have found that allowing for ¿¿0 does not lead to critical behavior either. Finite-size effects are rather severe, and in order to understand these effects, a finite-volume calculation for non-self-avoiding loops is presented. This model, like his 3D counterpart, exhibits very slow dynamics, but a careful analysis of dynamical observables reveals nonglassy evolution (unlike its 3D counterpart). We find, also in this ¿=0 case, the law that governs the long-time, low-temperature evolution of the system, through a dual description in terms of defects. A power, rather than logarithmic, law for the approach to equilibrium has been found.

Comparative performance of supertree algorithms in large data sets using the soapberry family (Sapindaceae) as a case study.

For the last 2 decades, supertree reconstruction has been an active field of research and has seen the development of a large number of major algorithms. Because of the growing popularity of the supertree methods, it has become necessary to evaluate the performance of these algorithms to determine which are the best options (especially with regard to the supermatrix approach that is widely used). In this study, seven of the most commonly used supertree methods are investigated by using a large empirical data set (in terms of number of taxa and molecular markers) from the worldwide flowering plant family Sapindaceae. Supertree methods were evaluated using several criteria: similarity of the supertrees with the input trees, similarity between the supertrees and the total evidence tree, level of resolution of the supertree and computational time required by the algorithm. Additional analyses were also conducted on a reduced data set to test if the performance levels were affected by the heuristic searches rather than the algorithms themselves. Based on our results, two main groups of supertree methods were identified: on one hand, the matrix representation with parsimony (MRP), MinFlip, and MinCut methods performed well according to our criteria, whereas the average consensus, split fit, and most similar supertree methods showed a poorer performance or at least did not behave the same way as the total evidence tree. Results for the super distance matrix, that is, the most recent approach tested here, were promising with at least one derived method performing as well as MRP, MinFlip, and MinCut. The output of each method was only slightly improved when applied to the reduced data set, suggesting a correct behavior of the heuristic searches and a relatively low sensitivity of the algorithms to data set sizes and missing data. Results also showed that the MRP analyses could reach a high level of quality even when using a simple heuristic search strategy, with the exception of MRP with Purvis coding scheme and reversible parsimony. The future of supertrees lies in the implementation of a standardized heuristic search for all methods and the increase in computing power to handle large data sets. The latter would prove to be particularly useful for promising approaches such as the maximum quartet fit method that yet requires substantial computing power.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.

The habenula encodes negative motivational value associated with primary punishment in humans.

Learning what to approach, and what to avoid, involves assigning value to environmental cues that predict positive and negative events. Studies in animals indicate that the lateral habenula encodes the previously learned negative motivational value of stimuli. However, involvement of the habenula in dynamic trial-by-trial aversive learning has not been assessed, and the functional role of this structure in humans remains poorly characterized, in part, due to its small size. Using high-resolution functional neuroimaging and computational modeling of reinforcement learning, we demonstrate positive habenula responses to the dynamically changing values of cues signaling painful electric shocks, which predict behavioral suppression of responses to those cues across individuals. By contrast, negative habenula responses to monetary reward cue values predict behavioral invigoration. Our findings show that the habenula plays a key role in an online aversive learning system and in generating associated motivated behavior in humans.

Stimulus statistics shape oscillations in nonlinear recurrent neural networks.

Rhythmic activity plays a central role in neural computations and brain functions ranging from homeostasis to attention, as well as in neurological and neuropsychiatric disorders. Despite this pervasiveness, little is known about the mechanisms whereby the frequency and power of oscillatory activity are modulated, and how they reflect the inputs received by neurons. Numerous studies have reported input-dependent fluctuations in peak frequency and power (as well as couplings across these features). However, it remains unresolved what mediates these spectral shifts among neural populations. Extending previous findings regarding stochastic nonlinear systems and experimental observations, we provide analytical insights regarding oscillatory responses of neural populations to stimulation from either endogenous or exogenous origins. Using a deceptively simple yet sparse and randomly connected network of neurons, we show how spiking inputs can reliably modulate the peak frequency and power expressed by synchronous neural populations without any changes in circuitry. Our results reveal that a generic, non-nonlinear and input-induced mechanism can robustly mediate these spectral fluctuations, and thus provide a framework in which inputs to the neurons bidirectionally regulate both the frequency and power expressed by synchronous populations. Theoretical and computational analysis of the ensuing spectral fluctuations was found to reflect the underlying dynamics of the input stimuli driving the neurons. Our results provide insights regarding a generic mechanism supporting spectral transitions observed across cortical networks and spanning multiple frequency bands.

Numerous potentially functional but non-genic conserved sequences on human chromosome 21.

The use of comparative genomics to infer genome function relies on the understanding of how different components of the genome change over evolutionary time. The aim of such comparative analysis is to identify conserved, functionally transcribed sequences such as protein-coding genes and non-coding RNA genes, and other functional sequences such as regulatory regions, as well as other genomic features. Here, we have compared the entire human chromosome 21 with syntenic regions of the mouse genome, and have identified a large number of conserved blocks of unknown function. Although previous studies have made similar observations, it is unknown whether these conserved sequences are genes or not. Here we present an extensive experimental and computational analysis of human chromosome 21 in an effort to assign function to sequences conserved between human chromosome 21 (ref. 8) and the syntenic mouse regions. Our data support the presence of a large number of potentially functional non-genic sequences, probably regulatory and structural. The integration of the properties of the conserved components of human chromosome 21 to the rapidly accumulating functional data for this chromosome will improve considerably our understanding of the role of sequence conservation in mammalian genomes.