Bayesian analysis and constraints on kinematic models from union SNIa

The kinematic expansion history of the universe is investigated by using the 307 supernovae type Ia from the Union Compilation set. Three simple model parameterizations for the deceleration parameter ( constant, linear and abrupt transition) and two different models that are explicitly parametrized by the cosmic jerk parameter ( constant and variable) are considered. Likelihood and Bayesian analyses are employed to find best fit parameters and compare models among themselves and with the flat Lambda CDM model. Analytical expressions and estimates for the deceleration and cosmic jerk parameters today (q(0) and j(0)) and for the transition redshift (z(t)) between a past phase of cosmic deceleration to a current phase of acceleration are given. All models characterize an accelerated expansion for the universe today and largely indicate that it was decelerating in the past, having a transition redshift around 0.5. The cosmic jerk is not strongly constrained by the present supernovae data. For the most realistic kinematic models the 1 sigma confidence limits imply the following ranges of values: q(0) is an element of [-0.96, -0.46], j(0) is an element of [-3.2,-0.3] and z(t) is an element of [0.36, 0.84], which are compatible with the Lambda CDM predictions, q(0) = -0.57 +/- 0.04, j(0) = -1 and z(t) = 0.71 +/- 0.08. We find that even very simple kinematic models are equally good to describe the data compared to the concordance Lambda CDM model, and that the current observations are not powerful enough to discriminate among all of them.

Programs to compute magnetization to density ratio and the magnetization inclination from 3-D gravity and magnetic anomalies

Vector field formulation based on the Poisson theorem allows an automatic determination of rock physical properties (magnetization to density ratio-MDR-and the magnetization inclination-MI) from combined processing of gravity and magnetic geophysical data. The basic assumptions (i.e., Poisson conditions) are: that gravity and magnetic fields share common sources, and that these sources have a uniform magnetization direction and MDR. In addition, the previously existing formulation was restricted to profile data, and assumed sufficiently elongated (2-D) sources. For sources that violate Poisson conditions or have a 3-D geometry, the apparent values of MDR and MI that are generated in this way have an unclear relationship to the actual properties in the subsurface. We present Fortran programs that estimate MDR and MI values for 3-D sources through processing of gridded gravity and magnetic data. Tests with simple geophysical models indicate that magnetization polarity can be successfully recovered by MDR-MI processing, even in cases where juxtaposed bodies cannot be clearly distinguished on the basis of anomaly data. These results may be useful in crustal studies, especially in mapping magnetization polarity from marine-based gravity and magnetic data. (c) 2007 Elsevier Ltd. All rights reserved.

Seed Cell Wall Storage Polysaccharides: Models to Understand Cell Wall Biosynthesis and Degradation

Cell wall storage polysaccharides (CWSPs) are found as the principal storage compounds in seeds of many taxonomically important groups of plants. These groups developed extremely efficient biochemical mechanisms to disassemble cell walls and use the products of hydrolysis for growth. To accumulate these storage polymers, developing seeds also contain relatively high activities of noncellulosic polysaccharide synthases and thus are interesting models to seek the discovery of genes and enzymes related to polysaccharide biosynthesis. CWSP systems offer opportunities to understand phenomena ranging from polysaccharide deposition during seed maturation to the control of source-sink relationship in developing seedlings. By studying polysaccharide biosynthesis and degradation and the consequences for cell and physiological behavior, we can use these models to develop future biotechnological applications.

Mitochondrial sequence data and Dipsacales phylogeny: Mixed models, partitioned Bayesian analyses, and model selection

Phylogenetic analyses of chloroplast DNA sequences, morphology, and combined data have provided consistent support for many of the major branches within the angiosperm, clade Dipsacales. Here we use sequences from three mitochondrial loci to test the existing broad scale phylogeny and in an attempt to resolve several relationships that have remained uncertain. Parsimony, maximum likelihood, and Bayesian analyses of a combined mitochondrial data set recover trees broadly consistent with previous studies, although resolution and support are lower than in the largest chloroplast analyses. Combining chloroplast and mitochondrial data results in a generally well-resolved and very strongly supported topology but the previously recognized problem areas remain. To investigate why these relationships have been difficult to resolve we conducted a series of experiments using different data partitions and heterogeneous substitution models. Usually more complex modeling schemes are favored regardless of the partitions recognized but model choice had little effect on topology or support values. In contrast there are consistent but weakly supported differences in the topologies recovered from coding and non-coding matrices. These conflicts directly correspond to relationships that were poorly resolved in analyses of the full combined chloroplast-mitochondrial data set. We suggest incongruent signal has contributed to our inability to confidently resolve these problem areas. (c) 2007 Elsevier Inc. All rights reserved.

Animal models for genetic neuromuscular diseases

The neuromuscular disorders are a heterogeneous group of genetic diseases, caused by mutations in genes coding sarcolemmal, sarcomeric, and citosolic muscle proteins. Deficiencies or loss of function of these proteins leads to variable degree of progressive loss of motor ability. Several animal models, manifesting phenotypes observed in neuromuscular diseases, have been identified in nature or generated in laboratory. These models generally present physiological alterations observed in human patients and can be used as important tools for genetic, clinic, and histopathological studies. The mdx mouse is the most widely used animal model for Duchenne muscular dystrophy (DMD). Although it is a good genetic and biochemical model, presenting total deficiency of the protein dystrophin in the muscle, this mouse is not useful for clinical trials because of its very mild phenotype. The canine golden retriever MD model represents a more clinically similar model of DMD due to its larger size and significant muscle weakness. Autosomal recessive limb-girdle MD forms models include the SJL/J mice, which develop a spontaneous myopathy resulting from a mutation in the Dysferlin gene, being a model for LGMD2B. For the human sarcoglycanopahties (SG), the BIO14.6 hamster is the spontaneous animal model for delta-SG deficiency, whereas some canine models with deficiency of SG proteins have also been identified. More recently, using the homologous recombination technique in embryonic stem cell, several mouse models have been developed with null mutations in each one of the four SG genes. All sarcoglycan-null animals display a progressive muscular dystrophy of variable severity and share the property of a significant secondary reduction in the expression of the other members of the sarcoglycan subcomplex and other components of the Dystrophin-glycoprotein complex. Mouse models for congenital MD include the dy/dy (dystrophia-muscularis) mouse and the allelic mutant dy(2J)/dy(2J) mouse, both presenting significant reduction of alpha 2-laminin in the muscle and a severe phenotype. The myodystrophy mouse (Large(myd)) harbors a mutation in the glycosyltransferase Large, which leads to altered glycosylation of alpha-DG, and also a severe phenotype. Other informative models for muscle proteins include the knockout mouse for myostatin, which demonstrated that this protein is a negative regulator of muscle growth. Additionally, the stress syndrome in pigs, caused by mutations in the porcine RYR1 gene, helped to localize the gene causing malignant hypertermia and Central Core myopathy in humans. The study of animal models for genetic diseases, in spite of the existence of differences in some phenotypes, can provide important clues to the understanding of the pathogenesis of these disorders and are also very valuable for testing strategies for therapeutic approaches.

The constrained compartmentalized knapsack problem: mathematical models and solution methods

The constrained compartmentalized knapsack problem can be seen as an extension of the constrained knapsack problem. However, the items are grouped into different classes so that the overall knapsack has to be divided into compartments, and each compartment is loaded with items from the same class. Moreover, building a compartment incurs a fixed cost and a fixed loss of the capacity in the original knapsack, and the compartments are lower and upper bounded. The objective is to maximize the total value of the items loaded in the overall knapsack minus the cost of the compartments. This problem has been formulated as an integer non-linear program, and in this paper, we reformulate the non-linear model as an integer linear master problem with a large number of variables. Some heuristics based on the solution of the restricted master problem are investigated. A new and more compact integer linear model is also presented, which can be solved by a branch-and-bound commercial solver that found most of the optimal solutions for the constrained compartmentalized knapsack problem. On the other hand, heuristics provide good solutions with low computational effort. (C) 2011 Elsevier BM. All rights reserved.

Destructive weighted Poisson cure rate models

In this paper, we develop a flexible cure rate survival model by assuming the number of competing causes of the event of interest to follow a compound weighted Poisson distribution. This model is more flexible in terms of dispersion than the promotion time cure model. Moreover, it gives an interesting and realistic interpretation of the biological mechanism of the occurrence of event of interest as it includes a destructive process of the initial risk factors in a competitive scenario. In other words, what is recorded is only from the undamaged portion of the original number of risk factors.

Continuum, Discrete, and Explicit Solvation Models for Describing the Low-Lying Absorption Spectrum of the Pterin Acid in Aqueous Environment

The absorption spectrum of the acid form of pterin in water was investigated theoretically. Different procedures using continuum, discrete, and explicit models were used to include the solvation effect on the absorption spectrum, characterized by two bands. The discrete and explicit models used Monte Carlo simulation to generate the liquid structure and time-dependent density functional theory (B3LYP/6-31G+(d)) to obtain the excitation energies. The discrete model failed to give the correct qualitative effect on the second absorption band. The continuum model, in turn, has given a correct qualitative picture and a semiquantitative description. The explicit use of 29 solvent molecules, forming a hydration shell of 6 angstrom, embedded in the electrostatic field of the remaining solvent molecules, gives absorption transitions at 3.67 and 4.59 eV in excellent agreement with the S(0)-S(1) and S(0)-S(2) absorption bands at of 3.66 and 4.59 eV, respectively, that characterize the experimental spectrum of pterin in water environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2371-2377, 2010

2-(4,5,6,7-tetrafluorobenzimidazol-2-yl)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazole-3-oxide-1-oxyl, a hydrogen-bonded organic quasi-1D ferromagnet

The title radical (F4BlmNN) is a stable nitronylnitroxide that forms hydrogen-bonded NH center dot center dot center dot ON chains in the solid state. The chains assemble the F4BlmNN molecules to form stacked contacts between the radical groups, in a geometry that is expected to exhibit ferromagnetic (FM) exchange based on spin polarization (SP) models. The experimental magnetic susceptibility of F4BlmNN confirms the expectation, showing 1-D Heisenberg chain FM exchange behavior over 1.8-300 K with an intrachain exchange constant Of J(chain)/k = +22 K. At lower temperatures, ac magnetic susceptibility and variable field heat capacity measurements show that F4BlmNN acts as a quasi-1-D ferromagnet. The dominant ferromagnetic exchange interaction is attributable to overlap between spin orbitals of molecules within the hydrogen-bonded chains, consistent with the SP model expectations. The chains appear to be antiferromagnetically exchange coupled, giving cusps in the ac susceptibility and zero field heat capacity at lower temperatures. The results indicate that the sample orders magnetically at about 0.7 K. The magnetic heat capacity ordering cusp shifts to lower temperatures as external magnetic field increases, consistent with forming a bulk antiferromagnetic phase below a Neel temperature of T-N(0) = 0.72 K, with a critical field of H-c approximate to 1800 Oe. The interchain exchange is estimated to be zJ/k congruent to (-)0.1 K.

Temperature Dependence of the Intergranular Critical Current Density in Uniaxially Pressed Bi(1.65)Pb(0.35)Sr(2)Ca(2)Cu(3)O(10+delta) Samples

We have studied the normal and superconducting transport properties of Bi(1.65)Pb(0.35)Sr(2)Ca(2)Cu(3)O(10+delta) (Bi-2223) ceramic samples. Four samples, from the same batch, were prepared by the solid-state reaction method and pressed uniaxially at different compacting pressures, ranging from 90 to 250 MPa before the last heat treatment. From the temperature dependence of the electrical resistivity, combined with current conduction models for cuprates, we were able to separate contributions arising from both the grain misalignment and microstructural defects. The behavior of the critical current density as a function of temperature at zero applied magnetic field, J (c) (T), was fitted to the relationship J (c) (T)ae(1-T/T (c) ) (n) , with na parts per thousand 2 in all samples. We have also investigated the behavior of the product J (c) rho (sr) , where rho (sr) is the specific resistance of the grain-boundary. The results were interpreted by considering the relation between these parameters and the grain-boundary angle, theta, with increasing the uniaxial compacting pressure. We have found that the above type of mechanical deformation improves the alignment of the grains. Consequently the samples exhibit an enhance in the intergranular properties, resulting in a decrease of the specific resistance of the grain-boundary and an increase in the critical current density.

Application of abrupt cut-off models in the analysis of the capacitance spectra of conjugated polymer devices

In this paper, a detailed study of the capacitance spectra obtained from Au/doped-polyaniline/Al structures in the frequency domain (0.05 Hz-10 MHz), and at different temperatures (150-340 K) is carried out. The capacitance spectra behavior in semiconductors can be appropriately described by using abrupt cut-off models, since they assume that the electronic gap states that can follow the ac modulation have response times varying rapidly with a certain abscissa, which is dependent on both temperature and frequency. Two models based on the abrupt cut-off concept, formerly developed to describe inorganic semiconductor devices, have been used to analyze the capacitance spectra of devices based on doped polyaniline (PANI), which is a well-known polymeric semiconductor with innumerous potential technological applications. The application of these models allowed the determination of significant parameters, such as Debye length (approximate to 20 nm), position of bulk Fermi level (approximate to 320 meV) and associated density of states (approximate to 2x10(18) eV(-1) cm(-3)), width of the space charge region (approximate to 70 nm), built-in potential (approximate to 780 meV), and the gap states` distribution.

Classical and hologram QSAR studies on a series of inhibitors of trypanosomatid Glyceraldehyde-3-Phosphate Dehydrogenase

Leishmaniasis and trypanosomiasis are major causes of morbidity and mortality in both tropical and subtropical regions of the world. The current available drugs are limited, ineffective, and require long treatment regimens. Due to the high dependence of trypanosomatids on glycolysis as a source of energy, some glycolytic enzymes have been identified as attractive targets for drug design. In the present work, classical Two-Dimensional Quantitative Structure -Activity Relationships (2D QSAR) and Hologram QSAR (HQSAR) studies were performed on a series of adenosine derivatives as inhibitors of Leishmania mexicana Glyceraldehyde-3-Phosphate Dehydrogenase (LmGAPDH). Significant correlation coefficients (classical QSAR, r(2)=0.83 and q(2) =0.81; HQSAR, r(2)=0.91 and q(2) =0.86) were obtained for the 56 training set compounds, indicating the potential of the models for untested compounds. The models were then externally validated using a test set of 14 structurally related compounds and the predicted values were in good agreement with the experimental results (classical QSAR, r(pred)(2) = 0.94; HQSAR, r(pred)(2) = 0.92).

Structural basis for selective inhibition of trypanosomatid glyceraldehyde-3-phosphate dehydrogenase: Molecular docking and 3D QSAR studies

The glycolytic enzyme glyceraldehyde-3 -phosphate dehydrogenase (GAPDH) is as an attractive target for the development of novel antitrypanosomatid agents. In the present work, comparative molecular field analysis and comparative molecular similarity index analysis were conducted on a large series of selective inhibitors of trypanosomatid GAPDH. Four statistically significant models were obtained (r(2) > 0.90 and q(2) > 0.70), indicating their predictive ability for untested compounds. The models were then used to predict the potency of an external test set, and the predicted values were in good agreement with the experimental results. Molecular modeling studies provided further insight into the structural basis for selective inhibition of trypanosomatid GAPDH.

Sampling, WLS, and Mixed Models

Mixed models may be defined with or without reference to sampling, and can be used to predict realized random effects, as when estimating the latent values of study subjects measured with response error. When the model is specified without reference to sampling, a simple mixed model includes two random variables, one stemming from an exchangeable distribution of latent values of study subjects and the other, from the study subjects` response error distributions. Positive probabilities are assigned to both potentially realizable responses and artificial responses that are not potentially realizable, resulting in artificial latent values. In contrast, finite population mixed models represent the two-stage process of sampling subjects and measuring their responses, where positive probabilities are only assigned to potentially realizable responses. A comparison of the estimators over the same potentially realizable responses indicates that the optimal linear mixed model estimator (the usual best linear unbiased predictor, BLUP) is often (but not always) more accurate than the comparable finite population mixed model estimator (the FPMM BLUP). We examine a simple example and provide the basis for a broader discussion of the role of conditioning, sampling, and model assumptions in developing inference.

Improved score tests in symmetric linear regression models

The class of symmetric linear regression models has the normal linear regression model as a special case and includes several models that assume that the errors follow a symmetric distribution with longer-than-normal tails. An important member of this class is the t linear regression model, which is commonly used as an alternative to the usual normal regression model when the data contain extreme or outlying observations. In this article, we develop second-order asymptotic theory for score tests in this class of models. We obtain Bartlett-corrected score statistics for testing hypotheses on the regression and the dispersion parameters. The corrected statistics have chi-squared distributions with errors of order O(n(-3/2)), n being the sample size. The corrections represent an improvement over the corresponding original Rao`s score statistics, which are chi-squared distributed up to errors of order O(n(-1)). Simulation results show that the corrected score tests perform much better than their uncorrected counterparts in samples of small or moderate size.