Prediction of mortality and conception rates of beef breeding cattle in northern Australia

In current simulation packages for the management of extensive beef-cattle enterprises, the relationships for the key biological rates (namely conception and mortality) are quite rudimentary. To better estimate these relationships, cohort-level data covering 17 100 cow-years from six sites across northern Australia were collated and analysed. Further validation data, from 7200 cow-years, were then used to test these relationships. Analytical problems included incomplete and non-standardised data, considerable levels of correlation among the 'independent' variables, and the close similarity of alternate possible models. In addition to formal statistical analyses of these data, the theoretical equations for predicting mortality and conception rates in the current simulation models were reviewed, and then reparameterised and recalibrated where appropriate. The final models explained up to 80% of the variation in the data. These are now proposed as more accurate and useful models to be used in the prediction of biological rates in simulation studies for northern Australia. © The State of Queensland (through the Department of Agriculture, Fisheries and Forestry) 2012. © CSIRO.

Inferential, non-parametric statistics to assess quality of probabilistic forecast systems

Many statistical forecast systems are available to interested users. In order to be useful for decision-making, these systems must be based on evidence of underlying mechanisms. Once causal connections between the mechanism and their statistical manifestation have been firmly established, the forecasts must also provide some quantitative evidence of `quality’. However, the quality of statistical climate forecast systems (forecast quality) is an ill-defined and frequently misunderstood property. Often, providers and users of such forecast systems are unclear about what ‘quality’ entails and how to measure it, leading to confusion and misinformation. Here we present a generic framework to quantify aspects of forecast quality using an inferential approach to calculate nominal significance levels (p-values) that can be obtained either by directly applying non-parametric statistical tests such as Kruskal-Wallis (KW) or Kolmogorov-Smirnov (KS) or by using Monte-Carlo methods (in the case of forecast skill scores). Once converted to p-values, these forecast quality measures provide a means to objectively evaluate and compare temporal and spatial patterns of forecast quality across datasets and forecast systems. Our analysis demonstrates the importance of providing p-values rather than adopting some arbitrarily chosen significance levels such as p < 0.05 or p < 0.01, which is still common practice. This is illustrated by applying non-parametric tests (such as KW and KS) and skill scoring methods (LEPS and RPSS) to the 5-phase Southern Oscillation Index classification system using historical rainfall data from Australia, The Republic of South Africa and India. The selection of quality measures is solely based on their common use and does not constitute endorsement. We found that non-parametric statistical tests can be adequate proxies for skill measures such as LEPS or RPSS. The framework can be implemented anywhere, regardless of dataset, forecast system or quality measure. Eventually such inferential evidence should be complimented by descriptive statistical methods in order to fully assist in operational risk management.

A new data source for fisheries resource assessment: genetic estimates of the effective number of spawners. Final Report to the Fisheries Research and Development Corporation.

The development of innovative methods of stock assessment is a priority for State and Commonwealth fisheries agencies. It is driven by the need to facilitate sustainable exploitation of naturally occurring fisheries resources for the current and future economic, social and environmental well being of Australia. This project was initiated in this context and took advantage of considerable recent achievements in genomics that are shaping our comprehension of the DNA of humans and animals. The basic idea behind this project was that genetic estimates of effective population size, which can be made from empirical measurements of genetic drift, were equivalent to estimates of the successful number of spawners that is an important parameter in process of fisheries stock assessment. The broad objectives of this study were to 1. Critically evaluate a variety of mathematical methods of calculating effective spawner numbers (Ne) by a. conducting comprehensive computer simulations, and by b. analysis of empirical data collected from the Moreton Bay population of tiger prawns (P. esculentus). 2. Lay the groundwork for the application of the technology in the northern prawn fishery (NPF). 3. Produce software for the calculation of Ne, and to make it widely available. The project pulled together a range of mathematical models for estimating current effective population size from diverse sources. Some of them had been recently implemented with the latest statistical methods (eg. Bayesian framework Berthier, Beaumont et al. 2002), while others had lower profiles (eg. Pudovkin, Zaykin et al. 1996; Rousset and Raymond 1995). Computer code and later software with a user-friendly interface (NeEstimator) was produced to implement the methods. This was used as a basis for simulation experiments to evaluate the performance of the methods with an individual-based model of a prawn population. Following the guidelines suggested by computer simulations, the tiger prawn population in Moreton Bay (south-east Queensland) was sampled for genetic analysis with eight microsatellite loci in three successive spring spawning seasons in 2001, 2002 and 2003. As predicted by the simulations, the estimates had non-infinite upper confidence limits, which is a major achievement for the application of the method to a naturally-occurring, short generation, highly fecund invertebrate species. The genetic estimate of the number of successful spawners was around 1000 individuals in two consecutive years. This contrasts with about 500,000 prawns participating in spawning. It is not possible to distinguish successful from non-successful spawners so we suggest a high level of protection for the entire spawning population. We interpret the difference in numbers between successful and non-successful spawners as a large variation in the number of offspring per family that survive – a large number of families have no surviving offspring, while a few have a large number. We explored various ways in which Ne can be useful in fisheries management. It can be a surrogate for spawning population size, assuming the ratio between Ne and spawning population size has been previously calculated for that species. Alternatively, it can be a surrogate for recruitment, again assuming that the ratio between Ne and recruitment has been previously determined. The number of species that can be analysed in this way, however, is likely to be small because of species-specific life history requirements that need to be satisfied for accuracy. The most universal approach would be to integrate Ne with spawning stock-recruitment models, so that these models are more accurate when applied to fisheries populations. A pathway to achieve this was established in this project, which we predict will significantly improve fisheries sustainability in the future. Regardless of the success of integrating Ne into spawning stock-recruitment models, Ne could be used as a fisheries monitoring tool. Declines in spawning stock size or increases in natural or harvest mortality would be reflected by a decline in Ne. This would be good for data-poor fisheries and provides fishery independent information, however, we suggest a species-by-species approach. Some species may be too numerous or experiencing too much migration for the method to work. During the project two important theoretical studies of the simultaneous estimation of effective population size and migration were published (Vitalis and Couvet 2001b; Wang and Whitlock 2003). These methods, combined with collection of preliminary genetic data from the tiger prawn population in southern Gulf of Carpentaria population and a computer simulation study that evaluated the effect of differing reproductive strategies on genetic estimates, suggest that this technology could make an important contribution to the stock assessment process in the northern prawn fishery (NPF). Advances in the genomics world are rapid and already a cheaper, more reliable substitute for microsatellite loci in this technology is available. Digital data from single nucleotide polymorphisms (SNPs) are likely to super cede ‘analogue’ microsatellite data, making it cheaper and easier to apply the method to species with large population sizes.

Non-destructive assessment of Arundo donax (Poaceae) leaf quality

Leaf carbon (C) content, leaf nitrogen (N) content, and C:N ratio are especially useful for understanding plant-herbivore interactions and may be important in developing control methods for the invasive riparian plant Arundo donax L. We measured C content, N content, C:N ratio, and chlorophyll index (SPAD 502 reading) for 768 leaves from A. donax collected over a five year period at several locations in California, Nevada, and Texas. Leaf N was more variable than leaf C, and thus we developed a linear regression equation for estimating A. donax leaf N from the leaf chlorophyll index (SPAD reading). When applied to two independent data sets, the equation (leaf N content % = -0.63 + 0.08 x SPAD) produced realistic estimates that matched seasonal and spatial trends reported from a natural A. donax population. Used in conjunction with the handheld SPAD 502 meter, the equation provides a rapid, non-destructive method for estimating A. donax leaf quality.

Annual forecasting of the Australian macadamia crop - integrating tree census data with statistical climate-adjustment models

To facilitate marketing and export, the Australian macadamia industry requires accurate crop forecasts. Each year, two levels of crop predictions are produced for this industry. The first is an overall longer-term forecast based on tree census data of growers in the Australian Macadamia Society (AMS). This data set currently accounts for around 70% of total production, and is supplemented by our best estimates of non-AMS orchards. Given these total tree numbers, average yields per tree are needed to complete the long-term forecasts. Yields from regional variety trials were initially used, but were found to be consistently higher than the average yields that growers were obtaining. Hence, a statistical model was developed using growers' historical yields, also taken from the AMS database. This model accounted for the effects of tree age, variety, year, region and tree spacing, and explained 65% of the total variation in the yield per tree data. The second level of crop prediction is an annual climate adjustment of these overall long-term estimates, taking into account the expected effects on production of the previous year's climate. This adjustment is based on relative historical yields, measured as the percentage deviance between expected and actual production. The dominant climatic variables are observed temperature, evaporation, solar radiation and modelled water stress. Initially, a number of alternate statistical models showed good agreement within the historical data, with jack-knife cross-validation R2 values of 96% or better. However, forecasts varied quite widely between these alternate models. Exploratory multivariate analyses and nearest-neighbour methods were used to investigate these differences. For 2001-2003, the overall forecasts were in the right direction (when compared with the long-term expected values), but were over-estimates. In 2004 the forecast was well under the observed production, and in 2005 the revised models produced a forecast within 5.1% of the actual production. Over the first five years of forecasting, the absolute deviance for the climate-adjustment models averaged 10.1%, just outside the targeted objective of 10%.

Non-specific conducting polymer-based array capable of monitoring odour emissions from a biofiltration system in a piggery building

A commercial non-specific gas sensor array system was evaluated in terms of its capability to monitor the odour abatement performance of a biofiltration system developed for treating emissions from a commercial piggery building. The biofiltration system was a modular system comprising an inlet ducting system, humidifier and closed-bed biofilter. It also included a gravimetric moisture monitoring and water application system for precise control of moisture content of an organic woodchip medium. Principal component analysis (PCA) of the sensor array measurements indicated that the biofilter outlet air was significantly different to both inlet air of the system and post-humidifier air. Data pre-processing techniques including normalising and outlier handling were applied to improve the odour discrimination performance of the non-specific gas sensor array. To develop an odour quantification model using the sensor array responses of the non-specific sensor array, PCA regression, artificial neural network (ANN) and partial least squares (PLS) modelling techniques were applied. The correlation coefficient (r(2)) values of the PCA, ANN, and PLS models were 0.44, 0.62 and 0.79, respectively.

Non-destructive prediction of 'Hass' avocado dry matter via FT-NIR spectroscopy.

BACKGROUND: The inability to consistently guarantee internal quality of horticulture produce is of major importance to the primary producer, marketers and ultimately the consumer. Currently, commercial avocado maturity estimation is based on the destructive assessment of percentage dry matter (%DM), and sometimes percentage oil, both of which are highly correlated with maturity. In this study the utility of Fourier transform (FT) near-infrared spectroscopy (NIRS) was investigated for the first time as a non-invasive technique for estimating %DM of whole intact 'Hass' avocado fruit. Partial least squares regression models were developed from the diffuse reflectance spectra to predict %DM, taking into account effects of intra-seasonal variation and orchard conditions. RESULTS: It was found that combining three harvests (early, mid and late) from a single farm in the major production district of central Queensland yielded a predictive model for %DM with a coefficient of determination for the validation set of 0.76 and a root mean square error of prediction of 1.53% for DM in the range 19.4-34.2%. CONCLUSION: The results of the study indicate the potential of FT-NIRS in diffuse reflectance mode to non-invasively predict %DM of whole 'Hass' avocado fruit. When the FT-NIRS system was assessed on whole avocados, the results compared favourably against data from other NIRS systems identified in the literature that have been used in research applications on avocados.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

Predicting Additive and Non-additive Genetic Effects from Trials Where Traits Are Affected by Interplot Competition

There are two key types of selection in a plant breeding program, namely selection of hybrids for potential commercial use and the selection of parents for use in future breeding. Oakey et al. (in Theoretical and Applied Genetics 113, 809-819, 2006) showed how both of these aims could be achieved using pedigree information in a mixed model analysis in order to partition genetic effects into additive and non-additive effects. Their approach was developed for field trial data subject to spatial variation. In this paper we extend the approach for data from trials subject to interplot competition. We show how the approach may be used to obtain predictions of pure stand additive and non-additive effects. We develop the methodology in the context of a single field trial using an example from an Australian sorghum breeding program.

The (Non)Effects of Lethal Population Control on the Diet of Australian Dingoes

Top-predators contribute to ecosystem resilience, yet individuals or populations are often subject to lethal control to protect livestock, managed game or humans from predation. Such management actions sometimes attract concern that lethal control might affect top-predator function in ways ultimately detrimental to biodiversity conservation. The primary function of a predator is predation, which is often investigated by assessing their diet. We therefore use data on prey remains found in 4,298 Australian dingo scats systematically collected from three arid sites over a four year period to experimentally assess the effects of repeated broad-scale poison-baiting programs on dingo diet. Indices of dingo dietary diversity and similarity were either identical or near-identical in baited and adjacent unbaited treatment areas in each case, demonstrating no control-induced change to dingo diets. Associated studies on dingoes' movement behaviour and interactions with sympatric mesopredators were similarly unaffected by poison-baiting. These results indicate that mid-sized top-predators with flexible and generalist diets (such as dingoes) may be resilient to ongoing and moderate levels of population control without substantial alteration of their diets and other related aspects of their ecological function.

Non-target species interaction with sodium fluoroacetate (1080) meat bait for controlling feral pigs (Sus scrofa)

Fresh meat baits containing sodium fluoroacetate (1080) are widely used for controlling feral pigs in Queensland, but there is a potential poisoning risk to non-target species. This study investigated the non-target species interactions with meat bait by comparing the time until first approach, investigation, sample and consumption, and whether dying bait green would reduce interactions. A trial assessing species interactions with undyed bait was completed at Culgoa Floodplain National Park, Queensland. Meat baits were monitored for 79 consecutive days with camera traps. Of 40 baits, 100% were approached, 35% investigated (moved) and 25% sampled, and 25% consumed. Monitors approached (P < 0.05) and investigated (P < 0.05) the bait more rapidly than pigs or birds, but the median time until first sampling was not significantly different (P > 0.05), and did not consume any entire bait. A trial was conducted at Whetstone State Forest, southern Queensland, with green-dyed and undyed baits monitored for eight consecutive days with cameras. Of 60 baits, 92% were approached and also investigated by one or more non-target species. Most (85%) were sampled and 57% were consumed, with monitors having slightly more interaction with undyed baits than with green-dyed baits. Mean time until first approach and sample differed significantly between species groups (P = 0.038 and 0.007 respectively) with birds approaching sooner (P < 0.05) and monitors sampling later (P < 0.05) than other (unknown) species (P > 0.05). Undyed bait was sampled earlier (mean 2.19 days) than green-dyed bait (2.7 days) (P = 0.003). Data from the two trials demonstrate that many non-target species regularly visit and sample baits. The use of green-dyed baits may help reduce non-target uptake, but testing is required to determine the effect on attractiveness to feral pigs. Further research is recommended to quantify the benefits of potential strategies to reduce the non-target uptake of meat baits to help improve the availability of bait to feral pigs.