Statistical tests for scaling in the inter-event times of earthquakes in California

We explore in depth the validity of a recently proposed scaling law for earthquake inter-event time distributions in the case of the Southern California, using the waveform cross-correlation catalog of Shearer et al. Two statistical tests are used: on the one hand, the standard two-sample Kolmogorov-Smirnov test is in agreement with the scaling of the distributions. On the other hand, the one-sample Kolmogorov-Smirnov statistic complemented with Monte Carlo simulation of the inter-event times, as done by Clauset et al., supports the validity of the gamma distribution as a simple model of the scaling function appearing on the scaling law, for rescaled inter-event times above 0.01, except for the largest data set (magnitude greater than 2). A discussion of these results is provided.

Some hypothesis tests for the covariance matrix when the dimension is large compared to the sample size

This paper analyzes whether standard covariance matrix tests work whendimensionality is large, and in particular larger than sample size. Inthe latter case, the singularity of the sample covariance matrix makeslikelihood ratio tests degenerate, but other tests based on quadraticforms of sample covariance matrix eigenvalues remain well-defined. Westudy the consistency property and limiting distribution of these testsas dimensionality and sample size go to infinity together, with theirratio converging to a finite non-zero limit. We find that the existingtest for sphericity is robust against high dimensionality, but not thetest for equality of the covariance matrix to a given matrix. For thelatter test, we develop a new correction to the existing test statisticthat makes it robust against high dimensionality.

R functions for quantifying nonindependence in standard dyadic and SRM designs

Interdependence is the main feature of dyadic relationships and, in recent years, various statistical procedures have been proposed for quantifying and testing this social attribute in different dyadic designs. The purpose of this paper is to develop several functions for this kind of statistical tests in an R package, known as nonindependence, for use by applied social researchers. A Graphical User Interface (GUI) is also developed to facilitate the use of the functions included in this package. Examples drawn from psychological research and simulated data are used to illustrate how the software works.

A priori parameterisation of the CERES soil-crop models and tests against several European data sets

Mechanistic soil-crop models have become indispensable tools to investigate the effect of management practices on the productivity or environmental impacts of arable crops. Ideally these models may claim to be universally applicable because they simulate the major processes governing the fate of inputs such as fertiliser nitrogen or pesticides. However, because they deal with complex systems and uncertain phenomena, site-specific calibration is usually a prerequisite to ensure their predictions are realistic. This statement implies that some experimental knowledge on the system to be simulated should be available prior to any modelling attempt, and raises a tremendous limitation to practical applications of models. Because the demand for more general simulation results is high, modellers have nevertheless taken the bold step of extrapolating a model tested within a limited sample of real conditions to a much larger domain. While methodological questions are often disregarded in this extrapolation process, they are specifically addressed in this paper, and in particular the issue of models a priori parameterisation. We thus implemented and tested a standard procedure to parameterize the soil components of a modified version of the CERES models. The procedure converts routinely-available soil properties into functional characteristics by means of pedo-transfer functions. The resulting predictions of soil water and nitrogen dynamics, as well as crop biomass, nitrogen content and leaf area index were compared to observations from trials conducted in five locations across Europe (southern Italy, northern Spain, northern France and northern Germany). In three cases, the model’s performance was judged acceptable when compared to experimental errors on the measurements, based on a test of the model’s root mean squared error (RMSE). Significant deviations between observations and model outputs were however noted in all sites, and could be ascribed to various model routines. In decreasing importance, these were: water balance, the turnover of soil organic matter, and crop N uptake. A better match to field observations could therefore be achieved by visually adjusting related parameters, such as field-capacity water content or the size of soil microbial biomass. As a result, model predictions fell within the measurement errors in all sites for most variables, and the model’s RMSE was within the range of published values for similar tests. We conclude that the proposed a priori method yields acceptable simulations with only a 50% probability, a figure which may be greatly increased through a posteriori calibration. Modellers should thus exercise caution when extrapolating their models to a large sample of pedo-climatic conditions for which they have only limited information.

Is response to fire influenced by dietary specialization and mobility? A comparative study with multiple animal assemblages

Fire is a major agent involved in landscape transformation and an indirect cause of changes in species composition. Responses to fire may vary greatly depending on life histories and functional traits of species. We have examined the taxonomic and functional responses to fire of eight taxonomic animal groups displaying a gradient of dietary and mobility patterns: Gastropoda, Heteroptera, Formicidae, Coleoptera, Araneae, Orthoptera, Reptilia and Aves. The fieldwork was conducted in a Mediterranean protected area on 3 sites (one unburnt and two burnt with different postfire management practices) with five replicates per site. We collected information from 4606 specimens from 274 animal species. Similarity in species composition and abundance between areas was measured by the Bray-Curtis index and ANOSIM, and comparisons between animal and plant responses by Mantel tests. We analyze whether groups with the highest percentage of omnivorous species, these species being more generalist in their dietary habits, show weak responses to fire (i.e. more similarity between burnt and unburnt areas), and independent responses to changes in vegetation. We also explore how mobility, i.e. dispersal ability, influences responses to fire. Our results demonstrate that differences in species composition and abundance between burnt and unburnt areas differed among groups. We found a tendency towards presenting lower differences between areas for groups with higher percentages of omnivorous species. Moreover, taxa with a higher percentage of omnivorous species had significantly more independent responses of changes in vegetation. High- (e.g. Aves) and low-mobility (e.g. Gastropoda) groups had the strongest responses to fire (higher R scores of the ANOSIM); however, we failed to find a significant general pattern with all the groups according to their mobility. Our results partially support the idea that functional traits underlie the response of organisms to environmental changes caused by fire.

The adequacy of different robust statistical tests in comparing two independent groups

In the current study, we evaluated various robust statistical methods for comparing two independent groups. Two scenarios for simulation were generated: one of equality and another of population mean differences. In each of the scenarios, 33 experimental conditions were used as a function of sample size, standard deviation and asymmetry. For each condition, 5000 replications per group were generated. The results obtained by this study show an adequate type error I rate but not a high power for the confidence intervals. In general, for the two scenarios studied (mean population differences and not mean population differences) in the different conditions analysed, the Mann-Whitney U-test demonstrated strong performance, and a little worse the t-test of Yuen-Welch.

Hausman tests for inefficient estimators: application to demand for health care service (revised)

The Hausman (1978) test is based on the vector of differences of two estimators. It is usually assumed that one of the estimators is fully efficient, since this simplifies calculation of the test statistic. However, this assumption limits the applicability of the test, since widely used estimators such as the generalized method of moments (GMM) or quasi maximum likelihood (QML) are often not fully efficient. This paper shows that the test may easily be implemented, using well-known methods, when neither estimator is efficient. To illustrate, we present both simulation results as well as empirical results for utilization of health care services.

Thermal tests on the LISA pathfinder engineering model

Report for the scientific sojourn carried out at Albert Einstein Institut in Germany, from April to July 2006.

Disseny d'itineraris d'educació ambiental en ambients de fons de vall

L’anàlisi de l’estat actual dels equipaments i activitats d’educació ambiental (EA) en l’àmbit del Parc Natural de l’Alt Pirineu (PNAP) ha detectat un grau de desenvolupament molt baix d’aquesta disciplina. Per aquest motiu, s’elabora una nova eina, amb l’objectiu d’establir unes pautes per a l’elaboració d’una activitat d’EA concreta: els itineraris d’EA; i així ampliar l’escassa oferta actual. El protocol per al disseny d’itineraris d’EA, dóna unes directrius per a estandarditzar la descripció i l’anàlisi de la zona proposada. La finalitat del qual és facilitar i agilitar la implantació d’itineraris que presentin unes condicions mínimes establertes. La eficiència d’aquesta eina, queda verificada amb els resultats obtinguts en la realització de dues proves pilot en dues zones amb possibilitats de desenvolupar-hi itineraris d’EA. Una de les zones ha quedat descartada, per l’incompliment d’un dels criteris que s’ha definit en el protocol com a imprescindibles pel bon desenvolupament de l’activitat que es planteja. Mentre que en l’altre zona, el compliment de totes aquestes condicions o criteris, fan viable la implantació de l’itinerari d’EA, duent a terme així el seu posterior disseny. Malgrat l’estudi es centri en un espai natural protegit concret, és una eina que pretén ser extrapolable en d’altres espais d’aquesta categoria.

Efficient knowledge retrieval to calibrate input variables in forest fire prediction

Forest fires are a serious threat to humans and nature from an ecological, social and economic point of view. Predicting their behaviour by simulation still delivers unreliable results and remains a challenging task. Latest approaches try to calibrate input variables, often tainted with imprecision, using optimisation techniques like Genetic Algorithms. To converge faster towards fitter solutions, the GA is guided with knowledge obtained from historical or synthetical fires. We developed a robust and efficient knowledge storage and retrieval method. Nearest neighbour search is applied to find the fire configuration from knowledge base most similar to the current configuration. Therefore, a distance measure was elaborated and implemented in several ways. Experiments show the performance of the different implementations regarding occupied storage and retrieval time with overly satisfactory results.

New evidence of the real interest rate parity for OECD countries using panel unit root tests with breaks

This paper tests for real interest parity (RIRP) among the nineteen major OECD countries over the period 1978:Q2-1998:Q4. The econometric methods applied consist of combining the use of several unit root or stationarity tests designed for panels valid under cross-section dependence and presence of multiple structural breaks. Our results strongly support the fulfillment of the weak version of the RIRP for the studied period once dependence and structural breaks are accounted for.

Memòria de recerca: paral·lelització dels diccionaris de la Standard Template Library de C++

Estudi elaborat a partir d’una estada a la Universität Karlsruhe entre gener i maig del 2007. Les biblioteques d’estructures de dades defineixen interfícies i implementen algorismes i estructures de dades fonamentals. Un exemple n’és la Satandard Template Library (STL ), que forma part del llenguatge de programació C++. En el marc d’una tesi, s’està treballant per obtenir implementacions més eficients i/o versàtils d’alguns components de la STL. Per a fer-ho s’utilitzen tècniques de la enginyeria d’algorismes. En particular, s’integra el coneixement de la comunitat algorítmica i es té en consideració la tecnologia existent. L’acció durant l’estada s’ha emmarcat en el desenvolupament la Multi Core STL (MCSTL ). La MCSTL és una implementació paral•lela de la STL per a màquines multi-core. Les màquines multi-core són actualment l’únic tipus de màquina disponible al mercat. Per tant, tot i que el paral•lelisme obtingut no sigui òptim, és preferible a tenir els processadors esperant, ja que , la tendència és que el nombre de processadors per computador augmenti.

Modelling of bird distribution under complex mediterranean landscape dynamics: open habitat birds and fire (MEDFIRE)

Report for the scientific sojourn at the Simon Fraser University, Canada, from July to September 2007. General context: landscape change during the last years is having significant impacts on biodiversity in many Mediterranean areas. Land abandonment, urbanisation and specially fire are profoundly transforming large areas in the Western Mediterranean basin and we know little on how these changes influence species distribution and in particular how these species will respond to further change in a context of global change including climate. General objectives: integrate landscape and population dynamics models in a platform allowing capturing species distribution responses to landscape changes and assessing impact on species distribution of different scenarios of further change. Specific objective 1: develop a landscape dynamic model capturing fire and forest succession dynamics in Catalonia and linked to a stochastic landscape occupancy (SLOM) (or spatially explicit population, SEPM) model for the Ortolan bunting, a species strongly linked to fire related habitat in the region. Predictions from the occupancy or spatially explicit population Ortolan bunting model (SEPM) should be evaluated using data from the DINDIS database. This database tracks bird colonisation of recently burnt big areas (&50 ha). Through a number of different SEPM scenarios with different values for a number of parameter, we should be able to assess different hypothesis in factors driving bird colonisation in new burnt patches. These factors to be mainly, landscape context (i.e. difficulty to reach the patch, and potential presence of coloniser sources), dispersal constraints, type of regenerating vegetation after fire, and species characteristics (niche breadth, etc).

Automatització de tests per a una aplicació J2EE

Aquest projecte consisteix en la creació d’una eina que permeti automatitzar els tests que s’han de fer a una aplicació comercial j2ee, amb el propòsit de facilitar i estalviar feina a les persones encarregades de testejar aquesta aplicació, ajudant-les així en la seva tasca de cerca d’errors. Concretament, s’ha creat una aplicació construïda per capes, modulable i fàcilment ampliable, la qual arriba més enllà de l’automatització dels tests més habituals, permetent executar un conjunt de tests per tal de validar si la versió de l’aplicació a testejar és vàlida o no.

Application of standard and refined heat balance integral methods to one-dimensional Stefan problems

The work in this paper concerns the study of conventional and refined heat balance integral methods for a number of phase change problems. These include standard test problems, both with one and two phase changes, which have exact solutions to enable us to test the accuracy of the approximate solutions. We also consider situations where no analytical solution is available and compare these to numerical solutions. It is popular to use a quadratic profile as an approximation of the temperature, but we show that a cubic profile, seldom considered in the literature, is far more accurate in most circumstances. In addition, the refined integral method can give greater improvement still and we develop a variation on this method which turns out to be optimal in some cases. We assess which integral method is better for various problems, showing that it is largely dependent on the specified boundary conditions.