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.

Dissolved organic carbon in rainwater from areas heavily impacted by sugar cane burning

This work reports on rainwater dissolved organic carbon (DOC) from Ribeirao Preto (RP) and Araraquara over a period of 3 years. The economies of these two cities, located in Sao Paulo state (Brazil), are based on agriculture and related industries, and the region is strongly impacted by the burning of sugar cane foliage before harvesting. Highest DOC concentrations were obtained when air masses traversed sugar cane fields burned on the same day as the rain event. Significant increases in the DOC volume weighted means (VWM) during the harvest period, for both sites, and a good linear correlation (r=0.83) between DOC and K (a biomass burning marker) suggest that regional scale organic carbon emissions prevail over long-range transport. The DOC VWMs and standard deviations were 272 +/- 22 mu mol L-1 (n=193) and 338 +/- 40 mu mol L-1 (n=80) for RP and Araraquara, respectively, values which are at least two times higher than those reported for other regions influenced by biomass burning, such as the Amazon. These high DOC levels are discussed in terms of agricultural activities, particularly the large usage of biogenic fuels in Brazil, as well as the analytical method used in this work, which includes volatile organic carbon when reporting DOC values. Taking into account rainfall volume, estimated annual rainwater DOC fluxes for RP (4.8 g C m(-2) yr(-1)) and Araraquara (5.4 g C m(-2) yr(-1)) were close to that previously found for the Amazon region (4.8 g C m(-2) yr(-1)). This work also discusses whether previous calculations of the global rainwater carbon flux may have been underestimated, since they did not consider large inputs from biomass combustion sources, and suffered from a possible analytical bias. (c) 2008 Elsevier Ltd. All rights reserved.

Loss of Ergodicity in the Transition from Annealed to Quenched Disorder in a Finite Kinetic Ising Model

We consider a kinetic Ising model which represents a generic agent-based model for various types of socio-economic systems. We study the case of a finite (and not necessarily large) number of agents N as well as the asymptotic case when the number of agents tends to infinity. The main ingredient are individual decision thresholds which are either fixed over time (corresponding to quenched disorder in the Ising model, leading to nonlinear deterministic dynamics which are generically non-ergodic) or which may change randomly over time (corresponding to annealed disorder, leading to ergodic dynamics). We address the question how increasing the strength of annealed disorder relative to quenched disorder drives the system from non-ergodic behavior to ergodicity. Mathematically rigorous analysis provides an explicit and detailed picture for arbitrary realizations of the quenched initial thresholds, revealing an intriguing ""jumpy"" transition from non-ergodicity with many absorbing sets to ergodicity. For large N we find a critical strength of annealed randomness, above which the system becomes asymptotically ergodic. Our theoretical results suggests how to drive a system from an undesired socio-economic equilibrium (e. g. high level of corruption) to a desirable one (low level of corruption).

Assessment of the differential linear coherent scattering coefficient of biological samples

New differential linear coherent scattering coefficient, mu(CS), data for four biological tissue types (fat pork, tendon chicken, adipose and fibroglandular human breast tissues) covering a large momentum transfer interval (0.07 <= q <= 70.5 nm(-1)), resulted from combining WAXS and SAXS data, are presented in order to emphasize the need to update the default data-base by including the molecular interference and the large-scale arrangements effect. The results showed that the differential linear coherent scattering coefficient demonstrates influence of the large-scale arrangement, mainly due to collagen fibrils for tendon chicken and fibroglandular breast samples, and triacylglycerides for fat pork and adipose breast samples at low momentum transfer region. While, at high momentum transfer, the mu(CS) reflects effects of molecular interference related to water for tendon chicken and fibroglandular samples and, fatty acids for fat pork and adipose samples. (C) 2009 Elsevier B.V. All rights reserved.

The association of prevalent medial meniscal pathology with cartilage loss in the medial tibiofemoral compartment over a 2-year period

Objective: To investigate the association of different types of magnetic resonance imaging (MRI)-detected medial meniscal pathology with subregional cartilage loss in the medial tibiofemoral compartment. Methods: A total of 152 women aged >= 40 years, with and without knee osteoarthritis (OA) were included in a longitudinal 24-month observational study. Spoiled gradient recalled acquisitions at steady state (SPGR) and T2-weighted fat-suppressed MRI sequences were acquired. Medial meniscal status of the anterior horn (AH), body, and posterior horn (PH) was graded at baseline: 0 (normal), 1 (intrasubstance meniscal signal changes), 2 (single tears), and 3 (complex tears/maceration). Cartilage segmentation was performed at baseline and 24-month follow-up in various tibiofemoral subregions using computation software. Multiple linear regression models were applied for the analysis with cartilage loss as the outcome. In a first model, the results were adjusted for age and body mass index (BMI). In a second model, the results were adjusted for age, BMI and medial meniscal extrusion. Results: After adjusting for age, BMI, and medial meniscal extrusion, cartilage loss in the total medial tibia (MT) (0.04 mm, P=0.04) and the external medial tibia (eMT) (0.068 mm, P=0.04) increased significantly for compartments with grade 3 lesions. Cartilage loss in the total central medial femoral condyle (cMF) (0.071 mm, P=0.03) also increased significantly for compartments with grade 2 lesions. Cartilage loss at the eMT was significantly related to tears of the PH (0.074 mm; P=0.03). Cartilage loss was not significantly increased for compartments with grade 1 lesions. Conclusion: The protective function of the meniscus appears to be preserved in the presence of intrasubstance meniscal signal changes. Prevalent single tears and meniscal maceration were found to be associated with increased cartilage loss in the same compartment, especially at the PH. (C) 2009 Osteoarthritis Research Society International. Published by Elsevier Ltd. All rights reserved.

Finite Element Analysis of Shear Versus Torsion Adhesive Strength Tests for Dental Resin Composites

Stress distributions in torsion and wire-loop shear tests were compared using three-dimensional (3-D) linear-elastic finite element method, in an attempt to predict the ideal conditions for testing adhesive strength of dental resin composites to dentin. The torsion test presented lower variability in stress concentration at the adhesive interface with changes in the proportion adhesive thickness/resin composite diameter, as well as lower variability with changes in the resin composite elastic modulus. Moreover, the torsion test eliminated variability from changes in loading distance, and reduced the cohesive fracture tendency in the dentin. The torsion test seems to be more appropriate than wire-loop shear test for testing the resin composite-tooth interface strength. (c) Koninklijke Brill NV, Leiden, 2009

Rossby wave propagation tracks in southern hemisphere mean basic flows associated to generalized frosts over southern South America

Based on previous observational studies on cold extreme events over southern South America, some recent studies suggest a possible relationship between Rossby wave propagation remotely triggered and the occurrence of frost. Using the concept of linear theory of Rossby wave propagation, this paper analyzes the propagation of such waves in two different basic states that correspond to austral winters with maximum and minimum generalized frost frequency of occurrence in the Wet Pampa (central-northwest Argentina). In order to determine the wave trajectories, the ray tracing technique is used in this study. Some theoretical discussion about this technique is also presented. The analysis of the basic state, from a theoretical point of view and based on the calculation of ray tracings, corroborates that remotely excited Rossby waves is the mechanism that favors the maximum occurrence of generalized frosts. The basic state in which the waves propagate is what conditions the places where they are excited. The Rossby waves are excited in determined places of the atmosphere, propagating towards South America along the jet streams that act as wave guides, favoring the generation of generalized frosts. In summary, this paper presents an overview of the ray tracing technique and how it can be used to investigate an important synoptic event, such as frost in a specific region, and its relationship with the propagation of large scale planetary waves.

Ground-based detection of sprites and their parent lightning flashes over Africa during the 2006 AMMA campaign

Sprites have been detected in video camera observations from Niger over mesoscale convective systems in Nigeria during the 2006 AMMA (African Monsoon Multidisciplinary Analysis) campaign The parent lightning flashes have been detected by multiple Extremely Low Frequency (ELF) receiving stations worldwide The recorded charge moments of the patent lightning flashes are often in excellent agreement between different receiving sites, and are furthermore consistent with conventional dielectric breakdown in the mesosphere as the origin of the sprites Analysis of the polarization of the horizontal magnetic field at the distant receivers provides evidence that the departure from linear magnetic polarization at ELF is caused primarily by the clay night asymmetry of the Earth-ionosphere cavity Copyright (C) 2009 Royal Meteorological Society

Are Reanalysis Data Useful for Calculating Climate Indices over South America?

Precipitation and temperature climate indices are calculated using the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis and validated against observational data from some stations over Brazil and other data sources. The spatial patterns of the climate indices trends are analyzed for the period 1961-1990 over South America. In addition, the correlation and linear regression coefficients for some specific stations were also obtained in order to compare with the reanalysis data. In general, the results suggest that NCEP/NCAR reanalysis can provide useful information about minimum temperature and consecutive dry days indices at individual grid cells in Brazil. However, some regional differences in the climate indices trends are observed when different data sets are compared. For instance, the NCEP/NCAR reanalysis shows a reversal signal for all rainfall annual indices and the cold night index over Argentina. Despite these differences, maps of the trends for most of the annual climate indices obtained from the NCEP/NCAR reanalysis and BRANT analysis are generally in good agreement with other available data sources and previous findings in the literature for large areas of southern South America. The pattern of trends for the precipitation annual indices over the 30 years analyzed indicates a change to wetter conditions over southern and southeastern parts of Brazil, Paraguay, Uruguay, central and northern Argentina, and parts of Chile and a decrease over southwestern South America. All over South America, the climate indices related to the minimum temperature (warm or cold nights) have clearly shown a warming tendency; however, no consistent changes in maximum temperature extremes (warm and cold days) have been observed. Therefore, one must be careful before suggesting an), trends for warm or cold days.

Are the magnetic fields of millisecond pulsars similar to 10(8) G?

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are similar to 10(8) G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The similar to 10(8)G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have < 10(11) G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit 10(11) G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict similar to 10(8) G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted similar to 10(8) G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields a parts per thousand << 10(8) G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields > 10(8) G with their magnetic dipole emission being balanced by low level accretion from their companions.

The analytic torsion of a cone over an odd dimensional manifold

We study the analytic torsion of a cone over an orientable odd dimensional compact connected Riemannian manifold W. We prove that the logarithm of the analytic torsion of the cone decomposes as the sum of the logarithm of the root of the analytic torsion of the boundary of the cone, plus a topological term, plus a further term that is a rational linear combination of local Riemannian invariants of the boundary. We show that this last term coincides with the anomaly boundary term appearing in the Cheeger Muller theorem [3, 2] for a manifold with boundary, according to Bruning and Ma (2006) [5]. We also prove Poincare duality for the analytic torsion of a cone. (C) 2010 Elsevier B.V. All rights reserved.

Finite-dimensional global attractors in Banach spaces

We provide bounds on the upper box-counting dimension of negatively invariant subsets of Banach spaces, a problem that is easily reduced to covering the image of the unit ball under a linear map by a collection of balls of smaller radius. As an application of the abstract theory we show that the global attractors of a very broad class of parabolic partial differential equations (semilinear equations in Banach spaces) are finite-dimensional. (C) 2010 Elsevier Inc. All rights reserved.

Invariant theory and reversible-equivariant vector fields

In this paper we present results for the systematic study of reversible-equivariant vector fields - namely, in the simultaneous presence of symmetries and reversing symmetries - by employing algebraic techniques from invariant theory for compact Lie groups. The Hilbert-Poincare series and their associated Molien formulae are introduced,and we prove the character formulae for the computation of dimensions of spaces of homogeneous anti-invariant polynomial functions and reversible-equivariant polynomial mappings. A symbolic algorithm is obtained for the computation of generators for the module of reversible-equivariant polynomial mappings over the ring of invariant polynomials. We show that this computation can be obtained directly from a well-known situation, namely from the generators of the ring of invariants and the module of the equivariants. (C) 2008 Elsevier B.V, All rights reserved.

A linear optimization approach to the combined production planning model

Two fundamental processes usually arise in the production planning of many industries. The first one consists of deciding how many final products of each type have to be produced in each period of a planning horizon, the well-known lot sizing problem. The other process consists of cutting raw materials in stock in order to produce smaller parts used in the assembly of final products, the well-studied cutting stock problem. In this paper the decision variables of these two problems are dependent of each other in order to obtain a global optimum solution. Setups that are typically present in lot sizing problems are relaxed together with integer frequencies of cutting patterns in the cutting problem. Therefore, a large scale linear optimizations problem arises, which is exactly solved by a column generated technique. It is worth noting that this new combined problem still takes the trade-off between storage costs (for final products and the parts) and trim losses (in the cutting process). We present some sets of computational tests, analyzed over three different scenarios. These results show that, by combining the problems and using an exact method, it is possible to obtain significant gains when compared to the usual industrial practice, which solve them in sequence. (C) 2010 The Franklin Institute. Published by Elsevier Ltd. All rights reserved.

A Direct Numerov Sixth-order Numerical Scheme to Accurately Solve the Unidimensional Poisson Equation with Dirichlet Boundary Conditions

In this article, we present an analytical direct method, based on a Numerov three-point scheme, which is sixth order accurate and has a linear execution time on the grid dimension, to solve the discrete one-dimensional Poisson equation with Dirichlet boundary conditions. Our results should improve numerical codes used mainly in self-consistent calculations in solid state physics.