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.

Long-run growth effects of taxation in a non-scale growth model with innovation

In previous studies, taxing income or consumption hinders long-run growth. Incorporating saving and leisure into the non-scale Schumpeterian model of [Journal of Political Economy 107 (1999) 715-730], we show that the usual growth effects of taxing consumption and labor income do not exist. (C) 2002 Elsevier Science B.V. All rights reserved.

The damped and coherently-driven Jaynes-Cummings model

We investigate the influence of a single-mode cavity on the Autler-Townes doublet that arises when a three-level atom is strongly driven by a laser field tuned to one of the atomic transitions and probed by a tunable, weak field coupled to the other transition. We assume that the cavity mode is coupled to the driven transition and the cavity and laser frequencies are equal to the atomic transition frequency. We find that the Autler-Townes spectrum can have one, two or three peaks depending on the relative magnitudes of the Rabi frequencies of the cavity and driving fields. We show that, in order to understand the three-peaked spectrum, it is necessary to go beyond the secular approximation, leading to interesting quantum interference effects. We find that the positions and relative intensities of the three spectral components are affected strongly by the atom-cavity coupling strength g and the cavity damping K. For an increasing g and/or decreasing K the triplet evolves into a single peak. This results in 'undressing' of the system such that the atom collapses into its ground state. We interpret the spectral features in terms of the semiclassical dressed-atom model, and also provide complementary views of the cavity effects in terms of quantum Langevin equations and the fully quantized, 'double -dressing' model.

Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR) -containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor. using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general. cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, interior/exterior orientations of side chains. and backbone conformation of the LRR structures can be predicted correctly. On the other hand. the analysis revealed that, in some cases. it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.

A compartmental model of hepatic disposition kinetics: 1. Model development and application to linear kinetics

The conventional convection-dispersion model is widely used to interrelate hepatic availability (F) and clearance (Cl) with the morphology and physiology of the liver and to predict effects such as changes in liver blood flow on F and Cl. The extension of this model to include nonlinear kinetics and zonal heterogeneity of the liver is not straightforward and requires numerical solution of partial differential equation, which is not available in standard nonlinear regression analysis software. In this paper, we describe an alternative compartmental model representation of hepatic disposition (including elimination). The model allows the use of standard software for data analysis and accurately describes the outflow concentration-time profile for a vascular marker after bolus injection into the liver. In an evaluation of a number of different compartmental models, the most accurate model required eight vascular compartments, two of them with back mixing. In addition, the model includes two adjacent secondary vascular compartments to describe the tail section of the concentration-time profile for a reference marker. The model has the added flexibility of being easy to modify to model various enzyme distributions and nonlinear elimination. Model predictions of F, MTT, CV2, and concentration-time profile as well as parameter estimates for experimental data of an eliminated solute (palmitate) are comparable to those for the extended convection-dispersion model.

Upstream porthole injection in a 2-D scramjet model

Injection from portholes upstream of the combustion chamber was investigated as a method of delivering fuel into a scramjet. This method reduces the viscous drag on a model by allowing a reduction in the length of the combustion chamber. At experimental enthalpies of 3.0 MJ/kg in the T4 shock tunnel, there was no evidence of combustion in the intake, either by shadowgraph or pressure measurements. Combustion was observed in the combustion chamber. A theoretical extension of these results is made to a hot wall scenario.

Drop penetration into porous powder beds

The kinetics of drop penetration were studied by filming single drops of several different fluids (water, PEG200, PEG600, and HPC solutions) as they penetrated into loosely packed beds of glass ballotini, lactose, zinc oxide, and titanium dioxide powders. Measured times ranged from 0.45 to 126 s and depended on the powder particle size,viscosity, surface tensions, and contact angle. The experimental drop penetration times were compared to existing theoretical predictions by M. Denesuk et al. (J. Colloid Interface Sci. 158, 114, 1993) and S. Middleman (Modeling Axisymmetric Flows: Dynamics of Films, Jets, and Drops, Academic Press, San Diego, 1995) but did not agree. Loosely packed powder beds tend to have a heterogeneous bed structure containing large macrovoids which do not participate in liquid flow but are included implicitly in the existing approach to estimating powder pore size. A new two-phase model was proposed where the total volume of the macrovoids was assumed to be the difference between the bed porosity and the tap porosity. A new parameter, the effective porosity (epsilon)eff, was defined as the tap porosity multiplied by the fraction of pores that terminate at a macrovoid and are effectively blocked pores. The improved drop penetration model was much more successful at estimating the drop penetration time on all powders and the predicted times were generally within an order of magnitude of the experimental results. (C) 2002 Elsevier Science (USA).

Mapping Gray's model of personality onto the Eysenck Personality Profiler (EPP)

The aim of this study was to determine how well Gray's model of personality [Gray, J.A. (1982). The neuropsychology of anxiety: an enquiry into the functions of the septo-hippocampal system. Oxfords: Oxford University Press, Gray, J.A. (1987). The psychology of fear and stress. Cambridge: Cambridge University Press], as measured by the Gray Wilson Personality Questionnaire (GWPQ), can provide a full description of personality as measured by the primary scales of the Eysenck Personality Profiler (EPP) and the type scales of the short version or the EPQ-R. Factor analysis of the GWPQ the Anxiety and linpulsivity scales of the EPP and the Learning Styles Questionnaire (LSQ) showed that the GWPQ seemed to measure general activation and inhibition factors, but not the finer features of Gray's theory. When the GWPQ scales were regressed against each scale of the EPP., it was round that they generally provide only a reasonable explanation of the EPP primary scales. It is concluded that the GWPQ measures general propel-ties of Gray's model, that the linpulsivity and Anxiety scales of the EPP also scent related to the GWPQ scales, and that Gray's model of personality provides only a partial explanation of personality in general. (C) 2002 Published by Elsevier Science Ltd. All rights reserved.

Dinoflagellate symbiosis: strategies and adaptations for the aquisition and fixation of inorganic carbon

Dinoflagellates exist in symbiosis with a number of marine invertebrates including giant clams, which are the largest of these symbiotic organisms. The dinoflagellates (Symbiodinium sp.) live intercellularly within tubules in the mantle of the host clam. The transport of inorganic carbon (Ci) from seawater to Symbiodinium (=zooxanthellae) is an essential function of hosts that derive the majority of their respiratory energy from the photosynthate exported by the zooxanthellae. Immunolocalisation studies show that the host has adapted its physiology to acquire, rather than remove CO2, from the haemolymph and clam tissues. Two carbonic anhydrase (CA) isoforms (32 and 70 kDa) play an essential part in this process. These have been localised to the mantle and gill tissues where they catalyse the interconversion of HCO3- to CO2, which then diffuses into the host tissues. The zooxanthellae exhibit a number of strategies to maximise Ci acquisition and utilisation. This is necessary as they express a form II Rubisco that has poor discrimination between CO2 and O-2. Evidence is presented for a carbon concentrating mechanism (CCM) to overcome. this disadvantage. The CCM incorporates the presence of a light-activated CA activity, a capacity to take up both HCO3- and CO2, an ability to accumulate an elevated concentration of Ci within the algal cell, and localisation of Rubisco to the pyrenoid. These algae also express both external and intracellular CAs, with the intracellular isoforms being localised to the thylakoid lumen and pyrenoid. These results have been incorporated into a model that explains the transport of Ci from seawater through the clam to the zooxanthellae.

Infusing the use of seasonal climate forecasting into crop management practice in North East Australia using discussion support software

Seasonal climate forecasting offers potential for improving management of crop production risks in the cropping systems of NE Australia. But how is this capability best connected to management practice? Over the past decade, we have pursued participative systems approaches involving simulation-aided discussion with advisers and decision-makers. This has led to the development of discussion support software as a key vehicle for facilitating infusion of forecasting capability into practice. In this paper, we set out the basis of our approach, its implementation and preliminary evaluation. We outline the development of the discussion support software Whopper Cropper, which was designed for, and in close consultation with, public and private advisers. Whopper Cropper consists of a database of simulation output and a graphical user interface to generate analyses of risks associated with crop management options. The charts produced provide conversation pieces for advisers to use with their farmer clients in relation to the significant decisions they face. An example application, detail of the software development process and an initial survey of user needs are presented. We suggest that discussion support software is about moving beyond traditional notions of supply-driven decision support systems. Discussion support software is largely demand-driven and can compliment participatory action research programs by providing cost-effective general delivery of simulation-aided discussions about relevant management actions. The critical role of farm management advisers and dialogue among key players is highlighted. We argue that the discussion support concept, as exemplified by the software tool Whopper Cropper and the group processes surrounding it, provides an effective means to infuse innovations, like seasonal climate forecasting, into farming practice. Crown Copyright (C) 2002 Published by Elsevier Science Ltd. All rights reserved.

Distribution kinetics of solutes in the isolated in-situ perfused rat head using the multiple indicator dilution technique and a physiological two-barrier model

The purpose of this study was to determine the pharmacokinetics of [C-14]diclofenac, [C-14]salicylate and [H-3]clonidine using a single pass rat head perfusion preparation. The head was perfused with 3-[N-morpholino] propane-sulfonic acid-buffered Ringer's solution. Tc-99m-red blood cells and a drug were injected in a bolus into the internal carotid artery and collected from the posterior facial vein over 28 min. A two-barrier stochastic organ model was used to estimate the statistical moments of the solutes. Plasma, interstitial and cellular distribution volumes for the solutes ranged from 1.0 mL (diclofenac) to 1.6 mL (salicylate), 2.0 mL (diclofenac) to 4.2 mL (water) and 3.9 mL (salicylate) to 20.9 mL (diclofenac), respectively. A comparison of these volumes to water indicated some exclusion of the drugs from the interstitial space and salicylate from the cellular space. Permeability-surface area (PS) products calculated from plasma to interstitial fluid permeation clearances (CLPI) (range 0.02-0.40 mL s(-1)) and fractions of solute unbound in the perfusate were in the order: diclofenac>salicylate >clonidine>sucrose (from 41.8 to 0.10 mL s(-1)). The slow efflux of diclofenac, compared with clonidine and salicylate, may be related to its low average unbound fraction in the cells. This work accounts for the tail of disposition curves in describing pharmacokinetics in the head.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

The Three-dimensional Solution Structure of NaD1, a New Floral Defensin from Nicotiana alata and its Application to a Homology Model of the Crop Defense Protein alfAFP

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of NaD1, a novel antifungal and insecticidal protein isolated from the flowers of Nicotiana alata. NaD1 is a basic, cysteine-rich protein of 47 residues and is the first example of a plant defensin from flowers to be characterized structurally. Its three-dimensional structure consists of an a-helix and a triple-stranded anti-parallel beta-sheet that are stabilized by four intramolecular disulfide bonds. NaD1 features all the characteristics of the cysteine-stabilized up motif that has been described for a variety of proteins of differing functions ranging from antibacterial insect defensins and ion channel-perturbing scorpion toxins to an elicitor of the sweet taste response. The protein is biologically active against insect pests, which makes it a potential candidate for use in crop protection. NaD1 shares 31% sequence identity with alfAFP, an antifungal protein from alfalfa that confers resistance to a fungal pathogen in transgenic potatoes. The structure of NaD1 was used to obtain a homology model of alfAFP, since NaD1 has the highest level of sequence identity with alfAFP of any structurally characterized antifungal defensin. The structures of NaD1 and alfAFP were used in conjunction with structure - activity data for the radish defensin Rs-AFP2 to provide an insight into structure-function relationships. In particular, a putative effector site was identified in the structure of NaD1 and in the corresponding homology model of alfAFP. (C) 2002 Elsevier Science Ltd. All rights reserved.

An investigation into commitment in non-Western industrial marketing relationships

This paper reports an investigation into the antecedents of commitment in non-Western industrial marketing relationships. The authors draw the antecedents from extant literature and posit that commitment is related to trust (integrity and reliability), communication quality, conflict, and similarity (social, ethnic, and economic). It is further argued that trust mediates the effects of communication, conflict, and similarity on commitment. As an extension, the authors examine the moderating effects of normative contracts (an implicit understanding of roles and responsibilities) on the construct interrelationships. The hypotheses are tested using data collected from approximately 150 industrial marketing relationships sampled from overseas Chinese firms. The results generally support the authors' framework; however, the mediating hypotheses are not supported. There is evidence of systematic differences in the effects of the studied antecedents on commitment and trust. Furthermore, a multigroup analysis provides evidence of significant moderating effects due to contracting mode. The study provides new insights into the theory and practice of industrial marketing. (C) 2003 Elsevier Science Inc. All rights reserved.

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.