Delineating environmental envelopes to improve mapping of species distributions, via a hurdle model with CART &/or MaxEnT

Species distribution modelling (SDM) typically analyses species’ presence together with some form of absence information. Ideally absences comprise observations or are inferred from comprehensive sampling. When such information is not available, then pseudo-absences are often generated from the background locations within the study region of interest containing the presences, or else absence is implied through the comparison of presences to the whole study region, e.g. as is the case in Maximum Entropy (MaxEnt) or Poisson point process modelling. However, the choice of which absence information to include can be both challenging and highly influential on SDM predictions (e.g. Oksanen and Minchin, 2002). In practice, the use of pseudo- or implied absences often leads to an imbalance where absences far outnumber presences. This leaves analysis highly susceptible to ‘naughty-noughts’: absences that occur beyond the envelope of the species, which can exert strong influence on the model and its predictions (Austin and Meyers, 1996). Also known as ‘excess zeros’, naughty noughts can be estimated via an overall proportion in simple hurdle or mixture models (Martin et al., 2005). However, absences, especially those that occur beyond the species envelope, can often be more diverse than presences. Here we consider an extension to excess zero models. The two-staged approach first exploits the compartmentalisation provided by classification trees (CTs) (as in O’Leary, 2008) to identify multiple sources of naughty noughts and simultaneously delineate several species envelopes. Then SDMs can be fit separately within each envelope, and for this stage, we examine both CTs (as in Falk et al., 2014) and the popular MaxEnt (Elith et al., 2006). We introduce a wider range of model performance measures to improve treatment of naughty noughts in SDM. We retain an overall measure of model performance, the area under the curve (AUC) of the Receiver-Operating Curve (ROC), but focus on its constituent measures of false negative rate (FNR) and false positive rate (FPR), and how these relate to the threshold in the predicted probability of presence that delimits predicted presence from absence. We also propose error rates more relevant to users of predictions: false omission rate (FOR), the chance that a predicted absence corresponds to (and hence wastes) an observed presence, and the false discovery rate (FDR), reflecting those predicted (or potential) presences that correspond to absence. A high FDR may be desirable since it could help target future search efforts, whereas zero or low FOR is desirable since it indicates none of the (often valuable) presences have been ignored in the SDM. For illustration, we chose Bradypus variegatus, a species that has previously been published as an exemplar species for MaxEnt, proposed by Phillips et al. (2006). We used CTs to increasingly refine the species envelope, starting with the whole study region (E0), eliminating more and more potential naughty noughts (E1–E3). When combined with an SDM fit within the species envelope, the best CT SDM had similar AUC and FPR to the best MaxEnt SDM, but otherwise performed better. The FNR and FOR were greatly reduced, suggesting that CTs handle absences better. Interestingly, MaxEnt predictions showed low discriminatory performance, with the most common predicted probability of presence being in the same range (0.00-0.20) for both true absences and presences. In summary, this example shows that SDMs can be improved by introducing an initial hurdle to identify naughty noughts and partition the envelope before applying SDMs. This improvement was barely detectable via AUC and FPR yet visible in FOR, FNR, and the comparison of predicted probability of presence distribution for pres/absence.

Ontogenetic depth partitioning by juvenile freshwater sawfish (Pristis microdon: Pristidae) in a riverine environment

The freshwater sawfish (Pristis microdon) is a critically endangered elasmobranch. Ontogenetic changes in the habitat use of juvenile P. microdon were studied using acoustic tracking in the Fitzroy River, Western Australia. Habitat partitioning was significant between 0+ (2007 year class) and larger 1+ (2006 year class) P. microdon. Smaller 0+ fish generally occupied shallower water (<0.6 m) compared with 1+ individuals, which mainly occurred in depths >0.6 m. Significant differences in hourly depth use were also revealed. The depth that 1+ P. microdon occupied was significantly influenced by lunar phase with these animals utilising a shallower and narrower depth range during the full moon compared with the new moon. This was not observed in 0+ individuals. Habitat partitioning was likely to be related to predator avoidance, foraging behaviours, and temperature and/or light regimes. The occurrence of 1+ P. microdon in deeper water may also result from a need for greater depths in which to manoeuvre. The present study demonstrates the utility of acoustic telemetry in monitoring P. microdon in a riverine environment. These results demonstrate the need to consider the habitat requirements of different P. microdon cohorts in the strategic planning of natural resources and will aid in the development of management strategies for this species.

Determine baseline dose responses of Alphitobius diaperinus to gamma cyhalothrin

Management of insecticide resistance.

Intracellular effects of atmospheric-pressure plasmas on melanoma cancer cells

Gas discharge plasmas formed at atmospheric pressure and near room temperature have recently been shown as a promising tool for cancer treatment. The mechanism of the plasma action is attributed to generation of reactive oxygen and nitrogen species, electric fields, charges, and photons. The relative importance of different modes of action of atmospheric-pressure plasmas depends on the process parameters and specific treatment objects. Hence, an in-depth understanding of biological mechanisms that underpin plasma-induced death in cancer cells is required to optimise plasma processing conditions. Here, the intracellular factors involved in the observed anti-cancer activity in melanoma Mel007 cells are studied, focusing on the effect of the plasma treatment dose on the expression of tumour suppressor protein TP73. Over-expression of TP73 causes cell growth arrest and/or apoptosis, and hence can potentially be targeted to enhance killing efficacy and selectivity of the plasma treatment. It is shown that the plasma treatment induces dose-dependent up-regulation of TP73 gene expression, resulting in significantly elevated levels of TP73 RNA and protein in plasma-treated melanoma cells. Silencing of TP73 expression by means of RNA interference inhibited the anticancer effects of the plasma, similar to the effect of caspase inhibitor z-VAD or ROS scavenger N-acetyl cysteine. These results confirm the role of TP73 protein in dose-dependent regulation of anticancer activity of atmospheric-pressure plasmas.

Distribution of ocular biometry in 7- and 14-year-old Chinese children

Purpose: To describe distributions of ocular biometry and their associations with refraction in 7- and 14-year-old children in urban areas of Anyang, central China. Methods: A total of 2271 grade 1 students aged 7.1 ± 0.4 years and 1786 grade 8 students aged 13.7 ± 0.5 years were measured with ocular biometry and cycloplegic refraction. A parental myopia questionnaire was administered to parents. Results: Mean axial length, anterior chamber depth, lens thickness, central corneal thickness, corneal diameter, corneal radius of curvature, axial length/corneal radius of curvature ratio, and spherical equivalent refraction were 22.72 ± 0.76 mm, 2.89 ± 0.24 mm, 3.61 ± 0.19 mm, 540.5 ± 31 μm, 12.06 ± 0.44 mm, 7.80 ± 0.25 mm, 2.91 ± 0.08, and +0.95 ± 1.05 diopters (D), respectively, in 7-year-old children. They were 24.39 ± 1.13 mm, 3.42 ± 0.41 mm, 3.18 ± 0.24 mm, 548.9 ± 33 μm, 12.03 ± 0.43 mm, 7.80 ± 0.26 mm, 3.13 ± 0.14, and −2.06 ± 2.20 D, respectively, in 14-year-old children. Compared with 7-year-old children, the older group had significantly more myopia (−3.0 D), longer axial length (1.7 mm), deeper anterior chamber depth (0.3 mm), thinner lens thickness (−0.2 mm), thicker central corneal thickness (10 μm), and greater axial length/corneal radius of curvature ratio (0.22) (all p < 0.001), as well as smaller corneal diameter (−0.03 mm, p = 0.02) and similar corneal radius of curvature. Sex differences were similar in both age groups, with boys having longer axial length (0.5 mm), deeper anterior chamber depth (0.1 mm), shorter lens thickness (0.03 mm), greater central corneal thickness (5 μm), greater corneal diameter (0.15 mm), and greater corneal radius of curvature (0.14 mm) than girls (all p < 0.01). The most important variables related to spherical equivalent refraction were vitreous length, corneal radius of curvature, and lens thickness. Conclusions: The 14-year-old group had larger parameter dimensions than the 7-year-old group except for corneal radius of curvature (unchanged) and lens thickness and corneal diameter (both smaller). Boys had large parameter dimensions than girls except for lens thickness (smaller). Axial length, corneal radius of curvature, and lens thickness were the most important determinants of refraction.

Design and assessment of weirs for fish passage under drowned conditions

In aquatic systems, in-stream structures such as dams, weirs and road crossings can act as barriers to fish movement along waterways. There is a growing array of technological fish-pass solutions for the movement of fish across large structures such as weirs and dams. However, most existing weir structures lack dedicated fishways, and fish often have to rely on drowned conditions to move upstream. In order to assess the adequacy of a given or proposed weir for upstream fish passage under drowned conditions, it is necessary to determine, firstly, the hydraulic properties of the drowned weir with respect to the requirements of the fish community and, secondly, the duration and timing of drowning flows with respect to the hydrograph for the site and the likely timing of fish movements. This paper primarily addresses the first issue. A computer program has been developed and incorporated in a simple-to-operate spreadsheet for the determination of the hydraulic characteristics of a drowned weir which are important to fish movement. The program is based on a theoretical analysis of drowned weirs and subsequent extensive verification in laboratory experiments. Inputs to the program include site information comprising channel cross-section data, channel slope, and channel roughness, and weir information comprising weir height and the required minimum drowned depth over the weir for migrating fish passage. The program then calculates the flow rate at which the required level of drowning occurs, the velocity characteristics above the weir (including transverse distributions), and flow depths and velocities upstream and downstream of the weir. The paper discusses (briefly) the theoretical background of the program and its experimental verification. A case study is then presented that illustrates the use of the program in the field to assess fish passage opportunities at an existing weir and to develop a case for retrofitting a fishway. Some discussion is also provided on the contribution of a modelled drownout volume to the assessment of how significant a barrier a weir is to fish passage. It is shown that the program is an important new additional tool in the assessment of the adequacy of weir structures in providing for fish movement and informing associated fish passage solutions. (C) 2011 Elsevier B.V. All rights reserved.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.

Effect of Increasing Low Dose Simplexin Exposure in Cattle consuming Pimelia trychostachya

Pimelea species (or desert riceflower) are small native plants endemic to the drier inland pastoral regions of Australia, which cause a unique syndrome in grazing cattle characterised by submandibular oedema and oedema in the brisket area as a result of right-sided heart failure attributed to the toxin simplexin. Field evidence suggests that poisoning can occur through minor, inadvertent consumption of Pimelea plant material, but the minimum simplexin intake required to induce Pimelea poisoning is not known. In this study, mild Pimelea poisoning was induced at a daily dose of 12.5 mg Pimelea/kg bodyweight per day, equivalent to 2.5 µg simplexin/kg bodyweight per day, demonstrating the high potential toxicity of these plant species. Effects in all animals diminished with prolonged low dose feeding and we postulate that these animals developed mechanisms for detoxifying simplexin, 1, possibly through rumen bacteria adaptation or activation of liver enzymes.

Dasatinib, nilotinib and standard-dose imatinib for the first line treatment of chronic myeloid leukaemia: Systematic reviews and economic analyses

- Background Nilotinib and dasatinib are now being considered as alternative treatments to imatinib as a first-line treatment of chronic myeloid leukaemia (CML). - Objective This technology assessment reviews the available evidence for the clinical effectiveness and cost-effectiveness of dasatinib, nilotinib and standard-dose imatinib for the first-line treatment of Philadelphia chromosome-positive CML. - Data sources Databases [including MEDLINE (Ovid), EMBASE, Current Controlled Trials, ClinicalTrials.gov, the US Food and Drug Administration website and the European Medicines Agency website] were searched from search end date of the last technology appraisal report on this topic in October 2002 to September 2011. - Review methods A systematic review of clinical effectiveness and cost-effectiveness studies; a review of surrogate relationships with survival; a review and critique of manufacturer submissions; and a model-based economic analysis. - Results Two clinical trials (dasatinib vs imatinib and nilotinib vs imatinib) were included in the effectiveness review. Survival was not significantly different for dasatinib or nilotinib compared with imatinib with the 24-month follow-up data available. The rates of complete cytogenetic response (CCyR) and major molecular response (MMR) were higher for patients receiving dasatinib than for those with imatinib for 12 months' follow-up (CCyR 83% vs 72%, p < 0.001; MMR 46% vs 28%, p < 0.0001). The rates of CCyR and MMR were higher for patients receiving nilotinib than for those receiving imatinib for 12 months' follow-up (CCyR 80% vs 65%, p < 0.001; MMR 44% vs 22%, p < 0.0001). An indirect comparison analysis showed no difference between dasatinib and nilotinib for CCyR or MMR rates for 12 months' follow-up (CCyR, odds ratio 1.09, 95% CI 0.61 to 1.92; MMR, odds ratio 1.28, 95% CI 0.77 to 2.16). There is observational association evidence from imatinib studies supporting the use of CCyR and MMR at 12 months as surrogates for overall all-cause survival and progression-free survival in patients with CML in chronic phase. In the cost-effectiveness modelling scenario, analyses were provided to reflect the extensive structural uncertainty and different approaches to estimating OS. First-line dasatinib is predicted to provide very poor value for money compared with first-line imatinib, with deterministic incremental cost-effectiveness ratios (ICERs) of between £256,000 and £450,000 per quality-adjusted life-year (QALY). Conversely, first-line nilotinib provided favourable ICERs at the willingness-to-pay threshold of £20,000-30,000 per QALY. - Limitations Immaturity of empirical trial data relative to life expectancy, forcing either reliance on surrogate relationships or cumulative survival/treatment duration assumptions. - Conclusions From the two trials available, dasatinib and nilotinib have a statistically significant advantage compared with imatinib as measured by MMR or CCyR. Taking into account the treatment pathways for patients with CML, i.e. assuming the use of second-line nilotinib, first-line nilotinib appears to be more cost-effective than first-line imatinib. Dasatinib was not cost-effective if decision thresholds of £20,000 per QALY or £30,000 per QALY were used, compared with imatinib and nilotinib. Uncertainty in the cost-effectiveness analysis would be substantially reduced with better and more UK-specific data on the incidence and cost of stem cell transplantation in patients with chronic CML. - Funding The Health Technology Assessment Programme of the National Institute for Health Research.

Focal brachytherapy treatment planning using multi-parametric MRI and biological dose optimisation

We describe a novel approach to treatment planning for focal brachytherapy utilizing a biologically based inverse optimization algorithm and biological imaging to target an ablative dose at known regions of significant tumour burden and a lower, therapeutic dose to low risk regions.

A new family of multivariate heavy-tailed distributions with variable marginal amounts of tailweight: Application to robust clustering

We propose a family of multivariate heavy-tailed distributions that allow variable marginal amounts of tailweight. The originality comes from introducing multidimensional instead of univariate scale variables for the mixture of scaled Gaussian family of distributions. In contrast to most existing approaches, the derived distributions can account for a variety of shapes and have a simple tractable form with a closed-form probability density function whatever the dimension. We examine a number of properties of these distributions and illustrate them in the particular case of Pearson type VII and t tails. For these latter cases, we provide maximum likelihood estimation of the parameters and illustrate their modelling flexibility on simulated and real data clustering examples.

Bayesian mixture model estimation of aerosol particle size distributions

In this paper, we examine approaches to estimate a Bayesian mixture model at both single and multiple time points for a sample of actual and simulated aerosol particle size distribution (PSD) data. For estimation of a mixture model at a single time point, we use Reversible Jump Markov Chain Monte Carlo (RJMCMC) to estimate mixture model parameters including the number of components which is assumed to be unknown. We compare the results of this approach to a commonly used estimation method in the aerosol physics literature. As PSD data is often measured over time, often at small time intervals, we also examine the use of an informative prior for estimation of the mixture parameters which takes into account the correlated nature of the parameters. The Bayesian mixture model offers a promising approach, providing advantages both in estimation and inference.

Optical and Physical Properties of Atmospheric Aerosols over the Bay of Bengal during ICARB

Simultaneous and collocated measurements of total and hemispherical backscattering coefficients (σ and β, respectively) at three wavelengths, mass size distributions, and columnar spectral aerosol optical depth (AOD) were made onboard an extensive cruise experiment covering, for the first time, the entire Bay of Bengal (BoB) and northern Indian Ocean. The results are synthesized to understand the optical properties of aerosols in the marine atmospheric boundary layer and their dependence on the size distribution. The observations revealed distinct spatial and spectral variations of all the aerosol parameters over the BoB and the presence of strong latitudinal gradients. The size distributions varied spatially, with the majority of accumulation modes decreasing from north to south. The scattering coefficient decreased from very high values (resembling those reported for continental/urban locations) in the northern BoB to very low values seen over near-pristine environments in the southeastern BoB. The average mass scattering efficiency of BoB aerosols was found to be 2.66 ± 0.1 m2 g−1 at 550 nm. The spectral dependence of columnar AOD deviated significantly from that of the scattering coefficients in the northern BoB, implying vertical heterogeneity in the aerosol type in that region. However, a more homogeneous scenario was observed in the southern BoB. Simultaneous lidar and in situ measurements onboard an aircraft over the ocean revealed the presence of elevated aerosol layers of enhanced extinction at altitudes of 1 to 3 km with an offshore extent of a few hundred kilometers. Back-trajectory analyses showed these layers to be associated with advection from west Asia and western India. The large spatial variations and vertical heterogeneity in aerosol properties, revealed by the present study, need to be included in the regional radiative forcing over the Bay of Bengal.

Low-dose methyl bromide fumigation as a quarantine disinfestation treatment for nectarines against Queensland fruit fly (Diptera:Tephritidae)

White nectarines (Prunus persica var. nucipersica) were fumigated with methyl bromide (MB) at a nominal treatment dose of 18 g m-3 at 18°C for 5 h and 30 min as a quarantine disinfestation treatment against Bactrocera tryoni, the Queensland fruit fly. Three large scale trials were conducted against each of the four immature lifestages, eggs and first, second and third instars. There were no survivors from the estimated 43,614 eggs, 41,873 first instars, 41,345 second instars and 33,549 third instars treated, thereby resulting in an efficacy of GROTERDAN99.99% mortality at the 95% confidence level for each lifestage. Of the 12 trials reported herein, the highest concentration of MB, sampled from the chamber headspace analysed by gas chromatography, was 18.7 g m-3. The maximum chamber temperature from 5 min readings was 19.7°C and the maximum fruit core temperature was 19.5°C. The treatment time for all trials was exactly 5.5 h. Thus the recommended treatment dose to disinfest nectarines from B. tryoni is 19.0 g m-3 MB at 20.0°C for 5.5 h. Fruit quality trials were conducted on white nectarines at three combinations of treatment parameters: 15 g m-3 MB at 19°C for 5.25 h; 18 g m-3 MB at 19°C for 5.5 h and 21 g m-3 MB at 19°C for 5.5 h. The fruit were stored at 0, 4 and 8 days at 4°C and 8 days at 4°C followed by 4 d at 22°C. They were then were assessed for skin colour, flesh colour, skin defects, flesh defects, fruit weight loss, flesh firmness, total soluble solids, titratable acidity and rots. There was no significant difference between untreated control and MB treated fruits in any of the parameters measured. Thus the treatments did not have adverse effects on fruit quality.