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.

Validation of a root water uptake model to estimate transpiration constraints

Experimental results obtained from a greenhouse trial with common bean (Phaseolus vulgaris L) plants performed to test model hypotheses regarding the onset of limiting hydraulic conditions and the shape of the transpiration reduction curve in the falling rate phase are presented. According to these hypotheses based on simulations with an upscaled single-root model, the matric flux potential at the onset of limiting hydraulic conditions is as a function of root length density and potential transpiration rate, while the relative transpiration in the falling rate phase equals the relative matric flux potential. Transpiration of bean plants in water stressed pots with four different soils was determined daily by weighing and compared to values obtained from non-stressed pots. This procedure allowed determining the onset of the falling rate phase and corresponding soil hydraulic conditions. At the onset of the falling rate phase, the value of matric flux potential M(I) showed to differ in order of magnitude from the model predicted value for three out of four soils. This difference between model and experiment can be explained by the heterogeneity of the root distribution which is not considered by the model. An empirical factor to deal with this heterogeneity should be included in the model to improve predictions. Comparing the predictions of relative transpiration in the falling rate phase using a linear shape with water content, pressure head or matric flux potential, the matric flux potential based reduction function, in agreement with the hypothesis, showed the best performance, while the pressure head based equation resulted in the highest deviations between observed and predicted values of relative transpiration rates. (C) 2010 Elsevier B.V. All rights reserved.

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.

A Log-Linear Regression Model for the Beta-Weibull Distribution

We introduce the log-beta Weibull regression model based on the beta Weibull distribution (Famoye et al., 2005; Lee et al., 2007). We derive expansions for the moment generating function which do not depend on complicated functions. The new regression model represents a parametric family of models that includes as sub-models several widely known regression models that can be applied to censored survival data. We employ a frequentist analysis, a jackknife estimator, and a parametric bootstrap for the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes, and censoring percentages, several simulations are performed. In addition, the empirical distribution of some modified residuals are displayed and compared with the standard normal distribution. These studies suggest that the residual analysis usually performed in normal linear regression models can be extended to a modified deviance residual in the proposed regression model applied to censored data. We define martingale and deviance residuals to evaluate the model assumptions. The extended regression model is very useful for the analysis of real data and could give more realistic fits than other special regression models.

A split-pot experiment with sorghum to test a root water uptake partitioning model

Correct modeling of root water uptake partitioning over depth is an important issue in hydrological and crop growth models. Recently a physically based model to describe root water uptake was developed at single root scale and upscaled to the root system scale considering a homogeneous distribution of roots per soil layer. Root water uptake partitioning is calculated over soil layers or compartments as a function of respective soil hydraulic conditions, specifically the soil matric flux potential, root characteristics and a root system efficiency factor to compensate for within-layer root system heterogeneities. The performance of this model was tested in an experiment performed in two-compartment split-pot lysimeters with sorghum plants. The compartments were submitted to different irrigation cycles resulting in contrasting water contents over time. The root system efficiency factor was determined to be about 0.05. Release of water from roots to soil was predicted and observed on several occasions during the experiment; however, model predictions suggested root water release to occur more often and at a higher rate than observed. This may be due to not considering internal root system resistances, thus overestimating the ease with which roots can act as conductors of water. Excluding these erroneous predictions from the dataset, statistical indices show model performance to be of good quality.

Study of lymphocyte subpopulations in bone marrow in a model of protein-energy malnutrition

Objective: Protein-energy malnutrition (PEM) is an important public health problem affecting millions of people worldwide. Hematopoietic tissue requires a high nutrient supply, and a reduction in leukocytes, especially lymphocytes, suggests that some nutritional deficiencies might be altering bone marrow function and decreasing its ability to produce lymphocytes. In this study, we evaluated the effect that PEM has on lymphocyte subtypes and the cell cycle of CD5(+) cells. Methods: Swiss mice were subjected to PEM using a low-protein diet containing 4% protein. When the experimental group had lost about 20% of their original body weight, we collected blood and bone marrow cells and evaluated the hemogram, the myelogram, bone marrow lymphoid markers using flow cytometry, and the cell cycle in CD5(+) bone marrow. Results: Malnourished animals presented anemia, reticulocytopenia, and leukopenia with lymphopenia. The bone marrow was hypocellular, and flow cytometric analyses of bone marrow cells showed cells that were CD45(+) (91.2%), CD2(+) (84.9%), CD5(+) (37.3%), CD3(+) (23.5%), CD19(+) (43.3%), CD22(+) (34.7%), CD19(+)/CD2(+) (51.2%), CD19(+)/CD3(+)(24.0%), CD19(+)/CD5(+) (13.2%), CD22(+)/CD2(+) (40.1%), CD22(+)/CD3(+) (30.3%), and CD22(+)/CD5(+) (1.1%) in malnourished animals and CD45(+) (97.5%), CD2(+) (42.9%), CD5(+) (91.5%), CD3(+) (92.0%), CD19(+) (52.0%), CD22(+) (75.6%), CD19(+)/CD2(+) (62.0%), CD19(+)/CD3(+) (55.4%), CD19(+)/CO5(+) (6.7%), CD22(+)/CD2(+) (70.3%), CD22(+)/CD3(+) (55.9%), and CD22(+)/ CD5(+) (8.4%) in control animals. Malnourished animals also presented more CD5(+) cells in the G0 phase of cell cycle development. Conclusion: Malnourished animals presented bone marrow hypoplasia, maturation interruption, prominent lymphopenia with depletion in the lymphoid lineage, and changes in cellular development. We suggest that these changes are some of the primary causes of lymphopenia in cases of PEM and partly explain the increase in susceptibility to infections found in malnourished individuals. Published by Elsevier Inc.

Stability and in vitro penetration study of rutin incorporated in a cosmetic emulsion through an alternative model biomembrane

Rutin is employed as antioxidant and to prevent the capillary fragility and, when incorporated in cosmetic emulsions, it must target the action site. In vitro cutaneous penetration studies through human skin is the ideal situation, however, there are difficulties to obtain and to maintain this tissue viability. Among the membrane models, shed snake skin presents itself as pure stratum corneum, providing barrier function similar to human and it is obtained without the animal sacrifice. The objectives of this research were the development and stability evaluation of a cosmetic emulsion containing rutin and propylene glycol (penetration enhancer) and the evaluation or rutin in vitro cutaneous penetration and retention from the emulsion, employing an alternative model biomembrane. Emulsion was developed with rutin and propylene glycol, both at 5.0% w/w. Active substance presented on the formulation was quantified by a validated spectrophotometric method at 361.0 nm. Rutin Rutin cutaneous penetration and retention was performed in vertical diffusion cells with shed snake skin of Crotalus durissus, as alternative model biomembrane, and distilled water and ethanol 99.5% (1:1), as receptor fluid. The experiment was conducted for six hours, at 37.0 +/- 0.5 degrees C with constant stirring of 300 rpm. Spectrophotometry at 410.0 nm, previously validated, determined the active substance after cutaneous penetration/ retention. Emulsion did not promote rutin cutaneous penetration through C. durissus skin, retaining 0.931 +/- 0.0391 mu g rutin/mg shed snake skin. The referred formulation was chemically stable for 30 days after stored at 25.0 +/- 2.0 degrees C, 5.0 +/- 0.5 degrees C and 45.0 +/- 0.5 degrees C. In conclusion, it has not been verified the active cutaneous penetration through the model biomembrane, but only its retention on the Crotalus durissus stratum corneum, condition considered stable for 30 days.

Excised Porcine Cornea Integrity Evaluation in an in vitro Model of Iontophoretic Ocular Research

Background/Aims: It is a challenge to adapt traditional in vitro diffusion experiments to ocular tissue. Thus, the aim of this work was to present experimental evidence on the integrity of the porcine cornea, barrier function and maintenance of electrical properties for 6 h of experiment when the tissue is mounted on an inexpensive and easy-to-use in vitro model for ocular iontophoresis. Methods: A modified Franz diffusion cell containing two ports for the insertion of the electrodes and a receiving compartment that does not need gassing with carbogen was used in the studies. Corneal electron transmission microscopy images were obtained, and diffusion experiments with fluorescent markers were performed to examine the integrity of the barrier function. The preservation of the negatively charged corneal epithelium was verified by the determination of the electro-osmotic flow of a hydrophilic and non-ionized molecule. Results: The diffusion cell was able to maintain the temperature, homogenization, porcine epithelial corneal structure integrity, barrier function and electrical characteristics throughout the 6 h of permeation experiment, without requiring CO(2) gassing when the receiving chamber was filled with 25 m M of HEPES buffer solution. Conclusion: The system described here is inexpensive, easy to handle and reliable as an in vitro model for iontophoretic ocular delivery studies. Copyright (C) 2010 S. Karger AG, Basel

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

Weibull model selection for reliability modelling

A large number of models have been derived from the two-parameter Weibull distribution and are referred to as Weibull models. They exhibit a wide range of shapes for the density and hazard functions, which makes them suitable for modelling complex failure data sets. The WPP and IWPP plot allows one to determine in a systematic manner if one or more of these models are suitable for modelling a given data set. This paper deals with this topic.

Integrability of the Russian doll BCS model

We show that integrability of the BCS model extends beyond Richardson's model (where all Cooper pair scatterings have equal coupling) to that of the Russian doll BCS model for which the couplings have a particular phase dependence that breaks time-reversal symmetry. This model is shown to be integrable using the quantum inverse scattering method, and the exact solution is obtained by means of the algebraic Bethe ansatz. The inverse problem of expressing local operators in terms of the global operators of the monodromy matrix is solved. This result is used to find a determinant formulation of a correlation function for fluctuations in the Cooper pair occupation numbers. These results are used to undertake exact numerical analysis for small systems at half-filling.

A model for assembly and activation of the GM-CSF, IL-3 and IL-5 receptors: Insights from activated mutants of the common beta subunit

Granulocyte-macrophage colony stimulating factor (GM-CSF), Interleukin-3 (IL-3) and Interleukin-5 (IL-5) have overlapping, pleiotropic effects on hematopoietic cells, including neutrophils, eosinophils, monocytes and early progenitor cells. The high-affinity receptors for human GM-CSF, IL-3, and IL-5 share a common beta-subunit (h beta(c)), which is essential for signalling and plays a major role in recruiting intracellular signalling molecules. While activation of the cytoplasmic tyrosine kinase JAK2 appears to be the initiating event for signalling, the immediate events that trigger this are still unclear. We have isolated a number of activated mutants of h beta(c), which can be grouped into classes defined by their state of receptor phosphorylation, their requirement for alpha subunit as a cofactor, and their activities in primary cells and cell lines. We discuss these findings with regard to the stoichiometry, activation, and signalling of the normal GM-CSF/IL-3/IL-5 receptor complexes. Specifically, this work has implications for the role of the ligand-specific alpha-subunits in initiating the signalling through the beta-subunit, the role of beta subunit dimerization as a receptor trigger, and the function of receptor tyrosine phosphorylation in generating growth and survival signals. Based on the properties of the activated mutants and the recent structures of erythropoietin receptor (Epo-R) complexes, we propose a model in which (1) activation of h beta(c) can occur via alternative states that differ with respect to stoichiometry and subunit assembly, but which all mediate proliferative responses, and (2) each of the different classes of activated mutants mimics one of these alternative states. (C) 2000 International Society for Experimental Hematology. Published by Elsevier Science Inc.

Model for gas-liquid flow through wet porous medium

A method involving bubbling of air through a fibrous filter immersed in water has recently been investigated (Agranovski et al. [1]). Experimental results showed that the removal efficiency for ultra-fine aerosols by such filters was greatly increased compared to dry filters. Nuclear Magnetic Resonance (NMR) imaging was used to examine the wet filter and to determine the nature of the gas flow inside the filter (Agranovski et al. [2]). It was found that tortuous preferential pathways (or flow tubes) develop within the filter through which the air flows and the distribution of air and water inside the porous medium has been investigated. The aim of this paper is to investigate the geometry of the pathways and to make estimates of the flow velocities and particle removal efficiency in such pathways. A mathematical model of the flow of air along the preferred pathways has been developed and verified experimentally. Even for the highest realistic gas velocity the flow field was essentially laminar (Re approximate to 250). We solved Laplace's equation for stream function to map trajectories of particles and gas molecules to investigate the possibility of their removal from the carrier.

CD40 ligation conditions dendritic cell antigen-presenting function through sustained activation of NF-kappa B

An understanding of the biochemical control of dendritic cell (DC) differentiation/activation is essential for improving T cell immunity by various immunotherapeutic approaches, including DC immunization. Ligation of CD40 enhances DC function, including conditioning for CTL priming. NF-kappaB, and particularly RelB, is an essential control pathway for myeloid DC differentiation. Furthermore, RelB regulates B cell Ag-presenting function. We hypothesized that CD40 ligand (CD40L) and TNF-alpha, which differ in their capacity to condition DC, would also differ in their capacity to activate NF-kappaB. DC differentiated for 2 days from monocytes in the presence of GM-CSF and IL-4 were used as a model, as NF-kappaB activity was constitutively low. The capacity of DC to activate T cells following CD40L treatment was enhanced compared with TNF-alpha treatment, and this was NF-kappaB dependent. Whereas RelB/p50 translocation induced by TNF-alpha was attenuated after 6 h, RelB/p50 nuclear translocation induced by CD40L was sustained for at least 24 h. The mechanism of this difference related to enhanced degradation of IkappaBalpha following CD40L stimulation. However, NF-kappaB activation induced by TNF-alpha could be sustained by blocking autocrine IL-10. These data indicate that NF-kappaB activation is essential for T cell activation by DC, and that this function is enhanced if DC NF-kappaB activation is prolonged. Because IL-10 moderates DC NF-kappaB activation by TNF-alpha, sustained NF-kappaB activation can be achieved by blocking IL-10 in the presence of stimuli that induce TNF-alpha.

Dormancy release in Lolium rigidum seeds is a function of thermal after-ripening time and seed water content

Dormancy release in seeds of Lolium rigidum Gaud. (annual ryegrass) was investigated in relation to temperature and seed water content. Freshly matured seeds were collected from cropping fields at Wongan Hills and Merredin, Western Australia. Seeds from Wongan Hills were equilibrated to water contents between 6 and 18% dry weight and after-ripened at constant temperatures between 9 and 50degreesC for up to 23 weeks. Wongan Hills and Merredin seeds at water contents between 7 and 17% were also after-ripened in full sun or shade conditions. Dormancy was tested at regular intervals during after-ripening by germinating seeds on agar at 12-h alternating 15degreesC (dark) and 25degreesC (light) periods. Rate of dormancy release for Wongan Hills seeds was a positive linear function of after-ripening temperature above a base temperature (T-b) of 5.4degreesC. A thermal after-ripening time model for dormancy loss accounting for seed moisture in the range 6-18% was developed using germination data for Wongan Hills seeds after-ripened at constant temperatures. The model accurately predicted dormancy release for Wongan Hills seeds after-ripened under naturally fluctuating temperatures. Seeds from Merredin responded similarly but had lower dormancy at collection and a faster rate of dormancy release in seeds below 9% water content.