Nested Sampling Research Articles

Inconvenient hyperdiversity - the traditional concept of "PheIdole pallIdula" includes four cryptic species (Hymenoptera: Formicidae)

Results of a pilot study are presented suggesting that Edward Wilson’s guess of some 1500 species in the hyperdiverse genus Pheidole Westwood, 1839 is likely to turn out as rather moderate estimate. Investigation of 112 nest samples containing 287 major workers of Westpalaearctic ants currently named Pheidole pallidula (Nylander, 1849) was performed by the explorative data analyses NC-Ward and NC-k-means clustering in combination with cross-validated linear discriminant analysis. Seventeen primary morphometric data were recorded by high-resolution stereomicroscopy. Allometric variance of shape variables was removed. Four cryptic species, forming the Ph. pallidula complex, were resolved: Ph. pallidula, Ph. balcanica nov. sp., Ph. koshewnikovi Ruzsky, 1905 and Ph. cicatricosa Stitz, 1917. The classification error varied between 0 and 2.8 % in NC clustering of nest samples and between 1.6 and 3.7 % in cross-validated linear discriminant analysis of individuals. Ph. cicatricosa has a North African distribution whereas Ph. pallidula, Ph. balcanica nov. sp. and Ph. koshewnikovi are Eurasian species with large sympatric ranges in the Balkans and Asia Minor. Colonization of urban regions north of the Alps and of offshore islands in the Mediterranean Sea indicates an invasive potential of the supercolonial social type of Ph. pallidula. There are no indications that any of these four taxa might represent an intraspecific polymorphism and signals for putative interspecific hybridization are not significant. The following synonymies were established: Ph. subdentata Mayr, 1853, Ph. pallidula var. obscura Santschi, 1936 and Xenoaphaenogaster inquilina Baroni Urbani, 1964 are junior synonyms of Ph. pallidula, Ph. pallidula var. arenarum Ruzsky, 1905 and Ph. pallidula ssp. orientalis Muller, 1923 are junior synonyms of Ph. koshewnikovi and Ph. pallidula var. recticeps Menozzi, 1932 is a junior synonym of Ph. cicatricosa. The senior synonymy of the following seven Westpalaearctic taxa with any member of the Ph. pallidula complex was excluded by type investigation and diagostic statements in the original descriptions: Ph. sinaitica Mayr, 1862, Ph. jordanica Saulcy, 1874, Ph. teneriffana Forel, 1893, Ph. laticeps Mayr, 1904, Ph. schmitzi Forel, 1911, Ph. obtusa Stitz, 1917 and Ph. pallidula selenia Ozdikmen, 2010. Species delimitation in the less differentiated caste of minor workers was not tested but seems possible when accessory morphological characters are included.

A comparative study of oribatid mite communities in red wood ant Formica polyctena nests and surrounding soil in a Finnish oak forest

Abstract Ant nests are also inhabited by other invertebrates, ant associates. The oribatid mite (Acari: Oribatida) fauna of nests is poorly known, both in terms of abundance and diversity. We compared the oribatid assemblages within and outside red wood ant Formica polyctena nests in a Finnish oak forest. Altogether ten mounds were investigated and comparable samples were collected from adjacent soil 4 m from each mound. A total of 10 529 specimens were recovered, 4959 of them in mounds and 5633 in soil. The specimens represented 74 species (and two higher taxa) belonging to 60 genera and 39 families. Significant differences in abundance or diversity were not observed between nests and soil habitats. For nest samples the abundance of specimens significantly correlated with diversity indicating that mounds form rich habitat patches with a diverse set of niches. Mound size had no effect on community composition. The nest and soil habitats were predominantly occupied by different species resulting in dissimilarities in community composition; 34% of the species showed a preference for nest habitat and 50% for soil habitats, whereas16% were indifferent with equal abundances in both habitats. Nonetheless, most of the species with a preference for nests are frequently found in various habitats in nature, and may be classified as facultative myrmecophiles that maintain populations also in ant nests. Ant nests host abundant and diverse oribatid assemblages, which differ from the fauna in adjacent soil, and wood ants are therefore important factors in maintaining biodiversity.

Exploiting molecular dynamics in Nested Sampling simulations of small peptides

Nested Sampling (NS) is a parameter space sampling algorithm which can be used for sampling the equilibrium thermodynamics of atomistic systems. NS has previously been used to explore the potential energy surface of a coarse-grained protein model and has significantly outperformed parallel tempering when calculating heat capacity curves of Lennard-Jones clusters. The original NS algorithm uses Monte Carlo (MC) moves; however, a variant, Galilean NS, has recently been introduced which allows NS to be incorporated into a molecular dynamics framework, so NS can be used for systems which lack efficient prescribed MC moves. In this work we demonstrate the applicability of Galilean NS to atomistic systems. We present an implementation of Galilean NS using the Amber molecular dynamics package and demonstrate its viability by sampling alanine dipeptide, both in vacuo and implicit solvent. Unlike previous studies of this system, we present the heat capacity curves of alanine dipeptide, whose calculation provides a stringent test for sampling algorithms. We also compare our results with those calculated using replica exchange molecular dynamics (REMD) and find good agreement. We show the computational effort required for accurate heat capacity estimation for small peptides. We also calculate the alanine dipeptide Ramachandran free energy surface for a range of temperatures and use it to compare the results using the latest Amber force field with previous theoretical and experimental results.

Shiny Cowbird (Molothrus bonariensis) Parasitism Records for Three Globally Threatened Species from the South American Pampas

ABSTRACT The Shiny Cowbird (Molothrus bonariensis) is a generalist brood parasite known to have parasitized more than 200 species, many of which are able to rear cowbird young. Here I provide new records of Shiny Cowbird parasitism on three globally threatened birds from the Pampas of South America: Black-and-white Monjita (Xolmis dominicanus; vulnerable), Pampas Meadowlark (Sturnella defilippii; vulnerable), and Yellow Cardinal (Gubernatrix cristata; endangered). Cowbirds parasitized 24 out of 32 monjita nests, four out of 17 cardinal nests and three out of 169 meadowlark nests. These are the largest nest samples reported for the monjita and for the meadowlark, documented detailed parasitism records were lacking until now. In the case of the cardinal, data reported here complement information recently obtained in northeastern Argentina. Although habitat loss and illegal trapping may represent more important threats than brood parasitism, more studies are needed in order to establish the impact of brood p...

Comparative Analysis of Microbial Diversity in Termite Gut and Termite Nest Using Ion Sequencing.

Termite gut and termite nest possess complex microbial communities. However, only limited information is available on the comparative investigation of termite gut- and nest-associated microbial communities. In the present study, we examined and compared the bacterial diversity of termite gut and their respective nest by high-throughput sequencing of V3 hypervariable region of 16S rDNA. A total of 14 barcoded libraries were generated from seven termite gut samples and their respective nest samples, and sequenced using Ion Torrent platform. The sequences of each group were pooled, which yielded 170,644 and 132,000 reads from termite gut and termite nest samples, respectively. Phylogenetic analysis revealed significant differences in the bacterial diversity and community structure between termite gut and termite nest samples. Phyla Verrucomicrobia and Acidobacteria were observed only in termite gut, whereas Synergistetes and Chlorobi were observed only in termite nest samples. These variations in microbial structure and composition could be attributed with the differences in physiological conditions prevailing in the termite gut (anoxic and alkaline) and termite nest (oxic, slightly acidic and rich in organic matter) environment. Overall, this study unmasked the complexity of bacterial population in the respective niche. Interestingly, majority of the sequence reads could be classified only up to the domain level indicating the presence of a huge number of uncultivable or unidentified novel bacterial species in both termite gut and nest samples. Whole metagenome sequencing and assessing the metabolic potential of these samples will be useful for biotechnological applications.

Point process-based Monte Carlo estimation

This paper addresses the issue of estimating the expectation of a real-valued random variable of the form $$X = g(\mathbf {U})$$X=g(U) where g is a deterministic function and $$\mathbf {U}$$U can be a random finite- or infinite-dimensional vector. Using recent results on rare event simulation, we propose a unified framework for dealing with both probability and mean estimation for such random variables, i.e. linking algorithms such as Tootsie Pop Algorithm or Last Particle Algorithm with nested sampling. Especially, it extends nested sampling as follows: first the random variable X does not need to be bounded any more: it gives the principle of an ideal estimator with an infinite number of terms that is unbiased and always better than a classical Monte Carlo estimator--in particular it has a finite variance as soon as there exists $$k \in \mathbb {R}> 1$$k?R>1 such that $${\text {E}}\left[ X^k \right] < \infty $$EXk<?. Moreover we address the issue of nested sampling termination and show that a random truncation of the sum can preserve unbiasedness while increasing the variance only by a factor up to 2 compared to the ideal case. We also build an unbiased estimator with fixed computational budget which supports a Central Limit Theorem and discuss parallel implementation of nested sampling, which can dramatically reduce its running time. Finally we extensively study the case where X is heavy-tailed.

Nested sampling of isobaric phase space for the direct evaluation of the isothermal-isobaric partition function of atomic systems.

Nested Sampling (NS) is a powerful athermal statistical mechanical sampling technique that directly calculates the partition function, and hence gives access to all thermodynamic quantities in absolute terms, including absolute free energies and absolute entropies. NS has been used predominately to compute the canonical (NVT) partition function. Although NS has recently been used to obtain the isothermal-isobaric (NPT) partition function of the hard sphere model, a general approach to the computation of the NPT partition function has yet to be developed. Here, we describe an isobaric NS (IBNS) method which allows for the computation of the NPT partition function of any atomic system. We demonstrate IBNS on two finite Lennard-Jones systems and confirm the results through comparison to parallel tempering Monte Carlo. Temperature-entropy plots are constructed as well as a simple pressure-temperature phase diagram for each system. We further demonstrate IBNS by computing part of the pressure-temperature phase diagram of a Lennard-Jones system under periodic boundary conditions.

Differentiation between house and cave edible bird's nests by chemometric analysis of amino acid composition data

A total of 60 pieces of house and cave edible bird's nest (EBN) samples were randomly collected from different locations in Malaysia and Indonesia. Amino acid compositions of the EBN samples were determined by gas chromatography-mass spectrometry (GC-MS). Data obtained was analyzed by Pearson Correlation analysis, Principal Component analysis (PCA) and Orthogonal Partial Least Square-Discriminant analysis (OPLS-DA). There were highly significant different correlations seen among amino acids in house- and cave-EBN samples. The model constructed by OPLS-DA was found to be a promising tool with high predictive power of 76.1%. Robustness of the model was validated and blind test samples were correctly assigned to their respective cluster. Tyrosine (TYR) and glutamic acid (GLU) were proposed as the promising markers for differentiating between house and cave EBN.

Security in big data

Security in big data

Bayesian peak bagging analysis of 19 low-mass low-luminosity red giants observed withKepler

The currently available Kepler light curves contain an outstanding amount of information but a detailed analysis of the individual oscillation modes in the observed power spectra, also known as peak bagging, is computationally demanding and challenging to perform on a large number of targets. Our intent is to perform for the first time a peak bagging analysis on a sample of 19 low-mass low-luminosity red giants observed by Kepler for more than four years. This allows us to provide high-quality asteroseismic measurements that can be exploited for an intensive testing of the physics used in stellar structure models, stellar evolution and pulsation codes, as well as for refining existing asteroseismic scaling relations in the red giant branch regime. For this purpose, powerful and sophisticated analysis tools are needed. We exploit the Bayesian code Diamonds, using an efficient nested sampling Monte Carlo algorithm, to perform both a fast fitting of the individual oscillation modes and a peak detection test based on the Bayesian evidence. We find good agreement for the parameters estimated in the background fitting phase with those given in the literature. We extract and characterize a total of 1618 oscillation modes, providing the largest set of detailed asteroseismic mode measurements ever published. We report on the evidence of a change in regime observed in the relation between linewidths and effective temperatures of the stars occurring at the bottom of the RGB. We show the presence of a linewidth depression or plateau around $\nu_\mathrm{max}$ for all the red giants of the sample. Lastly, we show a good agreement between our measurements of maximum mode amplitudes and existing maximum amplitudes from global analyses provided in the literature, useful as empirical tools to improve and simplify the future peak bagging analysis on a larger sample of evolved stars.

Spatial distribution of soil organic matter in two fields on tidal flat sediments (Zhejiang Province, China) differing in duration of paddy management

AbstractDespite of its ecological and economic importance, only limited information is available on spatial variability of soil organic C (OC) in paddy fields during soil development. To obtain deeper insights into the accumulation and spatial distribution of OC with ongoing paddy management, the present study investigates spatial heterogeneities in puddled topsoil layers, in the plough pan, as well as in the underlying subsoil of a 50‐ and 300‐y old paddy site at the area of Cixi (Hangzhou Bay, E China). Sampling of the plots was done on the basis of a regular, orthogonal grid. The grid combined a large‐scale design with distances of 5 m with a small‐scale nested sampling ranging within 0.5 to 1 m. Samples have been taken with an auger in three sampling depths: (1) puddled topsoil layer, (2) dense plough layer, (3) mixed subsoil horizon. All samples were analyzed for their OC and total N concentrations as well as for their soil color (lightness value) to assess relationships between soil lightness and OC amount and distribution. The spatial behavior of all analyzed parameters in the different soil layers was investigated by calculating variograms; the degree of spatial dependence was determined by nugget‐to‐sill ratios. The results show that for each analyzed soil layer, OC and N concentrations were significantly higher in the field with 300 y of paddy management than in the field with 50 y of paddy management. Topsoil OC concentrations changed from 13.8 to 23.6 mg g−1, OC in the plough layer changed from 3.3 to 9.6 mg g−1, and also the subsoil showed a slight increase from 2.9 to 3.3 mg g−1.Although only two points in time were compared, we assume that paddy soil formation is the most likely reason for differences in the spatial distribution of OC between the two sites. Comparing the spatial depth distribution of the 50‐ and 300‐y old paddy sites suggests that the development of a spatial pattern of OC proceeded with time from topsoil to the underlying soil horizons. Measuring soil lightness was insufficient to predict amount and distribution of soil OC in paddy soils, especially in soils with interfering dark‐colored burnt straw residues. The degree of spatial homogeneity was equal in the two sites, suggesting that cultivation practices (puddling, levelling, burning and incorporation of straw residues) most likely maintained a rather high random spatial variability of approx. 45%. The spatial pattern of OC in the topsoil was dissimilar to the plough layer, revealing that transport of OC in the plough pan occurred mainly in preferential flow paths along localized cracks and biopores.

Ellipsoidal nested sampling, expression of the model uncertainty and measurement

The measurand value, the conclusions, and the decisions inferred from measurements may depend on the models used to explain and to analyze the results. In this paper, the problems of identifying the most appropriate model and of assessing the model contribution to the uncertainty are formulated and solved in terms of Bayesian model selection and model averaging. As computational cost of this approach increases with the dimensionality of the problem, a numerical strategy, based on multimodal ellipsoidal nested sampling, to integrate over the nuisance parameters and to compute the measurand post-data distribution is outlined. In order to illustrate the numerical strategy, by use of MATHEMATICA an elementary example concerning a bimodal, two-dimensional distribution has also been studied.

Surface geochemical measurements applied to monitoring, verification, and accounting of leakage from sequestration projects

A detailed study of [Formula: see text] and [Formula: see text] surface fluxes, shallow soil gas composition, 10-m hole soil gas composition was made at Rangely, Colorado, Teapot Dome, Wyoming, USA, and the results of geochemical verification measurements at Weyburn, Saskatchewan, Canada. Summer and winter soil gas and gas flux measurements were made at Rangely, and winter only at Teapot Dome. The objectives were to determine if leakage of [Formula: see text] and/or [Formula: see text] could be detected in the overpressured Rangely and the underpressured Teapot Dome systems. At Weyburn, the objective was the determination of the presence or absence of an alleged leak. Seasonal surface [Formula: see text] fluxes at Rangely were similar at on-field locations and an off-field control area. Methane fluxes were much higher at on-field locations than in the control area, suggesting a reservoir source. Seasonal differences in [Formula: see text] fluxes suggested methanotrophic oxidation was occurring in the soils. Shallow soil gas concentration measurements of [Formula: see text] were similar at a 100-cm depth, on-field and in the control area. Methane soil gas concentrations were higher on-field than in the control area; summer and winter. These data were used to select locations with and without evidence for seepage allowing nested sampling up to a 10-m depth. More complex chemical and isotopic measurements were made in 10-m holes at Rangely and Teapot Dome. Methanotrophy was operational at the active Rangely and the passive Teapot Dome systems. A small methane leakage rate of 400–700 tonnes [Formula: see text] and [Formula: see text] of [Formula: see text] tonnes [Formula: see text] were determined at Rangely, and rates near zero at Teapot Dome. Inert gases, carbon-containing gases, as well as isotopic ratios supported the presence of pathways at geochemically anomalous locations at Teapot Dome. Measurements of soil gas and inert gas isotopes in shallow groundwater by other authors at Weyburn, Saskatchewan, found no evidence of leakage.

Abundance and genetic damage of barn swallows from Fukushima.

A number of studies have assessed or modeled the distribution of the radionuclides released by the accident at the Fukushima-Daiichi Nuclear Power Plant (FDNPP). Few studies however have investigated its consequences for the local biota. We tested whether exposure of barn swallow (Hirundo rustica) nestlings to low dose ionizing radiation increased genetic damage to their peripheral erythrocytes. We estimated external radiation exposure by using thermoluminescent dosimeters, and by measuring radioactivity of the nest material. We then assessed DNA damage by means of the neutral comet assay. In addition, we conducted standard point-count censuses of barn swallows across environmental radiation levels, and estimated their abundance and local age ratio. Radioactivity of nest samples was in the range 479–143,349 Bq kg−1, while external exposure varied between 0.15 and 4.9 mGy. Exposure to radioactive contamination did not correlate with higher genetic damage in nestlings. However, at higher levels of radioactive contamination the number of barn swallows declined and the fraction of juveniles decreased, indicating lower survival and lower reproduction and/or fledging rate. Thus, genetic damage to nestlings does not explain the decline of barn swallows in contaminated areas, and a proximate mechanism for the demographic effects documented here remains to be clarified.

Assessment of variation of nest survival for grassland birds due to method of nest discovery

Capsule Interpretation of nest survival estimates may be improved by incorporating the search method used to locate nests as a covariate.Aims To compare annual survival estimates for Dickcissel Spiza americana nests and determine if incorporating search method (structured, opportunistic, or behavioural searches) improved model fit.Methods Dickcissel nests were located using structured, opportunistic, or behavioural searches over three years (2011–2013) in Mississippi, USA. Models were used to estimate daily survival rates (DSRs) and to analyse factors influencing nest survival.Results DSRs for Dickcissels were best explained by quadratic date, nest age, age found, and year, but incorporating search method improved model fit. Daily survival was 1.51 times greater for nests located using opportunistic search methods relative to structured searches, but was not significantly different between structured and behavioural searches.Conclusions Survival estimates varied by search method, specifically between structured searches and opportunistically located nests. This might have arisen because heterogeneity in nest placement or parental behaviour may influence the sample of nests located with a given search method. Researchers may be able to account for this potential source of bias by including search method as a model covariate when using standard survey designs or modelling approaches.

Exploiting the potential energy landscape to sample free energy

We review a number of recently developed strategies for enhanced sampling of complex systems based on knowledge of the potential energy landscape. We describe four approaches, replica exchange, Kirkwood sampling, superposition‐enhanced nested sampling, and basin sampling, and show how each of them can exploit information for low‐lying potential energy minima obtained using basin‐hopping global optimization. Characterizing these minima is generally much faster than equilibrium thermodynamic sampling, because large steps in configuration space between local minima can be used without concern for maintaining detailed balance. WIREs Comput Mol Sci 2015, 5:273–289. doi: 10.1002/wcms.1217This article is categorized under: Structure and Mechanism &gt; Molecular Structures Structure and Mechanism &gt; Computational Biochemistry and Biophysics Structure and Mechanism &gt; Computational Materials Science

Temnothoraxcrasecundus sp. n. - a cryptic Eurocaucasian ant species (Hymenoptera, Formicidae) discovered by Nest Centroid Clustering.

The paper integrates two independent studies of numeric morphology-based alpha-taxonomy of the cryptic ant species Temnothorax crassispinus (Karavajev, 1926) and Temnothorax crasecundus sp. n. conducted by different investigators, using different equipment, considering different character combinations and evaluating different samples. Samples investigated included 603 individual workers from 203 nests – thereof 104 nest samples measured by Seifert and 99 by Csösz. The material originated from Europe, Asia Minor and Caucasia. There was a very strong interspecific overlap in any of the 29 shape characters recorded and subjective expert determination failed in many cases. Primary classification hypotheses were formed by the exploratory data analysis Nest Centroid (NC) clustering and corrected to final species hypotheses by an iterative linear discriminant analysis algorithm. The evaluation of Seifert’s and Csösz’s data sets arrived at fully congruent conclusions. NC-Ward and NC-K-means clustering disagreed from the final species hypothesis in only 1.9 and 1.9% of the samples in Seifert’s data set and by 1.1 and 2.1% in Csösz’s data set which is a strong argument for heterospecificity. The type series of Temnothorax crassispinus and Temnothorax crasecundus sp. n. were allocated to different clusters with p = 0.9851 and p = 0.9912 respectively. The type series of the junior synonym Temnothorax slavonicus (Seifert, 1995) was allocated to the Temnothorax crassispinus cluster with p = 0.9927. Temnothorax crasecundus sp. n. and Temnothorax crassispinus are parapatric species with a long contact zone stretching from the Peloponnisos peninsula across Bulgaria northeast to the southern Ukraine. There is no indication for occurrence of interspecifically mixed nests or intraspecific polymorphism. However, a significant reduction of interspecific morphological distance at sites with syntopic occurrence of both species indicates local hybridization. The results are discussed within the context of the Pragmatic Species Concept of Seifert (2014). The taxonomic description and a differential diagnosis of Temnothorax crasecundus sp. n. are given.

DIAMONDS: a new Bayesian nested sampling tool

In the context of high-quality asteroseismic data provided by the NASA Kepler mission, we developed a new code, termed Diamonds (high-DImensional And multi-MOdal NesteD Sampling), for fast Bayesian parameter estimation and model comparison by means of the Nested Sampling Monte Carlo (NSMC) algorithm, an efficient and powerful method very suitable for high-dimensional problems (like the peak bagging analysis of solar-like oscillations) and multi-modal problems (i.e. problems that show multiple solutions). We applied the code to the peak bagging analysis of solar-like oscillations observed in a challenging F-type star. By means of Diamonds one is able to detect the different backgrounds in the power spectrum of the star (e.g. stellar granulation and faculae activity) and to understand whether one or two oscillation peaks can be identified or not. In addition, we demonstrate a novel approach to peak bagging based on multimodality, which is able to reduce significantly the number of free parameters involved in the peak bagging model. This novel approach is therefore of great interest for possible future automatization of the entire analysis technique.

Peak Bagging of red giant stars observed byKepler: first results with a new method based on Bayesian nested sampling

The peak bagging analysis, namely the fitting and identification of single oscillation modes in stars' power spectra, coupled to the very high-quality light curves of red giant stars observed by Kepler, can play a crucial role for studying stellar oscillations of different flavor with an unprecedented level of detail. A thorough study of stellar oscillations would thus allow for deeper testing of stellar structure models and new insights in stellar evolution theory. However, peak bagging inferences are in general very challenging problems due to the large number of observed oscillation modes, hence of free parameters that can be involved in the fitting models. Effciency and robustness in performing the analysis is what may be needed to proceed further. For this purpose, we developed a new code implementing the Nested Sampling Monte Carlo (NSMC) algorithm, a powerful statistical method well suited for Bayesian analyses of complex problems. In this talk we show the peak bagging of a sample of high signal-to-noise red giant stars by exploiting recent Kepler datasets and a new criterion for the detection of an oscillation mode based on the computation of the Bayesian evidence. Preliminary results for frequencies and lifetimes for single oscillation modes, together with acoustic glitches, are therefore presented.

DIAMONDS: A new Bayesian nested sampling tool

To exploit the full potential of Kepler light curves, sophisticated and robust analysis tools are now required more than ever. Characterizing single stars with an unprecedented level of accuracy and subsequently analyzing stellar populations in detail are fundamental to further constrain stellar structure and evolutionary models. We developed a new code, termed Diamonds, for Bayesian parameter estimation and model comparison by means of the nested sampling Monte Carlo (NSMC) algorithm, an efficient and powerful method very suitable for high-dimensional and multi-modal problems. A detailed description of the features implemented in the code is given with a focus on the novelties and differences with respect to other existing methods based on NSMC. Diamonds is then tested on the bright F8 V star KIC~9139163, a challenging target for peak-bagging analysis due to its large number of oscillation peaks observed, which are coupled to the blending that occurs between $\ell=2,0$ peaks, and the strong stellar background signal. We further strain the performance of the approach by adopting a 1147.5 days-long Kepler light curve. The Diamonds code is able to provide robust results for the peak-bagging analysis of KIC~9139163. We test the detection of different astrophysical backgrounds in the star and provide a criterion based on the Bayesian evidence for assessing the peak significance of the detected oscillations in detail. We present results for 59 individual oscillation frequencies, amplitudes and linewidths and provide a detailed comparison to the existing values in the literature. Lastly, we successfully demonstrate an innovative approach to peak bagging that exploits the capability of Diamonds to sample multi-modal distributions, which is of great potential for possible future automatization of the analysis technique.