Disequilibrium Coefficient: A Bayesian Perspective

Hardy-Weinberg Equilibrium (HWE) is an important genetic property that populations should have whenever they are not observing adverse situations as complete lack of panmixia, excess of mutations, excess of selection pressure, etc. HWE for decades has been evaluated; both frequentist and Bayesian methods are in use today. While historically the HWE formula was developed to examine the transmission of alleles in a population from one generation to the next, use of HWE concepts has expanded in human diseases studies to detect genotyping error and disease susceptibility (association); Ryckman and Williams (2008). Most analyses focus on trying to answer the question of whether a population is in HWE. They do not try to quantify how far from the equilibrium the population is. In this paper, we propose the use of a simple disequilibrium coefficient to a locus with two alleles. Based on the posterior density of this disequilibrium coefficient, we show how one can conduct a Bayesian analysis to verify how far from HWE a population is. There are other coefficients introduced in the literature and the advantage of the one introduced in this paper is the fact that, just like the standard correlation coefficients, its range is bounded and it is symmetric around zero (equilibrium) when comparing the positive and the negative values. To test the hypothesis of equilibrium, we use a simple Bayesian significance test, the Full Bayesian Significance Test (FBST); see Pereira, Stern andWechsler (2008) for a complete review. The disequilibrium coefficient proposed provides an easy and efficient way to make the analyses, especially if one uses Bayesian statistics. A routine in R programs (R Development Core Team, 2009) that implements the calculations is provided for the readers.

Qualitative Analysis of the Interdisciplinary Interaction between Data Analysis Specialists and Novice Clinical Researchers

Background: The inherent complexity of statistical methods and clinical phenomena compel researchers with diverse domains of expertise to work in interdisciplinary teams, where none of them have a complete knowledge in their counterpart's field. As a result, knowledge exchange may often be characterized by miscommunication leading to misinterpretation, ultimately resulting in errors in research and even clinical practice. Though communication has a central role in interdisciplinary collaboration and since miscommunication can have a negative impact on research processes, to the best of our knowledge, no study has yet explored how data analysis specialists and clinical researchers communicate over time. Methods/Principal Findings: We conducted qualitative analysis of encounters between clinical researchers and data analysis specialists (epidemiologist, clinical epidemiologist, and data mining specialist). These encounters were recorded and systematically analyzed using a grounded theory methodology for extraction of emerging themes, followed by data triangulation and analysis of negative cases for validation. A policy analysis was then performed using a system dynamics methodology looking for potential interventions to improve this process. Four major emerging themes were found. Definitions using lay language were frequently employed as a way to bridge the language gap between the specialties. Thought experiments presented a series of ""what if'' situations that helped clarify how the method or information from the other field would behave, if exposed to alternative situations, ultimately aiding in explaining their main objective. Metaphors and analogies were used to translate concepts across fields, from the unfamiliar to the familiar. Prolepsis was used to anticipate study outcomes, thus helping specialists understand the current context based on an understanding of their final goal. Conclusion/Significance: The communication between clinical researchers and data analysis specialists presents multiple challenges that can lead to errors.

Selection of suitable housekeeping genes for expression analysis in glioblastoma using quantitative RT-PCR

Background: Considering the broad variation in the expression of housekeeping genes among tissues and experimental situations, studies using quantitative RT-PCR require strict definition of adequate endogenous controls. For glioblastoma, the most common type of tumor in the central nervous system, there was no previous report regarding this issue. Results: Here we show that amongst seven frequently used housekeeping genes TBP and HPRT1 are adequate references for glioblastoma gene expression analysis. Evaluation of the expression levels of 12 target genes utilizing different endogenous controls revealed that the normalization method applied might introduce errors in the estimation of relative quantities. Genes presenting expression levels which do not significantly differ between tumor and normal tissues can be considered either increased or decreased if unsuitable reference genes are applied. Most importantly, genes showing significant differences in expression levels between tumor and normal tissues can be missed. We also demonstrated that the Holliday Junction Recognizing Protein, a novel DNA repair protein over expressed in lung cancer, is extremely over-expressed in glioblastoma, with a median change of about 134 fold. Conclusion: Altogether, our data show the relevance of previous validation of candidate control genes for each experimental model and indicate TBP plus HPRT1 as suitable references for studies on glioblastoma gene expression.

Influence of Laser Photobiomodulation on Collagen IV During Skeletal Muscle Tissue Remodeling After Injury in Rats

Objective: The aim of the present study was to determine the effect of GaAlAs low-level laser therapy (LLLT) on collagen IV remodeling of the tibialis anterior (TA) muscle in rats after cryolesion. Background: Considerable interest exists in skeletal muscle regeneration in situations such as repair after exercise-induced muscle injury, after muscle transplantation, in muscular dystrophy, exercise-induced muscle injury, and the recovery of strength after atrophy due to disuse. A number of studies have demonstrated the potential of LLLT in facilitating the muscle-healing process; however, no consensus is found in the literature regarding the best laser-irradiation parameters. Methods: Adult male Wistar rats (n = 45) were used and randomly divided into three groups: control (n = 5); nontreated cryolesioned group (n = 20), and LLLT-cryolesioned group (n = 20). The cryolesioned groups were analyzed at 1, 7, 14, and 21 days after the injury procedure. Laser irradiation was performed 3 times per week on the injured region by using the GaAlAs laser (660 nm; beam spot of 0.04 cm(2), output power of 20 mW, power density of 500 mW/cm(2), and energy density of 5 J/cm(2), for 10 sec). The muscles were removed, frozen, cryosectioned, and then stained with hematoxylin-eosin for the visualization of general morphology or used for immunohistochemical analysis of collagen IV. Results: It was demonstrated that LLLT promotes an increase in collagen IV immunolabeling in skeletal muscle in the first 7 days after acute trauma caused by cryoinjury, but does not modify the duration of the tissue-repair process. Even with LLLT, the injured muscle tissue needs similar to 21 days to achieve the same state of organization as that in the noninjured muscle. Conclusion: The collagen IV content is modulated in regenerating skeletal muscle under LLLT, which might be associated with better tissue outcome, although the histologic analysis did not detect tissue improvement in the LLLT group.

Polarized Light (lambda 400-2000 nm) on Third-Degree Burns in Diabetic Rats: Immunohistochemical Study

Aim: The aim of this study was to evaluate with light microscopy the healing process of third-degree burns on diabetic rats treated with polarized light (lambda 400-2000 nm, 20 or 40 J/cm(2)/session, 40 mW/cm(2), 2.4 J/cm(2)/min, 5.5-cm beam diameter). Background: Uncontrolled diabetes mellitus causes severe disruption of the body's metabolism, including healing. Polarized light sources have been shown to be effective in improving healing in many situations. Animals and Methods: Diabetes mellitus was induced with streptozotocin (60 mg/kg) in 45 male Wistar albino rats, and a third-degree burn (1.5 by 1.5 cm) was created on the dorsum of each animal under general anesthesia. The animals were randomly distributed into three groups: control, 20 J/cm(2), and 40 J/cm(2). Each group was then divided into three subgroups based on time of death (7, 14, 21 d). Phototherapy (20 or 40 J/cm(2) per session) was carried out immediately after the burning and repeated daily until the day before death. Following animal death, specimens were removed, embedded in paraffin, sectioned, and stained with hematoxylin and eosin (HE) or Sirius Red or immunomarked with CK AE1/AE3 antibody. Qualitative and semiquantitative analyses were performed under light microscopy. The results were statistically analyzed. Results: The animals treated with 20 J/cm(2) showed significant differences with regard to revascularization and re-epithelialization. Although the 40 J/cm(2) group showed stimulation of fibroblastic proliferation as an isolated feature, no other difference from the control was observed. Conclusion: Our results suggest that the use of polarized light at 20 J/cm(2) effectively improves the healing of third-degree burns on diabetic animals at both early and late stages of repair.

The chemical evolution of IC 10

Context. Dwarf irregular galaxies are relatively simple unevolved objects where it is easy to test models of galactic chemical evolution. Aims. We attempt to determine the star formation and gas accretion history of IC 10, a local dwarf irregular for which abundance, gas, and mass determinations are available. Methods. We apply detailed chemical evolution models to predict the evolution of several chemical elements (He, O, N, S) and compared our predictions with the observational data. We consider additional constraints such as the present-time gas fraction, the star formation rate (SFR), and the total estimated mass of IC 10. We assume a dark matter halo for this galaxy and study the development of a galactic wind. We consider different star formation regimes: bursting and continuous. We explore different wind situations: i) normal wind, where all the gas is lost at the same rate and ii) metal-enhanced wind, where metals produced by supernovae are preferentially lost. We study a case without wind. We vary the star formation efficiency (SFE), the wind efficiency, and the time scale of the gas infall, which are the most important parameters in our models. Results. We find that only models with metal-enhanced galactic winds can reproduce the properties of IC 10. The star formation must have proceeded in bursts rather than continuously and the bursts must have been less numerous than similar to 10 over the whole galactic lifetime. Finally, IC 10 must have formed by a slow process of gas accretion with a timescale of the order of 8 Gyr.

Merging Resource Availability with Isotope Mixing Models: The Role of Neutral Interaction Assumptions

Background: Bayesian mixing models have allowed for the inclusion of uncertainty and prior information in the analysis of trophic interactions using stable isotopes. Formulating prior distributions is relatively straightforward when incorporating dietary data. However, the use of data that are related, but not directly proportional, to diet (such as prey availability data) is often problematic because such information is not necessarily predictive of diet, and the information required to build a reliable prior distribution for all prey species is often unavailable. Omitting prey availability data impacts the estimation of a predator's diet and introduces the strong assumption of consumer ultrageneralism (where all prey are consumed in equal proportions), particularly when multiple prey have similar isotope values. Methodology: We develop a procedure to incorporate prey availability data into Bayesian mixing models conditional on the similarity of isotope values between two prey. If a pair of prey have similar isotope values (resulting in highly uncertain mixing model results), our model increases the weight of availability data in estimating the contribution of prey to a predator's diet. We test the utility of this method in an intertidal community against independently measured feeding rates. Conclusions: Our results indicate that our weighting procedure increases the accuracy by which consumer diets can be inferred in situations where multiple prey have similar isotope values. This suggests that the exchange of formalism for predictive power is merited, particularly when the relationship between prey availability and a predator's diet cannot be assumed for all species in a system.

Particle competition for complex network community detection

In many real situations, randomness is considered to be uncertainty or even confusion which impedes human beings from making a correct decision. Here we study the combined role of randomness and determinism in particle dynamics for complex network community detection. In the proposed model, particles walk in the network and compete with each other in such a way that each of them tries to possess as many nodes as possible. Moreover, we introduce a rule to adjust the level of randomness of particle walking in the network, and we have found that a portion of randomness can largely improve the community detection rate. Computer simulations show that the model has good community detection performance and at the same time presents low computational complexity. (C) 2008 American Institute of Physics.

Nonaxisymmetric magnetorotational instability in ideal and viscous plasmas

The excitation of magnetorotational instability (MRI) in rotating laboratory plasmas is investigated. In contrast to astrophysical plasmas, in which gravitation plays an important role, in laboratory plasmas it can be neglected and the plasma rotation is equilibrated by the pressure gradient. The analysis is restricted to the simple model of a magnetic confinement configuration with cylindrical symmetry, in which nonaxisymmetric perturbations are investigated using the local approximation. Starting from the simplest case of an ideal plasma, the corresponding dispersion relations are derived for more complicated models including the physical effects of parallel and perpendicular viscosities. The Friemann-Rotenberg approach used for ideal plasmas is generalized for the viscous model and an analytical expression for the instability boundary is obtained. It is shown that, in addition to the standard effect of radial derivative of the rotation frequency (the Velikhov effect), which can be destabilizing or stabilizing depending on the sign of this derivative in the ideal plasma, there is a destabilizing effect proportional to the fourth power of the rotation frequency, or, what is the same, to the square of the plasma pressure gradient, and to the square of the azimuthal mode number of the perturbations. It is shown that the instability boundary also depends on the product of the plasma pressure and density gradients, which has a destabilizing effect when it is negative. In the case of parallel viscosity, the MRI looks like an ideal instability independent of viscosity, while, in the case of strong perpendicular viscosity, it is a dissipative instability with the growth rate inversely proportional to the characteristic viscous decay rate. We point out, however, that the modes of the continuous range of the magnetohydrodynamics spectrum are not taken into account in this paper, and they can be more dangerous than those that are considered. (c) 2008 American Institute of Physics.

Resummation of infrared divergences in the free energy of spin-two fields

We derive a closed form expression for the sum of all the infrared divergent contributions to the free energy of a gas of gravitons. An important ingredient of our calculation is the use of a gauge fixing procedure such that the graviton propagator becomes both traceless and transverse. This has been shown to be possible, in a previous work, using a general gauge fixing procedure, in the context of the lowest order expansion of the Einstein-Hilbert action, describing noninteracting spin-two fields. In order to encompass the problems involving thermal loops, such as the resummation of the free energy, in the present work, we have extended this procedure to the situations when the interactions are taken into account.

Bayesian online algorithms for learning in discrete Hidden Markov Models

We propose and analyze two different Bayesian online algorithms for learning in discrete Hidden Markov Models and compare their performance with the already known Baldi-Chauvin Algorithm. Using the Kullback-Leibler divergence as a measure of generalization we draw learning curves in simplified situations for these algorithms and compare their performances.

Aspects of a noncommutative scalar/tensor duality

We study the noncommutative massless Kalb-Ramond gauge field coupled to a dynamical U(1) gauge field in the adjoint representation together with a compensating vector field. We derive the Seiberg-Witten map and obtain the corresponding mapped action to first order in theta. The (emergent) gravity structure found in other situations is not present here. The off-shell dual scalar theory is derived and it does not coincide with the Seiberg-Witten mapped scalar theory. Dispersion relations are also discussed. The p-form generalization of the Seiberg-Witten map to order theta is also derived.

Dynamical mutation of dark energy

We discuss the intriguing possibility that dark energy may change its equation of state in situations where large dark energy fluctuations are present. We show indications of this dynamical mutation in some generic models of dark energy.

Gravitational wave recoil in Robinson-Trautman spacetimes

We consider the gravitational recoil due to nonreflection-symmetric gravitational wave emission in the context of axisymmetric Robinson-Trautman spacetimes. We show that regular initial data evolve generically into a final configuration corresponding to a Schwarzschild black hole moving with constant speed. For the case of (reflection-)symmetric initial configurations, the mass of the remnant black hole and the total energy radiated away are completely determined by the initial data, allowing us to obtain analytical expressions for some recent numerical results that have appeared in the literature. Moreover, by using the Galerkin spectral method to analyze the nonlinear regime of the Robinson-Trautman equations, we show that the recoil velocity can be estimated with good accuracy from some asymmetry measures (namely the first odd moments) of the initial data. The extension for the nonaxisymmetric case and the implications of our results for realistic situations involving head-on collision of two black holes are also discussed.

Conservation laws, integrability, and transport in one-dimensional quantum systems

In integrable one-dimensional quantum systems an infinite set of local conserved quantities exists which can prevent a current from decaying completely. For cases like the spin current in the XXZ model at zero magnetic field or the charge current in the attractive Hubbard model at half filling, however, the current operator does not have overlap with any of the local conserved quantities. We show that in these situations transport at finite temperatures is dominated by a diffusive contribution with the Drude weight being either small or even zero. For the XXZ model we discuss in detail the relation between our results, the phenomenological theory of spin diffusion, and measurements of the spin-lattice relaxation rate in spin chain compounds. Furthermore, we study the Haldane-Shastry model where a conserved spin current exists.