Generalized linear and generalized additive models in studies of species distributions: setting the scene

An important statistical development of the last 30 years has been the advance in regression analysis provided by generalized linear models (GLMs) and generalized additive models (GAMs). Here we introduce a series of papers prepared within the framework of an international workshop entitled: Advances in GLMs/GAMs modeling: from species distribution to environmental management, held in Riederalp, Switzerland, 6-11 August 2001.We first discuss some general uses of statistical models in ecology, as well as provide a short review of several key examples of the use of GLMs and GAMs in ecological modeling efforts. We next present an overview of GLMs and GAMs, and discuss some of their related statistics used for predictor selection, model diagnostics, and evaluation. Included is a discussion of several new approaches applicable to GLMs and GAMs, such as ridge regression, an alternative to stepwise selection of predictors, and methods for the identification of interactions by a combined use of regression trees and several other approaches. We close with an overview of the papers and how we feel they advance our understanding of their application to ecological modeling.

Bioclimatic constraints to Andean cat distribution: a modelling application for rare species

AimTo identify the bioclimatic niche of the endangered Andean cat (Leopardus jacobita), one of the rarest and least known felids in the world, by developing a species distribution model.LocationSouth America, High Andes and Patagonian steppe. Peru, Bolivia, Chile, Argentina.MethodsWe used 108 Andean cat records to build the models, and 27 to test them, applying the Maxent algorithm to sets of uncorrelated bioclimatic variables from global databases, including elevation. We based our biogeographical interpretations on the examination of the predicted geographic range, the modelled response curves and latitudinal variations in climatic variables associated with the locality data.ResultsSimple bioclimatic models for Andean cats were highly predictive with only 3-4 explanatory variables. The climatic niche of the species was defined by extreme diurnal variations in temperature, cold minimum and moderate maximum temperatures, and aridity, characteristic not only of the Andean highlands but also of the Patagonian steppe. Argentina had the highest representation of suitable climates, and Chile the lowest. The most favourable conditions were centrally located and spanned across international boundaries. Discontinuities in suitable climatic conditions coincided with three biogeographical barriers associated with climatic or topographic transitions.Main conclusionsSimple bioclimatic models can produce useful predictions of suitable climatic conditions for rare species, including major biogeographical constraints. In our study case, these constraints are also known to affect the distribution of other Andean species and the genetic structure of Andean cat populations. We recommend surveys of areas with suitable climates and no Andean cat records, including the corridor connecting two core populations. The inclusion of landscape variables at finer scales, crucially the distribution of Andean cat prey, would contribute to refine our predictions for conservation applications.

The distribution of tuberculosis in Porto Alegre: analysis of the magnitude and tuberculosis-HIV coinfection

Objective Analyzing the geographical distribution of the tuberculosis (TB), its incidence and prevalence and TB-HIV coinfection in the districts of Porto Alegre from 2007 to 2011. Method An ecological, descriptive study of time series that used descriptive and geoprocessing techniques. Results In total, were recorded 3,369 incident cases and 3,998 prevalent cases of pulmonary TB. In both contexts, there was predominance of cases in males and in Caucasians. Seventeen districts showed prevalence rates above 79.2 cases/100,000 inhabitants, considering that 15 of them had incidence rates above 73.7 cases/100,000 inhabitants. The TB-HIV coinfection rates reached 67% in some districts, which is above the city average value (30%). Conclusion The distribution analysis showed that the reformulation and restructuring of policies and health services in Porto Alegre are essential.


Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Nonlinear models and small sample performance of the generalized method of moments

In this paper I explore the issue of nonlinearity (both in the datageneration process and in the functional form that establishes therelationship between the parameters and the data) regarding the poorperformance of the Generalized Method of Moments (GMM) in small samples.To this purpose I build a sequence of models starting with a simple linearmodel and enlarging it progressively until I approximate a standard (nonlinear)neoclassical growth model. I then use simulation techniques to find the smallsample distribution of the GMM estimators in each of the models.

On information and competition in private value auctions

How much information does an auctioneer want bidders to have in a private value environment?We address this question using a novel approach to ordering information structures based on the property that in private value settings more information leads to a more disperse distribution of buyers updated expected valuations. We define the class of precision criteria following this approach and different notions of dispersion, and relate them to existing criteria of informativeness. Using supermodular precision, we obtain three results: (1) a more precise information structure yields a more efficient allocation; (2) the auctioneer provides less than the efficient level of information since more information increases bidder informational rents; (3) there is a strategic complementarity between information and competition, so that both the socially efficient and the auctioneer s optimal choice of precision increase with the number of bidders, and both converge as the number of bidders goes to infinity.

Constant fluctuating asymmetry but not directional asymmetry along the geographic distribution of Drosophila antonietae (Diptera, Drosophilidae)

ABSTRACT The population dynamics of a species tends to change from the core to the periphery of its distribution. Therefore, one could expect peripheral populations to be subject to a higher level of stress than more central populations (the center–periphery hypothesis) and consequently should present a higher level of fluctuating asymmetry. To test these predictions we study asymmetry in wing shape of five populations of Drosophila antonietae collected throughout the distribution of the species using fluctuating asymmetry as a proxy for developmental instability. More specifically, we addressed the following questions: (1) what types of asymmetry occur in populations of D. antonietae? (2) Does the level of fluctuating asymmetry vary among populations? (3) Does peripheral populations have a higher fluctuating asymmetry level than central populations? We used 12 anatomical landmarks to quantify patterns of asymmetry in wing shape in five populations of D. antonietae within the framework of geometric morphometrics. Net asymmetry – a composite measure of directional asymmetry + fluctuating asymmetry – varied significantly among populations. However, once net asymmetry of each population is decomposed into directional asymmetry and fluctuating asymmetry, most of the variation in asymmetry was explained by directional asymmetry alone, suggesting that populations of D. antonietae have the same magnitude of fluctuating asymmetry throughout the geographical distribution of the species. We hypothesize that larval development in rotting cladodes might play an important role in explaining our results. In addition, our study underscores the importance of understanding the interplay between the biology of a species and its geographical patterns of asymmetry.

Timeline and geographical distribution of Helicoverpa armigera (Hübner) (Lepidoptera, Noctuidae: Heliothinae) in Brazil

ABSTRACT This study presents registers of Helicoverpa armigera (Hübner) occurrence to assess its spatial and temporal distribution in Brazil. We used data from collections, especially from the Southern Region, systematic collections in Rio Grande do Sul, occasional collections of caterpillars and adults in different regions of Brazil, as well as literature registers. We conclude that the introduction of H. armigera in Brazil probably occurred before October 2008. We also register that in August 2012 H. armigera was already present from the extreme southern part (Rio Grande do Sul) to the extreme northern part (Amapá) of Brazil.

Comparison of the glomerular filtration rate in children by the new revised Schwartz formula and a new generalized formula.

The most widely used formula for estimating glomerular filtration rate (eGFR) in children is the Schwartz formula. It was revised in 2009 using iohexol clearances with measured GFR (mGFR) ranging between 15 and 75 ml/min × 1.73 m(2). Here we assessed the accuracy of the Schwartz formula using the inulin clearance (iGFR) method to evaluate its accuracy for children with less renal impairment comparing 551 iGFRs of 392 children with their Schwartz eGFRs. Serum creatinine was measured using the compensated Jaffe method. In order to find the best relationship between iGFR and eGFR, a linear quadratic regression model was fitted and a more accurate formula was derived. This quadratic formula was: 0.68 × (Height (cm)/serum creatinine (mg/dl))-0.0008 × (height (cm)/serum creatinine (mg/dl))(2)+0.48 × age (years)-(21.53 in males or 25.68 in females). This formula was validated using a split-half cross-validation technique and also externally validated with a new cohort of 127 children. Results show that the Schwartz formula is accurate until a height (Ht)/serum creatinine value of 251, corresponding to an iGFR of 103 ml/min × 1.73 m(2), but significantly unreliable for higher values. For an accuracy of 20 percent, the quadratic formula was significantly better than the Schwartz formula for all patients and for patients with a Ht/serum creatinine of 251 or greater. Thus, the new quadratic formula could replace the revised Schwartz formula, which is accurate for children with moderate renal failure but not for those with less renal impairment or hyperfiltration.

Integrated Analysis of High-Resolution Gene Expression and Copy Number Profiling Identified Biallelic Deletion of CDKN2A/2B Tumor Suppressor Locus As the Most Frequent and Unique Genomic Abnormality in Diffuse Large B-Cell Lymphoma (DLBCL) with Strong Prognostic Value in Both GCB and ABC Subtypes and Not Overcome by a Dose-Intensive Immunochemotherapy Regimen Plus Rituximab. Results of a Prospective GELA Clinical Trial Program

Background and aim of the study: Genomic gains and losses play a crucial role in the development and progression of DLBCL and are closely related to gene expression profiles (GEP), including the germinal center B-cell like (GCB) and activated B-cell like (ABC) cell of origin (COO) molecular signatures. To identify new oncogenes or tumor suppressor genes (TSG) involved in DLBCL pathogenesis and to determine their prognostic values, an integrated analysis of high-resolution gene expression and copy number profiling was performed. Patients and methods: Two hundred and eight adult patients with de novo CD20+ DLBCL enrolled in the prospective multicentric randomized LNH-03 GELA trials (LNH03-1B, -2B, -3B, 39B, -5B, -6B, -7B) with available frozen tumour samples, centralized reviewing and adequate DNA/RNA quality were selected. 116 patients were treated by Rituximab(R)-CHOP/R-miniCHOP and 92 patients were treated by the high dose (R)-ACVBP regimen dedicated to patients younger than 60 years (y) in frontline. Tumour samples were simultaneously analysed by high resolution comparative genomic hybridization (CGH, Agilent, 144K) and gene expression arrays (Affymetrix, U133+2). Minimal common regions (MCR), as defined by segments that affect the same chromosomal region in different cases, were delineated. Gene expression and MCR data sets were merged using Gene expression and dosage integrator algorithm (GEDI, Lenz et al. PNAS 2008) to identify new potential driver genes. Results: A total of 1363 recurrent (defined by a penetrance > 5%) MCRs within the DLBCL data set, ranging in size from 386 bp, affecting a single gene, to more than 24 Mb were identified by CGH. Of these MCRs, 756 (55%) showed a significant association with gene expression: 396 (59%) gains, 354 (52%) single-copy deletions, and 6 (67%) homozygous deletions. By this integrated approach, in addition to previously reported genes (CDKN2A/2B, PTEN, DLEU2, TNFAIP3, B2M, CD58, TNFRSF14, FOXP1, REL...), several genes targeted by gene copy abnormalities with a dosage effect and potential physiopathological impact were identified, including genes with TSG activity involved in cell cycle (HACE1, CDKN2C) immune response (CD68, CD177, CD70, TNFSF9, IRAK2), DNA integrity (XRCC2, BRCA1, NCOR1, NF1, FHIT) or oncogenic functions (CD79b, PTPRT, MALT1, AUTS2, MCL1, PTTG1...) with distinct distribution according to COO signature. The CDKN2A/2B tumor suppressor locus (9p21) was deleted homozygously in 27% of cases and hemizygously in 9% of cases. Biallelic loss was observed in 49% of ABC DLBCL and in 10% of GCB DLBCL. This deletion was strongly correlated to age and associated to a limited number of additional genetic abnormalities including trisomy 3, 18 and short gains/losses of Chr. 1, 2, 19 regions (FDR < 0.01), allowing to identify genes that may have synergistic effects with CDKN2A/2B inactivation. With a median follow-up of 42.9 months, only CDKN2A/2B biallelic deletion strongly correlates (FDR p.value < 0.01) to a poor outcome in the entire cohort (4y PFS = 44% [32-61] respectively vs. 74% [66-82] for patients in germline configuration; 4y OS = 53% [39-72] vs 83% [76-90]). In a Cox proportional hazard prediction of the PFS, CDKN2A/2B deletion remains predictive (HR = 1.9 [1.1-3.2], p = 0.02) when combined with IPI (HR = 2.4 [1.4-4.1], p = 0.001) and GCB status (HR = 1.3 [0.8-2.3], p = 0.31). This difference remains predictive in the subgroup of patients treated by R-CHOP (4y PFS = 43% [29-63] vs. 66% [55-78], p=0.02), in patients treated by R-ACVBP (4y PFS = 49% [28-84] vs. 83% [74-92], p=0.003), and in GCB (4y PFS = 50% [27-93] vs. 81% [73-90], p=0.02), or ABC/unclassified (5y PFS = 42% [28-61] vs. 67% [55-82] p = 0.009) molecular subtypes (Figure 1). Conclusion: We report for the first time an integrated genetic analysis of a large cohort of DLBCL patients included in a prospective multicentric clinical trial program allowing identifying new potential driver genes with pathogenic impact. However CDKN2A/2B deletion constitutes the strongest and unique prognostic factor of chemoresistance to R-CHOP, regardless the COO signature, which is not overcome by a more intensified immunochemotherapy. Patients displaying this frequent genomic abnormality warrant new and dedicated therapeutic approaches.

The induced 2-tuple linguistic generalized OWA operator and its application in linguistic decision making

[spa] Se presenta el operador de media ponderada ordenada generalizada lingüística de 2 tuplas inducida (2-TILGOWA). Es un nuevo operador de agregación que extiende los anteriores modelos a través de utilizar medias generalizadas, variables de ordenación inducidas e información lingüística representada mediante el modelo de las 2 tuplas lingüísticas. Su principal ventaja se encuentra en la posibilidad de incluir a un gran número de operadores de agregación lingüísticos como casos particulares. Por eso, el análisis puede ser visto desde diferentes perspectivas de forma que se obtiene una visión más completa del problema considerado y seleccionar la alternativa que parece estar en mayor concordancia con nuestros intereses o creencias. A continuación se desarrolla una generalización mayor a través de utilizar medias cuasi-aritméticas, obteniéndose el operador Quasi-2-TILOWA. El trabajo finaliza analizando la aplicabilidad del nuevo modelo en un problema de toma de decisiones sobre gestión de la producción.

Time of ruin in a risk model with generalized Erlang (n) interclaim times and a constant dividend barrier

In this paper we analyze the time of ruin in a risk process with the interclaim times being Erlang(n) distributed and a constant dividend barrier. We obtain an integro-differential equation for the Laplace Transform of the time of ruin. Explicit solutions for the moments of the time of ruin are presented when the individual claim amounts have a distribution with rational Laplace transform. Finally, some numerical results and a compare son with the classical risk model, with interclaim times following an exponential distribution, are given.

The induced 2-tuple linguistic generalized OWA operator and its application in linguistic decision making

[spa] Se presenta el operador de media ponderada ordenada generalizada lingüística de 2 tuplas inducida (2-TILGOWA). Es un nuevo operador de agregación que extiende los anteriores modelos a través de utilizar medias generalizadas, variables de ordenación inducidas e información lingüística representada mediante el modelo de las 2 tuplas lingüísticas. Su principal ventaja se encuentra en la posibilidad de incluir a un gran número de operadores de agregación lingüísticos como casos particulares. Por eso, el análisis puede ser visto desde diferentes perspectivas de forma que se obtiene una visión más completa del problema considerado y seleccionar la alternativa que parece estar en mayor concordancia con nuestros intereses o creencias. A continuación se desarrolla una generalización mayor a través de utilizar medias cuasi-aritméticas, obteniéndose el operador Quasi-2-TILOWA. El trabajo finaliza analizando la aplicabilidad del nuevo modelo en un problema de toma de decisiones sobre gestión de la producción.

Chemoendocrine compared with endocrine adjuvant therapies for node-negative breast cancer: predictive value of centrally reviewed expression of estrogen and progesterone receptors--International Breast Cancer Study Group.

PURPOSE: To centrally assess estrogen receptor (ER) and progesterone receptor (PgR) levels by immunohistochemistry and investigate their predictive value for benefit of chemo-endocrine compared with endocrine adjuvant therapy alone in two randomized clinical trials for node-negative breast cancer. PATIENTS AND METHODS: International Breast Cancer Study Group Trial VIII compared cyclophosphamide, methotrexate, and fluorouracil (CMF) chemotherapy for 6 cycles followed by endocrine therapy with goserelin with either modality alone in pre- and perimenopausal patients. Trial IX compared three cycles of CMF followed by tamoxifen for 5 years versus tamoxifen alone in postmenopausal patients. Central Pathology Office reviewed 883 (83%) of 1,063 patients on Trial VIII and 1,365 (82%) of 1,669 on Trial IX and determined ER and PgR by immunohistochemistry. Disease-free survival (DFS) was compared across the spectrum of expression of each receptor using the Subpopulation Treatment Effect Pattern Plot methodology. RESULTS: Both receptors displayed a bimodal distribution, with substantial proportions showing no staining (receptor absent) and most of the remainder showing a high percentage of stained cells. Chemo-endocrine therapy yielded DFS superior to endocrine therapy alone for patients with receptor-absent tumors, and in some cases also for those with low levels of receptor expression. Among patients with ER-expressing tumors, additional prediction of benefit was suggested in absent or low PgR in Trial VIII but not in Trial IX. CONCLUSION: Low levels of ER and PgR are predictive of the benefit of adding chemotherapy to endocrine therapy. Low PgR may add further prediction among pre- and perimenopausal but not postmenopausal patients whose tumors express ER.

Physiological traits affecting the distribution and wintering strategy of the bat Tadarida teniotis

The ability to enter torpor at low ambient temperature, which enables insectivorous bats to survive seasonal food shortage, is often seen as a prerequisite for colonizing cold environments. Free-tailed bats (Molossidae) show a distribution with a maximum latitudinal extension that appears to be intermediate between truly tropical and temperate-zone bat families. We therefore tested the hypothesis that Tadarida teniotis, the molossid species reaching the highest latitude worldwide (46 degrees N), lacks the extreme physiological adaptations to cold that enable other sympatric bats to enter further into the temperate zone. We studied the metabolism of individuals subjected to various ambient temperatures in the laboratory by respirometry, and we monitored the body temperature of free-ranging individuals in winter and early spring in the Swiss Alps using temperature-sensitive radio-tags. For comparison, metabolic data were obtained from Nyctalus noctula, a typically hibernating vespertilionid bat of similar body size and convergent foraging tactics. The metabolic data support the hypothesis that T. teniotis cannot experience such low ambient temperatures as sympatric temperate-zone vespertilionid bats without incurring much higher energetic costs for thermogenesis. The minimum rate of metabolism in torpor was obtained at 7.5 degrees-10 degrees C in T. teniotis, as compared to 2.5 degrees-5 degrees C in N. noctula. Field data showed that T. teniotis behaves as a classic thermo-conforming hibernator in the Alps, with torpor bouts lasting up to 8 d. This contradicts the widely accepted opinion that Molossidae are nonhibernating bars. However, average body temperature (10 degrees-13 degrees C) and mean arousal frequency (3.4 d in one bat in January) appear to be markedly higher than in other temperate-zone bat species. At the northern border of its range T. teniotis selects relatively warm roosts (crevices in tall, south-exposed limestone cliffs) in winter where temperatures oscillate around 10 degrees C. By this means, T. teniotis apparently avoids the risk of prolonged exposure to energetically critical ambient temperatures in torpor (<6.5 degrees-7.5 degrees C) during cold spells. Possibly shared by other Molossidae, the physiological pattern observed in T. teniotis may clearly be linked to the intermediate latitudinal extension of this bat family.