Robust Mixture Modelling Using the t Distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.

Water wave-driven seepage in marine sediments

The action of water waves moving over a porous seabed drives a seepage flux into and out of the marine sediments. The volume of fluid exchange per wave cycle may affect the rate of contaminant transport in the sediments. In this paper, the dynamic response of the seabed to ocean waves is treated analytically on the basis of pore-elastic theory applied to a porous seabed. The seabed is modelled as a semi-infinite, isotropic, homogeneous material. Most previous investigations on the wave-seabed interaction problem have assumed quasi-static conditions within the seabed, although dynamic behaviour often occurs in natural environments. Furthermore, wave pressures used in the previous approaches were obtained from conventional ocean wave theories: which are based on the assumption of an impermeable rigid seabed. By introducing a complex wave number, we derive a new wave dispersion equation, which includes the seabed characteristics (such as soil permeability, shear modulus, etc.). Based on the new closed-form analytical solution, the relative differences of the wave-induced seabed response under dynamic and quasi-static conditions are examined. The effects of wave and soil parameters on the seepage flux per wave cycle are also discussed in detail. (C) 2000 Elsevier Science Ltd. All rights reserved.

Quantum measurement and stochastic processes in mesoscopic conductors

In quantum measurement theory it is necessary to show how a, quantum source conditions a classical stochastic record of measured results. We discuss mesoscopic conductance using quantum stochastic calculus to elucidate the quantum nature of the measurement taking place in these systems. To illustrate the method we derive the current fluctuations in a two terminal mesoscopic circuit with two tunnel barriers containing a single quasi bound state on the well. The method enables us to focus on either the incoming/ outgoing Fermi fields in the leads, or on the irreversible dynamics of the well state itself. We show an equivalence between the approach of Buttiker and the Fermi quantum stochastic calculus for mesoscopic systems.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

The HI Parkes All Sky Survey: southern observations, calibration and robust imaging

The acquisition of HI Parkes All Shy Survey (HIPASS) southern sky data commenced at the Australia Telescope National Facility's Parkes 64-m telescope in 1997 February, and was completed in 2000 March. HIPASS is the deepest HI survey yet of the sky south of declination +2 degrees, and is sensitive to emission out to 170 h(75)(-1) Mpc. The characteristic root mean square noise in the survey images is 13.3 mJy. This paper describes the survey observations, which comprise 23 020 eight-degree scans of 9-min duration, and details the techniques used to calibrate and image the data. The processing algorithms are successfully designed to be statistically robust to the presence of interference signals, and are particular to imaging point (or nearly point) sources. Specifically, a major improvement in image quality is obtained by designing a median-gridding algorithm which uses the median estimator in place of the mean estimator.

Estimating HIV incidence rates from age prevalence data in epidemic situations

We present a method of estimating HIV incidence rates in epidemic situations from data on age-specific prevalence and changes in the overall prevalence over time. The method is applied to women attending antenatal clinics in Hlabisa, a rural district of KwaZulu/Natal, South Africa, where transmission of HIV is overwhelmingly through heterosexual contact. A model which gives age-specific prevalence rates in the presence of a progressing epidemic is fitted to prevalence data for 1998 using maximum likelihood methods and used to derive the age-specific incidence. Error estimates are obtained using a Monte Carlo procedure. Although the method is quite general some simplifying assumptions are made concerning the form of the risk function and sensitivity analyses are performed to explore the importance of these assumptions. The analysis shows that in 1998 the annual incidence of infection per susceptible woman increased from 5.4 per cent (3.3-8.5 per cent; here and elsewhere ranges give 95 per cent confidence limits) at age 15 years to 24.5 per cent (20.6-29.1 per cent) at age 22 years and declined to 1.3 per cent (0.5-2.9 per cent) at age 50 years; standardized to a uniform age distribution, the overall incidence per susceptible woman aged 15 to 59 was 11.4 per cent (10.0-13.1 per cent); per women in the population it was 8.4 per cent (7.3-9.5 per cent). Standardized to the age distribution of the female population the average incidence per woman was 9.6 per cent (8.4-11.0 per cent); standardized to the age distribution of women attending antenatal clinics, it was 11.3 per cent (9.8-13.3 per cent). The estimated incidence depends on the values used for the epidemic growth rate and the AIDS related mortality. To ensure that, for this population, errors in these two parameters change the age specific estimates of the annual incidence by less than the standard deviation of the estimates of the age specific incidence, the AIDS related mortality should be known to within +/-50 per cent and the epidemic growth rate to within +/-25 per cent, both of which conditions are met. In the absence of cohort studies to measure the incidence of HIV infection directly, useful estimates of the age-specific incidence can be obtained from cross-sectional, age-specific prevalence data and repeat cross-sectional data on the overall prevalence of HIV infection. Several assumptions were made because of the lack of data but sensitivity analyses show that they are unlikely to affect the overall estimates significantly. These estimates are important in assessing the magnitude of the public health problem, for designing vaccine trials and for evaluating the impact of interventions. Copyright (C) 2001 John Wiley & Sons, Ltd.

Modelling life history strategies with capture-recapture data: Evolutionary demography of the water skink Eulamprus tympanum

Matrix population models, elasticity analysis and loop analysis can potentially provide powerful techniques for the analysis of life histories. Data from a capture-recapture study on a population of southern highland water skinks (Eulamprus tympanum) were used to construct a matrix population model. Errors in elasticities were calculated by using the parametric bootstrap technique. Elasticity and loop analyses were then conducted to identify the life history stages most important to fitness. The same techniques were used to investigate the relative importance of fast versus slow growth, and rapid versus delayed reproduction. Mature water skinks were long-lived, but there was high immature mortality. The most sensitive life history stage was the subadult stage. It is suggested that life history evolution in E. tympanum may be strongly affected by predation, particularly by birds. Because our population declined over the study, slow growth and delayed reproduction were the optimal life history strategies over this period. Although the techniques of evolutionary demography provide a powerful approach for the analysis of life histories, there are formidable logistical obstacles in gathering enough high-quality data for robust estimates of the critical parameters.

The relationship between cannabis use and other substance use in the general population

This study examined if (1) there is an association in the general population between cannabis use, DSM-IV abuse and dependence, and other substance use and DSM-IV substance abuse/dependence; (2) if so, is it explained by demographic characteristics or levels of neuroticism? It used data from the Australian National Survey of Mental Health and Well-Being (NSMHWB), a stratified, multistage probability sample of 10641 adults, representative of the general population. DSM-IV diagnoses of substance abuse and dependence were derived using the Composite International Diagnostic Interview (CIDI). There was a strong bivariate association between involvement with cannabis use in the past 12 months and other substance use, abuse and dependence. In particular, cannabis abuse and dependence were highly associated with increased risks of other substance dependence. These associations remained after including other variables in multiple regression. Cannabis use without disorder was strongly related to other drug use, an association that was not explained by other variables considered here. The high likelihood of other substance use and substance use disorders needs to be considered among persons seeking treatment for cannabis use problems. (C) 2001 Elsevier Science Ireland Ltd. All rights reserved.

Patterns of co-morbidity between alcohol use and other substance use in the Australian population

The present study describes patterns of co-morbidity between alcohol use and other substance use problems in the Australian population using data from the 1997 National Survey of Mental Health and Well-Being. Multiple regression analyses examined whether the observed associations between alcohol and other drug use disorders were explained by other variables, including demographic characteristics and neuroticism. We also assessed whether the presence of co-morbid substance use disorders affected treatment seeking for a mental health problem. Alcohol use was related strongly to the use of other substances. Those who did not report alcohol use within the past 12 months were less likely to report using tobacco, cannabis, sedatives, stimulants or opiates. Higher rates again were observed among those with alcohol use disorders: half (51%) of those who were alcohol-dependent were regular tobacco smokers, one-third had used cannabis (32%); 15% reported other drug use; 15% met criteria for a cannabis use disorder and 7% met criteria for another drug use disorder. These associations were not accounted for by the demographic and other variables considered here. Co-morbid substance use disorders (sedatives, stimulants or opioids) predicted a high likelihood of seeking treatment for a mental health problem among alcohol-dependent people.

Cutting edge: Granzymes A and B are not essential for perforin-mediated tumor rejection

Controversy still exists regarding the biological function of granzyme serine proteases released with perforin from the cytotoxic granules of NK cells and CTLs. In particular, it is not clear whether the major granzymes, A and A play an essential role in tumor rejection mediated by the perforin pathway. We have now examined the relative importance of perforin and granzyme A and B clusters in five different tumor models that stringently distinguish their importance. We conclude that granzyme A and B clusters are not essential for CTL- and NK cell-mediated rejection of spontaneous and experimental tumors, raising the likelihood that either perforin alone or in combination with an additional granzyme or granule component(s) mediates cytotoxicity of tumor cells in vivo.

Factors associated with rural practice among Australian-trained general practitioners

Objective: To determine the factors associated with general practitioners' current practice location, with particular emphasis on rural location. Design: Observational, retrospective, case-control study using a self-administered questionnaire. Setting: Australian general practices in December 2000. Participants: 2414 Australian-trained rural and urban GPs. Main outcome measure: Current urban or rural practice location. Results: For Australia as a whole, rural GPs were more likely to be male (odds ratio [OR], 1.42; 95% CI, 1.17-1.73), Australian-born (OR, 1.95; 95% CI, 1.55-2.45), and to report attending a rural primary school for some (OR, 2.21; 95% CI, 1.69-2.89) or all (OR, 2.79; 95% CI, 1.94-4.00) of their primary schooling. Rural GPs' partners or spouses were also more likely to report some (OR, 2.75; 95% CI, 2.07-3.66) or all (OR, 2.86; 95% CI, 2.02-4.05) rural primary schooling. A rural background in both GP and partner produced the highest likelihood of rural practice (OR, 6.28; 95% CI, 4.26-9.25). For individual jurisdictions, a trend towards more rural GPs being men was only significant in Tasmania. In all jurisdictions except Tasmania and the Northern Territory, rural GPs were more likely to be Australian-born. Conclusions: GPs' and their partners' rural background (residence and primary and secondary schooling) influences choice of practice location, with partners' background appearing to exert more influence.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Phylogeny and evolution of anomalous roots in Daviesia (Fabaceae : Mirbelieae)

The phylogeny of the Australian legume genus Daviesia was estimated using sequences of the internal transcribed spacers of nuclear ribosomal DNA. Partial congruence was found with previous analyses using morphology, including strong support for monophyly of the genus and for a sister group relationship between the clade D. pachyloma and the rest of the genus. A previously unplaced bird-pollinated species, anceps + D. D. epiphyllum, was well supported as sister to the only other bird-pollinated species in the genus, D. speciosa, indicating a single origin of bird pollination in their common ancestor. Other morphological groups within Daviesia were not supported and require reassessment. A strong and previously unreported sister clade of Daviesia consists of the two monotypic genera Erichsenia and Viminaria. These share phyllode-like leaves and indehiscent fruits. The evolutionary history of cord roots, which have anomalous secondary thickening, was explored using parsimony. Cord roots are limited to three separate clades but have a complex history involving a small number of gains (most likely 0-3) and losses (0-5). The anomalous structure of cord roots ( adventitious vascular strands embedded in a parenchymatous matrix) may facilitate nutrient storage, and the roots may be contractile. Both functions may be related to a postfire resprouting adaptation. Alternatively, cord roots may be an adaptation to the low-nutrient lateritic soils of Western Australia. However, tests for association between root type, soil type, and growth habit were equivocal, depending on whether the variables were treated as phylogenetically dependent (insignificant) or independent ( significant).

Impact of undergraduate and postgraduate rural training, and medical school entry criteria on rural practice among Australian general practitioners: national study of 2414 doctors

Objective To determine the association between rural undergraduate training, rural postgraduate training and medical school entry criteria favouring rural students, on likelihood of working in rural Australian general practice. Methods National case - control study of 2414 rural and urban general practitioners (GPs) sampled from the Health Insurance Commission database. Participants completed a questionnaire providing information on demographics, current practice location and rural undergraduate and postgraduate experience. Results Rural GPs were more likely to report having had any rural undergraduate training [ odds ratio ( OR) 1.61, 95% confidence interval (CI) 1.32 - 1.95] than were urban GPs. Rural GPs were much more likely to report having had rural postgraduate training ( OR 3.14, 95% CI 2.57 - 3.83). As the duration of rural postgraduate training increased so did the likelihood of working as a rural GP: those reporting that more than half their postgraduate training was rural were most likely to be rural GPs ( OR 10.52, 95% CI 5.39 - 20.51). South Australians whose final high school year was rural were more likely to be rural GPs ( OR 3.18, 95% CI 0.99 - 10.22). Conclusions Undergraduate rural training, postgraduate training and medical school entry criteria favouring rural students, all are associated with an increased likelihood of being a rural GP. Longer rural postgraduate training is more strongly associated with rural practice. These findings argue for continuation of rural undergraduate training opportunities and rural entry schemes, and an expansion in postgraduate training opportunities for GPs.

Environmental influence on the stress corrosion cracking of rock bolts

In order to understand rock bolt Stress Corrosion Cracking (SCC), a series of experiments have been performed in Linearly Increasing Stress Test (LIST) apparatus. One series of experiments determined the threshold stress of various bolt metallurgies (900 MPa for Steel A, and 800 MPa for Steel B and C). The high values of threshold stress suggest that SCC begins in rock bolts when they are sheared by moving rock strata. Typical crack velocity values have been measured to be 2.5 x 10(-8) m s(-1), indicating that there is not much benefit for rock bolt steel of higher fracture toughness. Another series of experiments were performed to understand the environmental conditions causing SCC of steel A and galvanised Steel A rock bolt steel. SCC only occurred for environmental conditions for which produce hydrogen on the sample surface, leading to hydrogen embrittlement and SCC. Fracture surfaces of LIST samples failed by SCC were found to display the same fracture regions as fracture surfaces of rock bolts failed in service by SCC: Tearing Topography Surface (TTS), Corrugated Irregular Surface (CIS), quasi Micro Void Coalescence (qMVC) and Fast Fracture Surface (FFS). Water chemistry analysis were carried out on samples collected from various Australian mines in order to compare laboratory electrolyte conditions to those found in underground mines.