Detecting environmental impacts on metapopulations of mound spring invertebrates - Assessing an incidence function model

We use a stochastic patch occupancy model of invertebrates in the Mound Springs ecosystem of South Australia to assess the ability of incidence function models to detect environmental impacts on metapopulations. We assume that the probability of colonisation decreases with increasing isolation and the probability of extinction is constant across spring vents. We run the models to quasi-equilibrium, and then impose an impact by increasing the local extinction probability. We sample the output at various times pre- and postimpact, and examine the probability of detecting a significant change in population parameters. The incidence function model approach turns out to have little power to detect environmental impacts on metapopulations with small numbers of patches. (C) 2001 Elsevier Science Ltd. All rights reserved.

Influence of patient age and implantation technique on the probability of re-replacement of the homograft aortic valve

Background and aim of the study: Results of valve re-replacement (reoperation) in 898 patients undergoing aortic valve replacement with cryopreserved homograft valves between 1975 and 1998 are reported. The study aim was to provide estimates of unconditional probability of valve reoperation and cumulative incidence function (actual risk) of reoperation. Methods: Valves were implanted by subcoronary insertion (n = 500), inclusion cylinder (n = 46), and aortic root replacement (n = 352). Probability of reoperation was estimated by adopting a mixture model framework within which estimates were adjusted for two risk factors: patient age at initial replacement, and implantation technique. Results: For a patient aged 50 years, the probability of reoperation in his/her lifetime was estimated as 44% and 56% for non-root and root replacement techniques, respectively. For a patient aged 70 years, estimated probability of reoperation was 16% and 25%, respectively. Given that a reoperation is required, patients with non-root replacement have a higher hazard rate than those with root replacement (hazards ratio = 1.4), indicating that non-root replacement patients tend to undergo reoperation earlier before death than root replacement patients. Conclusion: Younger patient age and root versus non-root replacement are risk factors for reoperation. Valve durability is much less in younger patients, while root replacement patients appear more likely to live longer and hence are more likely to require reoperation.

A metafrontier production function for estimation of technical efficiencies and technology gaps for firms operating under different technologies

This paper presents a metafrontier production function model for firms in different groups having different technologies. The metafrontier model enables the calculation of comparable technical efficiencies for firms operating under different technologies. The model also enables the technology gaps to be estimated for firms under different technologies relative to the potential technology available to the industry as a whole. The metafrontier model is applied in the analysis of panel data on garment firms in five different regions of Indonesia, assuming that the regional stochastic frontier production function models have technical inefficiency effects with the time-varying structure proposed by Battese and Coelli ( 1992).

Lameness, metabolic and digestive disorders, and technical efficiency in Danish dairy herds: a stochastic frontier production function approach

The relationship between reported treatments of lameness, metabolic disorders (milk fever, ketosis), digestive disorders, and technical efficiency (TE) was investigated using neutral and non-neutral stochastic frontier analysis (SFA). TE is estimated relative to the stochastic frontier production function for a sample of 574 Danish dairy herds collected in 1997. Contrary to most published results, but in line with the expected negative impact of disorders on the average cow milk production, herds reporting higher frequencies of milk fever are less technically efficient. Unexpectedly, however, the opposite results were observed for lameness, ketosis, and digestive disorders. The non-neutral stochastic frontier indicated that the opposite results are due to the relative. high productivities of inputs. The productivity of the cows is also reflected by the direction of impact of herd management variables. Whereas efficient farms replace cows more frequently, enroll heifers in production at an earlier age, and have shorter calving intervals, they also report higher frequency of disorder treatments. The average estimated energy corrected milk loss per cow is 1036, 451 and 242 kg for low, medium and high efficient farms. The study demonstrates the benefit of the stochastic frontier production function involving the estimation of individual technical efficiencies to evaluate farm performance and investigate the source of inefficiency. (C) 2004 Elsevier B.V. All rights reserved.

Prediction of mortality using dobutamine echocardiography

OBJECTIVES We sought to find out whether dobutamine echocardiography (DbE) could provide independent prediction of total and cardiac mortality, incremental to clinical and angiographic variables. BACKGROUND Existing outcome studies with DbE have examined composite end points, rather than death, over a relatively short follow-up. METHODS Clinical and stress data were collected in 3,156 patients (age 63 +/- 12 years, 1,801 men) undergoing DbE. Significant stenoses (>50% diameter) were identified in 70% of 1,073 patients undergoing coronary angiography. Total and cardiac mortality were identified over nine years of follow-up (mean 3.8 +/- 1.9). Cox models were used to analyze the effect of ischemia and other variables, independent of other determinants of mortality. RESULTS The dobutamine echocardiogram was abnormal in 1,575 patients (50%). Death occurred in 716 patients (23%), 259 of whom (8%) were thought to have died from cardiac causes. Patients with normal DbE had a total mortality of 8% per year and a cardiac mortality of 1% per year over the first four years of follow-up. Ischemia and the extent of abnormal wall motion were independent predictors of cardiac death, together with age and heart failure. In sequential Cox models, the predictive power of clinical data alone (model chi-square 115) was strengthened by adding the resting left ventricular function (model chi-square 138) and the results of DbE (model chi-square 181). In the subgroup undergoing coronary angiography, the power of the model was increased to a minor degree by the addition of coronary anatomy data. CONCLUSIONS Dobutamine echocardiography is an independent predictor of death, incremental to other data. While a normal dobutamine echocardiogram predicts low risk of cardiac death ton the order of 1% per year), this risk increases with the extent of abnormal wall motion at rest and stress, (J Am Coil Cardiol 2001;37:754-60) (C) 2001 by the American College of Cardiology.

Modelling nucleation in wet granulation

Nucleation is the first stage in any granulation process where binder liquid first comes into contact with the powder. This paper investigates the nucleation process where binder liquid is added to a fine powder with a spray nozzle. The dimensionless spray flux approach of Hapgood et al. (Powder Technol. 141 (2004) 20) is extended to account for nonuniform spray patterns and allow for overlap of nuclei granules rather than spray drops. A dimensionless nuclei distribution function which describes the effects of the design and operating parameters of the nucleation process (binder spray characteristics, the nucleation area ratio between droplets and nuclei and the powder bed velocity) on the fractional surface area coverage of nuclei on a moving powder bed is developed. From this starting point, a Monte Carlo nucleation model that simulates full nuclei size distributions as a function of the design and operating parameters that were implemented in the dimensionless nuclei distribution function is developed. The nucleation model was then used to investigate the effects of the design and operating parameters on the formed nuclei size distributions and to correlate these effects to changes of the dimensionless nuclei distribution function. Model simulations also showed that it is possible to predict nuclei size distributions beyond the drop controlled nucleation regime in Hapgood's nucleation regime map. Qualitative comparison of model simulations and experimental nucleation data showed similar shapes of the nuclei size distributions. In its current form, the nucleation model can replace the nucleation term in one-dimensional population balance models describing wet granulation processes. Implementation of more sophisticated nucleation kinetics can make the model applicable to multi-dimensional population balance models.

Change in woolgrower profits from adopting staple strength-enhancing management strategies

Wool tenderness is a significant problem in Australia, especially in areas where sheep graze under highly seasonal conditions. In this study, a profit function model is specified, estimated and simulated to assess the economic impact of staple strength-enhancing research on the profits of Australian woolgrowers. The model is based on a number of fundamental characteristics of the Australian wool industry and the staple-strength enhancing technology being assessed. The model consists of a system of demand and supply equations that are specified in terms of effective, rather than actual, prices. The interrelationships between the inputs and outputs are allowed for in the model in a manner that is consistent with theoretical restrictions. The adoption of the new feed management strategy results in a 4.4% increase in the expected profits of Australian wool producers in the short-run, and a 2.2% increase in expected profits in the long-run.

A model for predicting the future incidence of coronary heart disease within percentiles of coronary heart disease risk

Background We present a method (The CHD Prevention Model) for modelling the incidence of fatal and nonfatal coronary heart disease (CHD) within various CHD risk percentiles of an adult population. The model provides a relatively simple tool for lifetime risk prediction for subgroups within a population. It allows an estimation of the absolute primary CHD risk in different populations and will help identify subgroups of the adult population where primary CHD prevention is most appropriate and cost-effective. Methods The CHD risk distribution within the Australian population was modelled, based on the prevalence of CHD risk, individual estimates of integrated CHD risk, and current CHD mortality rates. Predicted incidence of first fatal and nonfatal myocardial infarction within CHD risk strata of the Australian population was determined. Results Approximately 25% of CHD deaths were predicted to occur amongst those in the top 10 percentiles of integrated CHD risk, regardless of age group or gender. It was found that while all causes survival did not differ markedly between percentiles of CHD risk before the ages of around 50-60, event-free survival began visibly to differ about 5 years earlier. Conclusions The CHD Prevention Model provides a means of predicting future CHD incidence amongst various strata of integrated CHD risk within an adult population. It has significant application both in individual risk counselling and in the identification of subgroups of the population where drug therapy to reduce CHD risk is most cost-effective. J Cardiovasc Risk 8:31-37 (C) 2001 Lippincott Williams & Wilkins.

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.

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.

Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation

Modeling physiological processes using tracer kinetic methods requires knowledge of the time course of the tracer concentration in blood supplying the organ. For liver studies, however, inaccessibility of the portal vein makes direct measurement of the hepatic dual-input function impossible in humans. We want to develop a method to predict the portal venous time-activity curve from measurements of an arterial time-activity curve. An impulse-response function based on a continuous distribution of washout constants is developed and validated for the gut. Experiments with simultaneous blood sampling in aorta and portal vein were made in 13 anesthetized pigs following inhalation of intravascular [O-15] CO or injections of diffusible 3-O[ C-11] methylglucose (MG). The parameters of the impulse-response function have a physiological interpretation in terms of the distribution of washout constants and are mathematically equivalent to the mean transit time ( T) and standard deviation of transit times. The results include estimates of mean transit times from the aorta to the portal vein in pigs: (T) over bar = 0.35 +/- 0.05 min for CO and 1.7 +/- 0.1 min for MG. The prediction of the portal venous time-activity curve benefits from constraining the regression fits by parameters estimated independently. This is strong evidence for the physiological relevance of the impulse-response function, which includes asymptotically, and thereby justifies kinetically, a useful and simple power law. Similarity between our parameter estimates in pigs and parameter estimates in normal humans suggests that the proposed model can be adapted for use in humans.

Angiotensin II facilitates autoregulation in the perfused mouse kidney: An optimized in vitro model for assessment of renal vascular and tubular function

Background and Aims: We have optimized the isolated perfused mouse kidney (IPMK) model for studying renal vascular and tubular function in vitro using 24-28 g C57BL6J mice; the wild type controls for many transgenic mice. Methods and Results: Buffer composition was optimized for bovine serum albumin concentration (BSA). The effect of adding erythrocytes on renal function and morphology was assessed. Autoregulation was investigated during stepped increases in perfusion pressure. Perfusion for 60 min at 90-110 mmHg with Krebs bicarbonate buffer containing 5.5% BSA, and amino acids produced functional parameters within the in vivo range. Erythrocytes increased renal vascular resistance (3.8 +/- 0.2 vs 2.4 +/- 0.1 mL/min.mmHg, P < 0.05), enhanced sodium reabsorption (FENa = 0.3 +/- 0.08 vs 1.5 +/- 0.7%, P < 0.05), produced equivalent glomerular filtration rates (GFR; 364 +/- 38 vs 400 +/- 9 muL/min per gkw) and reduced distal tubular cell injury in the inner stripe (5.8 +/- 1.7 vs 23.7 +/- 3.1%, P < 0.001) compared to cell free perfusion. The IPMK was responsive to vasoconstrictor (angiotensin II, EC50 100 pM) and vasodilator (methacholine, EC50 75 nM) mediators and showed partial autoregulation of perfusate flow under control conditions over 65-85 mmHg; autoregulatory index (ARI) of 0.66 +/- 0.11. Angiotensin II (100 pM) extended this range (to 65-120 mmHg) and enhanced efficiency (ARI 0.21 +/- 0.02, P < 0.05). Angiotensin II facilitation was antagonized by methacholine (ARI 0.76 +/- 0.08) and papaverine (ARI 0.91 +/- 0.13). Conclusion: The IPMK model is useful for studying renal physiology and pathophysiology without systemic neurohormonal influences.

Functional-structural plant modelling using a combination of architectural analysis, L-systems and a canonical model of function

This paper presents a new method for producing a functional-structural plant model that simulates response to different growth conditions, yet does not require detailed knowledge of underlying physiology. The example used to present this method is the modelling of the mountain birch tree. This new functional-structural modelling approach is based on linking an L-system representation of the dynamic structure of the plant with a canonical mathematical model of plant function. Growth indicated by the canonical model is allocated to the structural model according to probabilistic growth rules, such as rules for the placement and length of new shoots, which were derived from an analysis of architectural data. The main advantage of the approach is that it is relatively simple compared to the prevalent process-based functional-structural plant models and does not require a detailed understanding of underlying physiological processes, yet it is able to capture important aspects of plant function and adaptability, unlike simple empirical models. This approach, combining canonical modelling, architectural analysis and L-systems, thus fills the important role of providing an intermediate level of abstraction between the two extremes of deeply mechanistic process-based modelling and purely empirical modelling. We also investigated the relative importance of various aspects of this integrated modelling approach by analysing the sensitivity of the standard birch model to a number of variations in its parameters, functions and algorithms. The results show that using light as the sole factor determining the structural location of new growth gives satisfactory results. Including the influence of additional regulating factors made little difference to global characteristics of the emergent architecture. Changing the form of the probability functions and using alternative methods for choosing the sites of new growth also had little effect. (c) 2004 Elsevier B.V. All rights reserved.

Entanglement between a qubit and the environment in the spin-boson model

The quantitative description of the quantum entanglement between a qubit and its environment is considered. Specifically, for the ground state of the spin-boson model, the entropy of entanglement of the spin is calculated as a function of α, the strength of the ohmic coupling to the environment, and ɛ, the level asymmetry. This is done by a numerical renormalization group treatment of the related anisotropic Kondo model. For ɛ=0, the entanglement increases monotonically with α, until it becomes maximal for α→1-. For fixed ɛ>0, the entanglement is a maximum as a function of α for a value, α=αM