Statistical methods for analysis of multi-harvest data from perennial pasture variety selection trials

Variety selection in perennial pasture crops involves identifying best varieties from data collected from multiple harvest times in field trials. For accurate selection, the statistical methods for analysing such data need to account for the spatial and temporal correlation typically present. This paper provides an approach for analysing multi-harvest data from variety selection trials in which there may be a large number of harvest times. Methods are presented for modelling the variety by harvest effects while accounting for the spatial and temporal correlation between observations. These methods provide an improvement in model fit compared to separate analyses for each harvest, and provide insight into variety by harvest interactions. The approach is illustrated using two traits from a lucerne variety selection trial. The proposed method provides variety predictions allowing for the natural sources of variation and correlation in multi-harvest data.

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.

How can warren destruction by ripping control European wild rabbits (Oryctolagus cuniculus) on large properties in the Australian arid zone?

Context: For over 100 years, control efforts have been unable to stop rabbits causing damage to cattle production and native plants and animals on large properties in arid parts of Australia. Warren destruction by ripping has shown promise, but doubts about long-term success and the perceived expense of treating vast areas have led to this technique not being commonly used. Aims: This study measured the long-term reduction in rabbit activity and calculated the potential cost saving associated with treating just the areas where rabbits are believed to survive drought. Wealso considered whether ripping should be used in a full-scale rabbit control program on a property where rabbits have been exceptionally resilient to the influence of biological and other control measures. Methods: Rabbits were counted along spotlight transects before warrens were ripped and during the two years after ripping, in treated and untreated plots. Rabbit activity was recorded to determine the immediate and long-term impact of ripping, up to seven years after treatment. The costs of ripping warrens within different distances from drought refuge areas were calculated. Key results: Destroying rabbit warrens by ripping caused an immediate reduction in rabbit activity and there were still 98% fewer rabbits counted by spotlight in ripped plots five months after ripping. Seven years after ripping no active warrens were found in ripped plots, whereas 57% of warrens in unripped plots showed signs of rabbit activity. The cost of ripping only the areas where rabbits were likely to seek refuge from drought was calculated to be less than 4%of the cost of ripping all warrens on the property. Conclusions: Destroying rabbit warrens by ripping is a very effective way of reducing rabbit numbers on large properties in arid Australia. Ripping should commence in areas used by rabbits to survive drought. It is possible that no further ripping will be required. Implications: Strategic destruction of warrens in drought refuge areas could provide an alternative to biological control for managing rabbits on large properties in the Australian arid zone.

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.

Characterization of north-eastern Australian environments using APSIM for increasing rainfed maize production

Recurring water stresses are a major risk factor for rainfed maize cropping across the highly diverse agro-ecological environments of Queensland (Qld) and northern New South Wales (NNSW). Enhanced understanding of such agro-ecological diversity is necessary to more consistently sample target production environments for testing and targeting release of improved germplasm, and to improve the efficiency of the maize pre-breeding and breeding programs of Qld and New South Wales. Here, we used the Agricultural Production Systems Simulator (APSIM) – a well validated maize crop model to characterize the key distinctive water stress patterns and risk to production across the main maize growing regions of Qld and NNSW located between 15.8° and 31.5°S, and 144.5° and 151.8°E. APSIM was configured to simulate daily water supply demand ratios (SDRs) around anthesis as an indicator of the degree of water stress, and the final grain yield. Simulations were performed using daily climatic records during the period between 1890 and 2010 for 32 sites-soils in the target production regions. The runs were made assuming adequate nitrogen supply for mid-season maize hybrid Pioneer 3153. Hierarchical complete linkage analyses of the simulated yield resulted in five major clusters showing distinct probability distribution of the expected yields and geographic patterns. The drought stress patterns and their frequencies using SDRs were quantified using multivariate statistical methods. The identified stress patterns included no stress, mid-season (flowering) stress, and three terminal stresses differing in terms of severity. The combined frequency of flowering and terminal stresses was highest (82.9%), mainly in sites-soils combinations in the west of Qld and NNSW. Yield variability across the different sites-soils was significantly related to the variability in frequencies of water stresses. Frequencies of water stresses within each yield cluster tended to be similar, but different across clusters. Sites-soils falling within each yield cluster therefore could be treated as distinct maize production environments for testing and targeting newly developed maize cultivars and hybrids for adaptation to water stress patterns most common to those environments.

Improving the effectiveness of age-abundance indicators in the management of fisheries in Queensland, Australia

The development of fishery indicators is a crucial undertaking as it ultimately provides evidence to stakeholders about the status of fished species such as population size and survival rates. In Queensland, as in many other parts of the world, age-abundance indicators (e.g. fish catch rate and/or age composition data) are traditionally used as the evidence basis because they provide information on species life history traits as well as on changes in fishing pressures and population sizes. Often, however, the accuracy of the information from age-abundance indicators can be limited due to missing or biased data. Consequently, improved statistical methods are required to enhance the accuracy, precision and decision-support value of age-abundance indicators.