A zero-augmented gamma mixed model for longitudinal data with many zeros

In many occupational safety interventions, the objective is to reduce the injury incidence as well as the mean claims cost once injury has occurred. The claims cost data within a period typically contain a large proportion of zero observations (no claim). The distribution thus comprises a point mass at 0 mixed with a non-degenerate parametric component. Essentially, the likelihood function can be factorized into two orthogonal components. These two components relate respectively to the effect of covariates on the incidence of claims and the magnitude of claims, given that claims are made. Furthermore, the longitudinal nature of the intervention inherently imposes some correlation among the observations. This paper introduces a zero-augmented gamma random effects model for analysing longitudinal data with many zeros. Adopting the generalized linear mixed model (GLMM) approach reduces the original problem to the fitting of two independent GLMMs. The method is applied to evaluate the effectiveness of a workplace risk assessment teams program, trialled within the cleaning services of a Western Australian public hospital.

A mixture model-based approach to the clustering of microarray expression data

Motivation: This paper introduces the software EMMIX-GENE that has been developed for the specific purpose of a model-based approach to the clustering of microarray expression data, in particular, of tissue samples on a very large number of genes. The latter is a nonstandard problem in parametric cluster analysis because the dimension of the feature space (the number of genes) is typically much greater than the number of tissues. A feasible approach is provided by first selecting a subset of the genes relevant for the clustering of the tissue samples by fitting mixtures of t distributions to rank the genes in order of increasing size of the likelihood ratio statistic for the test of one versus two components in the mixture model. The imposition of a threshold on the likelihood ratio statistic used in conjunction with a threshold on the size of a cluster allows the selection of a relevant set of genes. However, even this reduced set of genes will usually be too large for a normal mixture model to be fitted directly to the tissues, and so the use of mixtures of factor analyzers is exploited to reduce effectively the dimension of the feature space of genes. Results: The usefulness of the EMMIX-GENE approach for the clustering of tissue samples is demonstrated on two well-known data sets on colon and leukaemia tissues. For both data sets, relevant subsets of the genes are able to be selected that reveal interesting clusterings of the tissues that are either consistent with the external classification of the tissues or with background and biological knowledge of these sets.

A genomewide search for type 2 diabetes-susceptibility genes in indigenous Australians

The prevalence of type 2 diabetes among Australian residents is 7.5%; however, prevalence rates up to six times higher have been reported for indigenous Australian communities. Epidemiological evidence implicates genetic factors in the susceptibility of indigenous Australians to type 2 diabetes and supports the hypothesis of the thrifty genotype, but, to date, the nature of the genetic predisposition is unknown. We have ascertained clinical details from a community of indigenous Australian descent in North Stradbroke Island, Queensland. In this population, the phenotype is characterized by severe insulin resistance. We have conducted a genomewide scan, at an average resolution of 10 cM, for type 2 diabetes-susceptibility genes in a large multigeneration pedigree from this community. Parametric linkage analysis undertaken using FASTLINK version 4.1p yielded a maximum two-point LOD score of +2.97 at marker D2S2345. Multipoint analysis yielded a peak LOD score of +3.9

Positive-P and Wigner representations for quantum-optical systems with nonorthogonal modes

We generalize the basic concepts of the positive-P and Wigner representations to unstable quantum-optical systems that are based on nonorthogonal quasimodes. This lays the foundation for a quantum description of such systems, such as, for example an unstable cavity laser. We compare both representations by calculating the tunneling times for an unstable resonator optical parametric oscillator.

Quantum correlated twin atomic beams via photodissociation of a molecular Bose-Einstein condensate

We study the process of photodissociation of a molecular Bose-Einstein condensate as a potential source of strongly correlated twin atomic beams. We show that the two beams can possess nearly perfect quantum squeezing in their relative numbers.

Shock standoff on hypersonic blunt bodies in nonequilibrium gas flows

This paper presents a new theory of hypersonic blunt-nose shock standoff, based on a compressibility coordinate transformation for inviscid flow. It embraces a wide range of nonequilibrium shock-layer chemistry and gas mixtures including ionization and freestream dissociation. An extended binary scaling property of the analysis is also demonstrated. Specific application is made here to the family of arbitrarily diluted dissociating diatomic gases, with parametric study results presented for the scaled shock standoff distance as a function of an appropriate blunt-nose region Damkohler number. Comparisons with other theories and data in the case of nitrogen are also given and discussed.

Comparison of two population pharmacokinetic programs, NONMEM and P-PHARM, for tacrolimus

Objectives: To compare the population modelling programs NONMEM and P-PHARM during investigation of the pharmacokinetics of tacrolimus in paediatric liver-transplant recipients. Methods: Population pharmacokinetic analysis was performed using NONMEM and P-PHARM on retrospective data from 35 paediatric liver-transplant patients receiving tacrolimus therapy. The same data were presented to both programs. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F). Covariates screened for influence on these parameters were weight, age, gender, post-operative day, days of tacrolimus therapy, transplant type, biliary reconstructive procedure, liver function tests, creatinine clearance, haematocrit, corticosteroid dose, and potential interacting drugs. Results: A satisfactory model was developed in both programs with a single categorical covariate - transplant type - providing stable parameter estimates and small, normally distributed (weighted) residuals. In NONMEM, the continuous covariates - age and liver function tests - improved modelling further. Mean parameter estimates were CL/F (whole liver) = 16.3 1/h, CL/F (cut-down liver) = 8.5 1/h and V/F = 565 1 in NONMEM, and CL/F = 8.3 1/h and V/F = 155 1 in P-PHARM. Individual Bayesian parameter estimates were CL/F (whole liver) = 17.9 +/- 8.8 1/h, CL/F (cutdown liver) = 11.6 +/- 18.8 1/h and V/F = 712 792 1 in NONMEM, and CL/F (whole liver) = 12.8 +/- 3.5 1/h, CL/F (cut-down liver) = 8.2 +/- 3.4 1/h and V/F = 221 1641 in P-PHARM. Marked interindividual kinetic variability (38-108%) and residual random error (approximately 3 ng/ml) were observed. P-PHARM was more user friendly and readily provided informative graphical presentation of results. NONMEM allowed a wider choice of errors for statistical modelling and coped better with complex covariate data sets. Conclusion: Results from parametric modelling programs can vary due to different algorithms employed to estimate parameters, alternative methods of covariate analysis and variations and limitations in the software itself.

Exploring complex fitness surfaces: Multiple ornamentation and polymorphism in male guppies

The sexual ornamentation used by male guppies to attract females comprises many components, each of which varies considerably among males. Although natural and sexual selection have been shown to contribute to divergence among populations in male sexual ornaments, the role of sexual selection in maintaining polymorphism within populations is less clear. We used both parametric quadratic regression and nonparametric projection pursuit regression techniques to reveal the major axes of non-linear sexual selection on male ornaments. We visualized the fitness surfaces defined by these axes using thin-plate splines to allow a direct comparison of the two methodologies. Identification of the major axes of selection and their visualization was critical in determining the form and strength of nonlinear selection. Both types of analysis revealed fitness surfaces comprising three peaks, suggesting that there is more than one way to make an attractive guppy. Disruptive selection may be an important process underlying the presence of multiple sexual ornaments and may contribute to the maintenance of the high levels of polymorphism in male sexual ornaments found in guppy populations.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

A mechanistic model for carbon dioxide corrosion of mild steel in the presence of protective iron carbonate films - Part 3: Film growth model

A model of iron carbonate (FeCO3) film growth is proposed, which is an extension of the recent mechanistic model of carbon dioxide (CO2) corrosion by Nesic, et al. In the present model, the film growth occurs by precipitation of iron carbonate once saturation is exceeded. The kinetics of precipitation is dependent on temperature and local species concentrations that are calculated by solving the coupled species transport equations. Precipitation tends to build up a layer of FeCO3 on the surface of the steel and reduce the corrosion rate. On the other hand, the corrosion process induces voids under the precipitated film, thus increasing the porosity and leading to a higher corrosion rate. Depending on the environmental parameters such as temperature, pH, CO2 partial pressure, velocity, etc., the balance of the two processes can lead to a variety of outcomes. Very protective films and low corrosion rates are predicted at high pH, temperature, CO2 partial pressure, and Fe2+ ion concentration due to formation of dense protective films as expected. The model has been successfully calibrated against limited experimental data. Parametric testing of the model has been done to gain insight into the effect of various environmental parameters on iron carbonate film formation. The trends shown in the predictions agreed well with the general understanding of the CO2 corrosion process in the presence of iron carbonate films. The present model confirms that the concept of scaling tendency is a good tool for predicting the likelihood of protective iron carbonate film formation.

Localization of a novel melanoma susceptibility locus to 1p22

Over the past 20 years, the incidence of cutaneous malignant melanoma (CMM) has increased dramatically worldwide. A positive family history of the disease is among the most established risk factors for CMM; it is estimated that 10% of CMM cases result from an inherited predisposition. Although mutations in two genes, CDKN2A and CDK4, have been shown to confer an increased risk of CMM, they account for only 20%-25% of families with multiple cases of CMM. Therefore, to localize additional loci involved in melanoma susceptibility, we have performed a genomewide scan for linkage in 49 Australian pedigrees containing at least three CMM cases, in which CDKN2A and CDK4 involvement has been excluded. The highest two-point parametric LOD score (1.82; recombination fraction [theta] 0.2) was obtained at D1S2726, which maps to the short arm of chromosome 1 (1p22). A parametric LOD score of 4.65 (theta = 0) and a nonparametric LOD score of 4.19 were found at D1S2779 in nine families selected for early age at onset. Additional typing yielded seven adjacent markers with LOD scores 13 in this subset, with the highest parametric LOD score, 4.95 (theta = 0) ( nonparametric LOD score 5.37), at D1S2776. Analysis of 33 additional multiplex families with CMM from several continents provided further evidence for linkage to the 1p22 region, again strongest in families with the earliest mean age at diagnosis. A nonparametric ordered sequential analysis was used, based on the average age at diagnosis in each family. The highest LOD score, 6.43, was obtained at D1S2779 and occurred when the 15 families with the earliest ages at onset were included. These data provide significant evidence of a novel susceptibility gene for CMM located within chromosome band 1p22.

Rock fracture mean trace length estimation and confidence interval calculation using maximum likelihood methods

There has been a resurgence of interest in the mean trace length estimator of Pahl for window sampling of traces. The estimator has been dealt with by Mauldon and Zhang and Einstein in recent publications. The estimator is a very useful one in that it is non-parametric. However, despite some discussion regarding the statistical distribution of the estimator, none of the recent works or the original work by Pahl provide a rigorous basis for the determination a confidence interval for the estimator or a confidence region for the estimator and the corresponding estimator of trace spatial intensity in the sampling window. This paper shows, by consideration of a simplified version of the problem but without loss of generality, that the estimator is in fact the maximum likelihood estimator (MLE) and that it can be considered essentially unbiased. As the MLE, it possesses the least variance of all estimators and confidence intervals or regions should therefore be available through application of classical ML theory. It is shown that valid confidence intervals can in fact be determined. The results of the work and the calculations of the confidence intervals are illustrated by example. (C) 2003 Elsevier Science Ltd. All rights reserved.

A model to investigate hepatic extraction of oxygen during anaesthesia in the dog

Measurement of hepatic oxygen extraction was performed on six healthy Greyhound dogs over a two hour period. The Greyhounds were anaesthetised and a right subcostal surgical incision performed. Ultrasonic flow transducers were used to measure flow rate in the hepatic artery and the portal vein. The blood oxygen tensions in arterial blood and in the portal and hepatic veins were also measured. Hepatic oxygen extraction remained stable throughout the study, despite a steady decline in arterial blood pressure. The methodology described in this study provides a direct measure of oxygen uptake by the liver in the dog and could readily be used to investigate hepatic uptake of drugs. (C) 2003 Elsevier Ltd. All rights reserved.

Postbuckling of piezoelectric FGM plates subject to thermo-electro-mechanical loading

In this paper, we examine the postbuckling behavior of functionally graded material FGM rectangular plates that are integrated with surface-bonded piezoelectric actuators and are subjected to the combined action of uniform temperature change, in-plane forces, and constant applied actuator voltage. A Galerkin-differential quadrature iteration algorithm is proposed for solution of the non-linear partial differential governing equations. To account for the transverse shear strains, the Reddy higher-order shear deformation plate theory is employed. The bifurcation-type thermo-mechanical buckling of fully clamped plates, and the postbuckling behavior of plates with more general boundary conditions subject to various thermo-electro-mechanical loads, are discussed in detail. Parametric studies are also undertaken, and show the effects of applied actuator voltage, in-plane forces, volume fraction exponents, temperature change, and the character of boundary conditions on the buckling and postbuckling characteristics of the plates. (C) 2003 Elsevier Science Ltd. All rights reserved.