TESTING STATISTICAL HYPOTHESIS ON RANDOM TREES AND APPLICATIONS TO THE PROTEIN CLASSIFICATION PROBLEM

Efficient automatic protein classification is of central importance in genomic annotation. As an independent way to check the reliability of the classification, we propose a statistical approach to test if two sets of protein domain sequences coming from two families of the Pfam database are significantly different. We model protein sequences as realizations of Variable Length Markov Chains (VLMC) and we use the context trees as a signature of each protein family. Our approach is based on a Kolmogorov-Smirnov-type goodness-of-fit test proposed by Balding et at. [Limit theorems for sequences of random trees (2008), DOI: 10.1007/s11749-008-0092-z]. The test statistic is a supremum over the space of trees of a function of the two samples; its computation grows, in principle, exponentially fast with the maximal number of nodes of the potential trees. We show how to transform this problem into a max-flow over a related graph which can be solved using a Ford-Fulkerson algorithm in polynomial time on that number. We apply the test to 10 randomly chosen protein domain families from the seed of Pfam-A database (high quality, manually curated families). The test shows that the distributions of context trees coming from different families are significantly different. We emphasize that this is a novel mathematical approach to validate the automatic clustering of sequences in any context. We also study the performance of the test via simulations on Galton-Watson related processes.

FORAMINIFERAL RECORD OF CHANGES IN SUMMER MONSOON PRECIPITATION AT THE SOUTHEASTERN BRAZILIAN CONTINENTAL MARGIN SINCE THE LAST GLACIAL MAXIMUM

Changes in the oxygen isotopic composition of the planktonic foraminifer Globigerinoides ruber and in the foraminifera faunal composition in a core retrieved from the southeastern Brazilian continental margin were used to infer past changes in the hydrological balance and monsoon precipitation in the western South Atlantic since the Last Glacial Maximum (LGM). The results suggest a first-order orbital (precessional) control on the South American Monsoon precipitation. This agrees with previous studies based on continental proxies except for LGM estimates provided by pollen records. The causes for this disagreement are discussed.

Intermolecular vibrations and fast relaxations in supercooled ionic liquids

Short-time dynamics of ionic liquids has been investigated by low-frequency Raman spectroscopy (4 < omega < 100 cm(-1)) within the supercooled liquid range. Raman spectra are reported for ionic liquids with the same anion, bis(trifluoromethylsulfonyl)imide, and different cations: 1-butyl-3-methylimidazolium, 1-hexyl-3-methylimidazolium, 1-butyl-1-methylpiperidinium, trimethylbutylammonium, and tributylmethylammonium. It is shown that low-frequency Raman spectroscopy provides similar results as optical Kerr effect (OKE) spectroscopy, which has been used to study intermolecular vibrations in ionic liquids. The comparison of ionic liquids containing aromatic and non-aromatic cations identifies the characteristic feature in Raman spectra usually assigned to librational motion of the imidazolium ring. The strength of the fast relaxations (quasi-elastic scattering, QES) and the intermolecular vibrational contribution (boson peak) of ionic liquids with non-aromatic cations are significantly lower than imidazolium ionic liquids. A correlation length assigned to the boson peak vibrations was estimated from the frequency of the maximum of the boson peak and experimental data of sound velocity. The correlation length related to the boson peak (similar to 19 angstrom) does not change with the length of the alkyl chain in imidazolium cations, in contrast to the position of the first-sharp diffraction peak observed in neutron and X-ray scattering measurements of ionic liquids. The rate of change of the QES intensity in the supercooled liquid range is compared with data of excess entropy, free volume, and mean-squared displacement recently reported for ionic liquids. The temperature dependence of the QES intensity in ionic liquids illustrates relationships between short-time dynamics and long-time structural relaxation that have been proposed for glass-forming liquids. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3604533]

Last glacial maximum (LGM) vegetation changes in the Atlantic Forest, southeastern Brazil

This study was carried out in a continental Atlantic Forest located in the southern region of Sao Paulo State, southeastern Brazil. The aim of the study was to evaluate the vegetation dynamics in similar to 70 km forest ecosystem transect that occurred during the late Pleistocene and Holocene in this region, using the stable carbon isotopes (delta C-13) analysis on soil organic matter (SOM) and the C-14 dating of buried charcoal fragments and the humin fraction of SOM. The isotope data (delta C-13) of SOM in the deeper horizons, indicating the presence of more open vegetation than the present, with a probable mixture of C-3 and C-4 plants, suggesting the presence of a drier climate in the period of similar to 20 ka to similar to 16-14 ka BP. From similar to 16 to 14 ka BP to the present, a significant predominance Of C3 plants was observed, indicating an expansion of the forest, probably associated with the presence of a more humid climate than the previous period. The results indicated the presence of open vegetation during the late glacial, probably associated with a drier period, also observed in other regions of Brazil. The Atlantic Forest ecosystem seems to have developed at least since the early Holocene in Southeastern Brazil. (c) 2007 Elsevier Ltd and INQUA. All rights reserved.

A Markovian model market-Akerlof`s lemons and the asymmetry of information

In this work we study an agent based model to investigate the role of asymmetric information degrees for market evolution. This model is quite simple and may be treated analytically since the consumers evaluate the quality of a certain good taking into account only the quality of the last good purchased plus her perceptive capacity beta. As a consequence, the system evolves according to a stationary Markov chain. The value of a good offered by the firms increases along with quality according to an exponent alpha, which is a measure of the technology. It incorporates all the technological capacity of the production systems such as education, scientific development and techniques that change the productivity rates. The technological level plays an important role to explain how the asymmetry of information may affect the market evolution in this model. We observe that, for high technological levels, the market can detect adverse selection. The model allows us to compute the maximum asymmetric information degree before the market collapses. Below this critical point the market evolves during a limited period of time and then dies out completely. When beta is closer to 1 (symmetric information), the market becomes more profitable for high quality goods, although high and low quality markets coexist. The maximum asymmetric information level is a consequence of an ergodicity breakdown in the process of quality evaluation. (C) 2011 Elsevier B.V. All rights reserved.

The production of school numbers (1871-1931): contributions to a critical approach to statistical sources in History of Education

This article intends to contribute to the reflection on the Educational Statistics as being source for the researches on History of Education. The main concern was to reveal the way Educational Statistics related to the period from 1871 to 1931 were produced, in central government. Official reports - from the General Statistics Directory - and Statistics yearbooks released by that department were analyzed and, on this analysis, recommendations and definitions to perform the works were sought. By rending problematic to the documental issues on Educational Statistics and their usual interpretations, the intention was to reduce the ignorance about the origin of the school numbers, which are occasionally used in current researches without the convenient critical exam.

STATISTICAL EVALUATION OF MATURE LANDFILL LEACHATE TREATMENT BY HOMOGENEOUS CATALYTIC OZONATION

This study presents the results of a mature landfill leachate treated by a homogeneous catalytic ozonation process with ions Fe(2+) and Fe(3+) at acidic pH. Quality assessments were performed using Taguchi`s method (L(8) design). Strong synergism was observed statistically between molecular ozone and ferric ions, pointing to their catalytic effect on (center dot)OH generation. The achievement of better organic matter depollution rates requires an ozone flow of 5 L h(-1) (590 mg h(-1) O(3)) and a ferric ion concentration of 5 mg L(-1).

A study on the pretreatment of a sugarcane bagasse sample with dilute sulfuric acid

Experiments based on a 2(3) central composite full factorial design were carried out in 200-ml stainless-steel containers to study the pretreatment, with dilute sulfuric acid, of a sugarcane bagasse sample obtained from a local sugar-alcohol mill. The independent variables selected for study were temperature, varied from 112.5A degrees C to 157.5A degrees C, residence time, varied from 5.0 to 35.0 min, and sulfuric acid concentration, varied from 0.0% to 3.0% (w/v). Bagasse loading of 15% (w/w) was used in all experiments. Statistical analysis of the experimental results showed that all three independent variables significantly influenced the response variables, namely the bagasse solubilization, efficiency of xylose recovery in the hemicellulosic hydrolysate, efficiency of cellulose enzymatic saccharification, and percentages of cellulose, hemicellulose, and lignin in the pretreated solids. Temperature was the factor that influenced the response variables the most, followed by acid concentration and residence time, in that order. Although harsher pretreatment conditions promoted almost complete removal of the hemicellulosic fraction, the amount of xylose recovered in the hemicellulosic hydrolysate did not exceed 61.8% of the maximum theoretical value. Cellulose enzymatic saccharification was favored by more efficient removal of hemicellulose during the pretreatment. However, detoxification of the hemicellulosic hydrolysate was necessary for better bioconversion of the sugars to ethanol.

PVA-Hydrogel Entrapped Candida Guilliermondii for Xylitol Production from Sugarcane Hemicellulose Hydrolysate

Viable cells of Candida guilliermondii were immobilized by inclusion into polyvinyl alcohol (PVA) hydrogel using the freezing-thawing method. Entrapment experiments were planned according to a 2(3) full factorial design, using the PVA concentration (80, 100, and 120 g L(-1)), the freezing temperature (-10, -15, and -20 degrees C), and the number of freezing-thawing cycles (one, three, and five) as the independent variables, integrated with three additional tests to estimate the errors. The effectiveness of the immobilization procedure was checked in Erlenmeyer flasks as the pellet capability to catalyze the xylose-to-xylitol bioconversion of a medium based on sugarcane bagasse hemicellulosic hydrolysate. To this purpose, the yield of xylitol on consumed xylose, xylitol volumetric productivity, and cell retention yield were selected as the response variables. Cell pellets were then used to perform the same bioconversion in a stirred tank reactor operated at 400 rpm, 30 degrees C, and 1.04 vvm air flowrate. At the end of fermentation, a maximum xylitol concentration of 28.7 g L(-1), a xylitol yield on consumed xylose of 0.49 g g(-1) and a xylitol volumetric productivity of 0.24 g L(-1) h(-1) were obtained.

The influence of microstructure on the maximum load and fracture energy of refractory castables

Refractory castables are composed of fractions of fine to fairly coarse particles. The fine fraction is constituted primarily of raw materials and calcium aluminate cement, which becomes hydrated, forming chemical bonds that stiffen the concrete during the curing process. The present study focused on an evaluation of several characteristics of two refractory castables with similar chemical compositions but containing aggregates of different sizes. The features evaluated were the maximum load, the fracture energy, and the ""relative crack-propagation work"" of the two castables heat-treated at 110, 650, 1100 and 1550 degrees C. The results enabled us to draw the following conclusions: the heat treatment temperature exerts a significant influence on the matrix/aggregate interaction, different microstructures form in the castables with temperature, and a relationship was noted between the maximum load and the fracture energy. (C) 2009 Elsevier Ltd and Techna Group S.r.l. All rights reserved.

Effects of elastic indenter deformation on spherical instrumented indentation tests: the reduced elastic modulus

Although the Hertz theory is not applicable in the analysis of the indentation of elastic-plastic materials, it is common practice to incorporate the concept of indenter/specimen combined modulus to consider indenter deformation. The appropriateness was assessed of the use of reduced modulus to incorporate the effect of indenter deformation in the analysis of the indentation with spherical indenters. The analysis based on finite element simulations considered four values of the ratio of the indented material elastic modulus to that of the diamond indenter, E/E(i) (0, 0.04, 0.19, 0.39), four values of the ratio of the elastic reduced modulus to the initial yield strength, E(r)/Y (0, 10, 20, 100), and two values of the ratio of the indenter radius to maximum total displacement, R/delta(max) (3, 10). Indenter deformation effects are better accounted for by the reduced modulus if the indented material behaves entirely elastically. In this case, identical load-displacement (P - delta) curves are obtained with rigid and elastic spherical indenters for the same elastic reduced modulus. Changes in the ratio E/E(i), from 0 to 0.39, resulted in variations lower than 5% for the load dimensionless functions, lower than 3% in the contact area, A(c), and lower than 5% in the ratio H/E(r). However, deformations of the elastic indenter made the actual radius of contact change, even in the indentation of elastic materials. Even though the load dimensionless functions showed only a little increase with the ratio E/E(i), the hardening coefficient and the yield strength could be slightly overestimated when algorithms based on rigid indenters are used. For the unloading curves, the ratio delta(e)/delta(max), where delta(e) is the point corresponding to zero load of a straight line with slope S from the point (P(max), delta(max)), varied less than 5% with the ratio E/E(i). Similarly, the relationship between reduced modulus and the unloading indentation curve, expressed by Sneddon`s equation, did not reveal the necessity of correction with the ratio E/E(i). The most affected parameter in the indentation curve, as a consequence of the indentation deformation, was the ratio between the residual indentation depth after complete unloading and the maximum indenter displacement, delta(r)/delta(max) (up to 26%), but this variation did not significantly decrease the capability to estimate hardness and elastic modulus based on the ratio of the residual indentation depth to maximum indentation depth, h(r)/h(max). In general, the results confirm the convenience of the use of the reduced modulus in the spherical instrumented indentation tests.

Combining Mager and Steinmetz: The Effect of Grain Size and Maximum Induction on Hysteresis Energy Loss

Twelve samples with different grain sizes were prepared by normal grain growth and by primary recrystallization, and the hysteresis dissipated energy was measured by a quasi-static method. Results showed a linear relation between hysteresis energy loss and the inverse of grain size, which is here called Mager`s law, for maximum inductions from 0.6 to 1.5 T, and a Steinmetz power law relation between hysteresis loss and maximum induction for all samples. The combined effect is better described by a Mager`s law where the coefficients follow Steinmetz law.

Adaptive Carrier Tracking for Mars to Earth Communications During Entry, Descent, and Landing

We propose a robust and low complexity scheme to estimate and track carrier frequency from signals traveling under low signal-to-noise ratio (SNR) conditions in highly nonstationary channels. These scenarios arise in planetary exploration missions subject to high dynamics, such as the Mars exploration rover missions. The method comprises a bank of adaptive linear predictors (ALP) supervised by a convex combiner that dynamically aggregates the individual predictors. The adaptive combination is able to outperform the best individual estimator in the set, which leads to a universal scheme for frequency estimation and tracking. A simple technique for bias compensation considerably improves the ALP performance. It is also shown that retrieval of frequency content by a fast Fourier transform (FFT)-search method, instead of only inspecting the angle of a particular root of the error predictor filter, enhances performance, particularly at very low SNR levels. Simple techniques that enforce frequency continuity improve further the overall performance. In summary we illustrate by extensive simulations that adaptive linear prediction methods render a robust and competitive frequency tracking technique.

A method to estimate EMG crosstalk between two muscles based on the silent period following an H-reflex

The crosstalk phenomenon consists in recording the volume-conducted electromyographic activity of muscles other than that under study. This interference may impair the correct interpretation of the results in a variety of experiments. A new protocol is presented here for crosstalk assessment between two muscles based on changes in their electrical activity following a reflex discharge in one of the muscles in response to nerve stimulation. A reflex compound muscle action potential (H-reflex) was used to induce a silent period in the muscle that causes the crosstalk, called here the remote muscle. The rationale is that if the activity recorded in the target muscle is influenced by a distant source (the remote muscle) a silent period observed in the electromyogram (EMG) of the remote muscle would coincide with a decrease in the EMG activity of the target muscle. The new crosstalk index is evaluated based on the root mean square (RMS) values of the EMGs obtained in two distinct periods (background EMG and silent period) of both the remote and the target muscles. In the present work the application focused on the estimation of the degree of crosstalk from the soleus muscle to the tibialis anterior muscle during quiet stance. However, the technique may be extended to other pairs of muscles provided a silent period may be evoked in one of them. (C) 2009 IPEM. Published by Elsevier Ltd. All rights reserved.

Determining the effect of grain size and maximum induction upon coercive field of electrical steels

Although theoretical models have already been proposed, experimental data is still lacking to quantify the influence of grain size upon coercivity of electrical steels. Some authors consider a linear inverse proportionality, while others suggest a square root inverse proportionality. Results also differ with regard to the slope of the reciprocal of grain size-coercive field relation for a given material. This paper discusses two aspects of the problem: the maximum induction used for determining coercive force and the possible effect of lurking variables such as the grain size distribution breadth and crystallographic texture. Electrical steel sheets containing 0.7% Si, 0.3% Al and 24 ppm C were cold-rolled and annealed in order to produce different grain sizes (ranging from 20 to 150 mu m). Coercive field was measured along the rolling direction and found to depend linearly on reciprocal of grain size with a slope of approximately 0.9 (A/m)mm at 1.0 T induction. A general relation for coercive field as a function of grain size and maximum induction was established, yielding an average absolute error below 4%. Through measurement of B(50) and image analysis of micrographs, the effects of crystallographic texture and grain size distribution breadth were qualitatively discussed. (C) 2011 Elsevier B.V. All rights reserved.