Bootstrapping pairs in Distance-Based Regression

La regressió basada en distàncies és un mètode de predicció que consisteix en dos passos: a partir de les distàncies entre observacions obtenim les variables latents, les quals passen a ser els regressors en un model lineal de mínims quadrats ordinaris. Les distàncies les calculem a partir dels predictors originals fent us d'una funció de dissimilaritats adequada. Donat que, en general, els regressors estan relacionats de manera no lineal amb la resposta, la seva selecció amb el test F usual no és possible. En aquest treball proposem una solució a aquest problema de selecció de predictors definint tests estadístics generalitzats i adaptant un mètode de bootstrap no paramètric per a l'estimació dels p-valors. Incluim un exemple numèric amb dades de l'assegurança d'automòbils.

The empirics of spatial competition: evidence from European regions

The New Economic Geography literature allows detailed analysis of the factors that determine the location decisions of firms in integrated markets. However, the competitive process is modelled in a rather rudimentary way, and the empirical evidence has usually been obtained from reduced-form econometric specifications. This study describes a structural model that takes into account strategic interactions between firms. We investigate the relationship between the degree of perceived competition ¿ not only from local firms but from firms in other regions ¿ and geographic concentration. The preliminary results indicate that, in aggregate terms, local firms present stronger competition than firms in other regions. Moreover, it is confirmed that greater geographical concentration of production reduces market power, due to the intensification of local competition; however, its impact on production costs is unclear.

Spatial variability of soil chemical properties after coffee tree removal

Assessing the spatial variability of soil chemical properties has become an important aspect of soil management strategies with a view to higher crop yields with minimal environmental degradation. This study was carried out at the Centro Experimental of the Instituto Agronomico, in Campinas, São Paulo, Brazil. The aim was to characterize the spatial variability of chemical properties of a Rhodic Hapludox on a recently bulldozer-cleaned area after over 30 years of coffee cultivation. Soil samples were collected in a 20 x 20 m grid with 36 sampling points across a 1 ha area in the layers 0.0-0.2 and 0.2-0.4 m to measure the following chemical properties: pH, organic matter, K+, P, Ca2+, Mg2+, potential acidity, NH4-N, and NO3-N. Descriptive statistics were applied to assess the central tendency and dispersion moments. Geostatistical methods were applied to evaluate and to model the spatial variability of variables by calculating semivariograms and kriging interpolation. Spatial dependence patterns defined by spherical model adjusted semivariograms were made for all cited soil properties. Moderate to strong degrees of spatial dependence were found between 31 and 60 m. It was still possible to map soil spatial variability properties in the layers 0-20 cm and 20-40 cm after plant removal with bulldozers.

Plant species richness and environmental heterogeneity in a mountain landscape: effects of variability and spatial configuration

The loss of biodiversity has become a matter of urgent concern and a better understanding of local drivers is crucial for conservation. Although environmental heterogeneity is recognized as an important determinant of biodiversity, this has rarely been tested using field data at management scale. We propose and provide evidence for the simple hypothesis that local species diversity is related to spatial environmental heterogeneity. Species partition the environment into habitats. Biodiversity is therefore expected to be influenced by two aspects of spatial heterogeneity: 1) the variability of environmental conditions, which will affect the number of types of habitat, and 2) the spatial configuration of habitats, which will affect the rates of ecological processes, such as dispersal or competition. Earlier, simulation experiments predicted that both aspects of heterogeneity will influence plant species richness at a particular site. For the first time, these predictions were tested for plant communities using field data, which we collected in a wooded pasture in the Swiss Jura mountains using a four-level hierarchical sampling design. Richness generally increased with increasing environmental variability and "roughness" (i.e. decreasing spatial aggregation). Effects occurred at all scales, but the nature of the effect changed with scale, suggesting a change in the underlying mechanisms, which will need to be taken into account if scaling up to larger landscapes. Although we found significant effects of environmental heterogeneity, other factors such as history could also be important determinants. If a relationship between environmental heterogeneity and species richness can be shown to be general, recently available high-resolution environmental data can be used to complement the assessment of patterns of local richness and improve the prediction of the effects of land use change based on mean site conditions or land use history.

Spatial and temporal variability of crop yield and some Rhodic Hapludox properties under no-tillage

Soil properties are closely related with crop production and spite of the measures implemented, spatial variation has been repeatedly observed and described. Identifying and describing spatial variations of soil properties and their effects on crop yield can be a powerful decision-making tool in specific land management systems. The objective of this research was to characterize the spatial and temporal variations in crop yield and chemical and physical properties of a Rhodic Hapludox soil under no-tillage. The studied area of 3.42 ha had been cultivated since 1985 under no-tillage crop rotation in summer and winter. Yield and soil property were sampled in a regular 10 x 10 m grid, with 302 sample points. Yields of several crops were analyzed (soybean, maize, triticale, hyacinth bean and castor bean) as well as soil chemical (pH, Soil Organic Matter (SOM), P, Ca2+, Mg2+, H + Al, B, Fe, Mn, Zn, CEC, sum of bases (SB), and base saturation (V %)) and soil physical properties (saturated hydraulic conductivity, texture, density, total porosity, and mechanical penetration resistance). Data were analyzed using geostatistical analysis procedures and maps based on interpolation by kriging. Great variation in crop yields was observed in the years evaluated. The yield values in the Northern region of the study area were high in some years. Crop yields and some physical and soil chemical properties were spatially correlated.

Growth and externalities across economies : an empirical analysis using spatial econometrics

Recent theoretical models of economic growth have emphasised the role of external effects on the accumulation of factors of production. Although most of the literature has considered the externalities across firms within a region, in this paper we go a step further and consider the possibility that these externalities cross the barriers of regional economies. We assess the role of these external effects in explaining growth and economic convergence. We present a simple growth model, which includes externalities across economies, developing a methodology for testing their existence and estimating their strength. In our view, spatial econometrics is naturally suited to an empirical consideration of these externalities. We obtain evidence on the presence of significant externalities both across Spanish and European regions.

Spatial variability of physical attributes of an alfisol under different hillslope curvatures

The influence of relief forms has been studied by several authors and explains the variability in the soil attributes of a landscape. Soil physical attributes depend on relief forms, and their assessment is important in mechanized agricultural systems, such as of sugarcane. This study aimed to characterize the spatial variability in the physical soil attributes and their relationship to the hillslope curvatures in an Alfisol developed from sandstone and growing sugarcane. Grids of 100 x 100 m were delimited in a convex and a concave area. The grids had a regular spacing of 10 x 10 m, and the crossing points of this spacing determined a total of 121 georeferenced sampling points. Samples were collected to determine the physical attributes related to soil aggregates, porosity, bulk density, resistance to penetration and moisture within the 0-0.2 and 0.2-0.4 m depth. Statistical analyses, geostatistics and Student's t-tests were performed with the means of the areas. All attributes, except aggregates > 2 mm in the 0-0.2 m depth and macroporosity at both depths, showed significant differences between the hillslope curvatures. The convex area showed the highest values of the mean weighted diameter, mean geometric diameter, aggregates > 2 mm, 1-2 mm aggregates, total porosity and moisture and lower values of bulk density and resistance to penetration in both depth compared to the concave area. The number of soil attributes with greater spatial variability was higher in the concave area.

Pedotransfer functions related to spatial variability of water retention attributes for lowland soils

The estimation of non available soil variables through the knowledge of other related measured variables can be achieved through pedotransfer functions (PTF) mainly saving time and reducing cost. Great differences among soils, however, can yield non desirable results when applying this method. This study discusses the application of developed PTFs by several authors using a variety of soils of different characteristics, to evaluate soil water contents of two Brazilian lowland soils. Comparisons are made between PTF evaluated data and field measured data, using statistical and geostatistical tools, like mean error, root mean square error, semivariogram, cross-validation, and regression coefficient. The eight tested PTFs to evaluate gravimetric soil water contents (Ug) at the tensions of 33 kPa and 1,500 kPa presented a tendency to overestimate Ug 33 kPa and underestimate Ug1,500 kPa. The PTFs were ranked according to their performance and also with respect to their potential in describing the structure of the spatial variability of the set of measured values. Although none of the PTFs have changed the distribution pattern of the data, all resulted in mean and variance statistically different from those observed for all measured values. The PTFs that presented the best predictive values of Ug33 kPa and Ug1,500 kPa were not the same that had the best performance to reproduce the structure of spatial variability of these variables.

Triathlon event distance specialization: training and injury effects.

We conducted a preliminary, questionnaire-based, retrospective analysis of training and injury in British National Squad Olympic distance (OD) and Ironman distance (IR) triathletes. The main outcome measures were training duration and training frequency and injury frequency and severity. The number of overuse injuries sustained over a 5-year period did not differ between OD and IR. However, the proportions of OD and IR athletes who were affected by injury to particular anatomical sites differed (p < 0.05). Also, fewer OD athletes (16.7 vs. 36.8%, p < 0.05) reported that their injury recurred. Although OD sustained fewer running injuries than IR (1.6 +/- 0.5 vs. 1.9 +/- 0.3, p < 0.05), more subsequently stopped running (41.7 vs. 15.8%) and for longer (33.5 +/- 43.0 vs. 16.7 +/- 16.6 days, p < 0.01). In OD, the number of overuse injuries sustained inversely correlated with percentage training time, and number of sessions, doing bike hill repetitions (r = -0.44 and -0.39, respectively, both p < 0.05). The IR overuse injury number correlated with the amount of intensive sessions done (r = 0.67, p < 0.01 and r = 0.56, p < 0.05 for duration of "speed run" and "speed bike" sessions). Coaches should note that training differences between triathletes who specialize in OD or IR competition may lead to their exhibiting differential risk for injury to specific anatomical sites. It is also important to note that cycle and run training may have a "cumulative stress" influence on injury risk. Therefore, the tendency of some triathletes to modify rather than stop training when injured-usually by increasing load in another discipline from that in which the injury first occurred-may increase both their risk of injury recurrence and time to full rehabilitation.

Public health expenditure and spatial interactions in a decentralized national health system

One of the limitations of cross-country health expenditure analysis refers to the fact that the financing, the internal organization and political restraints of health care decision-making are country-specific and heterogeneous. Yet, a potential solution is to examine the influence of such effects in those countries that have undertaken decentralization processes. In such a setting, it is possible to examine potential expenditure spillovers across the geography of a country as well as the influence of the political ideology of regional incumbents on public health expenditure. This paper examines the determinants of public health expenditure within Spanish region-states (Autonomous Communities, ACs), most of them subject to similar financing structures although exhibiting significant heterogeneity as a result of the increasing decentralization, region-specific political factors along with different use of health care inputs, economic dimension and spatial interactions

Testing Spatial Autocorrelation in Weighted Networks : the Modes Permutation Test

In a weighted spatial network, as specified by an exchange matrix, the variances of the spatial values are inversely proportional to the size of the regions. Spatial values are no more exchangeable under independence, thus weakening the rationale for ordinary permutation and bootstrap tests of spatial autocorrelation. We propose an alternative permutation test for spatial autocorrelation, based upon exchangeable spatial modes, constructed as linear orthogonal combinations of spatial values. The coefficients obtain as eigenvectors of the standardised exchange matrix appearing in spectral clustering, and generalise to the weighted case the concept of spatial filtering for connectivity matrices. Also, two proposals aimed at transforming an acessibility matrix into a exchange matrix with with a priori fixed margins are presented. Two examples (inter-regional migratory flows and binary adjacency networks) illustrate the formalism, rooted in the theory of spectral decomposition for reversible Markov chains.

Spatial, temporal and subcellular localization of islet-brain 1 (IB1), a homologue of JIP-1, in mouse brain.

Islet-brain 1 (IB1) was recently identified as a DNA-binding protein of the GLUT2 gene promoter. The mouse IB1 is the rat and human homologue of the Jun-interacting protein 1 (JIP-1) which has been recognized as a key player in the regulation of c-Jun amino-terminal kinase (JNK) mitogen-activated protein kinase (MAPK) pathways. JIP-1 is involved in the control of apoptosis and may play a role in brain development and aging. Here, IB1 was studied in adult and developing mouse brain tissue by in situ hybridization, Northern and Western blot analysis at cellular and subcellular levels, as well as by immunocytochemistry in brain sections and cell cultures. IB1 expression was localized in the synaptic regions of the olfactory bulb, retina, cerebral and cerebellar cortex and hippocampus in the adult mouse brain. IB1 was also detected in a restricted number of axons, as in the mossy fibres from dentate gyrus in the hippocampus, and was found in soma, dendrites and axons of cerebellar Purkinje cells. After birth, IB1 expression peaks at postnatal day 15. IB1 was located in axonal and dendritic growth cones in primary telencephalon cells. By biochemical and subcellular fractionation of neuronal cells, IB1 was detected both in the cytosolic and membrane fractions. Taken together with previous data, the restricted neuronal expression of IB1 in developing and adult brain and its prominent localization in synapses suggest that the protein may be critical for cell signalling in developing and mature nerve terminals.

Characterization of gravitational rock slope deformations at different spatial scales based on field, remote sensing and numerical approaches

Rock slope instabilities such as rock slides, rock avalanche or deep-seated gravitational slope deformations are widespread in Alpine valleys. These phenomena represent at the same time a main factor that control the mountain belts erosion and also a significant natural hazard that creates important losses to the mountain communities. However, the potential geometrical and dynamic connections linking outcrop and slope-scale instabilities are often unknown. A more detailed definition of the potential links will be essential to improve the comprehension of the destabilization processes and to dispose of a more complete hazard characterization of the rock instabilities at different spatial scales. In order to propose an integrated approach in the study of the rock slope instabilities, three main themes were analysed in this PhD thesis: (1) the inventory and the spatial distribution of rock slope deformations at regional scale and their influence on the landscape evolution, (2) the influence of brittle and ductile tectonic structures on rock slope instabilities development and (3) the characterization of hazard posed by potential rock slope instabilities through the development of conceptual instability models. To prose and integrated approach for the analyses of these topics, several techniques were adopted. In particular, high resolution digital elevation models revealed to be fundamental tools that were employed during the different stages of the rock slope instability assessment. A special attention was spent in the application of digital elevation model for detailed geometrical modelling of past and potential instabilities and for the rock slope monitoring at different spatial scales. Detailed field analyses and numerical models were performed to complete and verify the remote sensing approach. In the first part of this thesis, large slope instabilities in Rhone valley (Switzerland) were mapped in order to dispose of a first overview of tectonic and climatic factors influencing their distribution and their characteristics. Our analyses demonstrate the key influence of neotectonic activity and the glacial conditioning on the spatial distribution of the rock slope deformations. Besides, the volumes of rock instabilities identified along the main Rhone valley, were then used to propose the first estimate of the postglacial denudation and filling of the Rhone valley associated to large gravitational movements. In the second part of the thesis, detailed structural analyses of the Frank slide and the Sierre rock avalanche were performed to characterize the influence of brittle and ductile tectonic structures on the geometry and on the failure mechanism of large instabilities. Our observations indicated that the geometric characteristics and the variation of the rock mass quality associated to ductile tectonic structures, that are often ignored landslide study, represent important factors that can drastically influence the extension and the failure mechanism of rock slope instabilities. In the last part of the thesis, the failure mechanisms and the hazard associated to five potential instabilities were analysed in detail. These case studies clearly highlighted the importance to incorporate different analyses and monitoring techniques to dispose of reliable and hazard scenarios. This information associated to the development of a conceptual instability model represents the primary data for an integrated risk management of rock slope instabilities. - Les mouvements de versant tels que les chutes de blocs, les éboulements ou encore les phénomènes plus lents comme les déformations gravitaires profondes de versant représentent des manifestations courantes en régions montagneuses. Les mouvements de versant sont à la fois un des facteurs principaux contrôlant la destruction progressive des chaines orogéniques mais aussi un danger naturel concret qui peut provoquer des dommages importants. Pourtant, les phénomènes gravitaires sont rarement analysés dans leur globalité et les rapports géométriques et mécaniques qui lient les instabilités à l'échelle du versant aux instabilités locales restent encore mal définis. Une meilleure caractérisation de ces liens pourrait pourtant représenter un apport substantiel dans la compréhension des processus de déstabilisation des versants et améliorer la caractérisation des dangers gravitaires à toutes les échelles spatiales. Dans le but de proposer un approche plus globale à la problématique des mouvements gravitaires, ce travail de thèse propose trois axes de recherche principaux: (1) l'inventaire et l'analyse de la distribution spatiale des grandes instabilités rocheuses à l'échelle régionale, (2) l'analyse des structures tectoniques cassantes et ductiles en relation avec les mécanismes de rupture des grandes instabilités rocheuses et (3) la caractérisation des aléas rocheux par une approche multidisciplinaire visant à développer un modèle conceptuel de l'instabilité et une meilleure appréciation du danger . Pour analyser les différentes problématiques traitées dans cette thèse, différentes techniques ont été utilisées. En particulier, le modèle numérique de terrain s'est révélé être un outil indispensable pour la majorité des analyses effectuées, en partant de l'identification de l'instabilité jusqu'au suivi des mouvements. Les analyses de terrain et des modélisations numériques ont ensuite permis de compléter les informations issues du modèle numérique de terrain. Dans la première partie de cette thèse, les mouvements gravitaires rocheux dans la vallée du Rhône (Suisse) ont été cartographiés pour étudier leur répartition en fonction des variables géologiques et morphologiques régionales. En particulier, les analyses ont mis en évidence l'influence de l'activité néotectonique et des phases glaciaires sur la distribution des zones à forte densité d'instabilités rocheuses. Les volumes des instabilités rocheuses identifiées le long de la vallée principale ont été ensuite utilisés pour estimer le taux de dénudations postglaciaire et le remplissage de la vallée du Rhône lié aux grands mouvements gravitaires. Dans la deuxième partie, l'étude de l'agencement structural des avalanches rocheuses de Sierre (Suisse) et de Frank (Canada) a permis de mieux caractériser l'influence passive des structures tectoniques sur la géométrie des instabilités. En particulier, les structures issues d'une tectonique ductile, souvent ignorées dans l'étude des instabilités gravitaires, ont été identifiées comme des structures très importantes qui contrôlent les mécanismes de rupture des instabilités à différentes échelles. Dans la dernière partie de la thèse, cinq instabilités rocheuses différentes ont été étudiées par une approche multidisciplinaire visant à mieux caractériser l'aléa et à développer un modèle conceptuel trois dimensionnel de ces instabilités. A l'aide de ces analyses on a pu mettre en évidence la nécessité d'incorporer différentes techniques d'analyses et de surveillance pour une gestion plus objective du risque associée aux grandes instabilités rocheuses.