Exact Search Directions for Optimization of Linear and Nonlinear Models Based on Generalized Orthonormal Functions

A novel technique for selecting the poles of orthonormal basis functions (OBF) in Volterra models of any order is presented. It is well-known that the usual large number of parameters required to describe the Volterra kernels can be significantly reduced by representing each kernel using an appropriate basis of orthonormal functions. Such a representation results in the so-called OBF Volterra model, which has a Wiener structure consisting of a linear dynamic generated by the orthonormal basis followed by a nonlinear static mapping given by the Volterra polynomial series. Aiming at optimizing the poles that fully parameterize the orthonormal bases, the exact gradients of the outputs of the orthonormal filters with respect to their poles are computed analytically by using a back-propagation-through-time technique. The expressions relative to the Kautz basis and to generalized orthonormal bases of functions (GOBF) are addressed; the ones related to the Laguerre basis follow straightforwardly as a particular case. The main innovation here is that the dynamic nature of the OBF filters is fully considered in the gradient computations. These gradients provide exact search directions for optimizing the poles of a given orthonormal basis. Such search directions can, in turn, be used as part of an optimization procedure to locate the minimum of a cost-function that takes into account the error of estimation of the system output. The Levenberg-Marquardt algorithm is adopted here as the optimization procedure. Unlike previous related work, the proposed approach relies solely on input-output data measured from the system to be modeled, i.e., no information about the Volterra kernels is required. Examples are presented to illustrate the application of this approach to the modeling of dynamic systems, including a real magnetic levitation system with nonlinear oscillatory behavior.

Size and power properties of some tests in the Birnbaum-Saunders regression model

The Birnbaum-Saunders distribution has been used quite effectively to model times to failure for materials subject to fatigue and for modeling lifetime data. In this paper we obtain asymptotic expansions, up to order n(-1/2) and under a sequence of Pitman alternatives, for the non-null distribution functions of the likelihood ratio, Wald, score and gradient test statistics in the Birnbaum-Saunders regression model. The asymptotic distributions of all four statistics are obtained for testing a subset of regression parameters and for testing the shape parameter. Monte Carlo simulation is presented in order to compare the finite-sample performance of these tests. We also present two empirical applications. (C) 2010 Elsevier B.V. All rights reserved.

Local power of some tests in exponential family nonlinear models

In this paper we obtain asymptotic expansions up to order n(-1/2) for the nonnull distribution functions of the likelihood ratio, Wald, score and gradient test statistics in exponential family nonlinear models (Cordeiro and Paula, 1989), under a sequence of Pitman alternatives. The asymptotic distributions of all four statistics are obtained for testing a subset of regression parameters and for testing the dispersion parameter, thus generalising the results given in Cordeiro et al. (1994) and Ferrari et al. (1997). We also present Monte Carlo simulations in order to compare the finite-sample performance of these tests. (C) 2010 Elsevier B.V. All rights reserved.

The Staircase Structure of the Southern Brazilian Continental Shelf

We show some evidences that the Southeastern Brazilian Continental Shelf (SBCS) has a devil's staircase structure, with a sequence of scarps and terraces with widths that obey fractal formation rules. Since the formation of these features is linked with the sea-level variations, we say that the sea level changes in an organized pulsating way. Although the proposed approach was applied in a particular region of the Earth, it is suitable to be applied in an integrated way to other shelves around the world, since the analyses favor the revelation of the global sea-level variations. Copyright (C) 2009 M. S. Baptista and L. A. Conti.

A simple procedure for reliability assessment of thin composite plates against buckling

In some circumstances ice floes may be modeled as beams. In general this modeling supposes constant thickness, which contradicts field observations. Action of currents, wind and the sequence of contacts, causes thickness to vary. Here this effect is taken into consideration on the modeling of the behavior of ice hitting inclined walls of offshore platforms. For this purpose, the boundary value problem is first equated. The set of equations so obtained is then transformed into a system of equations, that is then solved numerically. For this sake an implicit solution is developed, using a shooting method, with the accompanying Jacobian. In-plane coupling and the dependency of the boundary terms on deformation, make the problem non-linear and the development particular. Deformation and internal resultants are then computed for harmonic forms of beam profile. Forms of giving some additional generality to the problem are discussed.

Identification of a Highly Antigenic Linear B Cell Epitope within Plasmodium vivax Apical Membrane Antigen 1 (AMA-1)

Apical membrane antigen 1 (AMA-1) is considered to be a major candidate antigen for a malaria vaccine. Previous immunoepidemiological studies of naturally acquired immunity to Plasmodium vivax AMA-1 (PvAMA-1) have shown a higher prevalence of specific antibodies to domain II (DII) of AMA-1. In the present study, we confirmed that specific antibody responses from naturally infected individuals were highly reactive to both full-length AMA-1 and DII. Also, we demonstrated a strong association between AMA-1 and DII IgG and IgG subclass responses. We analyzed the primary sequence of PvAMA-1 for B cell linear epitopes co-occurring with intrinsically unstructured/ disordered regions (IURs). The B cell epitope comprising the amino acid sequence 290-307 of PvAMA-1 (SASDQPTQYEEEMTDYQK), with the highest prediction scores, was identified in domain II and further selected for chemical synthesis and immunological testing. The antigenicity of the synthetic peptide was identified by serological analysis using sera from P. vivax-infected individuals who were knowingly reactive to the PvAMA-1 ectodomain only, domain II only, or reactive to both antigens. Although the synthetic peptide was recognized by all serum samples specific to domain II, serum with reactivity only to the full-length protein presented 58.3% positivity. Moreover, IgG reactivity against PvAMA-1 and domain II after depletion of specific synthetic peptide antibodies was reduced by 18% and 33% (P = 0.0001 for both), respectively. These results suggest that the linear epitope SASDQPTQYEEEMTDYQK is highly antigenic during natural human infections and is an important antigenic region of the domain II of PvAMA-1, suggesting its possible future use in pre-clinical studies.

Unexpected Diversity of Cellular Immune Responses against Nef and Vif in HIV-1-Infected Patients Who Spontaneously Control Viral Replication

Background: HIV-1-infected individuals who spontaneously control viral replication represent an example of successful containment of the AIDS virus. Understanding the anti-viral immune responses in these individuals may help in vaccine design. However, immune responses against HIV-1 are normally analyzed using HIV-1 consensus B 15-mers that overlap by 11 amino acids. Unfortunately, this method may underestimate the real breadth of the cellular immune responses against the autologous sequence of the infecting virus. Methodology and Principal Findings: Here we compared cellular immune responses against nef and vif-encoded consensus B 15-mer peptides to responses against HLA class I-predicted minimal optimal epitopes from consensus B and autologous sequences in six patients who have controlled HIV-1 replication. Interestingly, our analysis revealed that three of our patients had broader cellular immune responses against HLA class I-predicted minimal optimal epitopes from either autologous viruses or from the HIV-1 consensus B sequence, when compared to responses against the 15-mer HIV-1 type B consensus peptides. Conclusion and Significance: This suggests that the cellular immune responses against HIV-1 in controller patients may be broader than we had previously anticipated.

Report of a chimeric origin of transposable elements in a bovine-coding gene

Despite the wide distribution of transposable elements (TEs) in mammalian genomes, part of their evolutionary significance remains to be discovered. Today there is a substantial amount of evidence showing that TEs are involved in the generation of new exons in different species. In the present study, we searched 22,805 genes and reported the occurrence of TE-cassettes in coding sequences of 542 cow genes using the RepeatMasker program. Despite the significant number (542) of genes with TE insertions in exons only 14 (2.6%) of them were translated into protein, which we characterized as chimeric genes. From these chimeric genes, only the FAST kinase domains 3 (FASTKD3) gene, present on chromosome BTA 20, is a functional gene and showed evidence of the exaptation event. The genome sequence analysis showed that the last exon coding sequence of bovine FASTKD3 is similar to 85% similar to the ART2A retrotransposon sequence. In addition, comparison among FASTKD3 proteins shows that the last exon is very divergent from those of Homo sapiens, Pan troglodytes and Canis familiares. We suggest that the gene structure of bovine FASTKD3 gene could have originated by several ectopic recombinations between TE copies. Additionally, the absence of TE sequences in all other species analyzed suggests that the TE insertion is clade-specific, mainly in the ruminant lineage.

EVOLUTION OF FLORAL MORPHOLOGY AND POLLINATION SYSTEM IN BIGNONIEAE (BIGNONIACEAE)

The radiation of angiosperms is associated with shifts among pollination modes that are thought to have driven the diversification of floral forms. However, the exact sequence of evolutionary events that led to such great diversity in floral traits is unknown for most plant groups. Here, we characterize the patterns of evolution of individual floral traits and overall floral morphologies in the tribe Bignonieae (Bignoniaceae). We identified 12 discrete traits that are associated with seven floral types previously described for the group and used a penalized likelihood tree of the tribe to reconstruct the ancestral states of those traits at all nodes of the phylogeny of Bignonieae. In addition, evolutionary correlations among traits were conducted using a maximum likelihood approach to test whether the evolution of individual floral traits followed the correlated patterns of evolution expected under the ""pollination syndrome"" concept. The ancestral Bignonieae flower presented an Anemopaegma-type morphology, which was followed by several parallel shifts in floral morphologies. Those shifts occurred through intermediate stages resulting in mixed floral morphologies as well as directly from the Anemopaegma-type morphology to other floral types. Positive and negative evolutionary correlations among traits fit patterns expected under the pollination syndrome perspective, suggesting that interactions between Bignonieae flowers and pollinators likely played important roles in the diversification of the group as a whole.

LIFE CYCLE AND POPULATION STRUCTURE OF AEGLA PAULENSIS (DECAPODA: ANOMURA: AEGLIDAE)

We describe growth, longevity, sex ratio, reproductive period, and recruitment of Aegla paulensis from Jaragua Stale Park, Sao Paulo, Brazil (23 degrees 27'27.9 '' S; 46 degrees 45'32.3 '' W). The population was sampled monthly (September 2007 through August 2009) with the aid of traps. Over five thousand individuals were captured, sexed, measured (carapace length = CL) and inspected for reproductive traits (females only), and then released back to the sampling site. The pattern of the reproductive cycle was strongly seasonal (austral mid autumn through late winter), with a single recruitment pulse per year. The obtained von Bertalanffy growth equations were CL = 21.25[1-e(-0.041(t + 1.250))] and CL = 16.52[1-e(-0.049(t + 1.823))] for males and females, respectively. Males (mean CL +/- SD = 11.86 +/- 2.79 mm) attain larger sizes than females (mean CL +/- SD = 10.84 +/- 2.36 mm). Aegla paulensis reproduces twice during an estimated life span of 40.2 months for females and 33.9 months for males. Temporal variation of sex ratio showed a distinctive pattern characterized by a sequence of three distinct periods that repeated from one year to another, and which suggested that a behavioral component influence the proportion of sex in adult specimens sampled with traps during reproductive and non-reproductive periods.

Musical genres: beating to the rhythms of different drums

Online music databases have increased significantly as a consequence of the rapid growth of the Internet and digital audio, requiring the development of faster and more efficient tools for music content analysis. Musical genres are widely used to organize music collections. In this paper, the problem of automatic single and multi-label music genre classification is addressed by exploring rhythm-based features obtained from a respective complex network representation. A Markov model is built in order to analyse the temporal sequence of rhythmic notation events. Feature analysis is performed by using two multi-variate statistical approaches: principal components analysis (unsupervised) and linear discriminant analysis (supervised). Similarly, two classifiers are applied in order to identify the category of rhythms: parametric Bayesian classifier under the Gaussian hypothesis (supervised) and agglomerative hierarchical clustering (unsupervised). Qualitative results obtained by using the kappa coefficient and the obtained clusters corroborated the effectiveness of the proposed method.

The mitochondrial genome of the phytopathogenic basidiomycete Moniliophthora perniciosa is 109 kb in size and contains a stable integrated plasmid

We present here the sequence of the mitochondrial genome of the basidiomycete phytopathogenic hemibiotrophic fungus Moniliophthora perniciosa, causal agent of the Witches` Broom Disease in Theobroma cacao. The DNA is a circular molecule of 109103 base pairs, with 31.9 % GC, and is the largest sequenced so far. This size is due essentially to the presence of numerous non-conserved hypothetical ORFs. It contains the 14 genes coding for proteins involved in the oxidative phosphorylation, the two rRNA genes, one ORF coding for a ribosomal protein (rps3), and a set of 26 tRNA genes that recognize codons for all amino acids. Seven homing endonucleases are located inside introns. Except atp8, all conserved known genes are in the same orientation. Phylogenetic analysis based on the cox genes agrees with the commonly accepted fungal taxonomy. An uncommon feature of this mitochondrial genome is the presence of a region that contains a set of four, relatively small, nested, inverted repeats enclosing two genes coding for polymerases with an invertron-type structure and three conserved hypothetical genes interpreted as the stable integration of a mitochondrial linear plasmid. The integration of this plasmid seems to be a recent evolutionary event that could have implications in fungal biology. This sequence is available under GenBank accession number AY376688. (c) 2008 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.

Nickel sorption capacity of ground xylem of Quercus ilex trees and effects of selected ligands present in the xylem sap

In this work the influence of four different ligands present in the xylem sap of Quercus ilex (histidine, citric, oxalic and aspartic acids) on Ni(II) adsorption by xylem was investigated. Grinded xylem was trapped in acrylic columns and solutions of Ni(II), in the absence and presence of the four ligands prepared in KNO(3) 0-1 mol L(-1) at pH 5.5, were percolated through the column. Aliquots of solutions were recovered in the column end for Ni determination by graphite furnace atomic absorption spectrometry (GFAAS). The experimental. data to describe Ni sorption by xylem in both the presence and absence of ligands was better explained by the Freundlich isotherm model. The decreasing affinity order of ligands for Ni was: oxalic acid > citric acid > histidine > aspartic acid. On the other hand, the Ni(II) adsorption by xylem increased following the inverse sequence of ligands. Potentiometric titrations of acidic groups were carried out to elucidate the sorption site groups available in Q. ilex xylem. The potentiometric titration has shown three sorption sites: pK(a) 2.6 (57.7% of the sorption sites), related to monobasic aliphatic carboxylic acids or nitrogen aromatic bases, pK(a) 8.1 (9.6%) and pK(a) 9.9 (32.7%), related to phenolic groups. (C) 2008 Elsevier GmbH. All rights reserved.

Effect of hemicellulose and lignin on enzymatic hydrolysis of cellulose from brewer`s spent grain

Enzymatic hydrolysis of brewer`s spent grain in three different forms: original (untreated), pretreated by dilute acid (cellulignin), and pretreated by a sequence of dilute acid and dilute alkali (cellulose pulp), was studied to verify the effect of hemicellulose and lignin on cellulose conversion into glucose. The hydrolysis was carried out using a commercial cellulase concentrate (Celluclast 1.5 L) in an enzyme/substrate ratio of 45 FPU/g, 2% (w/v) substrate concentration, 45 degrees C for 96 h. According to the results, the cellulose hydrolysis was affected by the presence of hemicellulose and/or lignin in the sample. The cellulose conversion ratio (defined as glucose yield + cellobiose yield) from cellulignin was 3.5-times higher than that from untreated sample, whereas from cellulose pulp such value was 4-times higher, correspondent to 91.8% (glucose yield of 85.6%). This best result was probably due to the strong modification in the material structure caused by the hemicellulose and lignin removal from the sample. As a consequence, the cellulose fibers were separated being more susceptible to the enzymatic attack. It was concluded that the lower the hemicellulose and lignin contents in the sample, the higher the efficiency of cellulose hydrolysis. (C) 2007 Elsevier Inc. All rights reserved.

A Constructive Genetic Algorithm for permutation flowshop scheduling

The general flowshop scheduling problem is a production problem where a set of n jobs have to be processed with identical flow pattern on in machines. In permutation flowshops the sequence of jobs is the same on all machines. A significant research effort has been devoted for sequencing jobs in a flowshop minimizing the makespan. This paper describes the application of a Constructive Genetic Algorithm (CGA) to makespan minimization on flowshop scheduling. The CGA was proposed recently as an alternative to traditional GA approaches, particularly, for evaluating schemata directly. The population initially formed only by schemata, evolves controlled by recombination to a population of well-adapted structures (schemata instantiation). The CGA implemented is based on the NEH classic heuristic and a local search heuristic used to define the fitness functions. The parameters of the CGA are calibrated using a Design of Experiments (DOE) approach. The computational results are compared against some other successful algorithms from the literature on Taillard`s well-known standard benchmark. The computational experience shows that this innovative CGA approach provides competitive results for flowshop scheduling; problems. (C) 2007 Elsevier Ltd. All rights reserved.