Nonlinear models and small sample performance of the generalized method of moments

In this paper I explore the issue of nonlinearity (both in the datageneration process and in the functional form that establishes therelationship between the parameters and the data) regarding the poorperformance of the Generalized Method of Moments (GMM) in small samples.To this purpose I build a sequence of models starting with a simple linearmodel and enlarging it progressively until I approximate a standard (nonlinear)neoclassical growth model. I then use simulation techniques to find the smallsample distribution of the GMM estimators in each of the models.

Adult sex ratio effects on male survivorship of Drosophila melanogaster Meigen (Diptera, Drosophilidae)

The behavioral biology has a central role in evolutionary biology mainly because the antagonistic relations that occur in the sexual reproduction. One involves the effect of reproduction on the future life expectation. In this scenario, changes in male operational sex ratio could lead to an increase in mortality due to costs associated with excessive courtship and mating displays. Thus, this work experimentally altered the male sex ratio of Drosophila melanogaster Meigen, 1830, to determine its impact on mortality. The results indicated that mortality increases as the sex ratio changes, including modifications in the survivorship curve type and in the curve concavity, measured by entropy.

Worst-case bounds for the logarithmic loss of predictors

We investigate on-line prediction of individual sequences. Given a class of predictors, the goal is to predict as well as the best predictor in the class, where the loss is measured by the self information (logarithmic) loss function. The excess loss (regret) is closely related to the redundancy of the associated lossless universal code. Using Shtarkov's theorem and tools from empirical process theory, we prove a general upper bound on the best possible (minimax) regret. The bound depends on certain metric properties of the class of predictors. We apply the bound to both parametric and nonparametric classes ofpredictors. Finally, we point out a suboptimal behavior of the popular Bayesian weighted average algorithm.

Comparative phylogeography of mutualists and the effect of the host on the genetic structure of its partners.

Whether or not species participating in specialized and obligate interactions display similar and simultaneous demographic variations at the intraspecific level remains an open question in phylogeography. In the present study, we used the mutualistic nursery pollination occurring between the European globeflower Trollius europaeus and its specialized pollinators in the genus Chiastocheta as a case study. Explicitly, we investigated if the phylogeographies of the pollinating flies are significantly different from the expectation under a scenario of plant-insect congruence. Based on a large-scale sampling, we first used mitochondrial data to infer the phylogeographical histories of each fly species. Then, we defined phylogeographical scenarios of congruence with the plant history, and used maximum likelihood and Bayesian approaches to test for plant-insect phylogeographical congruence for the three Chiastocheta species. We show that the phylogeographical histories of the three fly species differ. Only Chiastocheta lophota and Chiastocheta dentifera display strong spatial genetic structures, which do not appear to be statistically different from those expected under scenarios of phylogeographical congruence with the plant. The results of the present study indicate that the fly species responded in independent and different ways to shared evolutionary forces, displaying varying levels of congruence with the plant genetic structure

Model selection and error estimation

We study model selection strategies based on penalized empirical loss minimization. We point out a tight relationship between error estimation and data-based complexity penalization: any good error estimate may be converted into a data-based penalty function and the performance of the estimate is governed by the quality of the error estimate. We consider several penalty functions, involving error estimates on independent test data, empirical {\sc vc} dimension, empirical {\sc vc} entropy, andmargin-based quantities. We also consider the maximal difference between the error on the first half of the training data and the second half, and the expected maximal discrepancy, a closely related capacity estimate that can be calculated by Monte Carlo integration. Maximal discrepancy penalty functions are appealing for pattern classification problems, since their computation is equivalent to empirical risk minimization over the training data with some labels flipped.

Stable isotope and trace element stratigraphy across the Permian-Triassic transition: A redefinition of the boundary in the Velebit Mountain, Croatia

Stable isotopes of carbonates (delta(13)C(carb), delta(18)O(carb)), organic matter (delta(13)C(org), delta(15)N(org)) and major, trace and rare earth element (REE) compositions of marine carbonate rocks of Late Permian to Early Triassic age were used to establish the position of the Permian-Triassic boundary (PTB) at two continuous sections in the Velebit Mountain, Croatia. The chosen sections - Rizvanusa and Brezimenjaca - are composed of two lithostratigraphic units, the Upper Permian Transitional Dolomite and the overlying Sandy Dolomite. The contact between these units, characterized by the erosional features and sudden occurrence of ooids and siliciclastic grains, was previously considered as the chronostratigraphic PTB. The Sandy Dolomite is characterized by high content of non-carbonate material (up to similar to 30 wt.% insoluble residue), originated from erosion of the uplifted hinterland. A relatively rich assemblage of Permian fossils (including Geinitzina, Globivalvulina, Hemigordius, bioclasts of gastropods, ostracods and brachiopods) was found for the first time in Sandy Dolomite, 5 m above the lithologic boundary in the Rizvanusa section. A rather abrupt negative delta(13)C(carb) excursion in both sections appears in rocks showing no recognizable facies change within the Sandy Dolomite, -2 parts per thousand at Rizvanusa and -1.2 parts per thousand at Brezimenjaca, 11 m and 0.2 m above the lithologic contact, respectively. This level within the lower part of the Sandy Dolomite is proposed as the chemostratigraphic PTB. In the Rizvanusa section, the delta(13)C(org) values decline gradually from similar to-25 parts per thousand in the Upper Permian to similar to-29 parts per thousand in the Lower Triassic. The first negative delta(13)C(org) excursion occurs above the lithologic contact, within the uppermost Permian deposits, and appears to be related to the input of terrigenous material. The release of isotopically light microbial soil-biomass into the shallow-marine water may explain this sudden decrease of delta(13)C(org) values below the PTB. This would support the hypothesis that in the western Tethyan realm the land extinction, triggering a sudden drop of woody vegetation and related land erosion, preceded the marine extinction. The relatively low delta(15)N(org) values at the Permian-Triassic (P-Tr) transition level, close to approximate to 0 parts per thousand, and a secondary negative delta(13)C(org) excursion of -0.5 parts per thousand point to significant terrestrial input and primary contribution of cyanobacteria. The profiles of the concentrations of redox-sensitive elements (Ce, Mn, Fe, V), biogenic or biogenic-scavenged elements (P, Ba, Zn, V), Ce/Ce* values, and normalized trace elements, including Ba/Al, Ba/Fe, Ti/Al, Al/(Al + Fe + Mn) and Mn/Ti show clear excursions at the Transitional Dolomite-Sandy Dolomite lithologic boundary and the chemostratigraphic P-Tr boundary. The stratigraphic variations indicate a major regression phase marking the lithologic boundary, transgressive phases in the latest Permian and a gradual change into shallow/stagnant anoxic marine environment towards the P-Tr boundary level and during the earliest Triassic. (C) 2010 Elsevier B.V. All rights reserved.

Characterization of a major V3 duplication in the NS5A region of genotype 1b HCV by quasispecies analysis in a French multicenter study

Background and Aims: The NS5A protein of the HCV is known tobe involved in viral replication and assembly and probably in theresistance to Interferon based-therapy. Previous studies identifiedinsertions or deletions from 1 to 12 nucleotides in several genomicregions. In a multicenter study (17 French and 1 Swiss laboratoriesof virology), we identified for the first time a 31 amino acidsinsertion leading to a duplication of the V3 domain in the NS5Aregion with a high prevalence. Quasispecies of each strain withduplication were characterized and the inserted V3 domain wasidentified.Methods: Between 2006 and 2008, 1067 patients chronicallyinfected with a 1b HCV were consecutively included in the study.We first amplified the V3 region by RT-PCR to detect duplication(919 samples successfully amplified). The entire NS5A region wasthen amplified, cloned and sequenced in strains bearing theduplication. V3 sequences (called R1 and R2) from each clonewere analyzed with BioEdit and compared to a V3 consensussequence (C) built from the Database Los Alamos Hepatitis C.Entropy was determined at each position.Results: V3 duplications were identified in 25 patients representinga prevalence of 2.72%. We sequenced 2043 clones from which776 had a complete coding NS5A sequence (corresponding toa mean of 30 clones per patient). At the intra-individual level,6 to 17 variants were identified per V3 region, with a maximum of3 different amino acids. At the inter-individual level, a differenceof 7 and 2 amino acids was observed between C and R1 and R2sequences, respectively. Moreover few positions presented entropyhigher than 1 (4 for the R1, 2 for the R2 and 2 for the C). Among allthe sequenced clones, more than 60% were defective virus (partialfragment of NS5A or stop codon).Conclusions: We identified a duplication of the V3 domain ingenotype 1b HCV with a high prevalence. The R2 domain, which wasthe most similar to the C region, might probably be the "original"domain, whereas R1 should be the inserted domain. Phylogeneticanalyses are under process to confirm this hypothesis.

Scenic Byways of Iowa, June 12, 2007

Bypass traffic and experience a “scenic” change of pace by traveling along Iowa’s scenic byways. Iowa’s eight state-designated and two nationally-designated scenic byways are a great way to experience Iowa’s natural beauty, history and culture. Stop to smell the wildflowers or listen to the songbirds, or follow an impulse to take a side trip to one of the many attractions and countryside hamlets. A camera is a must for these postcard-perfect vistas. You never know when you will encounter a bald eagle along the Mississippi River, rare plants and animals in the Loess Hills, or the exceptional architecture of unique barns, churches and other buildings along the routes. This brochure identifies each scenic byway route and the approximate mileage in terms of hard-surfaced and gravel roadways. Estimated driving time ranges from one and one-half hour to three and one-half hours, depending on your speed and the number of stops. These routes are offered for those of you who want to relax and stop often to enjoy the sights.

La relació entre el valor N30 de L'SPT i la resistència a la compressió simple

The unconfined compression strength is a geotechnical parameter frequently used to determine the bearing capacity of foundations on clayey soils. Its value can be obtained from the unconfined compression test and, among several in situ tests, the standard penetration test (SPT). However, a comparison between the values of such a strength determined by both tests in CL, CL-ML and CH soils shows that it does not exist an approximate relationship. In any case it is remarkable that values of the test conducted in the laboratory fit between the lower and the higher ones calculated from the standard penetration number. This is very important for those geotechnical reports in which the unconfined compression strength is only based on the SPT. In this way, if the lower value is taken to evaluate the ultimate soilbearing capacity the result will not be very accurate but it will probably be by the side of safety

Clinicopathologic and molecular analysis of a choroidal pigmented schwannoma in the context of a PTEN hamartoma tumor syndrome.

PURPOSE: To report the first case of choroidal schwannoma in a patient affected by PTEN hamartoma tumor syndrome (PHTS) and investigate the molecular involvement of the phosphatase and tensin homolog (PTEN) and neurofibromin 2 (NF2) genes in this rare intraocular tumor. DESIGN: Observational case report. PARTICIPANT: A 10-year-old girl diagnosed with PHTS. METHODS: The enucleated specimen underwent histologic, immunohistochemical, and transmission electronic microscopy. The expression of PTEN and NF2 and their protein products were evaluated by reverse transcription-polymerase chain reaction and immunohistochemistry. Somatic mutations of PTEN and NF2, as well as allelic loss, were investigated by direct sequencing of DNA extracted from the tumor. PTEN epigenetic silencing was investigated by pyrosequencing. MAIN OUTCOME MEASURES: Histopathologic and molecular characterization of a choroidal pigmented schwannoma. RESULTS: Histopathologic, immunohistochemical, and electron microscopic analysis demonstrated features consistent with a pigmented cellular schwannoma of the choroid. We found no loss of heterozygosity at the genomic level for the PTEN germline mutation and no promoter hypermethylation or other somatic intragenic mutations. However, we observed an approximate 40% reduction of PTEN expression at both the mRNA and the protein level, indicating that the tumor was nonetheless functionally deficient for PTEN. Although DNA sequencing of NF2 failed to identify any pathologic variants, its expression was abolished within the tumor. CONCLUSIONS: We report the first description of a pigmented choroidal schwannoma in the context of a PHTS. This rare tumor showed a unique combination of reduction of PTEN and absence of NF2 expression.

Active learning methods for remote sensing image classification

In this paper, we propose two active learning algorithms for semiautomatic definition of training samples in remote sensing image classification. Based on predefined heuristics, the classifier ranks the unlabeled pixels and automatically chooses those that are considered the most valuable for its improvement. Once the pixels have been selected, the analyst labels them manually and the process is iterated. Starting with a small and nonoptimal training set, the model itself builds the optimal set of samples which minimizes the classification error. We have applied the proposed algorithms to a variety of remote sensing data, including very high resolution and hyperspectral images, using support vector machines. Experimental results confirm the consistency of the methods. The required number of training samples can be reduced to 10% using the methods proposed, reaching the same level of accuracy as larger data sets. A comparison with a state-of-the-art active learning method, margin sampling, is provided, highlighting advantages of the methods proposed. The effect of spatial resolution and separability of the classes on the quality of the selection of pixels is also discussed.

Brownian dynamics simulation of DNA condensation.

DNA condensation observed in vitro with the addition of polyvalent counterions is due to intermolecular attractive forces. We introduce a quantitative model of these forces in a Brownian dynamics simulation in addition to a standard mean-field Poisson-Boltzmann repulsion. The comparison of a theoretical value of the effective diameter calculated from the second virial coefficient in cylindrical geometry with some experimental results allows a quantitative evaluation of the one-parameter attractive potential. We show afterward that with a sufficient concentration of divalent salt (typically approximately 20 mM MgCl(2)), supercoiled DNA adopts a collapsed form where opposing segments of interwound regions present zones of lateral contact. However, under the same conditions the same plasmid without torsional stress does not collapse. The condensed molecules present coexisting open and collapsed plectonemic regions. Furthermore, simulations show that circular DNA in 50% methanol solutions with 20 mM MgCl(2) aggregates without the requirement of torsional energy. This confirms known experimental results. Finally, a simulated DNA molecule confined in a box of variable size also presents some local collapsed zones in 20 mM MgCl(2) above a critical concentration of the DNA. Conformational entropy reduction obtained either by supercoiling or by confinement seems thus to play a crucial role in all forms of condensation of DNA.

Determinación de sistemas de flujo regionales y locales en las depresiones tectónicas del Baix Empordà y La Selva (NE de España) en base a datos hidroquímicos e isotópicos

The Baix Empordà-Selva-Gavarres aquifer system is related to the fault set that created the tectonic basins of Empordà and Selva areas (NE Spain) during the Neogene. In this work, we describe groundwater hydrogeological, hydrochemical and isotopical (3H, δD, δ18O, and the 87Sr/86Sr ratio) characteristics of this system in order to illustrate the relevance of fault zones in groundwater flow-paths and the recharge. In that way, we identify two flow systems, with distinct hydrochemistry and isotopes. A local flow system originates at the Gavarres Range, and it flows towards the basins of the Baix Empordà and Selva, with an approximate residence time of 20 years. Additionally, a regional flow system has only been identified in the Selva basin. This one is related to the main fault zones, as preferential flow paths. Its recharge is located in mountain ranges with higher altitudes, namely the Transversal and Guilleries Ranges, with residence times larger than 50 years. Isotopical data has also shown mixing processes between both flow systems and rainfall recharge while multivariate statistical analysis of principal components has shown the main processes that control hydrochemistry of each flow systems

Extremal climatic states simulated by a 2-dimensional model Part I: Sensitivity of the model and present state

The criterion, based on the thermodynamics theory, that the climatic system tends to extremizesome function has suggested several studies. In particular, special attention has been devoted to the possibility that the climate reaches an extremal rate of planetary entropy production.Due to both radiative and material effects contribute to total planetary entropy production,climatic simulations obtained at the extremal rates of total, radiative or material entropy production appear to be of interest in order to elucidate which of the three extremal assumptions behaves more similar to current data. In the present paper, these results have been obtainedby applying a 2-dimensional (2-Dim) horizontal energy balance box-model, with a few independent variables (surface temperature, cloud-cover and material heat fluxes). In addition, climatic simulations for current conditions by assuming a fixed cloud-cover have been obtained. Finally,sensitivity analyses for both variable and fixed cloud models have been carried out

Extremal climatic states simulated by a 2-dimensional model Part II: Different climatic scenarios

Different climatic simulations have been obtained by using a 2-Dim horizontal energy balancemodel (EBM), which has been constrained to satisfy several extremal principles on dissipationand convection. Moreover, 2 different versions of the model with fixed and variable cloud-coverhave been used. The assumption of an extremal type of behaviour for the climatic system canacquire additional support depending on the similarities found with measured data for pastconditions as well as with usual projections for possible future scenarios