Double generalized linear model for tissue culture proportion data: a Bayesian perspective

Joint generalized linear models and double generalized linear models (DGLMs) were designed to model outcomes for which the variability can be explained using factors and/or covariates. When such factors operate, the usual normal regression models, which inherently exhibit constant variance, will under-represent variation in the data and hence may lead to erroneous inferences. For count and proportion data, such noise factors can generate a so-called overdispersion effect, and the use of binomial and Poisson models underestimates the variability and, consequently, incorrectly indicate significant effects. In this manuscript, we propose a DGLM from a Bayesian perspective, focusing on the case of proportion data, where the overdispersion can be modeled using a random effect that depends on some noise factors. The posterior joint density function was sampled using Monte Carlo Markov Chain algorithms, allowing inferences over the model parameters. An application to a data set on apple tissue culture is presented, for which it is shown that the Bayesian approach is quite feasible, even when limited prior information is available, thereby generating valuable insight for the researcher about its experimental results.

Parametric correlation functions to model the structure of permanent environmental (co)variances in milk yield random regression models

The objective of the present study was to estimate milk yield genetic parameters applying random regression models and parametric correlation functions combined with a variance function to model animal permanent environmental effects. A total of 152,145 test-day milk yields from 7,317 first lactations of Holstein cows belonging to herds located in the southeastern region of Brazil were analyzed. Test-day milk yields were divided into 44 weekly classes of days in milk. Contemporary groups were defined by herd-test-day comprising a total of 2,539 classes. The model included direct additive genetic, permanent environmental, and residual random effects. The following fixed effects were considered: contemporary group, age of cow at calving (linear and quadratic regressions), and the population average lactation curve modeled by fourth-order orthogonal Legendre polynomial. Additive genetic effects were modeled by random regression on orthogonal Legendre polynomials of days in milk, whereas permanent environmental effects were estimated using a stationary or nonstationary parametric correlation function combined with a variance function of different orders. The structure of residual variances was modeled using a step function containing 6 variance classes. The genetic parameter estimates obtained with the model using a stationary correlation function associated with a variance function to model permanent environmental effects were similar to those obtained with models employing orthogonal Legendre polynomials for the same effect. A model using a sixth-order polynomial for additive effects and a stationary parametric correlation function associated with a seventh-order variance function to model permanent environmental effects would be sufficient for data fitting.

Adsorption in mesopores: A molecular-continuum model with application to MCM-41

The classical model of surface layering followed by capillary condensation during adsorption in mesopores, is modified here by consideration of the adsorbate solid interaction potential. The new theory accurately predicts the capillary coexistence curve as well as pore criticality, matching that predicted by density functional theory. The model also satisfactorily predicts the isotherm for nitrogen adsorption at 77.4 K on MCM-41 material of various pore sizes, synthesized and characterized in our laboratory, including the multilayer region, using only data on the variation of condensation pressures with pore diameter. The results indicate a minimum mesopore diameter for the surface layering model to hold as 14.1 Å, below which size micropore filling must occur, and a minimum pore diameter for mechanical stability of the hemispherical meniscus during desorption as 34.2 Å. For pores in-between these two sizes reversible condensation is predicted to occur, in accord with the experimental data for nitrogen adsorption on MCM-41 at 77.4 K.

The Coefficient of Variance as an index of L2 lexical processing skill

The Coefficient of Variance (mean standard deviation/mean Response time) is a measure of response time variability that corrects for differences in mean Response time (RT) (Segalowitz & Segalowitz, 1993). A positive correlation between decreasing mean RTs and CVs (rCV-RT) has been proposed as an indicator of L2 automaticity and more generally as an index of processing efficiency. The current study evaluates this claim by examining lexical decision performance by individuals from three levels of English proficiency (Intermediate ESL, Advanced ESL and L1 controls) on stimuli from four levels of item familiarity, as defined by frequency of occurrence. A three-phase model of skill development defined by changing rCV-RT.values was tested. Results showed that RTs and CVs systematically decreased as a function of increasing proficiency and frequency levels, with the rCV-RT serving as a stable indicator of individual differences in lexical decision performance. The rCV-RT and automaticity/restructuring account is discussed in light of the findings. The CV is also evaluated as a more general quantitative index of processing efficiency in the L2.

A mathematical model of HIV transmission in NSW prisons

A mathematical model was developed to estimate HIV incidence in NSW prisons. Data included: duration of imprisonment; number of inmates using each needle; lower and higher number of shared injections per IDU per week; proportion of IDUs using bleach; efficacy of bleach; HIV prevalence and probability of infection. HIV prevalence in IDUs in prison was estimated to have risen from 0.8 to 5.7% (12.2%) over 180 weeks when using lower (and higher) values for frequency of shared injections. The estimated minimum (and maximum) number of IDU inmates infected with HIV in NSW prisons was 38 (and 152) in 1993 according to the model. These figures require confirmation by seroincidence studies. (C) 1998 Published by Elsevier Science Ireland Ltd. All rights reserved.

Modeling of hepatic elimination and organ distribution kinetics with the extended convection-dispersion model

The conventional convection-dispersion (also called axial 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. An extended form of the convection-dispersion model has been developed to adequately describe the outflow concentration-time profiles for vascular markers at both short and long times after bolus injections into perfused livers. The model, based on flux concentration and a convolution of catheters and large vessels, assumes that solute elimination in hepatocytes follows either fast distribution into or radial diffusion in hepatocytes. The model includes a secondary vascular compartment, postulated to be interconnecting sinusoids. Analysis of the mean hepatic transit time (MTT) and normalized variance (CV2) of solutes with extraction showed that the discrepancy between the predictions of MTT and CV2 for the extended and conventional models are essentially identical irrespective of the magnitude of rate constants representing permeability, volume, and clearance parameters, providing that there is significant hepatic extraction. In conclusion, the application of a newly developed extended convection-dispersion model has shown that the unweighted conventional convection-dispersion model can be used to describe the disposition of extracted solutes and, in particular, to estimate hepatic availability and clearance in booth experimental and clinical situations.

A Comparison Between Pulse Pressure Variation and Right End Diastolic Volume Index as Guides to Resuscitation in a Model of Hemorrhagic Shock in Pigs

Background: Different hemodynamic parameters including static indicators of cardiac preload as right ventricular end-diastolic volume index (RVEDVI) and dynamic parameters as pulse pressure variation (PPV) have been used in the decision-making process regarding volume expansion in critically ill patients. The objective of this study was to compare fluid resuscitation guided by either PPV or RVEDVI after experimentally induced hemorrhagic shock. Methods: Twenty-six anesthetized and mechanically ventilated pigs were allocated into control (group I), PPV (group II), or RVEDVI (group III) group. Hemorrhagic shock was induced by blood withdrawal to target mean arterial pressure of 40 mm Hg, maintained for 60 minutes. Parameters were measured at baseline, time of shock, 60 minutes after shock, immediately after resuscitation with hydroxyethyl starch 6% (130/0.4), 1 hour and 2 hours thereafter. The endpoint of fluid resuscitation was determined as the baseline values of PPV and RVEDVI. Statistical analysis of data was based on analysis of variance for repeated measures followed by the Bonferroni test (p < 0.05). Results: Volume and time to resuscitation were higher in group III than in group II (group III = 1,305 +/- 331 mL and group II = 965 +/- 245 mL, p < 0.05; and group III = 24.8 +/- 4.7 minutes and group II = 8.8 +/- 1.3 minutes, p < 0.05, respectively). All static and dynamic parameters and biomarkers of tissue oxygenation were affected by hemorrhagic shock and nearly all parameters were restored after resuscitation in both groups. Conclusion: In the proposed model of hemorrhagic shock, resuscitation to the established endpoints was achieved within a smaller amount of time and with less volume when guided by PPV than when guided by pulmonary artery catheter-derived RVEDVI.

Comparison Between Perfadex and Locally Manufactured Low-Potassium Dextran Solution for Pulmonary Preservation in an Ex Vivo Isolated Lung Perfusion Model

Introduction. Lung tranplantation, a consolidated treatment for end-stage lung disease, utilizes preservation solutions, such as low potassium dextran (LPD), to mitigate ischemia reperfusion injury. We sought the local development of LPD solutions in an attempt to facilitate access and enhance usage. We also sought to evaluate the effectiveness of a locally manufactured LPD solution in a rat model of ex vivo lung perfusion. Methods. We randomized the following groups \?\adult of male Wistar rats (n = 25 each): Perfadex (LPD; Vitro life, Sweden); locally manufactured LPD-glucose (LPDnac) (Farmoterapica, Brazil), and normal saline solution (SAL) with 3 ischemic times (6, 12, and 24 hours). The harvested heart lung blocks were flushed with solution at 4 C. After storage, the blocks were connected to an IL-2 Isolated Perfused Rat or Guinea Pig Lung System (Harvard Apparatus) and reperfused with homologous blood for 60 minutes. Respiratory mechanics, pulmonary artery pressure, perfusate blood gas analysis, and lung weight were measured at 10-minute intervals. Comparisons between groups and among ischemic times were performed using analysis of variance with a 5% level of significance. Results. Lungs preserved for 24 hours were nonviable and therefore excluded from the analysis. Those preserved for 6 hours showed better ventilatory mechanics when compared with 12 hours. The oxygenation capacity was not different between lungs flushed with LPD or LPDnac, regardless of the ischemic time. SAL lungs showed higher PCO(2) values than the other solutions. Lung weight increased over time during perfusion; however, there were no significant differences among the tested solutions (LPD, P = .23; LPDnac, P = .41; SAL, P = .26). We concluded that the LPDnac solution results in gas exchange were comparable to the original LPD (Perfadex); however ventilatory mechanics and edema formation were better with LPD, particularly among lungs undergoing 6 hours of cold ischemia.

A mathematical model for Escherichia coli debris size reduction during high pressure homogenisation based on grinding theory

Experimental data for E. coli debris size reduction during high-pressure homogenisation at 55 MPa are presented. A mathematical model based on grinding theory is developed to describe the data. The model is based on first-order breakage and compensation conditions. It does not require any assumption of a specified distribution for debris size and can be used given information on the initial size distribution of whole cells and the disruption efficiency during homogenisation. The number of homogeniser passes is incorporated into the model and used to describe the size reduction of non-induced stationary and induced E. coil cells during homogenisation. Regressing the results to the model equations gave an excellent fit to experimental data ( > 98.7% of variance explained for both fermentations), confirming the model's potential for predicting size reduction during high-pressure homogenisation. This study provides a means to optimise both homogenisation and disc-stack centrifugation conditions for recombinant product recovery. (C) 1997 Elsevier Science Ltd.

Model for Differential Nursing Diagnosis of Alterations in Urinary Elimination Based on Fuzzy Logic

Nursing diagnoses associated with alterations of urinary elimination require different interventions, Nurses, who are not specialists, require support to diagnose and manage patients with disturbances of urine elimination. The aim of this study was to present a model based on fuzzy logic for differential diagnosis of alterations in urinary elimination, considering nursing diagnosis approved by the North American Nursing Diagnosis Association, 2001-2002. Fuzzy relations and the maximum-minimum composition approach were used to develop the system. The model performance was evaluated with 195 cases from the database of a previous study, resulting in 79.0% of total concordance and 19.5% of partial concordance, when compared with the panel of experts. Total discordance was observed in only three cases (1.5%). The agreement between model and experts was excellent (kappa = 0.98, P < .0001) or substantial (kappa = 0.69, P < .0001) when considering the overestimative accordance (accordance was considered when at least one diagnosis was equal) and the underestimative discordance (discordance was considered when at least one diagnosis was different), respectively. The model herein presented showed good performance and a simple theoretical structure, therefore demanding few computational resources.

Comparative study of pressure- and volume-controlled ventilation on pulse pressure variation in a model of hypovolaemia in rabbits

Background and objective: Dynamic indices represented by systolic pressure variation and pulse pressure variation have been demonstrated to be more accurate than filling pressures in predicting fluid responsiveness. However, the literature is scarce concerning the impact of different ventilatory modes on these indices. We hypothesized that systolic pressure variation or pulse pressure variation could be affected differently by volume-controlled ventilation and pressure-controlled ventilation in an experimental model, during normovolaemia and hypovolaemia. Method: Thirty-two anaesthetized rabbits were randomly allocated into four groups according to ventilatory modality and volaemic status where G1-ConPCV was the pressure-controlled ventilation control group, G2-HemPCV was associated with haemorrhage, G3-ConVCV was the volume-controlled ventilation control group and G4-HemVCV was associated with haemorrhage. In the haemorrhage groups, blood was removed in two stages: 15% of the estimated blood volume withdrawal at M1, and, 30 min later, an additional 15% at M2. Data were submitted to analysis of variance for repeated measures; a value of P < 0.05 was considered to be statistically significant. Results: At MO (baseline), no significant differences were observed among groups. At M1, dynamic parameters differed significantly among the control and hypovolaemic groups (P < 0.05) but not between ventilation modes. However, when 30% of the estimated blood volume was removed (M2), dynamic parameters became significantly higher in animals under volume-controlled ventilation when compared with those under pressure-controlled ventilation. Conclusions: Under normovolaemia and moderate haemorrhage, dynamic parameters were not influenced by either ventilatory modalities. However, in the second stage of haemorrhage (30%), animals in volume-controlled ventilation presented higher values of systolic pressure variation and pulse pressure variation when compared with those submitted to pressure-controlled ventilation.

Antenatal steroid and tracheal occlusion restore vascular endothelial growth factor receptors in congenital diaphragmatic hernia rat model

OBJECTIVE: Investigate the effects of antenatal steroids and tracheal occlusion on pulmonary expression of vascular endothelial growth factor receptors in rats with nitrofen-induced congenital diaphragmatic hernia. STUDY DESIGN: Fetuses were exposed to nitrofen at embryonic day 9.5. Subgroups received dexamethasone or were operated on for tracheal occlusion, or received combined treatment. Morphologic variables were recorded. To analyze vascular endothelial growth factor receptor 1 and vascular endothelial growth factor receptor 2 expression, we performed Western blotting and immunohistochemistry. Morphologic variables were analyzed by analysis of variance and immunohistochemistry by Kruskal-Wallis test. RESULTS: Congenital diaphragmatic hernia decreased body weight, total lung weight, and lung-to-body weight ratio. Tracheal occlusion increased total lung weight and lung-to-body weight ratio (P < .05). Fetuses with congenital diaphragmatic hernia had reduced vascular endothelial growth factor receptor 1 and vascular endothelial growth factor receptor 2 expression, whereas steroids and tracheal occlusion increased their expression. Combined treatment increased expression of receptors, but had no additive effect. CONCLUSION: Vascular endothelial growth factor signaling disruption may be associated with pulmonary hypertension in congenital diaphragmatic hernia. Tracheal occlusion and steroids provide a pathway for restoring expression of vascular endothelial growth factor receptors.

Heart rate and arterial pressure variability in the experimental renovascular hypertension model in rats

This study was conducted in one kidney, one clip (1K1C) Goldblatt hypertensive rats to evaluate vascular and cardiac autonomic control using different approaches: 1) evaluation of the autonomic modulation of heart rate (HR) and systolic arterial pressure (SAP) by means of autoregressive power spectral analysis 2) assessment of the cardiac baroreflex sensitivity; and 3) double blockade with methylatropine and propranolol. The 1K1C group developed hypertension and tachycardia. The 1K1C group also presented reduction in variance as well as in LF (0.23 +/- 0.1 vs. 1.32 +/- 0.2 ms(2)) and HF (6.6 +/- 0.49 vs. 15.1 +/- 0.61 ms(2)) oscillations of pulse interval. Autoregressive spectral analysis of SAP showed that 1K1C rats had an increase in variance and LF band (13.3 +/- 2.7 vs. 7.4 +/- 1.01 mmHg(2)) in comparison with the sham group. The baroreflex gain was attenuated in the hypertensive 1K1C (- 1.83 +/- 0.05 bpm/mmHg) rats in comparison with normotensive sham (-3.23 +/- 0.06 bpm/MmHg) rats. The autonomic blockade caused an increase in the intrinsic HR and sympathetic predominance on the basal HR of 1K1C rats. Overall, these data indicate that the tachycardia observed in the 1K1C group may be attributed to intrinsic cardiac mechanisms (increased intrinsic heart rate) and to a shift in the sympathovagal balance towards cardiac sympathetic over-activity and vagal suppression associated to depressed baroreflex sensitivity. Finally, the increase in the LF components of SAP also suggests an increase in sympathetic activity to peripheral vessels. (c) 2008 Elsevier B.V. All rights reserved.

Habitat complexity, spatial interference, and "minimum risk distribution": a framework for population stability

In the past century, the debate over whether or not density-dependent factors regulate populations has generally focused on changes in mean population density, ignoring the spatial variance around the mean as unimportant noise. In an attempt to provide a different framework for understanding population dynamics based on individual fitness, this paper discusses the crucial role of spatial variability itself on the stability of insect populations. The advantages of this method are the following: (1) it is founded on evolutionary principles rather than post hoc assumptions; (2) it erects hypotheses that can be tested; and (3) it links disparate ecological schools, including spatial dynamics, behavioral ecology, preference-performance, and plant apparency into an overall framework. At the core of this framework, habitat complexity governs insect spatial variance. which in turn determines population stability. First, the minimum risk distribution (MRD) is defined as the spatial distribution of individuals that results in the minimum number of premature deaths in a population given the distribution of mortality risk in the habitat (and, therefore, leading to maximized population growth). The greater the divergence of actual spatial patterns of individuals from the MRD, the greater the reduction of population growth and size from high, unstable levels. Then, based on extensive data from 29 populations of the processionary caterpillar, Ochrogaster lunifer, four steps are used to test the effect of habitat interference on population growth rates. (1) The costs (increasing the risk of scramble competition) and benefits (decreasing the risk of inverse density-dependent predation) of egg and larval aggregation are quantified. (2) These costs and benefits, along with the distribution of resources, are used to construct the MRD for each habitat. (3) The MRD is used as a benchmark against which the actual spatial pattern of individuals is compared. The degree of divergence of the actual spatial pattern from the MRD is quantified for each of the 29 habitats. (4) Finally, indices of habitat complexity are used to provide highly accurate predictions of spatial divergence from the MRD, showing that habitat interference reduces population growth rates from high, unstable levels. The reason for the divergence appears to be that high levels of background vegetation (vegetation other than host plants) interfere with female host-searching behavior. This leads to a spatial distribution of egg batches with high mortality risk, and therefore lower population growth. Knowledge of the MRD in other species should be a highly effective means of predicting trends in population dynamics. Species with high divergence between their actual spatial distribution and their MRD may display relatively stable dynamics at low population levels. In contrast, species with low divergence should experience high levels of intragenerational population growth leading to frequent habitat-wide outbreaks and unstable dynamics in the long term. Six hypotheses, erected under the framework of spatial interference, are discussed, and future tests are suggested.

Patchy populations in stochastic environments: Critical number of patches for persistence

We introduce a model for the dynamics of a patchy population in a stochastic environment and derive a criterion for its persistence. This criterion is based on the geometric mean (GM) through time of the spatial-arithmetic mean of growth rates. For the population to persist, the GM has to be greater than or equal to1. The GM increases with the number of patches (because the sampling error is reduced) and decreases with both the variance and the spatial covariance of growth rates. We derive analytical expressions for the minimum number of patches (and the maximum harvesting rate) required for the persistence of the population. As the magnitude of environmental fluctuations increases, the number of patches required for persistence increases, and the fraction of individuals that can be harvested decreases. The novelty of our approach is that we focus on Malthusian local population dynamics with high dispersal and strong environmental variability from year to year. Unlike previous models of patchy populations that assume an infinite number of patches, we focus specifically on the effect that the number of patches has on population persistence. Our work is therefore directly relevant to patchily distributed organisms that are restricted to a small number of habitat patches.