Do stacked species distribution models reflect altitudinal diversity patterns?

The objective of this study was to evaluate the performance of stacked species distribution models in predicting the alpha and gamma species diversity patterns of two important plant clades along elevation in the Andes. We modelled the distribution of the species in the Anthurium genus (53 species) and the Bromeliaceae family (89 species) using six modelling techniques. We combined all of the predictions for the same species in ensemble models based on two different criteria: the average of the rescaled predictions by all techniques and the average of the best techniques. The rescaled predictions were then reclassified into binary predictions (presence/absence). By stacking either the original predictions or binary predictions for both ensemble procedures, we obtained four different species richness models per taxa. The gamma and alpha diversity per elevation band (500 m) was also computed. To evaluate the prediction abilities for the four predictions of species richness and gamma diversity, the models were compared with the real data along an elevation gradient that was independently compiled by specialists. Finally, we also tested whether our richness models performed better than a null model of altitudinal changes of diversity based on the literature. Stacking of the ensemble prediction of the individual species models generated richness models that proved to be well correlated with the observed alpha diversity richness patterns along elevation and with the gamma diversity derived from the literature. Overall, these models tend to overpredict species richness. The use of the ensemble predictions from the species models built with different techniques seems very promising for modelling of species assemblages. Stacking of the binary models reduced the over-prediction, although more research is needed. The randomisation test proved to be a promising method for testing the performance of the stacked models, but other implementations may still be developed.

Optimization of particle bombardment parameters for the genetic transformation of Brazilian maize inbred lines

The objective of this work was to develop a genetic transformation system for tropical maize genotypes via particle bombardment of immature zygotic embryos. Particle bombardment was carried out using a genetic construct with bar and uidA genes under control of CaMV35S promoter. The best conditions to transform maize tropical inbred lines L3 and L1345 were obtained when immature embryos were cultivated, prior to the bombardment, in higher osmolarity during 4 hours and bombarded at an acceleration helium gas pressure of 1,100 psi, two shots per plate, and a microcarrier flying distance of 6.6 cm. Transformation frequencies obtained using these conditions ranged from 0.9 to 2.31%. Integration of foreign genes into the genome of maize plants was confirmed by Southern blot analysis as well as bar and uidA gene expressions. The maize genetic transformation protocol developed in this work will possibly improve the efficiency to produce new transgenic tropical maize lines expressing desirable agronomic characteristics.

Investigation of Negative Moment Reinforcing in Bridge Decks, TR-660, 2015

Multi-span pre-tensioned pre-stressed concrete beam (PPCB) bridges made continuous usually experience a negative live load moment region over the intermediate supports. Conventional thinking dictates that sufficient reinforcement must be provided in this region to satisfy the strength and serviceability requirements associated with the tensile stresses in the deck. The American Association of State Highway and Transportation Officials (AASHTO) Load and Resistance Factor Design (LRFD) Bridge Design Specifications recommend the negative moment reinforcement (b2 reinforcement) be extended beyond the inflection point. Based upon satisfactory previous performance and judgment, the Iowa Department of Transportation (DOT) Office of Bridges and Structures (OBS) currently terminates b2 reinforcement at 1/8 of the span length. Although the Iowa DOT policy results in approximately 50% shorter b2 reinforcement than the AASHTO LRFD specifications, the Iowa DOT has not experienced any significant deck cracking over the intermediate supports. The primary objective of this project was to investigate the Iowa DOT OBS policy regarding the required amount of b2 reinforcement to provide the continuity over bridge decks. Other parameters, such as termination length, termination pattern, and effects of the secondary moments, were also studied. Live load tests were carried out on five bridges. The data were used to calibrate three-dimensional finite element models of two bridges. Parametric studies were conducted on the bridges with an uncracked deck, a cracked deck, and a cracked deck with a cracked pier diaphragm for live load and shrinkage load. The general conclusions were as follows: -- The parametric study results show that an increased area of the b2 reinforcement slightly reduces the strain over the pier, whereas an increased length and staggered reinforcement pattern slightly reduce the strains of the deck at 1/8 of the span length. -- Finite element modeling results suggest that the transverse field cracks over the pier and at 1/8 of the span length are mainly due to deck shrinkage. -- Bridges with larger skew angles have lower strains over the intermediate supports. -- Secondary moments affect the behavior in the negative moment region. The impact may be significant enough such that no tensile stresses in the deck may be experienced.

Cement Stabilization of Embankment Materials, 2015

Embankment subgrade soils in Iowa are generally rated as fair to poor as construction materials. These soils can exhibit low bearing strength, high volumetric instability, and freeze/thaw or wet/dry durability problems. Cement stabilization offers opportunities to improve these soils conditions. The objective of this study was to develop relationships between soil index properties, unconfined compressive strength and cement content. To achieve this objective, a laboratory study was conducted on 28 granular and non-granular materials obtained from 9 active construction sites in Iowa. The materials consisted of glacial till, loess, and alluvium sand. Type I/II portland cement was used for stabilization. Stabilized and unstabilized specimens were prepared using Iowa State University 2 in. by 2 in. compaction apparatus. Specimens were prepared, cured, and tested for unconfined compressive strength (UCS) with and without vacuum saturation. Percent fines content (F200), AASHTO group index (GI), and Atterberg limits were tested before and after stabilization. The results were analyzed using multi-variate statistical analysis to assess influence of the various soil index properties on post-stabilization material properties. Results indicated that F200, liquid limit, plasticity index, and GI of the materials generally decreased with increasing cement content. The UCS of the stabilized specimens increased with increasing cement content, as expected. The average saturated UCS of the unstabilized materials varied between 0 and 57 psi. The average saturated UCS of stabilized materials varied between 44 and 287 psi at 4% cement content, 108 and 528 psi at t 8% cement content, and 162 and 709 psi at 12% cement content. The UCS of the vacuum saturated specimens was on average 1.5 times lower than that of the unsaturated specimens. Multi-variate statistical regression models are provided in this report to predict F200, plasticity index, GI, and UCS after treatment, as a function of cement content and soil index properties.

Estimating parameters for generalized mass action models using constraint propagation.

As modern molecular biology moves towards the analysis of biological systems as opposed to their individual components, the need for appropriate mathematical and computational techniques for understanding the dynamics and structure of such systems is becoming more pressing. For example, the modeling of biochemical systems using ordinary differential equations (ODEs) based on high-throughput, time-dense profiles is becoming more common-place, which is necessitating the development of improved techniques to estimate model parameters from such data. Due to the high dimensionality of this estimation problem, straight-forward optimization strategies rarely produce correct parameter values, and hence current methods tend to utilize genetic/evolutionary algorithms to perform non-linear parameter fitting. Here, we describe a completely deterministic approach, which is based on interval analysis. This allows us to examine entire sets of parameters, and thus to exhaust the global search within a finite number of steps. In particular, we show how our method may be applied to a generic class of ODEs used for modeling biochemical systems called Generalized Mass Action Models (GMAs). In addition, we show that for GMAs our method is amenable to the technique in interval arithmetic called constraint propagation, which allows great improvement of its efficiency. To illustrate the applicability of our method we apply it to some networks of biochemical reactions appearing in the literature, showing in particular that, in addition to estimating system parameters in the absence of noise, our method may also be used to recover the topology of these networks.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Degree of multicollinearity and variables involved in linear dependence in additive-dominant models

The objective of this work was to assess the degree of multicollinearity and to identify the variables involved in linear dependence relations in additive-dominant models. Data of birth weight (n=141,567), yearling weight (n=58,124), and scrotal circumference (n=20,371) of Montana Tropical composite cattle were used. Diagnosis of multicollinearity was based on the variance inflation factor (VIF) and on the evaluation of the condition indexes and eigenvalues from the correlation matrix among explanatory variables. The first model studied (RM) included the fixed effect of dam age class at calving and the covariates associated to the direct and maternal additive and non-additive effects. The second model (R) included all the effects of the RM model except the maternal additive effects. Multicollinearity was detected in both models for all traits considered, with VIF values of 1.03 - 70.20 for RM and 1.03 - 60.70 for R. Collinearity increased with the increase of variables in the model and the decrease in the number of observations, and it was classified as weak, with condition index values between 10.00 and 26.77. In general, the variables associated with additive and non-additive effects were involved in multicollinearity, partially due to the natural connection between these covariables as fractions of the biological types in breed composition.

Predictors for the emergence of the 2 multi-nucleoside/nucleotide resistance mutations 69 insertion and Q151M and their impact on clinical outcome in the Swiss HIV cohort study.

The 69 insertion and Q151M mutations are multi-nucleoside/nucleotide resistance mutations (MNR). The prevalence among 4078 antiretroviral therapy (ART)-experienced individuals was <1.3%. Combined ART fully prevented MNR in subtype B infections. Case-control studies were performed to identify risk factors. Control subjects were patients with ≥ 3 thymidine-analogue mutations. The 69 insertion study (27 control subjects, 14 case patients) identified didanosine exposure as a risk (odds ratio, 5.0 per year; P = .019), whereas the Q151M study (which included 44 control subjects and 25 case patients) detected no associations. Following detection, individuals with Q151M tended to have lower suppression rates and higher mortality rates, relative to control subjects. Additional studies are needed to verify these findings in non-subtype B infections.

A multi-temporal approach to model endangered species distribution in Europe: the case of the Eurasian otter in Italy

We analysed the relationship between changes in land cover patterns and the Eurasian otter occurrence over the course of about 20 years (1985-2006) using multi-temporal Species Distribution Models (SDMs). The study area includes five river catchments covering most of the otter's Italian range. Land cover and topographic data were used as proxies of the ecological requirements of the otter within a 300-m buffer around river courses. We used species presence, pseudo-absence data, and environmental predictors to build past (1985) and current (2006) SDMs by applying an ensemble procedure through the BIOMOD modelling package. The performance of each model was evaluated by measuring the area under the curve (AUC) of the receiver-operating characteristic (ROC). Multi-temporal analyses of species distribution and land cover maps were performed by comparing the maps produced for 1985 and 2006. The ensemble procedure provided a good overall modelling accuracy, revealing that elevation and slope affected the otter's distribution in the past; in contrast, land cover predictors, such as cultivations and forests, were more important in the present period. During the transition period, 20.5% of the area became suitable, with 76% of the new otter presence data being located in these newly available areas. The multi-temporal analysis suggested that the quality of otter habitat improved in the last 20 years owing to the expansion of forests and to the reduction of cultivated fields in riparian belts. The evidence presented here stresses the great potential of riverine habitat restoration and environmental management for the future expansion of the otter in Italy

Comments on "On the relationship between fractal dimension and the performance of multi-resonant dipole antennas using Koch curves"

Comentaris referits a l'article següent: K. J. Vinoy, J. K. Abraham, and V. K. Varadan, “On the relationshipbetween fractal dimension and the performance of multi-resonant dipoleantennas using Koch curves,” IEEE Transactions on Antennas and Propagation, 2003, vol. 51, p. 2296–2303.

The IST METRA Project

This article summarizes the main achievementsof the Multi-Element Transmit andReceive Antennas (METRA) Project, an ISTresearch and technological development project carried out between January 2000 and June 2001 by Universitat Politècnica de Catalunya, the Center for Personkommunikation of Aalborg University, Nokia Networks, Nokia Mobile Phones, and Vodafone Group Research and Development.The main objective of METRA was the performanceevaluation of multi-antenna terminals incombination with adaptive antennas at the basestation in UMTS communication systems. 1 AMIMO channel sounder was developed that providedrealistic multi-antenna channel measurements.Using these measured data, stochasticchannel models were developed and properly validated.These models were also evaluated inorder to estimate their corresponding channelcapacity. Different MIMO configurations andprocessing schemes were developed for both theFDD and TDD modes of UTRA, and their linkperformance was assessed. Performance evaluationwas completed by system simulations thatillustrated the benefits of MIMO configurationsto the network operator. Implementation cost vs.performance improvement was also covered bythe project, including the base station and terminalmanufacturer and network operator viewpoints.Finally, significant standards contributionswere generated by the project and presented to the pertinent 3GPP working groups.

Analysis of Industrial Granular Flow Applications by Using Advanced Collision Models

Granular flow phenomena are frequently encountered in the design of process and industrial plants in the traditional fields of the chemical, nuclear and oil industries as well as in other activities such as food and materials handling. Multi-phase flow is one important branch of the granular flow. Granular materials have unusual kinds of behavior compared to normal materials, either solids or fluids. Although some of the characteristics are still not well-known yet, one thing is confirmed: the particle-particle interaction plays a key role in the dynamics of granular materials, especially for dense granular materials. At the beginning of this thesis, detailed illustration of developing two models for describing the interaction based on the results of finite-element simulation, dimension analysis and numerical simulation is presented. The first model is used to describing the normal collision of viscoelastic particles. Based on some existent models, more parameters are added to this model, which make the model predict the experimental results more accurately. The second model is used for oblique collision, which include the effects from tangential velocity, angular velocity and surface friction based on Coulomb's law. The theoretical predictions of this model are in agreement with those by finite-element simulation. I n the latter chapters of this thesis, the models are used to predict industrial granular flow and the agreement between the simulations and experiments also shows the validation of the new model. The first case presents the simulation of granular flow passing over a circular obstacle. The simulations successfully predict the existence of a parabolic steady layer and show how the characteristics of the particles, such as coefficients of restitution and surface friction affect the separation results. The second case is a spinning container filled with granular material. Employing the previous models, the simulation could also reproduce experimentally observed phenomena, such as a depression in the center of a high frequency rotation. The third application is about gas-solid mixed flow in a vertically vibrated device. Gas phase motion is added to coherence with the particle motion. The governing equations of the gas phase are solved by using the Large eddy simulation (LES) and particle motion is predicted by using the Lagrangian method. The simulation predicted some pattern formation reported by experiment.

Random regression models to estimate genetic parameters for milk production of Guzerat cows using orthogonal Legendre polynomials

The objective of this work was to compare random regression models for the estimation of genetic parameters for Guzerat milk production, using orthogonal Legendre polynomials. Records (20,524) of test-day milk yield (TDMY) from 2,816 first-lactation Guzerat cows were used. TDMY grouped into 10-monthly classes were analyzed for additive genetic effect and for environmental and residual permanent effects (random effects), whereas the contemporary group, calving age (linear and quadratic effects) and mean lactation curve were analized as fixed effects. Trajectories for the additive genetic and permanent environmental effects were modeled by means of a covariance function employing orthogonal Legendre polynomials ranging from the second to the fifth order. Residual variances were considered in one, four, six, or ten variance classes. The best model had six residual variance classes. The heritability estimates for the TDMY records varied from 0.19 to 0.32. The random regression model that used a second-order Legendre polynomial for the additive genetic effect, and a fifth-order polynomial for the permanent environmental effect is adequate for comparison by the main employed criteria. The model with a second-order Legendre polynomial for the additive genetic effect, and that with a fourth-order for the permanent environmental effect could also be employed in these analyses.

Semivariogram models for estimating fig fly population density throughout the year

The objective of this work was to select semivariogram models to estimate the population density of fig fly (Zaprionus indianus; Diptera: Drosophilidae) throughout the year, using ordinary kriging. Nineteen monitoring sites were demarcated in an area of 8,200 m2, cropped with six fruit tree species: persimmon, citrus, fig, guava, apple, and peach. During a 24 month period, 106 weekly evaluations were done in these sites. The average number of adult fig flies captured weekly per trap, during each month, was subjected to the circular, spherical, pentaspherical, exponential, Gaussian, rational quadratic, hole effect, K-Bessel, J-Bessel, and stable semivariogram models, using ordinary kriging interpolation. The models with the best fit were selected by cross-validation. Each data set (months) has a particular spatial dependence structure, which makes it necessary to define specific models of semivariograms in order to enhance the adjustment to the experimental semivariogram. Therefore, it was not possible to determine a standard semivariogram model; instead, six theoretical models were selected: circular, Gaussian, hole effect, K-Bessel, J-Bessel, and stable.

Optimointijärjestelmän kehittäminen sähköä ja lämpöä tuottavaan voimalaitokseen

Työn tavoitteena oli kehittää automaattinen optimointijärjestelmä energiayhtiön omistamaan pieneen sähkön- ja lämmöntuotantolaitokseen (CHP-laitos). Optimointitarve perustuu energiayhtiön sähkön hankintaan sähköpörssistä, kaasun hankintahintaan, kohteen paikallisiin sähkö- ja lämpökuormituksiin ja muihin laitoksen talouteen vaikuttaviin tekijöihin. Kehitettävällä optimointijärjestelmällä ontarkoitus tulevaisuudessa hallita useita hajautetun energiantuotannon yksiköitäkeskitetysti. Työssä kehitettiin algoritmi, joka optimoi voimalaitoksen taloutta sähkötehoa säätävillä ajomalleilla ja suoralla sähköteho-ohjeella. Työssä kehitetyn algoritmin tuottamia hyötyjä selvitettiin Harjun oppimiskeskuksen CHP-laitoksen mittaushistoriatiedoilla. CHP-laitosten käytön optimointiin luotiin keskitettyyn laskentaan ja hajautettuun ohjaukseen perustuva järjestelmä. Se ohjaa CHP-laitoksia reaaliaikaisesti ja ennustaa historiatietoihin perustuvalla aikasarjamallilla laitoksen tulevaa käyttöä. Optimointijärjestelmän toimivuus ja saatu hyöty selvitettiin Harjun oppimiskeskuksen CHP-laitoksella vertaamalla mittauksista laskettua toteutunutta hyötyä optimointijärjestelmän laskemaan ennustettuun hyötyyn.