abctools: an R package for tuning approximate Bayesian computation analyses

Approximate Bayesian computation (ABC) is a popular family of algorithms which perform approximate parameter inference when numerical evaluation of the likelihood function is not possible but data can be simulated from the model. They return a sample of parameter values which produce simulations close to the observed dataset. A standard approach is to reduce the simulated and observed datasets to vectors of summary statistics and accept when the difference between these is below a specified threshold. ABC can also be adapted to perform model choice. In this article, we present a new software package for R, abctools which provides methods for tuning ABC algorithms. This includes recent dimension reduction algorithms to tune the choice of summary statistics, and coverage methods to tune the choice of threshold. We provide several illustrations of these routines on applications taken from the ABC literature.

Rapid selective separation of americium/curium from simulated nuclear forensic matrices using triazine ligands

In analysis of complex nuclear forensic samples containing lanthanides, actinides and matrix elements, rapid selective extraction of Am/Cm for quantification is challenging, in particular due the difficult separation of Am/Cm from lanthanides. Here we present a separation process for Am/Cm(III) which is achieved using a combination of AG1-X8 chromatography followed by Am/Cm extraction with a triazine ligand. The ligands tested in our process were CyMe4-BTPhen, CyMe4- BTBP, CA-BTP and CA-BTPhen. Our process allows for purification and quantification of Am and Cm (recoveries 80%–100%) and other major actinides in < 2d without the use of multiple columns or thiocyanate. The process is unaffected by high level Ca(II)/Fe(III)/Al(III) (10mg mL−1) and thus requires little pre-treatment of samples.

Proper evaluation of the external convective heat transfer for the thermal analysis of cool roofs

Cool materials are characterized by high solar reflectance and high thermal emittance; when applied to the external surface of a roof, they make it possible to limit the amount of solar irradiance absorbed by the roof, and to increase the rate of heat flux emitted by irradiation to the environment, especially during nighttime. However, a roof also releases heat by convection on its external surface; this mechanism is not negligible, and an incorrect evaluation of its entity might introduce significant inaccuracy in the assessment of the thermal performance of a cool roof, in terms of surface temperature and rate of heat flux transferred to the indoors. This issue is particularly relevant in numerical simulations, which are essential in the design stage, therefore it deserves adequate attention. In the present paper, a review of the most common algorithms used for the calculation of the convective heat transfer coefficient due to wind on horizontal building surfaces is presented. Then, with reference to a case study in Italy, the simulated results are compared to the outcomes of a measurement campaign. Hence, the most appropriate algorithms for the convective coefficient are identified, and the errors deriving by an incorrect selection of this coefficient are discussed.

The effect of simulated flooding post-anthesis on pre-harvest sprouting and subsequent seed longevity in contrasting rice cultivars

Unpredictable flooding is a major constraint to rice production. It can occur at any growth stage. The effect of simulated flooding post-anthesis on yield and subsequent seed quality of pot-grown rice (Oryza sativa L.) plants was investigated in glasshouses and controlled-environment growth cabinets. Submergence post-anthesis (9-40 DAA) for 3 or 5 days reduced seed weight of japonica rice cv. Gleva, with considerable pre-harvest sprouting (up to 53%). The latter was greater the later in seed development and maturation that flooding occurred. Sprouted seed had poor ability to survive desiccation or germinate normally upon rehydration, whereas the effects of flooding on the subsequent air-dry seed storage longevity (p50) of the non-sprouted seed fraction was negligible. The indica rice cvs IR64 and IR64Sub1 (introgression of submergence tolerance gene Submergence1A-1) were both far more tolerant to flooding post-anthesis than cv. Gleva: four days’ submergence of these two near-isogenic cultivars at 10-40 DAA resulted less than 1% sprouted seeds. The presence of the Sub1A-1 allele in cv. IR64Sub1 was verified by gel electrophoresis and DNA sequencing. It had no harmful effect on loss in seed viability during storage compared with IR64 in both control and flooded environments. Moreover, the germinability and changes in dormancy during seed development and maturation were very similar to IR64. The efficiency of using chemical spray to increase seed dormancy was investigated in the pre-harvest sprouting susceptible rice cv. Gleva. Foliar application of molybdenum at 100 mg L-1 reduced sprouted seeds by 15-21% following 4 days’ submergence at 20-30 DAA. Analyses confirmed that the treatment did result in molybdenum uptake by the plants, and also tended to increase seed abscisic acid concentration. The latter was reduced by submergence and declined exponentially during grain ripening. The selection of submergence-tolerant varieties was more successful than application of molybdenum in reducing pre-harvest sprouting.

Rainfastness of poly(vinyl alcohol) deposits on Vicia faba leaf surfaces: from laboratory-scale washing to simulated rain

Rainfastness is the ability of agrochemical deposits to resist wash-off by rain and other related environmental phenomena. This work reports laboratory-scale and raintower studies of the rainfastness of fluorescently labeled poly(vinyl alcohol) (PVA) using fluorescent microscopy combined with image analysis. Samples of hydrolyzed PVA exhibit improved rainfastness over a threshold molecular weight, which correlates with PVA film dissolution, swelling, and crystalline properties. It was also established that the rainfastness of PVA scaled with the molecular weight over this threshold. These PVA samples were further characterized in order to determine the effect of the crystallinity on rainfastness. The quantification of rainfastness is of great interest to the field of agrochemical formulation development in order to improve the efficacy of pesticides and their adjuvants.

Zero attracting recursive least squares algorithms

The l1-norm sparsity constraint is a widely used technique for constructing sparse models. In this contribution, two zero-attracting recursive least squares algorithms, referred to as ZA-RLS-I and ZA-RLS-II, are derived by employing the l1-norm of parameter vector constraint to facilitate the model sparsity. In order to achieve a closed-form solution, the l1-norm of the parameter vector is approximated by an adaptively weighted l2-norm, in which the weighting factors are set as the inversion of the associated l1-norm of parameter estimates that are readily available in the adaptive learning environment. ZA-RLS-II is computationally more efficient than ZA-RLS-I by exploiting the known results from linear algebra as well as the sparsity of the system. The proposed algorithms are proven to converge, and adaptive sparse channel estimation is used to demonstrate the effectiveness of the proposed approach.

Effects of rain shelter or simulated rain during grain filling and maturation on subsequent wheat grain quality in the UK

The effects of simulated additional rain (ear wetting, 25 mm) or of rain shelter imposed at different periods after anthesis on grain quality at maturity and the dynamics of grain filling and desiccation were investigated in UK field-grown crops of wheat (Triticum aestivum L., cvar Tybalt) in 2011 and in 2012 when June–August rainfall was 255.0 and 214.6 mm, respectively, and above the decadal mean (157.4 mm). Grain filling and desiccation were quantified well by broken-stick regressions and Gompertz curves, respectively. Rain shelter for 56 (2011) or 70 d (2012) after anthesis, and to a lesser extent during late maturation only, resulted in more rapid desiccation and hence progress to harvest maturity whereas ear wetting had negligible effects, even when applied four times. Grain-filling duration was also affected as above in 2011, but with no significant effect in 2012. In both years, there were strong positive associations between final grain dry weight and duration of filling. The treatments affected all grain quality traits in 2011: nitrogen (N) and sulphur (S) concentrations, N:S ratio, sodium dodecyl sulphate (SDS) sedimentation volume, Hagberg Falling Number (HFN), and the incidence of blackpoint. Only N concentration and blackpoint were affected significantly by treatments in 2012. Rain shelter throughout grain filling reduced N concentration, whereas rain shelter reduced the incidence of blackpoint and ear wetting increased it. In 2011, rain shelter throughout reduced S concentration, increased N:S ratio and reduced SDS. Treatment effects on HFN were not consistent within or between years. Nevertheless, a comparison between the extreme treatment means in 2012 indicated damage from late rain combined with ear wetting resulted in a reduction of c. 0.7 s in HFN/mm August rainfall, whilst that between samples taken immediately after ear wetting at harvest maturity or 7 d later suggested recovery from damage to HFN upon re-drying in planta. Hence, the incidence of blackpoint was the only grain quality trait affected consistently by the diverse treatments. The remaining aspects of grain quality were comparatively resilient to rain incident upon developing and maturing ears of cvar Tybalt. No consistent temporal patterns of sensitivity to shelter or ear wetting were detected for any aspect of grain quality.