On-line estimation of Allan variance parameters

A new online method is presented for estimation of the angular random walk and rate random walk coefficients of IMU (inertial measurement unit) gyros and accelerometers. The online method proposes a state space model and proposes parameter estimators for quantities previously measured from off-line data techniques such as the Allan variance graph. Allan variance graphs have large off-line computational effort and data storage requirements. The technique proposed here requires no data storage and computational effort of O(100) calculations per data sample.

Seasonal occurrence and distribution of Bryde's whales in the Hauraki Gulf, New Zealand

The Hauraki Gulf is a large, shallow embayment located north of Auckland City (36°51′S, 174°46′E), New Zealand. Bryde's whales (Balaenoptera edeni) are the most frequently observed balaenopterid in these waters. To assess the use of the Hauraki Gulf for this species, we examined the occurrence and distribution in relation to environmental parameters. Data were collected from a platform of opportunity during 674 daily surveys between March 2003 and February 2006. A total of 760 observations of Bryde's whales were recorded throughout the study period during 371 surveys. The number of Bryde's whales sighted/day was highest in winter, coinciding with the coolest median sea-surface temperature (14.6°C). Bryde's whales were recorded throughout the Hauraki Gulf in water depths ranging from 12.1–59.8 m (mean = 42.3, SD = 5.1). Cow–calf pairs were most frequently observed during the austral autumn in water depths of 29.9–53.9 m (mean = 40.8, SD = 5.2). Data from this study suggest Bryde's whales in the Hauraki Gulf exhibit a mix of both “inshore” and “offshore” characteristics from the Bryde's whales examined off the coast of South Africa. Based on complete mitochondrial DNA sequences, Sasaki et al. (2006) recognized two sister species of Bryde's whales: Balaenoptera brydei and B. edeni, with the latter including small-type, more coastal Bryde's whales from Japan, Hong Kong, and Australia. Their samples and samples in previous analyses of small-type whales, all originated from eastern and southeastern Asia. These authors did not include the forms of Bryde's whales that occur in other regions, e.g., in the Pacific off Peru (Valdivia et al. 1981), in the Atlantic off Brazil (Best 1977) and in the western Indian Ocean off South Africa (Best 1977). Recent genetic analysis using mtDNA from the “inshore” and “offshore” forms from South Africa confirms the offshore form is B. brydei, and establishes that the inshore form is more closely related to B. brydei than to B. edeni (Penry 2010). These different forms do vary considerably in their habitat use and ecology (refer to Table 1 for a detailed comparison between the South African inshore and offshore forms, as described by Best (1967, 1977) and the Bryde's whales from New Zealand (Wiseman 2008). Recent genetic analysis on the Bryde's whales in the Hauraki Gulf suggests they are B. brydei (Wiseman 2008). However, pending resolution of the uncertainty within and between species of this genus, we follow the Society of Marine Mammal's committee on taxonomy, who state that B. edeni applies to all Bryde's whales.

Identification of New Zealand bats in flight from analysis of echolocation calls by artificial neural networks.

Time-expanded and heterodyned echolocation calls of the New Zealand long-tailed Chalinolobus tuberculatus and lesser short-tailed bat Mystacina tuberculata were recorded and digitally analysed. Temporal and spectral parameters were measured from time-expanded calls and power spectra generated for both time-expanded and heterodyned calls. Artificial neural networks were trained to classify the calls of both species using temporal and spectral parameters and power spectra as input data. Networks were then tested using data not previously seen. Calls could be unambiguously identified using parameters and power spectra from time-expanded calls. A neural network, trained and tested using power spectra of calls from both species recorded using a heterodyne detector set to 40 kHz (the frequency with the most energy of the fundamental of C. tuberculatus call), could identify 99% and 84% of calls of C. tuberculatus and M. tuberculata, respectively. A second network, trained and tested using power spectra of calls from both species recorded using a heterodyne detector set to 27 kHz (the frequency with the most energy of the fundamental of M. tuberculata call), could identify 34% and 100% of calls of C. tuberculatus and M. tuberculata, respectively. This study represents the first use of neural networks for the identification of bats from their echolocation calls. It is also the first study to use power spectra of time-expanded and heterodyned calls for identification of chiropteran species. The ability of neural networks to identify bats from their echolocation calls is discussed, as is the ecology of both species in relation to the design of their echolocation calls.

A new approach to estimating trends in chlamydia incidence

Objectives Directly measuring disease incidence in a population is difficult and not feasible to do routinely. We describe the development and application of a new method of estimating at a population level the number of incident genital chlamydia infections, and the corresponding incidence rates, by age and sex using routine surveillance data. Methods A Bayesian statistical approach was developed to calibrate the parameters of a decision-pathway tree against national data on numbers of notifications and tests conducted (2001-2013). Independent beta probability density functions were adopted for priors on the time-independent parameters; the shape parameters of these beta distributions were chosen to match prior estimates sourced from peer-reviewed literature or expert opinion. To best facilitate the calibration, multivariate Gaussian priors on (the logistic transforms of) the time-dependent parameters were adopted, using the Matérn covariance function to favour changes over consecutive years and across adjacent age cohorts. The model outcomes were validated by comparing them with other independent empirical epidemiological measures i.e. prevalence and incidence as reported by other studies. Results Model-based estimates suggest that the total number of people acquiring chlamydia per year in Australia has increased by ~120% over 12 years. Nationally, an estimated 356,000 people acquired chlamydia in 2013, which is 4.3 times the number of reported diagnoses. This corresponded to a chlamydia annual incidence estimate of 1.54% in 2013, increased from 0.81% in 2001 (~90% increase). Conclusions We developed a statistical method which uses routine surveillance (notifications and testing) data to produce estimates of the extent and trends in chlamydia incidence.

Investigation into chemistry of new particle formation and growth in subtropical urban environment

The role of different chemical compounds, particularly organics, involved in the new particle formation (NPF) and its consequent growth are not fully understood. Therefore, this study was conducted to investigate the chemistry of aerosol particles during NPF events in an urban subtropical environment. Aerosol chemical composition was measured along with particle number size distribution (PNSD) and several other air quality parameters at five sites across an urban subtropical environment. An Aerodyne compact Time-of-Flight Aerosol Mass Spectrometer (c-TOF-AMS) and a TSI Scanning Mobility Particle Sizer (SMPS) measured aerosol chemical composition and PNSD, respectively. Five NPF events, with growth rates in the range 3.3-4.6 nm, were detected at two sites. The NPF events happened on relatively warmer days with lower humidity and higher solar radiation. Temporal percent fractions of nitrate, sulphate, ammonium and organics were modelled using the Generalised Additive Model (GAM), with a basis of penalised spline. Percent fractions of organics increased after the NPF events, while the mass fraction of ammonium and sulphate decreased. This uncovered the important role of organics in the growth of newly formed particles. Three organic markers, factors f43, f44 and f57, were calculated and the f44 vs f43 trends were compared between nucleation and non-nucleation days. f44 vs f43 followed a different pattern on nucleation days compared to non-nucleation days, whereby f43 decreased for vehicle emission generated particles, while both f44 and f43 decreased for NPF generated particles. It was found for the first time that vehicle generated and newly formed particles cluster in different locations on f44 vs f43 plot and this finding can be used as a tool for source apportionment of measured particles.

Melanoma cell colony expansion parameters revealed by approximate Bayesian computation

In vitro studies and mathematical models are now being widely used to study the underlying mechanisms driving the expansion of cell colonies. This can improve our understanding of cancer formation and progression. Although much progress has been made in terms of developing and analysing mathematical models, far less progress has been made in terms of understanding how to estimate model parameters using experimental in vitro image-based data. To address this issue, a new approximate Bayesian computation (ABC) algorithm is proposed to estimate key parameters governing the expansion of melanoma cell (MM127) colonies, including cell diffusivity, D, cell proliferation rate, λ, and cell-to-cell adhesion, q, in two experimental scenarios, namely with and without a chemical treatment to suppress cell proliferation. Even when little prior biological knowledge about the parameters is assumed, all parameters are precisely inferred with a small posterior coefficient of variation, approximately 2–12%. The ABC analyses reveal that the posterior distributions of D and q depend on the experimental elapsed time, whereas the posterior distribution of λ does not. The posterior mean values of D and q are in the ranges 226–268 µm2h−1, 311–351 µm2h−1 and 0.23–0.39, 0.32–0.61 for the experimental periods of 0–24 h and 24–48 h, respectively. Furthermore, we found that the posterior distribution of q also depends on the initial cell density, whereas the posterior distributions of D and λ do not. The ABC approach also enables information from the two experiments to be combined, resulting in greater precision for all estimates of D and λ.

Colonoscopy requirements of population screening for colorectal cancer in New Zealand

Aim: To estimate the colonoscopy burden of introducing population screening for colorectal cancer in New Zealand. Methods: Screening for colorectal cancer using biennial immunochemical faecal occult blood tests offered to people aged 50-74 years of age was modelled using population estimates from Statistics New Zealand for 2011-2031. Modelling to determine colonoscopy requirements was based on participation and test positivity rates from published results of screening programmes. Estimates of the number of procedures required for ongoing adenoma surveillance were calculated using screening literature results of adenoma yield, and New Zealand Guidelines for Adenoma Surveillance. Sensitivity analysis was undertaken on key parameters. Results: For a test positivity of 6.4%, biennial screening using immunochemical faecal occult blood testing with a 60% participation rate, would require 18,000 colonoscopies nationally, increasing to 28,000 by 2031. The majority of procedures are direct referrals from a positive FOBT, with surveillance colonoscopy numbers building over time. Conclusion: Colonoscopy requirements for immunochemical faecal occult blood based population screening for colorectal cancer are high. Significant expansion of services is required and careful management of surveillance procedures to ensure timely delivery of initial colonoscopies whilst maintaining symptomatic services. A model re-run informed by data from the screening pilot will allow improved estimates for the New Zealand setting.

Development of a new design expression for the in-plane shear capacity of the partially grouted concrete masonry walls

Partially grouted masonry walls subjected to in-plane shear exhibit a complex behaviour because of the influence of the aspect ratio, the pre-compression, the grouting pattern, the ratios of the horizontal and the vertical reinforcements, the boundary conditions and the characteristics of the constituent materials. The existing in-plane shear expressions for the partially grouted masonry are formulated as sum of strength of three parameters, namely, the masonry, the reinforcement and the axial force. The parameter ‘masonry’ includes the wall aspect ratio and the masonry compressive strength; the aspect ratio of the unreinforced panel inscribed into the grouted cores and bond beams are not considered, although failure is often dominated by these unreinforced masonry panels. This paper describes the dominance of these panels, particularly those that are squat, to the shear capacity of whole of shear walls. Further, the current design formulae are shown highly un-conservative by many researchers; this paper provides a potential reason for this un-conservativeness. It is shown that by including an additional term of the unreinforced panel aspect ratio a rational design formula could be established. This new expression is validated with independent test results reported in the literature – both Australian and overseas; the predictions are shown to be conservative.

A simple method for estimating growth parameters from multiple length-frequency data in presence of continuous recruitment

The extended recruitment season for short-lived species such as prawns biases the estimation of growth parameters from length-frequency data when conventional methods are used. We propose a simple method for overcoming this bias given a time series of length-frequency data. The difficulties arising from extended recruitment are eliminated by predicting the growth of the succeeding samples and the length increments of the recruits in previous samples. This method requires that some maximum size at recruitment can be specified. The advantages of this multiple length-frequency method are: it is simple to use; it requires only three parameters; no specific distributions need to be assumed; and the actual seasonal recruitment pattern does not have to be specified. We illustrate the new method with length-frequency data on the tiger prawn Penaeus esculentus from the north-western Gulf of Carpentaria, Australia.

A generalized estimating equations approach for analysis of the impact of new technology on a trawl fishery

The article describes a generalized estimating equations approach that was used to investigate the impact of technology on vessel performance in a trawl fishery during 1988-96, while accounting for spatial and temporal correlations in the catch-effort data. Robust estimation of parameters in the presence of several levels of clustering depended more on the choice of cluster definition than on the choice of correlation structure within the cluster. Models with smaller cluster sizes produced stable results, while models with larger cluster sizes, that may have had complex within-cluster correlation structures and that had within-cluster covariates, produced estimates sensitive to the correlation structure. The preferred model arising from this dataset assumed that catches from a vessel were correlated in the same years and the same areas, but independent in different years and areas. The model that assumed catches from a vessel were correlated in all years and areas, equivalent to a random effects term for vessel, produced spurious results. This was an unexpected finding that highlighted the need to adopt a systematic strategy for modelling. The article proposes a modelling strategy of selecting the best cluster definition first, and the working correlation structure (within clusters) second. The article discusses the selection and interpretation of the model in the light of background knowledge of the data and utility of the model, and the potential for this modelling approach to apply in similar statistical situations.

Atmospheric visibility and PM10 as indicators of new particle formation in an urban environment

It is well-known that new particle formation (NPF) in the atmosphere is inhibited by pre-existing particles in the air that act as condensation sinks to decrease the concentration and, thus, the supersaturation of precursor gases. In this study, we investigate the effects of two parameters - atmospheric visibility, expressed as the particle back-scatter coefficient (BSP), and PM10 particulate mass concentration, on the occurrences of NPF events in an urban environment where the majority of precursor gases originate from motor vehicle and industrial sources. This is the first attempt to derive direct relationships between each of these two parameters and the occurrence of NPF. NPF events were identified from data obtained with a neutral cluster and air ion spectrometer over 245 days within a calendar year. Bayesian logistic regression was used to determine the probability of observing NPF as functions of BSP and PM10. We show that the BSP at 08 h on a given day is a reliable indicator of an NPF event later that day. The posterior median probability of observing an NPF event was greater than 0.5 (95%) when the BSP at 08 h was less than 6.8 Mm-1.

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.

Virtual Special Issue: “The impact of temperature on human health: new insights”

Ambient temperature is one of the basic parameters characterising human comfort: are we too hot, too cold, or just right? The impact of temperature goes beyond comfort: inadequate temperature and temperature variations have consequences on human health, as the increasing numbers of studies have demonstrated. The topic is of particular significance at the times when climate change shifts the traditional – as we know them- temperature zones, and brings much wider temperature variations. For these reasons the impact of temperature on health has been one of the most popular topics among the articles submitted and published in Science of the Total Environment over the last few years. This Virtual Special Issue compiles 18 articles published in our journal on this topic since 2012. It is worth briefly summarizing the rich scientific insights brought by these articles, as well as broader considerations, particularly those extending to management, discussed by the authors of the articles.