Classification using distance nearest neighbours

This paper proposes a new probabilistic classification algorithm using a Markov random field approach. The joint distribution of class labels is explicitly modelled using the distances between feature vectors. Intuitively, a class label should depend more on class labels which are closer in the feature space, than those which are further away. Our approach builds on previous work by Holmes and Adams (2002, 2003) and Cucala et al. (2008). Our work shares many of the advantages of these approaches in providing a probabilistic basis for the statistical inference. In comparison to previous work, we present a more efficient computational algorithm to overcome the intractability of the Markov random field model. The results of our algorithm are encouraging in comparison to the k-nearest neighbour algorithm.Comment: 12 pages, 2 figures. To appear in Statistics and Computin

Satisfaction survey with DNA cards method to collect genetic samples for pharmacogenetics studies

BACKGROUND: Pharmacogenetic studies are essential in understanding the interindividual variability of drug responses. DNA sample collection for genotyping is a critical step in genetic studies. A method using dried blood samples from finger-puncture, collected on DNA-cards, has been described as an alternative to the usual venepuncture technique. The purpose of this study is to evaluate the implementation of the DNA cards method in a multicentre clinical trial, and to assess the degree of investigators' satisfaction and the acceptance of the patients perceived by the investigators. METHODS: Blood samples were collected on DNA-cards. The quality and quantity of DNA recovered were analyzed. Investigators were questioned regarding their general interest, previous experience, safety issues, preferences and perceived patient satisfaction. RESULTS: 151 patients' blood samples were collected. Genotyping of GST polymorphisms was achieved in all samples (100%). 28 investigators completed the survey. Investigators perceived patient satisfaction as very good (60.7%) or good (39.3%), without reluctance to finger puncture. Investigators preferred this method, which was considered safer and better than the usual methods. All investigators would recommend using it in future genetic studies. CONCLUSION: Within the clinical trial setting, the DNA-cards method was very well accepted by investigators and patients (in perception of investigators), and was preferred to conventional methods due to its ease of use and safety

Resilience of seagrass populations to thermal stress does not reflect regional differences in ocean climate

The prevalence of local adaptation and phenotypic plasticity among populations is critical to accurately predicting when and where climate change impacts will occur. Currently, comparisons of thermal performance between populations are untested for most marine species or overlooked by models predicting the thermal sensitivity of species to extirpation. Here we compared the ecological response and recovery of seagrass populations (Posidonia oceanica) to thermal stress throughout a year-long translocation experiment across a 2800-km gradient in ocean climate. Transplants in central and warm-edge locations experienced temperatures > 29°C, representing thermal anomalies > 5°C above long-term maxima for cool-edge populations, 1.5°C for central and < 1°C for warm-edge populations. Cool-edge, central and warm-edge populations differed in thermal performance when grown under common conditions, but patterns contrasted with expectations based on thermal geography. Cool-edge populations did not differ from warm-edge populations under common conditions and performed significantly better than central populations in growth and survival. Our findings reveal that thermal performance does not necessarily reflect the thermal geography of a species. We demonstrate that warm-edge populations can be less sensitive to thermal stress than cooler, central populations suggesting that Mediterranean seagrasses have greater resilience to warming than current paradigms suggest

Ghost admixture in eastern gorillas

Archaic admixture has had a substantial impact on human evolution with multiple events across different clades, including from extinct hominins such as Neanderthals and Denisovans into modern humans. In great apes, archaic admixture has been identified in chimpanzees and bonobos but the possibility of such events has not been explored in other species. Here, we address this question using high-coverage whole-genome sequences from all four extant gorilla subspecies, including six newly sequenced eastern gorillas from previously unsampled geographic regions. Using approximate Bayesian computation with neural networks to model the demographic history of gorillas, we find a signature of admixture from an archaic ‘ghost’ lineage into the common ancestor of eastern gorillas but not western gorillas. We infer that up to 3% of the genome of these individuals is introgressed from an archaic lineage that diverged more than 3 million years ago from the common ancestor of all extant gorillas. This introgression event took place before the split of mountain and eastern lowland gorillas, probably more than 40 thousand years ago and may have influenced perception of bitter taste in eastern gorillas. When comparing the introgression landscapes of gorillas, humans and bonobos, we find a consistent depletion of introgressed fragments on the X chromosome across these species. However, depletion in protein-coding content is not detectable in eastern gorillas, possibly as a consequence of stronger genetic drift in this species

Ghost admixture in eastern gorillas

Archaic admixture has had a substantial impact on human evolution with multiple events across different clades, including from extinct hominins such as Neanderthals and Denisovans into modern humans. In great apes, archaic admixture has been identified in chimpanzees and bonobos but the possibility of such events has not been explored in other species. Here, we address this question using high-coverage whole-genome sequences from all four extant gorilla subspecies, including six newly sequenced eastern gorillas from previously unsampled geographic regions. Using approximate Bayesian computation with neural networks to model the demographic history of gorillas, we find a signature of admixture from an archaic 'ghost' lineage into the common ancestor of eastern gorillas but not western gorillas. We infer that up to 3% of the genome of these individuals is introgressed from an archaic lineage that diverged more than 3 million years ago from the common ancestor of all extant gorillas. This introgression event took place before the split of mountain and eastern lowland gorillas, probably more than 40 thousand years ago and may have influenced perception of bitter taste in eastern gorillas. When comparing the introgression landscapes of gorillas, humans and bonobos, we find a consistent depletion of introgressed fragments on the X chromosome across these species. However, depletion in protein-coding content is not detectable in eastern gorillas, possibly as a consequence of stronger genetic drift in this species

Bayesian Computation with Intractable Likelihoods

This article surveys computational methods for posterior inference with intractable likelihoods, that is where the likelihood function is unavailable in closed form, or where evaluation of the likelihood is infeasible. We review recent developments in pseudo-marginal methods, approximate Bayesian computation (ABC), the exchange algorithm, thermodynamic integration, and composite likelihood, paying particular attention to advancements in scalability for large datasets. We also mention R and MATLAB source code for implementations of these algorithms, where they are available.Comment: arXiv admin note: text overlap with arXiv:1503.0806

Fuzzy maximum capacity and occupancy time rate measurements in urban railway lines

Abstract Nowadays, railway transport capacity is an important bottleneck for many railway operators that face its ever-increasing demand. This poses a challenge to existing lines as operation under conditions close to saturation tends be unstable. Capacity in urban railway systems depends largely on dwell times at platforms, but capacity measures proposed in literature rarely include the uncertainty associated to these times. In this paper this uncertainty is modeled as fuzzy numbers and two new capacity measures are proposed: the Fuzzy Maximum Capacity and the Fuzzy Occupancy Time Rate. The proposed model makes use of a railway simulator that enables route compression to obtain the conflict-free compressed time of the section under study. Three practical capacity problems from the perspective of the railway traffic operator have been presented and solved. The new measures provide more information to the railway operator than the standard UIC method that does not include uncertainty regarding dwell times. Finally, the model has been applied to the section Gràcia-Sarrià, belonging to the Spanish railway operator FGC