A genetic and spatial Bayesian analysis of mastitis resistance

A nationwide health card recording system for dairy cattle was introduced in Norway in 1975 (the Norwegian Cattle Health Services). The data base holds information on mastitis occurrences on an individual cow basis. A reduction in mastitis frequency across the population is desired, and for this purpose risk factors are investigated. In this paper a Bayesian proportional hazards model is used for modelling the time to first veterinary treatment of clinical mastitis, including both genetic and environmental covariates. Sire effects were modelled as shared random components, and veterinary district was included as an environmental effect with prior spatial smoothing. A non-informative smoothing prior was assumed for the baseline hazard, and Markov chain Monte Carlo methods (MCMC) were used for inference. We propose a new measure of quality for sires, in terms of their posterior probability of being among the, say 10% best sires. The probability is an easily interpretable measure that can be directly used to rank sires. Estimating these complex probabilities is straightforward in an MCMC setting. The results indicate considerable differences between sires with regards to their daughters disease resistance. A regional effect was also discovered with the lowest risk of disease in the south-eastern parts of Norway.     

Novel Monte Carlo dose calculation algorithm for robotic radiosurgery with multi leaf collimator: Dosimetric evaluation

PurposeAt introduction in 2014, dose calculation for the first MLC on a robotic SRS/SBRT platform was limited to a correction-based Finite-Size Pencil Beam (FSPB) algorithm. We report on the dosimetric accuracy of a novel Monte Carlo (MC) dose calculation algorithm for this MLC, included in the Precision™ treatment planning system.MethodsA phantom was built of one slab (5.0 cm) of lung-equivalent material (0.09…0.29 g/cc) enclosed by 3.5 cm (above) and 5 cm (below) slabs of solid water (1.045 g/cc). This was irradiated using rectangular (15.4 × 15.4 mm² to 53.8 × 53.7 mm²) and two irregular MLC-fields. Radiochromic film (EBT3) was positioned perpendicular to the slabs and parallel to the beam. Calculated dose distributions were compared to film measurements using line scans and 2D gamma analysis.ResultsMeasured and MC calculated percent depth dose curves showed a characteristic dose drop within the low-density region, which was not correctly reproduced by FSPB. Superior average gamma pass rates (2%/1 mm) were found for MC (91.2 ± 1.5%) compared to FSPB (55.4 ± 2.7%). However, MC calculations exhibited localized anomalies at mass density transitions around 0.15 g/cc, which were traced to a simplification in electron transport. Absence of these anomalies was confirmed in a modified build of the MC engine, which increased gamma pass rates to 96.6 ± 1.2%.ConclusionsThe novel MC algorithm greatly improves dosimetric accuracy in heterogeneous tissue, potentially expanding the clinical use of robotic radiosurgery with MLC. In-depth, independent validation is paramount to identify and reduce the residual uncertainties in any software solution.     

Robust Selection of the Most Discriminative Variables in the Dichotomous Problem with Application to Some Respiratory Disease Data

A robust method for selection of variables with the greatest discriminatory power is presented in the paper. The method deals with the two groups of data problem. An application of the method to some respiratory disease data and comparisons with classical procedures are given, also.     

The effect of dielectric polarization of the electrode on anomalous temperature effects in the electrical double layer

Monte Carlo simulations and a modified Poisson–Boltzmann (MPB) theory are used to investigate the temperature dependence of the capacitance (around the potential of zero charge) of an electric double layer in the presence of surface polarization due to a dielectric boundary. Within the context of the restricted primitive model planar double layer, whose solvent dielectric constant is ε₂, the cases when the electrode is an insulator (ε₁ = 1), when the electrode and the electrolyte have the same permittivity (ε₁ = ε₂, no polarization), and when the electrode is a conductor (ε₁ → ∞) are studied for the case where the electrolyte concentration is 0.1 M. The simulations reveal a capacitance anomaly, that is, a positive temperature dependence of the capacitance at low temperatures for the former two situations. The MPB theory also shows this effect for these two situations and is in qualitative or better agreement with the simulation data. In these two cases, both the simulations and theory show a dramatic increase of the diffuse layer potential in the temperature regime where capacitance anomaly occurs. However, in the latter situation, where the electrode is metallic, the capacitance always has a negative temperature derivative for the MPB theory and probably also for the simulation data.     

On computational efficiency of the hybrid Monte Carlo method applied to the multicanonical ensemble

Abstract

In the present paper, computational efficiency of the hybrid Monte Carlo (HMC) method applied to the multicanonical ensemble is studied; the HMC is an equation of motion guided Monte Carlo method. As in the standard HMC for the canonical ensemble, the multicanonical HMC calculations with high acceptance ratio show better efficiency; about 60% acceptance yields the best performance for the system examined.     

Reference values for feeding parameters of isopods (Porcellio scaber,Isopoda, Crustacea)

The advantage of using terrestrial isopods in toxicity studies is that a battery of parameters can be tested at different levels of biological complexity. Feeding parameters for example link organism level response to potential ecological consequences but a problem with using feeding parameters in toxicity tests with terrestrial isopods is their high variability. The aim of our study was to set benchmark values for feeding parameters for isopod Porcellio scaber (Isopoda, Crustacea) in laboratory-controlled experiments. In the work presented here, the daily feeding rate of the central 50% of the control population of Porcellio scaber and a correlation between feeding rate and isopod weight were set. Values outside these ranges need additional evaluation to increase the relevance of test outcomes. We suggest using benchmark values for feeding parameters as well as the coefficient of variation (a) to identify animals with altered feeding parameters with respect to controls, and (b) to assess the data quality in each experiment.     

Anatomically-based skeletal coordinate systems for use with impact biomechanics data intended for anthropomorphic test device development

Post-mortem human subjects (PMHS) are frequently used to characterize biomechanical response and injury tolerance of humans to various types of loading by means of instrumentation installed directly on the skeleton. Data extracted from such tests are often used to develop and validate anthropomorphic test devices (ATDs), which function as human surrogates in tests for injury assessment. Given that the location and orientation of installed instrumentation differs between subjects, nominally similar measurements made on different PMHS must be transformed to standardized, skeletal-based local coordinate systems (LCS) before appropriate data comparisons can be made. Standardized PMHS LCS that correspond to ATD instrumentation locations and orientations have not previously been published. This paper introduces anatomically-defined PMHS LCS for body regions in which kinematic measurements are made using ATDs. These LCS include the head, sternum, single vertebrae, pelvis, femurs (distal and proximal), and tibiae (distal and proximal) based upon skeletal landmarks extracted from whole body CT scans. The proposed LCS provide a means to standardize the reporting of PMHS data, and facilitate both the comparison of PMHS impact data across institutions and the application of PMHS data to the development and validation of ATDs.     

Progeny testing for recessive genes: procedures and interpretations

Summary Methods for calculating the probability of detecting a carrier of a recessive gene by utilizing matings among related individuals are presented for single and litter bearing species. The confidence level for detection of heterozygosity depends upon: (1) the genetic relationship between mates, (2) the number of mates per male and the number of offspring per mate, (3) whether an estimate of recessive gene frequency before selection is available and (4) the magnitude of that frequency. Methods of computing probability of heterozygosity vs homozygosity utilizing Bayes theorem also are presented. In the conventional progeny test method, a sire initially is assumed heterozygous before calculations are made, but no prior information concerning his probable genotype is utilized. In the method using Bayes theorem, prior sources of information from relatives or from estimates of population allele frequency are utilized. This method gives the exact probability that a sire is not a carrier, given prior information and that he produces all normal offspring. These methods could be used in any sexually reproducing species to identify not only detrimental genes but beneficial genes as well.     

Studies on evolutionary and selective properties of hypercycles using a Monte Carlo method

Summary The most relevant properties of hypercycles were previously studied mainly from a theoretical point of view. We have developed a Monte Carlo method simulating hypercyclic organization to obtain information about the dynamics of this prebiotic organization. Nucleation, growth, and selective properties have been tested and the results obtained are in good agreement with those of the theoretical predictions. The influence of hypercyclic organization of the error threshold has also been studied. As a consequence of the emergence of a hypercycle, the value of this threshold decreases. The amount of this decrease depends on the population size. Moreover, for some interval of quality factor values, either the hypercycle organization or an error catastrophe can be produced, depending on the initial conditions. The influence of these phenomena on both the dynamic behavior and evolutionary advantages of the hypercycle, as well as their decisive roles on genome size, are discussed.Presented at the FEBS Symposium on Genome Organization and Evolution, held in Crete, Greece, September 1–5, 1986     

Volume estimation of biological objects by systematic sections

The absolute volume of biological objects is often estimated stereologically from an exhaustive set of systematic sections. The usual volume estimator is the sum of the section contents times the distance between sections. For systematic sectioning with a random start, it has been recently shown that is unbiased when m, the ratio between projected object length and section distance, is an integer number (Cruz-Orive 1985). As this quantity is no integer in the real world, we have explored the properties of in the general and realistic situation m . The unbiasedness of under appropriate sampling conditions is demonstrated for the arbitrary compact set in 3 dimensions by a rigorous proof. Exploration of further properties of for the general triaxial ellipsoid leads to a new class of non-elementary real functions with common formal structure which we denote as np-functions. The relative mean square error (CE ²) of in ellipsoids is an oscillating differentiable np-function, which reduces to the known result CE ²= 1/(5m ⁴) for integer m. As a biological example the absolute volumes of 10 left cardiac ventricles and their internal cavities were estimated from systematic sections. Monte Carlo simulation of replicated systematic sectioning is shown to be improved by using m instead of m . In agreement with the geometric model of ellipsoids with some added shape irregularities, mean empirical CE was proportional to m ^–1.36 and m^–1.73 in the cardiac ventricle and its cavity. The considerable variance reduction by systematic sectioning is shown to be a geometric realization of the principle of antithetic variates.