Modeling spatial distribution of cancer incidence in Cordoba, Argentina

Background: Cancer is the second leading cause of death in Argentina, and there is little knowledge about its incidence. The first study based on population-based cancer registry described spatial incidence and indicated that there existed at least county-level aggregation. The aim of the present work is to model the incidence patterns for the most incidence cancer in Córdoba Province, Argentina, using information from the Córdoba Cancer Registry by performing multilevel mixed model approach to deal with dependence and unobserved heterogeneity coming from the geo-reference cancer occurrence. Methods: Standardized incidence rates (world standard population) (SIR) by sex based on 5-year age groups were calculated for 109 districts nested on 26 counties for the most incidence cancers in Cordoba using 2004 database. A Poisson twolevel random effect model representing unobserved heterogeneity between first level-districts and second level-counties was fitted to assess the spatial distribution of the overall and site specific cancer incidence rates. Results: SIR cancer at Córdoba province shown an average of 263.53±138.34 and 200.45±98.30 for men and women, respectively. Considering the ratio site specific mean SIR to the total mean, breast cancer ratio was 0.25±0.19, prostate cancer ratio was 0.12±0.10 and lower values for lung and colon cancer for both sexes. The Poisson two-level random intercepts model fitted for SIR data distributed with overdispersion shown significant hierarchical structure for the cancer incidence distribution. Conclusions: a strong spatial-nested effect for the cancer incidence in Córdoba was observed and will help to begin the study of the factors associated with it.

A General Framework for the Analysis of Animal Resource Selection from Telemetry Data

We propose a general framework for the analysis of animal telemetry data through the use of weighted distributions. It is shown that several interpretations of resource selection functions arise when constructed from the ratio of a use and availability distribution. Through the proposed general framework, several popular resource selection models are shown to be special cases of the general model by making assumptions about animal movement and behavior. The weighted distribution framework is shown to be easily extended to readily account for telemetry data that are highly auto-correlated; as is typical with use of new technology such as global positioning systems animal relocations. An analysis of simulated data using several models constructed within the proposed framework is also presented to illustrate the possible gains from the flexible modeling framework. The proposed model is applied to a brown bear data set from southeast Alaska.

Distance software: design and analysis of distance sampling surveys for estimating population size

1. Distance sampling is a widely used technique for estimating the size or density of biological populations. Many distance sampling designs and most analyses use the software Distance. 2. We briefly review distance sampling and its assumptions, outline the history, structure and capabilities of Distance, and provide hints on its use. 3. Good survey design is a crucial prerequisite for obtaining reliable results. Distance has a survey design engine, with a built-in geographic information system, that allows properties of different proposed designs to be examined via simulation, and survey plans to be generated. 4. A first step in analysis of distance sampling data is modeling the probability of detection. Distance contains three increasingly sophisticated analysis engines for this: conventional distance sampling, which models detection probability as a function of distance from the transect and assumes all objects at zero distance are detected; multiple-covariate distance sampling, which allows covariates in addition to distance; and mark–recapture distance sampling, which relaxes the assumption of certain detection at zero distance. 5. All three engines allow estimation of density or abundance, stratified if required, with associated measures of precision calculated either analytically or via the bootstrap. 6. Advanced analysis topics covered include the use of multipliers to allow analysis of indirect surveys (such as dung or nest surveys), the density surface modeling analysis engine for spatial and habitat-modeling, and information about accessing the analysis engines directly from other software. 7. Synthesis and applications. Distance sampling is a key method for producing abundance and density estimates in challenging field conditions. The theory underlying the methods continues to expand to cope with realistic estimation situations. In step with theoretical developments, state-of- the-art software that implements these methods is described that makes the methods accessible to practicing ecologists.

Temporal and Spatial Variations of Ions, Isotopes and Agricultural Contaminants in Surface Waters and Groundwater of Nebraska's Rainwater Basin Wetland Region

The wetlands of south-central Nebraska’s Rainwater Basin region are considered of international importance as a habitat for millions of migratory birds, but are being endangered by agricultural practices. The Rainwater Basin extends across 17 counties and covers 4,000 square miles. The purpose of this study was to assemble baseline chemical data for several representative wetlands across the Rainwater Basin region, and determine the use of these chemical data for investigating groundwater recharge. Eight representative wetlands were chosen across the Rainwater Basin to monitor surface and groundwater chemistry. At each site, a shallow well and deep well were installed and sampled once in the summer of 2009 and again in the spring of 2010. Wetland surface water was sampled monthly from April, 2009 to May, 2010. Waters were analyzed for major ions, nutrients, pesticides and oxygen-18 and deuterium isotopes at the University of Nebraska Water Sciences Laboratory. Geochemical analysis of surface waters presents a range of temporal and spatial variations. Wetlands had variable water volumes, isotopic compositions, ion chemistries and agricultural contaminant levels throughout the year and, except for a few trends, theses variations cannot be predicted with certainty year-to-year or wetland-to-wetland. Isotopic compositions showed evaporation was a contributor to water loss, and thus, did impact water chemistry. Surface water nitrate concentrations ranged from <0.10 to 4.04 mg/L. The nitrate levels are much higher in the groundwater, ranging from <0.10 to 18.4 mg/L, and are of concern because they are found above the maximum contaminant level (MCL) of 10 mg/L. Atrazine concentrations in surface waters ranged from <0.05 to 10.3 ppb. Groundwater atrazine concentrations ranged from <0.05 to 0.28 ppb. The high atrazine concentrations in surface waters are of concern as they are above the MCL of 3 ppb, and the highest levels occur during the spring bird migration. Most sampled groundwaters had detectable tritium indicating a mix of modern (<5 to 10 years old) and submodern (older than 1950s) recharge. The groundwater also had differences in chemical and isotope composition, and in some cases, increased nitrate concentrations, between the two sampling periods. Modern groundwater tritium ages and changes in groundwater chemical and isotopic compositions may indicate connections with surface waters in the Rainwater Basin.

Contributions of demography and dispersal parameters to the spatial spread of a stage-structured insect invasion

Stage-structured models that integrate demography and dispersal can be used to identify points in the life cycle with large effects on rates of population spatial spread, information that is vital in the development of containment strategies for invasive species. Current challenges in the application of these tools include: (1) accounting for large uncertainty in model parameters, which may violate assumptions of ‘‘local’’ perturbation metrics such as sensitivities and elasticities, and (2) forecasting not only asymptotic rates of spatial spread, as is usually done, but also transient spatial dynamics in the early stages of invasion. We developed an invasion model for the Diaprepes root weevil (DRW; Diaprepes abbreviatus [Coleoptera: Curculionidae]), a generalist herbivore that has invaded citrus-growing regions of the United States. We synthesized data on DRW demography and dispersal and generated predictions for asymptotic and transient peak invasion speeds, accounting for parameter uncertainty. We quantified the contributions of each parameter toward invasion speed using a ‘‘global’’ perturbation analysis, and we contrasted parameter contributions during the transient and asymptotic phases. We found that the asymptotic invasion speed was 0.02–0.028 km/week, although the transient peak invasion speed (0.03– 0.045 km/week) was significantly greater. Both asymptotic and transient invasions speeds were most responsive to weevil dispersal distances. However, demographic parameters that had large effects on asymptotic speed (e.g., survival of early-instar larvae) had little effect on transient speed. Comparison of the global analysis with lower-level elasticities indicated that local perturbation analysis would have generated unreliable predictions for the responsiveness of invasion speed to underlying parameters. Observed range expansion in southern Florida (1992–2006) was significantly lower than the invasion speed predicted by the model. Possible causes of this mismatch include overestimation of dispersal distances, demographic rates, and spatiotemporal variation in parameter values. This study demonstrates that, when parameter uncertainty is large, as is often the case, global perturbation analyses are needed to identify which points in the life cycle should be targets of management. Our results also suggest that effective strategies for reducing spread during the asymptotic phase may have little effect during the transient phase. Includes Appendix.

Numerical analysis of different methods to calculate residual stresses in thin films based on instrumented indentation data

In this work, different methods to estimate the value of thin film residual stresses using instrumented indentation data were analyzed. This study considered procedures proposed in the literature, as well as a modification on one of these methods and a new approach based on the effect of residual stress on the value of hardness calculated via the Oliver and Pharr method. The analysis of these methods was centered on an axisymmetric two-dimensional finite element model, which was developed to simulate instrumented indentation testing of thin ceramic films deposited onto hard steel substrates. Simulations were conducted varying the level of film residual stress, film strain hardening exponent, film yield strength, and film Poisson's ratio. Different ratios of maximum penetration depth h(max) over film thickness t were also considered, including h/t = 0.04, for which the contribution of the substrate in the mechanical response of the system is not significant. Residual stresses were then calculated following the procedures mentioned above and compared with the values used as input in the numerical simulations. In general, results indicate the difference that each method provides with respect to the input values depends on the conditions studied. The method by Suresh and Giannakopoulos consistently overestimated the values when stresses were compressive. The method provided by Wang et al. has shown less dependence on h/t than the others.

Clusters of leprosy transmission and of late diagnosis in a highly endemic area in Brazil: focus on different spatial analysis approaches

Objective The Brazilian National Hansens Disease Control Program recently identified clusters with high disease transmission. Herein, we present different spatial analytical approaches to define highly vulnerable areas in one of these clusters. Method The study area included 373 municipalities in the four Brazilian states Maranha o, Para ', Tocantins and Piaui '. Spatial analysis was based on municipalities as the observation unit, considering the following disease indicators: (i) rate of new cases / 100 000 population, (ii) rate of cases < 15 years / 100 000 population, (iii) new cases with grade-2 disability / 100 000 population and (iv) proportion of new cases with grade-2 disabilities. We performed descriptive spatial analysis, local empirical Bayesian analysis and spatial scan statistic. Results A total of 254 (68.0%) municipalities were classified as hyperendemic (mean annual detection rates > 40 cases / 100 000 inhabitants). There was a concentration of municipalities with higher detection rates in Para ' and in the center of Maranha o. Spatial scan statistic identified 23 likely clusters of new leprosy case detection rates, most of them localized in these two states. These clusters included only 32% of the total population, but 55.4% of new leprosy cases. We also identified 16 significant clusters for the detection rate < 15 years and 11 likely clusters of new cases with grade-2. Several clusters of new cases with grade-2 / population overlap with those of new cases detection and detection of children < 15 years of age. The proportion of new cases with grade-2 did not reveal any significant clusters. Conclusions Several municipality clusters for high leprosy transmission and late diagnosis were identified in an endemic area using different statistical approaches. Spatial scan statistic is adequate to validate and confirm high-risk leprosy areas for transmission and late diagnosis, identified using descriptive spatial analysis and using local empirical Bayesian method. National and State leprosy control programs urgently need to intensify control actions in these highly vulnerable municipalities.

Bayesian analysis of weak gravitational lensing and Sunyaev-Zel'dovich data for six galaxy clusters

We present an analysis of observations made with the Arcminute Microkelvin Imager (AMI) and the CanadaFranceHawaii Telescope (CFHT) of six galaxy clusters in a redshift range of 0.160.41. The cluster gas is modelled using the SunyaevZeldovich (SZ) data provided by AMI, while the total mass is modelled using the lensing data from the CFHT. In this paper, we (i) find very good agreement between SZ measurements (assuming large-scale virialization and a gas-fraction prior) and lensing measurements of the total cluster masses out to r200; (ii) perform the first multiple-component weak-lensing analysis of A115; (iii) confirm the unusual separation between the gas and mass components in A1914 and (iv) jointly analyse the SZ and lensing data for the relaxed cluster A611, confirming our use of a simulation-derived masstemperature relation for parametrizing measurements of the SZ effect.

Long-term effects of oil pollution in mangrove forests (Baixada Santista, Southeast Brazil) detected using a GIS-based multitemporal analysis of aerial photographs

Oil spills are potential threats to the integrity of highly productive coastal wetlands, such as mangrove forests. In October 1983, a mangrove area of nearly 300 ha located on the southeastern coast of Brazil was impacted by a 3.5 million liter crude oil spill released by a broken pipeline. In order to assess the long-term effects of oil pollution on mangrove vegetation, we carried out a GIS-based multitemporal analysis of aerial photographs of the years 1962, 1994, 2000 and 2003. Photointerpretation, visual classification, class quantification, ground-truth and vegetation structure data were combined to evaluate the oil impact. Before the spill, the mangroves exhibited a homogeneous canopy and well-developed stands. More than ten years after the spill, the mangrove vegetation exhibited three distinct zones reflecting the long-term effects of the oil pollution. The most impacted zone (10.5 ha) presented dead trees, exposed substrate and recovering stands with reduced structural development. We suggest that the distinct impact and recovery zones reflect the spatial variability of oil removal rates in the mangrove forest. This study identifies the multitemporal analysis of aerial photographs as a useful tool for assessing a system's capacity for recovery and monitoring the long-term residual effects of pollutants on vegetation dynamics, thus giving support to mangrove forest management and conservation.

SPATIAL PATTERN OF FOUR FOREST COMMUNITIES BY RIPLEY'S K FUNCTION IN SAO PAULO STATE

Investigating tree's spatial patterns according to their size classes and according to their more abundant species can provide evidences about the structure of the vegetal community, since the spatial pattern is a key question for forestry ecology studies. The tree spatial organization patterns on the environment depend on several ecological processes and on the specific characteristics of each environment, so that the best understanding of this frame provides important elements for the knowledge on forestry formation. This paper aimed to study tree spatial patterns, according to the diameter classes and from four most abundant species in different forests, in order to provide evidences regarding to the ecology of each vegetal community. The spatial pattern description in each forestry formation was developed using Ripley's K function. The studied forestry formations were: Ombrophilous Forest, Cerradao, Seasonal Forest and Restinga Forest. In this work, a 10.24 ha plot was installed in each forestry formation, and every tree, with a circumference at breast height (CBH) larger than 15 cm were measured, georeferenced and identified. The obtained data highlights the aggregated character in tropical forests, as observed in every studied forest. The 'Cerraddo' and 'Restinga' forest trees showed close aggregate patterns. In the Ombrophilous forest, for all distance scales, the aggregate pattern was meaningful. In the seasonal forest, a random tendency was observed, although a meaningful aggregation was observed in all short distances. The spatial pattern by diameter classes was generally aggregated for trees smaller than 10 cm of diameter and between 10 and 20 cm and random for the others, proving the existence of a tendency which in young trees is more aggregated than in old ones. The spatial pattern of the dominant species is always strongly similar to the general pattern of each forestry formation. The differences between the spatial patterns of two or three coincident species, among the forestry formations, indicate that its pattern is influenced by each different environment. This stands out the importance of its self-ecology and of its ecological processes, intrinsic of each community that can explain the observed patterns.

Water balance and spatial distribution of an anuran community from Brazil

Terrestrial amphibians may dehydrate when exposed to low humidity, representing an important factor affecting spatial distribution and community composition. In this study we investigated whether rates of dehydration and rehydration are able to explain the spatial distribution of an anuran community in a Restinga environment at the northern coast of the State of Bahia, Brazil, represented by 11 species distributed in 27 sample units. The environmental data set containing 20 variables was reduced to a few synthetic axes by principal component analysis (PCA). Physiological variables measured were rates of dehydration, rehydration from water, and rehydration from a neutral substrate. Multiple regression analyses were used to test the null hypothesis of no association between the environmental data set (synthetic axes of PCA) and each axis representative of a physiological variable, which was rejected (P < 0.001). Of 15 possible partial regressions only rehydration rate from neutral substrate vs. PC1. and PC2, rehydration rate from water vs. PC1, and dehydration rate vs. PC2 were significant. Our analysis was influenced by a gradient between two different groups of sample units: a beach area with high density of bromeliads and an environment without bodies of water with low density of bromeliads. Species of very specific natural history and morphological characters occur in these environments: Phyllodytes melanomystax and Scinax auratus, species frequently occurring in terrestrial bromeliads, and Ischnocnema paulodutrai, common along the northern coast of Bahia and usually found in forest remnants within environments with low number of bodies of water. In dry environments species with lower rates of dehydration were dominant, whereas species showing greater rates of dehydration were found predominantly in microhabitats with greater moisture or abundance of bodies of water.

Hyperinsulinemic hypoglycemia of the infancy. Analysis of clinical data from a Brazilian sample

Objective: To review the presentation of hyperinsulinemic hypoglycemia of the infancy (HHI), its treatment and histology in Brazilian pediatric endocrinology sections. Materials and method: The protocol analyzed data of birth, laboratory results, treatment, surgery, and pancreas histology. Results: Twenty-five cases of HHI from six centers were analyzed: 15 male, 3/25 born by vaginal delivery. The average age at diagnosis was 10.3 days. Glucose and insulin levels in the critical sample showed an average of 24.7 mg/dL and 26.3 UI/dL. Intravenous infusion of the glucose was greater than 10 mg/kg/min in all cases (M:19,1). Diazoxide was used in 15/25 of the cases, octreotide in 10, glucocorticoid in 8, growth hormone in 3, nifedipine in 2 and glucagon in 1. Ten of the cases underwent pancreatectomy and histology results showed the diffuse form of disease. Conclusion: This is the first critic review of a Brazilian sample with congenital HHI. Arq Bras Endocrinol Metab. 2012; 56(9): 666-71

Semi-supervised dimensionality reduction based on partial least squares for visual analysis of high dimensional data

Dimensionality reduction is employed for visual data analysis as a way to obtaining reduced spaces for high dimensional data or to mapping data directly into 2D or 3D spaces. Although techniques have evolved to improve data segregation on reduced or visual spaces, they have limited capabilities for adjusting the results according to user's knowledge. In this paper, we propose a novel approach to handling both dimensionality reduction and visualization of high dimensional data, taking into account user's input. It employs Partial Least Squares (PLS), a statistical tool to perform retrieval of latent spaces focusing on the discriminability of the data. The method employs a training set for building a highly precise model that can then be applied to a much larger data set very effectively. The reduced data set can be exhibited using various existing visualization techniques. The training data is important to code user's knowledge into the loop. However, this work also devises a strategy for calculating PLS reduced spaces when no training data is available. The approach produces increasingly precise visual mappings as the user feeds back his or her knowledge and is capable of working with small and unbalanced training sets.

How far can we go in simplifying biomonitoring assessments? An integrated analysis of taxonomic surrogacy, taxonomic sufficiency and numerical resolution in a megadiverse region

The need for biodiversity conservation is increasing at a rate much faster than the acquisition of knowledge of biodiversity, such as descriptions of new species and mapping species distributions. As global changes are winning the race against the acquisition of knowledge, many researchers resort to the use of surrogate groups to aid in conservation decisions. Reductions in taxonomic and numerical resolution are also desirable, because they could allow more rapid the acquisition of knowledge while requiring less effort, if little important information is lost. In this study, we evaluated the congruence among 22 taxonomic groups sampled in a tropical forest in the Amazon basin. Our aim was to evaluate if any of these groups could be used as surrogates for the others in monitoring programs. We also evaluated if the taxonomic or numerical resolution of possible surrogates could be reduced without greatly reducing the overall congruence. Congruence among plant groups was high, whereas the congruence among most animal groups was very low, except for anurans in which congruence values were only slightly lower than for plants. Liana (Bignoniaceae) was the group with highest congruence, even using genera presence-absence data. The congruence among groups was related to environmental factors, specifically the clay and phosphorous contents of soil. Several groups showed strong spatial clumping, but this was unrelated to the congruence among groups. The high degree of congruence of lianas with the other groups suggests that it may be a reasonable surrogate group, mainly for the other plant groups analyzed, if soil data are not available. Although lianas are difficult to count and identify, the number of studies on the ecology of lianas is increasing. Most of these studies have concluded that lianas are increasing in abundance in tropical forests. In addition to the high congruence, lianas are worth monitoring in their own right because they are sensitive to global warming and the increasing frequency and severity of droughts in tropical regions. Our findings suggest that the use of data on surrogate groups with relatively low taxonomic and numerical resolutions can be a reliable shortcut for biodiversity assessments, especially in megadiverse areas with high rates of habitat conversion, where the lack of biodiversity knowledge is pervasive. (c) 2012 Elsevier Ltd. All rights reserved.

BayGO: Bayesian analysis of ontology term enrichment in microarray data

Abstract Background The search for enriched (aka over-represented or enhanced) ontology terms in a list of genes obtained from microarray experiments is becoming a standard procedure for a system-level analysis. This procedure tries to summarize the information focussing on classification designs such as Gene Ontology, KEGG pathways, and so on, instead of focussing on individual genes. Although it is well known in statistics that association and significance are distinct concepts, only the former approach has been used to deal with the ontology term enrichment problem. Results BayGO implements a Bayesian approach to search for enriched terms from microarray data. The R source-code is freely available at http://blasto.iq.usp.br/~tkoide/BayGO in three versions: Linux, which can be easily incorporated into pre-existent pipelines; Windows, to be controlled interactively; and as a web-tool. The software was validated using a bacterial heat shock response dataset, since this stress triggers known system-level responses. Conclusion The Bayesian model accounts for the fact that, eventually, not all the genes from a given category are observable in microarray data due to low intensity signal, quality filters, genes that were not spotted and so on. Moreover, BayGO allows one to measure the statistical association between generic ontology terms and differential expression, instead of working only with the common significance analysis.