Testing for poverty dominance: an application to Canada

The paper proposes and applies statistical tests for poverty dominance that check for whether poverty comparisons can be made robustly over ranges of poverty lines and classes of poverty indices. This helps provide both normative and statistical confidence in establishing poverty rankings across distributions. The tests, which can take into account the complex sampling procedures that are typically used by statistical agencies to generate household-level surveys, are implemented using the Canadian Survey of Labour and Income Dynamics (SLID) for 1996, 1999 and 2002. Although the yearly cumulative distribution functions cross at the lower tails of the distributions, the more recent years tend to dominate earlier years for a relatively wide range of poverty lines. Failing to take into account SLID's sampling variability (as is sometimes done) can inflate significantly one's confidence in ranking poverty. Taking into account SLID's complex sampling design (as has not been done before) can also decrease substantially the range of poverty lines over which a poverty ranking can be inferred.

On the existence of hysteresis in the Kuramoto model with bimodal frequency distributions

We investigate the transition to synchronization in the Kuramoto model with bimodal distributions of the natural frequencies. Previous studies have concluded that the model exhibits a hysteretic phase transition if the bimodal distribution is close to a unimodal one, due to the shallowness the central dip. Here we show that proximity to the unimodal-bimodal border does not necessarily imply hysteresis when the width, but not the depth, of the central dip tends to zero. We draw this conclusion from a detailed study of the Kuramoto model with a suitable family of bimodal distributions.

Knowing the past to predict the future: land-use change and the distribution of invasive bullfrogs

Biological invasions and land-use changes are two major causes of the global modifications of biodiversity. Habitat suitability models are the tools of choice to predict potential distributions of invasive species. Although land-use is a key driver of alien species invasions, it is often assumed that land-use is constant in time. Here we combine historical and present day information, to evaluate whether land-use changes could explain the dynamic of invasion of the American bullfrog Rana catesbeiana (=Lithobathes catesbeianus) in Northern Italy, from the 1950s to present-day. We used maxent to build habitat suitability models, on the basis of past (1960s, 1980s) and present-day data on land-uses and species distribution. For example, we used models built using the 1960s data to predict distribution in the 1980s, and so on. Furthermore, we used land-use scenarios to project suitability in the future. Habitat suitability models predicted well the spread of bullfrogs in the subsequent temporal step. Models considering land-use changes predicted invasion dynamics better than models assuming constant land-use over the last 50 years. Scenarios of future land-use suggest that suitability will remain similar in the next years. Habitat suitability models can help to understand and predict the dynamics of invasions; however, land-use is not constant in time: land-use modifications can strongly affect invasions; furthermore, both land management and the suitability of a given land-use class may vary in time. An integration of land-use changes in studies of biological invasions can help to improve management strategies.

How to set up and apply reference levels in fluoroscopy at a national level.

A nationwide survey was launched to investigate the use of fluoroscopy and establish national reference levels (RL) for dose-intensive procedures. The 2-year investigation covered five radiology and nine cardiology departments in public hospitals and private clinics, and focused on 12 examination types: 6 diagnostic and 6 interventional. A total of 1,000 examinations was registered. Information including the fluoroscopy time (T), the number of frames (N) and the dose-area product (DAP) was provided. The data set was used to establish the distributions of T, N and the DAP and the associated RL values. The examinations were pooled to improve the statistics. A wide variation in dose and image quality in fixed geometry was observed. As an example, the skin dose rate for abdominal examinations varied in the range of 10 to 45 mGy/min for comparable image quality. A wide variability was found for several types of examinations, mainly complex ones. DAP RLs of 210, 125, 80, 240, 440 and 110 Gy cm2 were established for lower limb and iliac angiography, cerebral angiography, coronary angiography, biliary drainage and stenting, cerebral embolization and PTCA, respectively. The RL values established are compared to the data published in the literature.

Modelling the distribution of bats in relation to landscape structure in a temperate mountain environment

1. Landscape modification is often considered the principal cause of population decline in many bat species. Thus, schemes for bat conservation rely heavily on knowledge about species-landscape relationships. So far, however, few studies have quantified the possible influence of landscape structure on large-scale spatial patterns in bat communities. 2. This study presents quantitative models that use landscape structure to predict (i) spatial patterns in overall community composition and (ii) individual species' distributions through canonical correspondence analysis and generalized linear models, respectively. A geographical information system (GIS) was then used to draw up maps of (i) overall community patterns and (ii) distribution of potential species' habitats. These models relied on field data from the Swiss Jura mountains. 3. Fight descriptors of landscape structure accounted for 30% of the variation in bat community composition. For some species, more than 60% of the variance in distribution could be explained by landscape structure. Elevation, forest or woodland cover, lakes and suburbs, were the most frequent predictors. 4. This study shows that community composition in bats is related to landscape structure through species-specific relationships to resources. Due to their nocturnal activities and the difficulties of remote identification, a comprehensive bat census is rarely possible, and we suggest that predictive modelling of the type described here provides an indispensable conservation tool.

Understanding the low-temperature limitations to forest growth through calibration of a forest dynamics model with tree-ring data

The sensitivity of altitudinal and latitudinal tree-line ecotones to climate change, particularly that of temperature, has received much attention. To improve our understanding of the factors affecting tree-line position, we used the spatially explicit dynamic forest model TreeMig. Although well-suited because of its landscape dynamics functions, TreeMig features a parabolic temperature growth response curve, which has recently been questioned. and the species parameters are not specifically calibrated for cold temperatures. Our main goals were to improve the theoretical basis of the temperature growth response curve in the model and develop a method for deriving that curve's parameters from tree-ring data. We replaced the parabola with an asymptotic curve, calibrated for the main species at the subalpine (Swiss Alps: Pinus cembra, Larix decidua, Picea abies) and boreal (Fennoscandia: Pinus sylvestris, Betula pubescens, P. abies) tree-lines. After fitting new parameters, the growth curve matched observed tree-ring widths better. For the subalpine species, the minimum degree-day sum allowing, growth (kDDMin) was lowered by around 100 degree-days; in the case of Larix, the maximum potential ring-width was increased to 5.19 mm. At the boreal tree-line, the kDDMin for P. sylvestris was lowered by 210 degree-days and its maximum ring-width increased to 2.943 mm; for Betula (new in the model) kDDMin was set to 325 degree-days and the maximum ring-width to 2.51 mm; the values from the only boreal sample site for Picea were similar to the subalpine ones, so the same parameters were used. However, adjusting the growth response alone did not improve the model's output concerning species' distributions and their relative importance at tree-line. Minimum winter temperature (MinWiT, mean of the coldest winter month), which controls seedling establishment in TreeMig, proved more important for determining distribution. Picea, P. sylvestris and Betula did not previously have minimum winter temperature limits, so these values were set to the 95th percentile of each species' coldest MinWiT site (respectively -7, -11, -13). In a case study for the Alps, the original and newly calibrated versions of TreeMig were compared with biomass data from the National Forest Inventor), (NFI). Both models gave similar, reasonably realistic results. In conclusion, this method of deriving temperature responses from tree-rings works well. However, regeneration and its underlying factors seem more important for controlling species' distributions than previously thought. More research on regeneration ecology, especially at the upper limit of forests. is needed to improve predictions of tree-line responses to climate change further.

The Potential for Dispersal of Onchocerciasis in Ecuador in Relation to the Distribution of the Vector Simulium exiguum (Diptera:Simuliidae)

The future dispersal of onchocerciasis in Ecuador is dependent on the distribution of cytotypes of the vector species complex Simulium exiguum. Over the last 14 years, collections of larvae have been made from over 25 rivers, between 80-1600 m altitude, from various sites on both sides of the Andes. Analysis of larval polytene chromosomes was used to determine the distributions of each cytotype. On the western side of the Andes, the Cayapa cytotype (the only cytotype directly incriminated as a vector) has a distribution from Santo Domingo de los Colorados northwards. The Quevedo and Bucay cytotypes occur from Santo Domingo de los Colorados southwards. On the eastern side of the Andes, the Aguarico cytotype occurs in the Rio Aguarico and a new cytotype is present in the tributaries of the Rio Napo. Whether the disease will spread south of Santo Domingo and on the eastern side of the Andes depends on vector capacity of the cytotypes and the dispersal patterns of individuals infected with onchocerciasis. At present the Aguarico, Bucay and Quevedo cytotypes are known to be efficient hosts, but their biting preferences and biting densities have not yet been evaluated

Estimation of parametric and nonparametric models for univariate claim severity distributions - an approach using R

This paper presents an analysis of motor vehicle insurance claims relating to vehicle damage and to associated medical expenses. We use univariate severity distributions estimated with parametric and non-parametric methods. The methods are implemented using the statistical package R. Parametric analysis is limited to estimation of normal and lognormal distributions for each of the two claim types. The nonparametric analysis presented involves kernel density estimation. We illustrate the benefits of applying transformations to data prior to employing kernel based methods. We use a log-transformation and an optimal transformation amongst a class of transformations that produces symmetry in the data. The central aim of this paper is to provide educators with material that can be used in the classroom to teach statistical estimation methods, goodness of fit analysis and importantly statistical computing in the context of insurance and risk management. To this end, we have included in the Appendix of this paper all the R code that has been used in the analysis so that readers, both students and educators, can fully explore the techniques described

Universality of rain event size distributions

We compare rain event size distributions derived from measurements in climatically different regions, which we find to be well approximated by power laws of similar exponents over broad ranges. Differences can be seen in the large-scale cutoffs of the distributions. Event duration distributions suggest that the scale-free aspects are related to the absence of characteristic scales in the meteorological mesoscale.

Chemical and carbon isotopic evolution of hydrocarbons during prograde metamorphism from 100°C to 550°C: Case study in the Liassic black shale formation of Central Swiss Alps

Hydrocarbon distributions and stable isotope ratios of carbonates (delta(13)C(car), delta(18)O(car)), kerogen (delta(13)C(ker)), extractable organic matter (delta(13)C(EOM)) and individual hydrocarbons of Liassic black shale samples from a prograde metamorphic sequence in the Swiss Alps were used to identify the major organic reactions with increasing metamorphic grade. The studied samples range from the diagenetic zone (< 100 degrees C) to amphibolite facies (similar to 550 degrees C). The samples within the diagenetic zones (< 100 and 150 degrees C) are characterized by the dominance of C-< 20 n-alkanes, suggesting an origin related with marine and/or bacterial inputs. The metamorphic samples (200 to 550 degrees C) have distributions significantly dominated by C-12 and C-13 n-alkanes, C-14, C-16 and C-18 n-alkylcyclopentanes and to a lesser extend C-15, C-17 and C-21 n-alkylcyclohexanes. The progressive C-13-enrichment (up to 3.9 parts per thousand) with metamorphism of the C-> 17 n-alkanes suggests the occurrence of cracking reactions of high molecular weight compounds. The isotopically heavier (up to 5.6 parts per thousand) C-< 17 n-alkanes in metamorphic samples are likely originated by thermal degradation of long-chain homologous with preferential release of isotopically light C-1 and C-2 radicals. The dominance of specific even C-number n-alkylcyclopentanes suggests an origin related to direct cyclization mechanism (without decarboxylation step) of algal or bacterial fatty acids occurring in reducing aqueous metamorphic fluid conditions. The regular increase of the concentrations of n-alkylcycloalkanes vs. C-> 13 n-alkanes with metamorphism suggests progressive thermal release of kerogen-linked fatty acid precursors and degradation of n-alkanes. Changes of the steroid and terpenoid distributions are clearly related to increasing metamorphic temperatures. The absence of 18 alpha(H)-22,29,30-trisnorneohopane (Ts), the occurrence of 17 beta(H)-trisnorhopane, 17 beta(H), 21 alpha(H)-hopanes in the C-29 to C-31 range and 5 alpha(H),14 alpha(H),17 alpha(H)-20R C-27, C-29 steranes in the low diagenetic samples (< 100 degrees C) are characteristic of immature bitumens. The higher thermal stress within the upper diagenetic zone (150 degrees C) is marked by the presence of Ts, the disappearance of 17 beta(H)-trisnorhopane and thermodynamic equilibrium of the 22S/(22S + 22R) homohopane ratios. The increase of the alpha alpha alpha-sterane 20S/(20S + 20R) and 20R beta beta/(beta beta + alpha alpha) ratios (from 0.0 to 0.55 and from 0.0 to 0.40, respectively) in the upper diagenetic zone indicates the occurrence of isomerization reactions already at < 150 degrees C. However, the isomerization at C-20 (R -> S) reaches thermodynamic equilibrium values already at the upper diagenesis (similar to 150 degrees C) whereas the epimerisation at C-14 and C-17 (alpha alpha ->beta beta) arrives to constant values in the lower anchizone (similar to 200 degrees C). The ratios Ts vs. 17 alpha(H)-22,29,30-trisnorneohopane [(Ts/(Ts + Tm)] and 18 alpha(H)-30-norneohopane (C29Ts) vs. 17 alpha(H),21 beta(H)-30-norhopane [C29Ts/(C29Ts + C-29)] increase until the medium anchizone (200 to 250 degrees C) from 0.0 to 0.96 and from 0.0 to 0.44, respectively. An opposite trend owards lower values is observed in the higher metamorphic samples. The occurrence of specific hydrocarbons (e.g., n-alkylcyclopentanes, cadalene, hydrogenated aromatic compounds) in metamorphic samples points to kerogen degradation reactions most probably occurring in the presence of water and under reducing conditions. The changes of hydrocarbon distributions and carbon isotopic compositions of n-alkanes related to metamorphism suggest that the organic geochemistry may help to evaluate the lowest grades of prograde metamorphism. Copyright (c) 2005 Elsevier Ltd.

Break zones in the distributions of alleles and species in alpine plants

Aim  We test for the congruence between allele-based range boundaries (break zones) in silicicolous alpine plants and species-based break zones in the silicicolous flora of the European Alps. We also ask whether such break zones coincide with areas of large elevational variation.Location  The European Alps.Methods  On a regular grid laid across the entire Alps, we determined areas of allele- and species-based break zones using respective clustering algorithms, identifying discontinuities in cluster distributions (breaks), and quantifying integrated break densities (break zones). Discontinuities were identified based on the intra-specific genetic variation of 12 species and on the floristic distribution data from 239 species, respectively. Coincidence between the two types of break zones was tested using Spearman's correlation. Break zone densities were also regressed on topographical complexity to test for the effect of elevational variation.Results  We found that two main break zones in the distribution of alleles and species were significantly correlated. Furthermore, we show that these break zones are in topographically complex regions, characterized by massive elevational ranges owing to high mountains and deep glacial valleys. We detected a third break zone in the distribution of species in the eastern Alps, which is not correlated with topographic complexity, and which is also not evident from allelic distribution patterns. Species with the potential for long-distance dispersal tended to show larger distribution ranges than short-distance dispersers.Main conclusions  We suggest that the history of Pleistocene glaciations is the main driver of the congruence between allele-based and species-based distribution patterns, because occurrences of both species and alleles were subject to the same processes (such as extinction, migration and drift) that shaped the distributions of species and genetic lineages. Large elevational ranges have had a profound effect as a dispersal barrier for alleles during post-glacial immigration. Because plant species, unlike alleles, cannot spread via pollen but only via seed, and thus disperse less effectively, we conclude that species break zones are maintained over longer time spans and reflect more ancient patterns than allele break zones.Conny Thiel-Egenter and Nadir Alvarez contributed equally to this paper and are considered joint first authors.

Morphometric comparison of Simulium perflavum larvae (Diptera: Simuliidae) in relation to season and gender in Central Amazônia, Brazil

Number of larval instars, age structure and environmental effects on these parameters represent basic information in the study of insect population biology. When species have economic importance, this information is essential in order to choose the best period to apply different control methods and to determine the stages of the life cycle of the insect that are most susceptible to each treatment. The family Simuliidae has many species of medical/veterinary importance in the world, and some studies in the temperate region have suggested that the number of larval instars and the larval size can vary according to the season, gender and some environmental factors, such as temperature and diet. This study, with the zoophilic species Simulium perflavum Roubaud, is the first in the Neotropics observing some of these factors and will serve as a template for other species of medical importance in the region. S. perflavum larvae were collected in five streams in Central Amazônia (Manaus and Presidente Figueiredo counties, State of Amazonas), in Sept./Oct. 1996 (dry season) and Feb./Mar. 1997 (rainy season). These larvae were measured (lateral length of head capsule and width of cephalic apodema) to determine the number of larval instars (n=3985), to compare the larval size between seasons and genders (last and penultimate larval instars, n=200). Seven larval instars were determined for this species using frequency distributions, t-tests and Crosby´s growth rule. Significant differences were not detected (t-test, p>0.05) in larval size between seasons and genders. Our results differ from some found in temperate regions suggesting that in the Neotropical region the larval size in different seasons and different genders remains constant, although some environmental parameters, such as diet, change depending on the season.

Sensitivity of predictive species distribution models to change in grain size

Predictive species distribution modelling (SDM) has become an essential tool in biodiversity conservation and management. The choice of grain size (resolution) of environmental layers used in modelling is one important factor that may affect predictions. We applied 10 distinct modelling techniques to presence-only data for 50 species in five different regions, to test whether: (1) a 10-fold coarsening of resolution affects predictive performance of SDMs, and (2) any observed effects are dependent on the type of region, modelling technique, or species considered. Results show that a 10 times change in grain size does not severely affect predictions from species distribution models. The overall trend is towards degradation of model performance, but improvement can also be observed. Changing grain size does not equally affect models across regions, techniques, and species types. The strongest effect is on regions and species types, with tree species in the data sets (regions) with highest locational accuracy being most affected. Changing grain size had little influence on the ranking of techniques: boosted regression trees remain best at both resolutions. The number of occurrences used for model training had an important effect, with larger sample sizes resulting in better models, which tended to be more sensitive to grain. Effect of grain change was only noticeable for models reaching sufficient performance and/or with initial data that have an intrinsic error smaller than the coarser grain size.

An algebraic operator approach to the analysis of Gerber-Shiu functions

We introduce an algebraic operator framework to study discounted penalty functions in renewal risk models. For inter-arrival and claim size distributions with rational Laplace transform, the usual integral equation is transformed into a boundary value problem, which is solved by symbolic techniques. The factorization of the differential operator can be lifted to the level of boundary value problems, amounting to iteratively solving first-order problems. This leads to an explicit expression for the Gerber-Shiu function in terms of the penalty function.

Bats of Borneo: diversity, distributions and representation in protected areas

Protected areas are valuable in conserving tropical biodiversity, but an insufficient understanding of species diversity and distributions makes it difficult to evaluate their effectiveness. This is especially true on Borneo, a species rich island shared by three countries, and is particularly concerning for bats, a poorly known component of mammal diversity that may be highly susceptible to landscape changes. We reviewed the diversity, distributions and conservation status of 54 bat species to determine the representation of these taxa in Borneo's protected areas, and whether these reserves complement each other in terms of bat diversity. Lower and upper bound estimates of bat species composition were characterised in 23 protected areas and the proposed boundaries of the Heart of Borneo conservation area. We used lower and upper bound estimates of species composition. By using actual inventories, species representation was highly irregular, and even if some reserves were included in the Heart of Borneo, the protected area network would still exhibit low complementarity. By inferring species presence from distributions, composition between most reserves was similar, and complementarity was much higher. Predicting species richness using abundance information suggested that bat species representation in reserves may lie between these two extremes. We recommend that researchers better sample biodiversity over the island and address the conservation threats faced in Borneo both within and outside protected areas. While the Heart of Borneo Initiative is commendable, it should not divert attention from other conservation areas.