Gibbs point process approximation: Total variation bounds using Stein's method

We obtain upper bounds for the total variation distance between the distributions of two Gibbs point processes in a very general setting. Applications are provided to various well-known processes and settings from spatial statistics and statistical physics, including the comparison of two Lennard-Jones processes, hard core approximation of an area interaction process and the approximation of lattice processes by a continuous Gibbs process. Our proof of the main results is based on Stein's method. We construct an explicit coupling between two spatial birth-death processes to obtain Stein factors, and employ the Georgii-Nguyen-Zessin equation for the total bound.

Total Variation Blind Deconvolution: The Devil Is in the Details

In this paper we study the problem of blind deconvolution. Our analysis is based on the algorithm of Chan and Wong [2] which popularized the use of sparse gradient priors via total variation. We use this algorithm because many methods in the literature are essentially adaptations of this framework. Such algorithm is an iterative alternating energy minimization where at each step either the sharp image or the blur function are reconstructed. Recent work of Levin et al. [14] showed that any algorithm that tries to minimize that same energy would fail, as the desired solution has a higher energy than the no-blur solution, where the sharp image is the blurry input and the blur is a Dirac delta. However, experimentally one can observe that Chan and Wong's algorithm converges to the desired solution even when initialized with the no-blur one. We provide both analysis and experiments to resolve this paradoxical conundrum. We find that both claims are right. The key to understanding how this is possible lies in the details of Chan and Wong's implementation and in how seemingly harmless choices result in dramatic effects. Our analysis reveals that the delayed scaling (normalization) in the iterative step of the blur kernel is fundamental to the convergence of the algorithm. This then results in a procedure that eludes the no-blur solution, despite it being a global minimum of the original energy. We introduce an adaptation of this algorithm and show that, in spite of its extreme simplicity, it is very robust and achieves a performance comparable to the state of the art.

A total variation discontinuous Galerkin approach for image restoration

The focal point of this paper is to propose and analyze a P 0 discontinuous Galerkin (DG) formulation for image denoising. The scheme is based on a total variation approach which has been applied successfully in previous papers on image processing. The main idea of the new scheme is to model the restoration process in terms of a discrete energy minimization problem and to derive a corresponding DG variational formulation. Furthermore, we will prove that the method exhibits a unique solution and that a natural maximum principle holds. In addition, a number of examples illustrate the effectiveness of the method.

A Clearer Picture of Total Variation Blind Deconvolution

Blind deconvolution is the problem of recovering a sharp image and a blur kernel from a noisy blurry image. Recently, there has been a significant effort on understanding the basic mechanisms to solve blind deconvolution. While this effort resulted in the deployment of effective algorithms, the theoretical findings generated contrasting views on why these approaches worked. On the one hand, one could observe experimentally that alternating energy minimization algorithms converge to the desired solution. On the other hand, it has been shown that such alternating minimization algorithms should fail to converge and one should instead use a so-called Variational Bayes approach. To clarify this conundrum, recent work showed that a good image and blur prior is instead what makes a blind deconvolution algorithm work. Unfortunately, this analysis did not apply to algorithms based on total variation regularization. In this manuscript, we provide both analysis and experiments to get a clearer picture of blind deconvolution. Our analysis reveals the very reason why an algorithm based on total variation works. We also introduce an implementation of this algorithm and show that, in spite of its extreme simplicity, it is very robust and achieves a performance comparable to the top performing algorithms.

Grain yield variation and association of major traits in brown-seeded genotypes of tef [Eragrostis tef (Zucc.)Trotter]

Background Tef [Eragrostis tef (Zucc.) Trotter] is the major cereal crop of Ethiopia where it is annually cultivated on more than three million hectares of land by over six million small-scale farmers. It is broadly grouped into white and brown-seeded type depending on grain color, although some intermediate color grains also exist. Earlier breeding experiments focused on white-seeded tef, and a number of improved varieties were released to the farming community. Thirty-six brown-seeded tef genotypes were evaluated using a 6 × 6 simple lattice design at three locations in the central highlands of Ethiopia to assess the productivity, heritability, and association among major pheno-morphic traits. Results The mean square due to genotypes, locations, and genotype by locations were significant (P < 0.01) for all traits studied. Genotypic and phenotypic coefficients of variations ranged from 2.5 to 20.3 % and from 4.3 to 21.7 %, respectively. Grain yield showed significant (P < 0.01) genotypic correlation with shoot biomass and harvest index, while it had highly significant (P < 0.01) phenotypic correlation with all the traits evaluated. Besides, association of lodging index with biomass and grain yield was negative and significant at phenotypic level while it was not significant at genotypic level. Cluster analysis grouped the 36 test genotypes into seven distinct classes. Furthermore, the first three principal components with eigenvalues greater than unity extracted 78.3 % of the total variation. Conclusion The current study, generally, revealed the identification of genotypes with superior grain yield and other desirable traits for further evaluation and eventual release to the farming community.

Cosmogenic 180W variations in meteorites and re-assessment of a possible 184Os–180W decay system

We measured tungsten (W) isotopes in 23 iron meteorites and the metal phase of the CB chondrite Gujba in order to ascertain if there is evidence for a large-scale nucleosynthetic heterogeneity in the p-process isotope 180W in the solar nebula as recently suggested by Schulz et al. (2013). We observed large excesses in 180W (up to ≈ 6 ε) in some irons. However, significant within-group variations in magmatic IIAB and IVB irons are not consistent with a nucleosynthetic origin, and the collateral effects on 180W from an s-deficit in IVB irons cannot explain the total variation. We present a new model for the combined effects of spallation and neutron capture reactions on 180W in iron meteorites and show that at least some of the observed within-group variability is explained by cosmic ray effects. Neutron capture causes burnout of 180W, whereas spallation reactions lead to positive shifts in 180W. These effects depend on the target composition and cosmic-ray exposure duration; spallation effects increase with Re/W and Os/W ratios in the target and with exposure age. The correlation of 180W/184W with Os/W ratios in iron meteorites results in part from spallogenic production of 180W rather than from 184Os decay, contrary to a recent study by Peters et al. (2014). Residual ε180W excesses after correction for an s-deficit and for cosmic ray effects may be due to ingrowth of 180W from 184Os decay, but the magnitude of this ingrowth is at least a factor of ≈2 smaller than previously suggested. These much smaller effects strongly limit the applicability of the putative 184Os-180W system to investigate geological problems.

Phylogenetic relationships of German heavy draught horse breeds inferred from mitochondrial DNA D-loop variation

We analysed a 610-bp mitochondrial (mt)DNA D-loop fragment in a sample of German draught horse breeds and compared the polymorphic sites with sequences from Arabian, Hanoverian, Exmoor, Icelandic, Sorraia and Przewalski's Horses as well as with Suffolk, Shire and Belgian horses. In a total of 65 horses, 70 polymorphic sites representing 47 haplotypes were observed. The average percentage of polymorphic sites was 11.5% for the mtDNA fragment analysed. In the nine different draught horse breeds including South German, Mecklenburg, Saxon Thuringa coldblood, Rhenisch German, Schleswig Draught Horse, Black Forest Horse, Shire, Suffolk and Belgian, 61 polymorphic sites and 24 haplotypes were found. The phylogenetic analysis failed to show monophyletic groups for the draught horses. The analysis indicated that the draught horse populations investigated consist of diverse genetic groups with respect to their maternal lineage.

Variation in human platelet glycoprotein VI content modulates glycoprotein VI-specific prothrombinase activity

- Glycoprotein VI (GPVI) is a platelet-specific receptor for collagen that figures prominently in signal transduction. An addition to binding to type I and III collagens, GPVI is also bound specifically by collagen-related peptide and convulxin (CVX), a snake venom protein. We developed a quantitative assay of platelet GPVI in which biotin-conjugated CVX binds selectively to GPVI in separated total platelet proteins by a ligand blot procedure. Using this approach, we have documented a 5-fold range in platelet GPVI content among 23 normal healthy subjects. In addition, we have determined that CVX-induced or collagen-related peptide-induced prothrombinase activity is directly proportional to the platelet content of GPVI. A statistically significant correlation was observed at 2 CVX concentrations: 14.7 ng/mL (R(2)=0.854 and P<0.001, n=11) and 22 ng/mL (R(2)=0.776 and P<0.001, n=12). In previous studies, we established a similar range of expression of the integrin collagen receptor alpha(2)beta(1) on platelets of normal subjects. Among 15 donors, there is a direct correlation between platelet alpha(2)beta(1) density and GPVI content (R(2)=0.475 and P=0.004). In view of the well-documented association of GPVI with platelet procoagulant activity, this study suggests that the variation in GPVI content is a potential risk factor that may predispose individuals to hemorrhagic or thromboembolic disorders.

Variation of a test's sensitivity and specificity with disease prevalence

BACKGROUND Anecdotal evidence suggests that the sensitivity and specificity of a diagnostic test may vary with disease prevalence. Our objective was to investigate the associations between disease prevalence and test sensitivity and specificity using studies of diagnostic accuracy. METHODS We used data from 23 meta-analyses, each of which included 10-39 studies (416 total). The median prevalence per review ranged from 1% to 77%. We evaluated the effects of prevalence on sensitivity and specificity using a bivariate random-effects model for each meta-analysis, with prevalence as a covariate. We estimated the overall effect of prevalence by pooling the effects using the inverse variance method. RESULTS Within a given review, a change in prevalence from the lowest to highest value resulted in a corresponding change in sensitivity or specificity from 0 to 40 percentage points. This effect was statistically significant (p < 0.05) for either sensitivity or specificity in 8 meta-analyses (35%). Overall, specificity tended to be lower with higher disease prevalence; there was no such systematic effect for sensitivity. INTERPRETATION The sensitivity and specificity of a test often vary with disease prevalence; this effect is likely to be the result of mechanisms, such as patient spectrum, that affect prevalence, sensitivity and specificity. Because it may be difficult to identify such mechanisms, clinicians should use prevalence as a guide when selecting studies that most closely match their situation.

Interannual variation in land-use intensity enhances grassland multidiversity

Although temporal heterogeneity is a well-accepted driver of biodiversity, effects of interannual variation in land-use intensity (LUI) have not been addressed yet. Additionally, responses to land use can differ greatly among different organisms; therefore, overall effects of land-use on total local biodiversity are hardly known. To test for effects of LUI (quantified as the combined intensity of fertilization, grazing, and mowing) and interannual variation in LUI (SD in LUI across time), we introduce a unique measure of whole-ecosystem biodiversity, multidiversity. This synthesizes individual diversity measures across up to 49 taxonomic groups of plants, animals, fungi, and bacteria from 150 grasslands. Multidiversity declined with increasing LUI among grasslands, particularly for rarer species and aboveground organisms, whereas common species and belowground groups were less sensitive. However, a high level of interannual variation in LUI increased overall multidiversity at low LUI and was even more beneficial for rarer species because it slowed the rate at which the multidiversity of rare species declined with increasing LUI. In more intensively managed grasslands, the diversity of rarer species was, on average, 18% of the maximum diversity across all grasslands when LUI was static over time but increased to 31% of the maximum when LUI changed maximally over time. In addition to decreasing overall LUI, we suggest varying LUI across years as a complementary strategy to promote biodiversity conservation.

A genome-to-genome analysis of associations between human genetic variation, HIV-1 sequence diversity, and viral control

HIV-1 sequence diversity is affected by selection pressures arising from host genomic factors. Using paired human and viral data from 1071 individuals, we ran >3000 genome-wide scans, testing for associations between host DNA polymorphisms, HIV-1 sequence variation and plasma viral load (VL), while considering human and viral population structure. We observed significant human SNP associations to a total of 48 HIV-1 amino acid variants (p<2.4 × 10−12). All associated SNPs mapped to the HLA class I region. Clinical relevance of host and pathogen variation was assessed using VL results. We identified two critical advantages to the use of viral variation for identifying host factors: (1) association signals are much stronger for HIV-1 sequence variants than VL, reflecting the ‘intermediate phenotype’ nature of viral variation; (2) association testing can be run without any clinical data. The proposed genome-to-genome approach highlights sites of genomic conflict and is a strategy generally applicable to studies of host–pathogen interaction.

Variation in fat mobilization during early lactation differently affects feed intake, body condition, and lipid and glucose metabolism in high-yielding dairy cows.

Fat mobilization to meet energy requirements during early lactation is inevitable because of insufficient feed intake, but differs greatly among high-yielding dairy cows. Therefore, we studied milk production, feed intake, and body condition as well as metabolic and endocrine changes in high-yielding dairy cows to identify variable strategies in metabolic and endocrine adaptation to overcome postpartum metabolic load attributable to milk production. Cows used in this study varied in fat mobilization around calving, as classified by mean total liver fat concentrations (LFC) postpartum. German Holstein cows (n=27) were studied from dry off until d 63 postpartum in their third lactation. All cows were fed the same total mixed rations ad libitum during the dry period and lactation. Plasma concentrations of metabolites and hormones were measured in blood samples taken at d 56, 28, 15, and 5 before expected calving and at d 1 and once weekly up to d 63 postpartum. Liver biopsies were taken on d 56 and 15 before calving, and on d 1, 14, 28, and 49 postpartum to measure LFC and glycogen concentrations. Cows were grouped accordingly to mean total LFC on d 1, 14, and 28 in high, medium, and low fat-mobilizing cows. Mean LFC (±SEM) differed among groups and were 351±14, 250±10, and 159±9 mg/g of dry matter for high, medium, and low fat-mobilizing cows, respectively, whereas hepatic glycogen concentrations postpartum were the highest in low fat-mobilizing cows. Cows in the low group showed the highest dry matter intake and the least negative energy balance postpartum, but energy-corrected milk yield was similar among groups. The decrease in body weight postpartum was greatest in high fat-mobilizing cows, but the decrease in backfat thickness was greatest in medium fat-mobilizing cows. Plasma concentrations of nonesterified fatty acids and β-hydroxybutyrate were highest around calving in high fat-mobilizing cows. Plasma triglycerides were highest in the medium group and plasma cholesterol concentrations were lowest in the high group at calving. During early lactation, the decrease in plasma glucose concentrations was greatest in the high group, and plasma insulin concentrations postpartum were highest in the low group. The revised quantitative insulin sensitivity check index values decreased during the transition period and postpartum, and were highest in the medium group. Plasma cortisol concentrations during the transition period and postpartum period and plasma leptin concentrations were highest in the medium group. In conclusion, cows adapted differently to the metabolic load and used variable strategies for homeorhetic regulation of milk production. Differences in fat mobilization were part of these strategies and contributed to the individual adaptation of energy metabolism to milk production.

Pelvic tilt is minimally changed by total hip arthroplasty

BACKGROUND While surgical navigation offers the opportunity to accurately place an acetabular component, questions remain as to the best goal for acetabular component positioning in individual patients. Overall functional orientation of the pelvis after surgery is one of the most important variables for the surgeon to consider when determining the proper goal for acetabular component orientation. QUESTIONS/PURPOSES We measured the variation in pelvic tilt in 30 patients before THA and the effect of THA on pelvic tilt in the same patients more than a year after THA. METHODS Each patient had a CT study for CT-based surgical navigation and standing and supine radiographs before and after surgery. Pelvic tilt was calculated for each of the radiographs using a novel and validated two-dimensional/three-dimensional matching technique. RESULTS Mean supine pelvic tilt changed less than 2°, from 4.4° ± 6.4° (range, -7.7° to 20.8°) before THA to 6.3° ± 6.6° (range, -5.7° to 19.6°) after THA. Mean standing pelvic tilt changed less than 1°, from 1.5° ± 7.2° (range, -13.1° to 12.8°) before THA to 2.0° ± 8.3° (range, -12.3° to 16.8°) after THA. Preoperative pelvic tilt correlated with postoperative tilt in both the supine (r(2) = 0.75) and standing (r(2) = 0.87) positions. CONCLUSIONS In this population, pelvic tilt had a small and predictable change after surgery. However, intersubject variability of pelvic tilt was high, suggesting preoperative pelvic tilt should be considered when determining desired acetabular component positioning on a patient-specific basis.

“Sweating meteorites”-Water-soluble salts and temperature variation in ordinary chondrites and soil from the hot desert of Oman

The common appearance of hygroscopic brine (“sweating”) on ordinary chondrites (OCs) from Oman during storage under room conditions initiated a study on the role of water-soluble salts on the weathering of OCs. Analyses of leachates from OCs and soils, combined with petrography of alteration features and a 11-month record of in situ meteorite and soil temperatures, are used to evaluate the role of salts in OC weathering. Main soluble ions in soils are Ca2+, SO42−, HCO3−, Na+, and Cl−, while OC leachates are dominated by Mg2+ (from meteoritic olivine), Ca2+ (from soil), Cl− (from soil), SO42− (from meteoritic troilite and soil), and iron (meteoritic). “Sweating meteorites” mainly contain Mg2+ and Cl−. The median Na/Cl mass ratio of leachates changes from 0.65 in soils to 0.07 in meteorites, indicating the precipitation of a Na-rich phase or loss of an efflorescent Na-salt. The total concentrations of water-soluble ions in bulk OCs ranges from 600 to 9000 μg g−1 (median 2500 μg g−1) as compared to 187–14140 μg g−1 in soils (median 1148 μg g−1). Soil salts dissolved by rain water are soaked up by meteorites by capillary forces. Daily heating (up to 66.3 °C) and cooling of the meteorites cause a pumping effect, resulting in a strong concentration of soluble ions in meteorites over time. The concentrations of water-soluble ions in meteorites, which are complex mixtures of ions from the soil and from oxidation and hydrolysis of meteoritic material, depend on the degree of weathering and are highest at W3. Input of soil contaminants generally dominates over the ions mobilized from meteorites. Silicate hydrolysis preferentially affects olivine and is enhanced by sulfide oxidation, producing local acidic conditions as evidenced by jarosite. Plagioclase weathering is negligible. After completion of troilite oxidation, the rate of chemical weathering slows down with continuing Ca-sulfate contamination.