A new linear initialization in SOM for biomolecular data

In the past decade, the amount of data in biological field has become larger and larger; Bio-techniques for analysis of biological data have been developed and new tools have been introduced. Several computational methods are based on unsupervised neural network algorithms that are widely used for multiple purposes including clustering and visualization, i.e. the Self Organizing Maps (SOM). Unfortunately, even though this method is unsupervised, the performances in terms of quality of result and learning speed are strongly dependent from the neuron weights initialization. In this paper we present a new initialization technique based on a totally connected undirected graph, that report relations among some intersting features of data input. Result of experimental tests, where the proposed algorithm is compared to the original initialization techniques, shows that our technique assures faster learning and better performance in terms of quantization error.

Sequence-selective assembly of tweezer-molecules on linear templates enables frameshift-reading of sequence information

Monomer-sequence information in synthetic copolyimides can be recognised by tweezer-type molecules binding to adjacent triplet-sequences on the polymer chains. In the present paper different tweezer-molecules are found to have different sequence-selectivities, as demonstrated in solution by 1H NMR spectroscopy and in the solid state by single crystal X-ray analyses of tweezer-complexes with linear and macrocyclic oligo-imides. This work provides clear-cut confirmation of polyimide chain-folding and adjacent-tweezer-binding. It also reveals a new and entirely unexpected mechanism for sequence-recognition which, by analogy with a related process in biomolecular information processing, may be termed "frameshift-reading". The ability of one particular tweezer-molecule to detect, with exceptionally high sensitivity, long-range sequence-information in chain-folding aromatic copolyimides, is readily explained by this novel process.

An old problem revisited - and solved: Upton and Dalton's linear production response: A note

This note presents a robust method for estimating response surfaces that consist of linear response regimes and a linear plateau. The linear response-and-plateau model has fascinated production scientists since von Liebig (1855) and, as Upton and Dalton indicated, some years ago in this Journal, the response-and-plateau model seems to fit the data in many empirical studies. The estimation algorithm evolves from Bayesian implementation of a switching-regression (finite mixtures) model and demonstrates routine application of Gibbs sampling and data augmentation-techniques that are now in widespread application in other disciplines.

Alternatives to linear analysis of energy balance data from lactating dairy cows

The current energy requirements system used in the United Kingdom for lactating dairy cows utilizes key parameters such as metabolizable energy intake (MEI) at maintenance (MEm), the efficiency of utilization of MEI for 1) maintenance, 2) milk production (k(l)), 3) growth (k(g)), and the efficiency of utilization of body stores for milk production (k(t)). Traditionally, these have been determined using linear regression methods to analyze energy balance data from calorimetry experiments. Many studies have highlighted a number of concerns over current energy feeding systems particularly in relation to these key parameters, and the linear models used for analyzing. Therefore, a database containing 652 dairy cow observations was assembled from calorimetry studies in the United Kingdom. Five functions for analyzing energy balance data were considered: straight line, two diminishing returns functions, (the Mitscherlich and the rectangular hyperbola), and two sigmoidal functions (the logistic and the Gompertz). Meta-analysis of the data was conducted to estimate k(g) and k(t). Values of 0.83 to 0.86 and 0.66 to 0.69 were obtained for k(g) and k(t) using all the functions (with standard errors of 0.028 and 0.027), respectively, which were considerably different from previous reports of 0.60 to 0.75 for k(g) and 0.82 to 0.84 for k(t). Using the estimated values of k(g) and k(t), the data were corrected to allow for body tissue changes. Based on the definition of k(l) as the derivative of the ratio of milk energy derived from MEI to MEI directed towards milk production, MEm and k(l) were determined. Meta-analysis of the pooled data showed that the average k(l) ranged from 0.50 to 0.58 and MEm ranged between 0.34 and 0.64 MJ/kg of BW0.75 per day. Although the constrained Mitscherlich fitted the data as good as the straight line, more observations at high energy intakes (above 2.4 MJ/kg of BW0.75 per day) are required to determine conclusively whether milk energy is related to MEI linearly or not.

Effects of volatile fatty acid supply on their absorption and on water kinetics in the rumen of sheep sustained by intragastric infusions

Three sheep fitted with a ruminal cannula and an abomasal catheter were used to study water kinetics and absorption of VFA infused continuously into the rumen. The effects of changing VFA concentrations in the rumen by shifting VFA infusion rates were investigated in an experiment with a 3 x 3 Latin square design. On experimental days, the animals received the basal infusion rate of VFA (271 mmol/h) during the first 2 h. Each animal then received VFA at a different rate (135, 394, or 511 mmol/h) for the next 7.5 h. Using soluble markers (polyethylene glycol and Cr-EDTA), ruminal volume, liquid outflow, apparent water absorption, and VFA absorption rates were estimated. There were no significant effects of VFA infusion rate on ruminal volume and water kinetics. As the VFA infusion rate was increased, VFA concentration and osmolality in the rumen were increased and pH was decreased. There was a biphasic response of liquid outflow to changes in the total VFA concentration in the rumen, as both variables increased together up to a total VFA concentration of 80.1 mM, whereas, beyond that concentration, liquid outflow remained stable at an average rate of 407 mL/h. There were significant linear (P = 0.003) and quadratic (P = 0.001) effects of VFA infusion rate on the VFA absorption rate, confirming that VFA absorption in the rumen is mainly a concentration-dependent process. The proportion of total VFA supplied that was absorbed in the rumen was 0.845 (0.822, 0.877, and 0.910 for acetate, propionate, and butyrate, respectively). The molar proportions of acetate, propionate, and butyrate absorbed were affected by the level of VFA infusion in the rumen, indicating that this level affected to a different extent the absorption of the different acids.

Inclusion of nonlinear demand-supply relationships within large-scale partial equilibrium linear programming models

This paper presents a new method for the inclusion of nonlinear demand and supply relationships within a linear programming model. An existing method for this purpose is described first and its shortcomings are pointed out before showing how the new approach overcomes those difficulties and how it provides a more accurate and 'smooth' (rather than a kinked) approximation of the nonlinear functions as well as dealing with equilibrium under perfect competition instead of handling just the monopolistic situation. The workings of the proposed method are illustrated by extending a previously available sectoral model for the UK agriculture.

The effect of wilting and soaking Eupatorium adenophorum on its digestibility in vitro and voluntary intake by goats

An experiment was conducted to determine what effect simple treatments might have on the voluntary intake by goats in Nepal of Eupatorium adenophorum, an invasive weed that is usually only consumed by goats to a very limited extent. Samples of E. adenophorum were collected and either untreated, soaked for 2 h or wilted for 2 h before being oven dried (60 degrees C) and ground. Soaking and wilting had little effect on the chemical composition of E. adenophorum, but did increase (P=0.036) its in vitro organic matter degradability, by approximately 8%. The short-term intake rate (STIR) of treated and untreated E. adenophorum was then estimated with eight goats. Soaking time (from 2 to 24 h) was not related to STIR (r = -0.111, P=0.198), but the time E. adenophorum was left to wilt (from 2 to 48h), was positively related to STIR (r=0.521, P<0.001), with values of STIR (g dry matter/min kg goat liveweight(0.75)) being 0.405, 0.649,1.058, S.E.M. 0.088 for E. adenophorum, that had been wilted for 0, 24 and 48 h respectively (P<0.001). Liveweight change of goats and voluntary intake of E. adenophorum by goats was then estimated with 24 goats. E. adenophorum was fed either unwilted, or wilted for 24 or 48 h. It was fed as the sole forage or as a 3:1 mixture (dry matter basis) with Ficus cunia. There was a linear (P<0.001) and quadratic (P<0.01) increase in the intake of total forage and E. adenophorum with wilting time of E. adenophorum. Offering Ficus cunia increased total forage intake, but decreased E. adenophorum intake (P<0.05). After four weeks, there was virtually no change in goat liveweight and no significant difference between treatments. The results suggest that wilting E adenophorum for 24 h could increase its intake by goats, and thereby increase its usefulness, as a potential source of forage in the dry season of Nepal, when forage scarcity is a common constraint to livestock production. (C) 2007 Elsevier B.V. All rights reserved.

Analysis of feeding strategies for small-scale dairy systems in central Mexico using linear programming

A limitation of small-scale dairy systems in central Mexico is that traditional feeding strategies are less effective when nutrient availability varies through the year. In the present work, a linear programming (LP) model that maximizes income over feed cost was developed, and used to evaluate two strategies: the traditional one used by the small-scale dairy producers in Michoacan State, based on fresh lucerne, maize grain and maize straw; and an alternative strategy proposed by the LIP model, based on ryegrass hay, maize silage and maize grain. Biological and economic efficiency for both strategies were evaluated. Results obtained with the traditional strategy agree with previously published work. The alternative strategy did not improve upon the performance of the traditional strategy because of low metabolizable protein content of the maize silage considered by the model. However, the Study recommends improvement of forage quality to increase the efficiency of small-scale dairy systems, rather than looking for concentrate supplementation.

Forage production and nutrient availability in small-scale dairy systems in central Mexico using linear programming and partial budgeting

Small-scale dairy systems play an important role in the Mexican dairy sector and farm planning activities related to resource allocation have a significant impact on the profitability of such enterprises. Linear programming is a technique widely used for planning and ration formulation, and partial budgeting is a technique for assessing the impact of changes on the profitability of an enterprise. This study used both methods to optimise land use for forage production and nutrient availability, and to evaluate the economic impact of such changes in small-scale Mexican dairy systems. The model showed satisfactory performance when optimal solutions were compared with the traditional strategy. The strategy using fresh ryegrass, maize silage and oat hay, and the strategy using a combination of alfalfa hay, maize silage, fresh ryegrass and oat hay appeared attractive options for providing a better nutrient supply and maintaining a higher stocking rate throughout the year than the traditional strategy.

A Bayesian approach for dose-escalation in a phase I clinical trial incorporating pharmacodynamic endpoints

Bayesian decision procedures have already been proposed for and implemented in Phase I dose-escalation studies in healthy volunteers. The procedures have been based on pharmacokinetic responses reflecting the concentration of the drug in blood plasma and are conducted to learn about the dose-response relationship while avoiding excessive concentrations. However, in many dose-escalation studies, pharmacodynamic endpoints such as heart rate or blood pressure are observed, and it is these that should be used to control dose-escalation. These endpoints introduce additional complexity into the modeling of the problem relative to pharmacokinetic responses. Firstly, there are responses available following placebo administrations. Secondly, the pharmacodynamic responses are related directly to measurable plasma concentrations, which in turn are related to dose. Motivated by experience of data from a real study conducted in a conventional manner, this paper presents and evaluates a Bayesian procedure devised for the simultaneous monitoring of pharmacodynamic and pharmacokinetic responses. Account is also taken of the incidence of adverse events. Following logarithmic transformations, a linear model is used to relate dose to the pharmacokinetic endpoint and a quadratic model to relate the latter to the pharmacodynamic endpoint. A logistic model is used to relate the pharmacokinetic endpoint to the risk of an adverse event.

Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.