Model Evaluation Based on the Distribution of Estimated Absolute Prediction Error

The construction of a reliable, practically useful prediction rule for future response is heavily dependent on the "adequacy" of the fitted regression model. In this article, we consider the absolute prediction error, the expected value of the absolute difference between the future and predicted responses, as the model evaluation criterion. This prediction error is easier to interpret than the average squared error and is equivalent to the mis-classification error for the binary outcome. We show that the distributions of the apparent error and its cross-validation counterparts are approximately normal even under a misspecified fitted model. When the prediction rule is "unsmooth", the variance of the above normal distribution can be estimated well via a perturbation-resampling method. We also show how to approximate the distribution of the difference of the estimated prediction errors from two competing models. With two real examples, we demonstrate that the resulting interval estimates for prediction errors provide much more information about model adequacy than the point estimates alone.

A Diagnostic Test for the Mixing Distribution in a Generalised Linear Mixed Model

We introduce a diagnostic test for the mixing distribution in a generalised linear mixed model. The test is based on the difference between the marginal maximum likelihood and conditional maximum likelihood estimates of a subset of the fixed effects in the model. We derive the asymptotic variance of this difference, and propose a test statistic that has a limiting chi-square distribution under the null hypothesis that the mixing distribution is correctly specified. For the important special case of the logistic regression model with random intercepts, we evaluate via simulation the power of the test in finite samples under several alternative distributional forms for the mixing distribution. We illustrate the method by applying it to data from a clinical trial investigating the effects of hormonal contraceptives in women.

Distribution and genetic variability among Campylobacter spp. isolates from different animal species and humans in Switzerland

In Switzerland, a national database with 1028 Campylobacter isolates from poultry, pigs, cats, dogs, cattle, humans, zoo animals and water has been created. The database contains the genetic fingerprint and background information of each Campylobacter isolate. Dominant species could be identified in the different sources with a majority of Campylobacter jejuni in poultry (73%), humans (79%), cattle (95%), zoo animals (40%) and water (100%), of Campylobacter coli in pigs (72%), and of Campylobacter upsaliensis/helveticus in cats and dogs (55%). The comparison of three genotyping methods, amplified fragment length polymorphism (AFLP), pulsed field gel electrophoresis and restriction fragment length polymorphism, revealed that AFLP allows discrimination between the different Campylobacter species and is the most appropriate method to distinguish specific strains within the same species. Genotyping analysis demonstrated that the Campylobacter population is heterogeneous among the different sources and that no dominant clone is spread in the country. Genotyping and the resulting database are useful tools to trace back future Campylobacter infections.

Evolutionary demography of long-lived monocarpic perennials: a time-lagged integral projection model

1. The evolution of flowering strategies (when and at what size to flower) in monocarpic perennials is determined by balancing current reproduction with expected future reproduction, and these are largely determined by size-specific patterns of growth and survival. However, because of the difficulty in following long-lived individuals throughout their lives, this theory has largely been tested using short-lived species (< 5 years). 2. Here, we tested this theory using the long-lived monocarpic perennial Campanula thyrsoides which can live up to 16 years. We used a novel approach that combined permanent plot and herb chronology data from a 3-year field study to parameterize and validate integral projection models (IPMs). 3. Similar to other monocarpic species, the rosette leaves of C. thyrsoides wither over winter and so size cannot be measured in the year of flowering. We therefore extended the existing IPM framework to incorporate an additional time delay that arises because flowering demography must be predicted from rosette size in the year before flowering. 4. We found that all main demographic functions (growth, survival probability, flowering probability and fecundity) were strongly size-dependent and there was a pronounced threshold size of flowering. There was good agreement between the predicted distribution of flowering ages obtained from the IPMs and that estimated in the field. Mostly, there was good agreement between the IPM predictions and the direct quantitative field measurements regarding the demographic parameters lambda, R-0 and T. We therefore conclude that the model captures the main demographic features of the field populations. 5. Elasticity analysis indicated that changes in the survival and growth function had the largest effect (c. 80%) on lambda and this was considerably larger than in short-lived monocarps. We found only weak selection pressure operating on the observed flowering strategy which was close to the predicted evolutionary stable strategy. 6. Synthesis. The extended IPM accurately described the demography of a long-lived monocarpic perennial using data collected over a relatively short period. We could show that the evolution of flowering strategies in short- and long-lived monocarps seem to follow the same general rules but with a longevity-related emphasis on survival over fecundity.

Phylogenetic relationship of equine Actinobacillus species and distribution of RTX toxin genes among clusters

Equine Actinobacillus species were analysed phylogenetically by 16S rRNA gene (rrs) sequencing focusing on the species Actinobacillus equuli, which has recently been subdivided into the non-haemolytic A. equuli subsp. equuli and the haemolytic A. equuli subsp. haemolyticus. In parallel we determined the profile for RTX toxin genes of the sample of strains by PCR testing for the presence of the A. equuli haemolysin gene aqx, and the toxin genes apxI, apxII, apxIII and apxIV, which are known in porcine pathogens such as Actinobacillus pleuropneumoniae and Actinobacillus suis. The rrs-based phylogenetic analysis revealed two distinct subclusters containing both A. equuli subsp. equuli and A. equuli subsp. haemolyticus distributed through both subclusters with no correlation to taxonomic classification. Within one of the rrs-based subclusters containing the A. equuli subsp. equuli type strain, clustered as well the porcine Actinobacillus suis strains. This latter is known to be also phenotypically closely related to A. equuli. The toxin gene analysis revealed that all A. equuli subsp. haemolyticus strains from both rrs subclusters specifically contained the aqx gene while the A. suis strains harboured the genes apxI and apxII. The aqx gene was found to be specific for A. equuli subsp. haemolyticus, since A. equuli subsp. equuli contained no aqx nor any of the other RTX genes tested. The specificity of aqx for the haemolytic equine A. equuli and ApxI and ApxII for the porcine A. suis indicates a role of these RTX toxins in host species predilection of the two closely related species of bacterial pathogens and allows PCR based diagnostic differentiation of the two.

Effects of two traits of the ecological state equation on our understanding of species coexistence and ecosystem services

Species coexistence has been a fundamental issue to understand ecosystem functioning since the beginnings of ecology as a science. The search of a reliable and all-encompassing explanation for this issue has become a complex goal with several apparently opposing trends. On the other side, seemingly unconnected with species coexistence, an ecological state equation based on the inverse correlation between an indicator of dispersal that fits gamma distribution and species diversity has been recently developed. This article explores two factors, whose effects are inconspicuous in such an equation at the first sight, that are used to develop an alternative general theoretical background in order to provide a better understanding of species coexistence. Our main outcomes are: (i) the fit of dispersal and diversity values to gamma distribution is an important factor that promotes species coexistence mainly due to the right-skewed character of gamma distribution; (ii) the opposite correlation between species diversity and dispersal implies that any increase of diversity is equivalent to a route of “ecological cooling” whose maximum limit should be constrained by the influence of the third law of thermodynamics; this is in agreement with the well-known asymptotic trend of diversity values in space and time; (iii) there are plausible empirical and theoretical ways to apply physical principles to explain important ecological processes; (iv) the gap between theoretical and empirical ecology in those cases where species diversity is paradoxically high could be narrowed by a wave model of species coexistence based on the concurrency of local equilibrium states. In such a model, competitive exclusion has a limited but indispensable role in harmonious coexistence with functional redundancy. We analyze several literature references as well as ecological and evolutionary examples that support our approach, reinforcing the meaning equivalence between important physical and ecological principles.

Spatial Distribution of Cryptic Species Diversity in European Freshwater Amphipods (Gammarus fossarum) as Revealed by Pyrosequencing

In order to understand and protect ecosystems, local gene pools need to be evaluated with respect to their uniqueness. Cryptic species present a challenge in this context because their presence, if unrecognized, may lead to serious misjudgement of the distribution of evolutionarily distinct genetic entities. In this study, we describe the current geographical distribution of cryptic species of the ecologically important stream amphipod Gammarus fossarum (types A, B and C). We use a novel pyrosequencing assay for molecular species identification and survey 62 populations in Switzerland, plus several populations in Germany and eastern France. In addition, we compile data from previous publications (mainly Germany). A clear transition is observed from type A in the east (Danube and Po drainages) to types B and, more rarely, C in the west (Meuse, Rhone, and four smaller French river systems). Within the Rhine drainage, the cryptic species meet in a contact zone which spans the entire G. fossarum distribution range from north to south. This large-scale geographical sorting indicates that types A and B persisted in separate refugia during Pleistocene glaciations. Within the contact zone, the species rarely co-occur at the same site, suggesting that ecological processes may preclude long-term coexistence. The clear phylogeographical signal observed in this study implies that, in many parts of Europe, only one of the cryptic species is present.

Are cryptic host species also cryptic to parasites? Host specificity and geographical distribution of acanthocephalan parasites infecting freshwater Gammarus

Many parasites infect multiple host species. In coevolving host–parasite interactions, theory predicts that parasites should be adapted to locally common hosts, which could lead to regional shifts in host preferences. We studied the interaction between freshwater Gammarus (Crustacea, Amphipoda) and their acanthocephalan parasites using a large-scale field survey and experiments, combined with molecular identification of cryptic host and parasite species. Gammarus pulex is a common host for multiple species of Acanthocephala in Europe but, in Switzerland, is less common than two cryptic members of the Gammarus fossarum species complex (type A and type B). We found that natural populations of these cryptic species were frequently infected by Pomphorhynchus tereticollis and Polymorphus minutus. Four additional parasite species occurred only locally. Parasites were more common in G. fossarum type B than in type A. Infection experiments using several host and parasite sources confirmed consistently lower infection rates in G. pulex than in G. fossarum type A, suggesting a general difference in susceptibility between the two species. In conclusion, we could show that cryptic host species differ in their interactions with parasites, but that these differences were much less dramatic than differences between G. fossarum (type A) and G. pulex. Our data suggest that the acanthocephalans in Switzerland have adapted to the two most common Gammarus species in this region where host species frequencies differ from near-by regions in Europe.

Model-Based Estimates of Right Whale Habitat Use in the Gulf of Maine

Balancing human uses of the marine environment with the recovery of protected species requires accurate information on when and where species of interest are likely to be present. Here, we describe a system that can produce useful estimates of right whale Eubalaena glacialis presence and abundance on their feeding grounds in the Gulf of Maine. The foundation of our system is a coupled physical-biological model of the copepod Calan us finmarchicus, the preferred prey of right whales. From the modeled prey densities, we can estimate when whales will appear in the Great South Channel feeding ground. Based on our experience with the system, we consider how the relationship between right whales and copepods changes across spatial scales. The scale-dependent relationship between whales and copepods provides insight into how to improve future estimates of the distribution of right whales and other pelagic predators.

Temperature Prediction Model for Bone Drilling Based on Density Distribution and In Vivo Experiments for Minimally Invasive Robotic Cochlear Implantation

Surgical robots have been proposed ex vivo to drill precise holes in the temporal bone for minimally invasive cochlear implantation. The main risk of the procedure is damage of the facial nerve due to mechanical interaction or due to temperature elevation during the drilling process. To evaluate the thermal risk of the drilling process, a simplified model is proposed which aims to enable an assessment of risk posed to the facial nerve for a given set of constant process parameters for different mastoid bone densities. The model uses the bone density distribution along the drilling trajectory in the mastoid bone to calculate a time dependent heat production function at the tip of the drill bit. Using a time dependent moving point source Green's function, the heat equation can be solved at a certain point in space so that the resulting temperatures can be calculated over time. The model was calibrated and initially verified with in vivo temperature data. The data was collected in minimally invasive robotic drilling of 12 holes in four different sheep. The sheep were anesthetized and the temperature elevations were measured with a thermocouple which was inserted in a previously drilled hole next to the planned drilling trajectory. Bone density distributions were extracted from pre-operative CT data by averaging Hounsfield values over the drill bit diameter. Post-operative [Formula: see text]CT data was used to verify the drilling accuracy of the trajectories. The comparison of measured and calculated temperatures shows a very good match for both heating and cooling phases. The average prediction error of the maximum temperature was less than 0.7 °C and the average root mean square error was approximately 0.5 °C. To analyze potential thermal damage, the model was used to calculate temperature profiles and cumulative equivalent minutes at 43 °C at a minimal distance to the facial nerve. For the selected drilling parameters, temperature elevation profiles and cumulative equivalent minutes suggest that thermal elevation of this minimally invasive cochlear implantation surgery may pose a risk to the facial nerve, especially in sclerotic or high density mastoid bones. Optimized drilling parameters need to be evaluated and the model could be used for future risk evaluation.

Past and future evolution of Abies alba forests in Europe - comparison of a dynamic vegetation model with palaeo data and observations

Information on how species distributions and ecosystem services are impacted by anthropogenic climate change is important for adaptation planning. Palaeo data suggest that Abies alba formed forests under significantly warmer-than-present conditions in Europe and might be a native substitute for widespread drought-sensitive temperate and boreal tree species such as beech (Fagus sylvatica) and spruce (Picea abies) under future global warming conditions. Here, we combine pollen and macrofossil data, modern observations, and results from transient simulations with the LPX-Bern dynamic global vegetation model to assess past and future distributions of A. alba in Europe. LPX-Bern is forced with climate anomalies from a run over the past 21 000 years with the Community Earth System Model, modern climatology, and with 21st-century multimodel ensemble results for the high-emission RCP8.5 and the stringent mitigation RCP2.6 pathway. The simulated distribution for present climate encompasses the modern range of A. alba, with the model exceeding the present distribution in north-western and southern Europe. Mid-Holocene pollen data and model results agree for southern Europe, suggesting that at present, human impacts suppress the distribution in southern Europe. Pollen and model results both show range expansion starting during the Bølling–Allerød warm period, interrupted by the Younger Dryas cold, and resuming during the Holocene. The distribution of A. alba expands to the north-east in all future scenarios, whereas the potential (currently unrealized) range would be substantially reduced in southern Europe under RCP8.5. A. alba maintains its current range in central Europe despite competition by other thermophilous tree species. Our combined palaeoecological and model evidence suggest that A. alba may ensure important ecosystem services including stand and slope stability, infrastructure protection, and carbon sequestration under significantly warmer-than-present conditions in central Europe.

Spatial Distribution of Stem Cell-Like Keratinocytes in Dissected Compound Hair Follicles of the Dog

Hair cycle disturbances are common in dogs and comparable to some alopecic disorders in humans. A normal hair cycle is maintained by follicular stem cells which are predominately found in an area known as the bulge. Due to similar morphological characteristics of the bulge area in humans and dogs, the shared particularity of compound hair follicles as well as similarities in follicular biomarker expression, the dog is a promising model to study human hair cycle and stem cell disorders. To gain insight into the spatial distribution of follicular keratinocytes with stem cell potential in canine compound follicles, we microdissected hair follicles in anagen and telogen from skin samples of freshly euthanized dogs. The keratinocytes isolated from different locations were investigated for their colony forming efficiency, growth and differentiation potential as well as clonal growth. Our results indicate that i) compound and single hair follicles exhibit a comparable spatial distribution pattern with respect to cells with high growth potential and stem cell-like characteristics, ii) the lower isthmus (comprising the bulge) harbors most cells with high growth potential in both, the anagen and the telogen hair cycle stage, iii) unlike in other species, colonies with highest growth potential are rather small with an irregular perimeter and iv) the keratinocytes derived from the bulbar region exhibit characteristics of actively dividing transit amplifying cells. Our results now provide the basis to conduct comparative studies of normal dogs and those with hair cycle disorders with the possibility to extend relevant findings to human patients.