A systematic method of parameterisation estimation using data assimilation

In numerical weather prediction, parameterisations are used to simulate missing physics in the model. These can be due to a lack of scientific understanding or a lack of computing power available to address all the known physical processes. Parameterisations are sources of large uncertainty in a model as parameter values used in these parameterisations cannot be measured directly and hence are often not well known; and the parameterisations themselves are also approximations of the processes present in the true atmosphere. Whilst there are many efficient and effective methods for combined state/parameter estimation in data assimilation (DA), such as state augmentation, these are not effective at estimating the structure of parameterisations. A new method of parameterisation estimation is proposed that uses sequential DA methods to estimate errors in the numerical models at each space-time point for each model equation. These errors are then fitted to pre-determined functional forms of missing physics or parameterisations that are based upon prior information. We applied the method to a one-dimensional advection model with additive model error, and it is shown that the method can accurately estimate parameterisations, with consistent error estimates. Furthermore, it is shown how the method depends on the quality of the DA results. The results indicate that this new method is a powerful tool in systematic model improvement.

Environmental controls on the distribution and diversity of lentic Chironomidae (Insecta: Diptera) across an altitudinal gradient in tropical South America

To predict the response of aquatic ecosystems to future global climate change, data on the ecology and distribution of keystone groups in freshwater ecosystems are needed. In contrast to mid- and high-latitude zones, such data are scarce across tropical South America (Neotropics). We present the distribution and diversity of chironomid species using surface sediments of 59 lakes from the Andes to the Amazon (0.1–17°S and 64–78°W) within the Neotropics. We assess the spatial variation in community assemblages and identify the key variables influencing the distributional patterns. The relationships between environmental variables (pH, conductivity, depth, and sediment organic content), climatic data, and chironomid assemblages were assessed using multivariate statistics (detrended correspondence analysis and canonical correspondence analysis). Climatic parameters (temperature and precipitation) were most significant in describing the variance in chironomid assemblages. Temperature and precipitation are both predicted to change under future climate change scenarios in the tropical Andes. Our findings suggest taxa of Orthocladiinae, which show a preference to cold high-elevation oligotrophic lakes, will likely see range contraction under future anthropogenic-induced climate change. Taxa abundant in areas of high precipitation, such as Micropsectra and Phaenopsectra, will likely become restricted to the inner tropical Andes, as the outer tropical Andes become drier. The sensitivity of chironomids to climate parameters makes them important bio-indicators of regional climate change in the Neotropics. Furthermore, the distribution of chironomid taxa presented here is a vital first step toward providing urgently needed autecological data for interpreting fossil chironomid records of past ecological and climate change from the tropical Andes.

Sparse density estimation on multinomial manifold combining local component analysis

A new sparse kernel density estimator is introduced based on the minimum integrated square error criterion combining local component analysis for the finite mixture model. We start with a Parzen window estimator which has the Gaussian kernels with a common covariance matrix, the local component analysis is initially applied to find the covariance matrix using expectation maximization algorithm. Since the constraint on the mixing coefficients of a finite mixture model is on the multinomial manifold, we then use the well-known Riemannian trust-region algorithm to find the set of sparse mixing coefficients. The first and second order Riemannian geometry of the multinomial manifold are utilized in the Riemannian trust-region algorithm. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with competitive accuracy to existing kernel density estimators.

Gaze stabilization testing: The efficiency and test-retest reliability of user operated control versus the parameter estimation (PEST) algorithm

The primary objective of this research study is to determine which form of testing, the PEST algorithm or an operator-controlled condition is most accurate and time efficient for administration of the gaze stabilization test

Sex and socioeconomic inequalities of lung cancer mortality in Barcelona, Spain and Sao Paulo, Brazil

The objective of this paper was to assess sex and socioeconomic inequalities in lung cancer mortality in two major cities of Europe and South America. Official information on mortality and population allowed the estimation of sex- and age-specific death rates for Barcelona, Spain and Sao Paulo, Brazil (1995-2003). Mortality trends and levels were independently assessed for each city and subsequently compared. Rate ratios assessed by Poisson regression analysis addressed hypotheses of association between the outcome and socioeconomic covariates (human development index, unemployment and schooling) at the inner-city area level. Barcelona had a higher mortality in men (76.9/100000 inhabitants) than Sao Paulo (38.2/100 000 inhabitants); although rates were decreasing for the former (-2%/year) and levelled-off for the [after. Mortality in women ranked similarly (9.1 for Barcelona, 11.5 for Sao Paulo); with an increasing trend for women aged 35-64 years (+ 7.7%/year in Barcelona and + 2.4%/year in Sao Paulo). The socioeconomic gradient of mortality in men was negative for Barcelona and positive for Sao Paulo; for women, the socioeconomic gradient was positive in both cities. Negative gradients indicate that deprived areas suffer a higher burden of disease; positive gradients suggest that prosmoking lifestyles may have been more prevalent in more affluent areas during the last decades. Sex and socioeconomic inequalities of lung cancer mortality reinforce the hypothesis that the epidemiologic profile of cancer can be improved by an expanded access to existing technology of healthcare and prevention. The continuous monitoring of inequalities in health may contribute to the concurrent promotion of well-being and social justice.

DETECTABILITY AND ERROR ESTIMATION IN ORBITAL FITS OF RESONANT EXTRASOLAR PLANETS

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than similar to 4. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are currently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically overestimated, although the inner planet`s eccentricity suffers a much smaller effect. If the initial conditions correspond to small-amplitude oscillations around stable apsidal corotation resonances, the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

Integrated marine gravity field in the Brazilian coast from altimeter-derived sea surface gradient and shipborne gravity

Sea surface gradients derived from the Geosat and ERS-1 satellite altimetry geodetic missions were integrated with marine gravity data from the National Geophysical Data Center and Brazilian national surveys. Using the least squares collocation method, models of free-air gravity anomaly and geoid height were calculated for the coast of Brazil with a resolution of 2` x 2`. The integration of satellite and shipborne data showed better statistical results in regions near the coast than using satellite data only, suggesting an improvement when compared to the state-of-the-art global gravity models. Furthermore, these results were obtained with considerably less input information than was used by those reference models. The least squares collocation presented a very low content of high-frequency noise in the predicted gravity anomalies. This may be considered essential to improve the high resolution representation of the gravity field in regions of ocean-continent transition. (C) 2010 Elsevier Ltd. All rights reserved.

Germination rate and sporeling development of Chondracanthus chamissoi (Rhodophyta, Gigartinales) varies along a latitudinal gradient on the coast of Chile

Chondracanthus chamissoi (C. Agardh) Kutzing is an economically important red seaweed with an extended latitudinal distribution along the south-east Pacific. Here we report on the seasonal in vitro germination of carpospores and tetraspores from four populations distributed from 27 to 41 degrees S on the Chilean coast. Our results show that both types of spores exhibited a different physiological behavior related to the geographic origin of the specimens. Germination occurred throughout the year for both spore types in the four populations. However, for the northern locations (Calderilla, La Herradura and Puerto Aldea) germination was higher in spring, while for the southern location (Lechagua), germination was higher in summer. The growth rate of carposporelings and tetrasporelings varied seasonally in ail locations studied, with higher growth in spring. Among all, carposporelings from Lechagua specimens reached the highest growth rates (9.3 +/- 0.2% d(-1)). However, spores from Herradura and P. Aldea had a good germination and SGR in all seasons and would be good candidates to start spores-based cultivation of this valuable resource in Chile. (C) 2009 Elsevier B.V. All rights reserved

Seed dispersal and predation in the endemic Atlantic rainforest palm Astrocaryum aculeatissimum across a gradient of seed disperser abundance

Tropical forests have been subject to intense hunting of medium and large frugivores that are important in dispersing large-seeded species. It has been hypothesized that in areas with extinction or low abundance of medium and large-bodied animals the density of small rodents may increase. Therefore, this increment in the density of small rodents may compensate for the absence or low abundance of medium and large frugivores on seed removal and seed dispersal. Here, we fill up this gap in the literature by determining if seed removal, seed dispersal, and seed predation by small rodents (spiny rats, Trinomys inheringi and squirrels, Sciurus ingrami) are maintained in defaunated areas. We accessed seed removal, seed dispersal, seed predation, and seedling recruitment of an endemic Atlantic rainforest palm, Astrocaryum aculeatissimum, in a gradient of abundance of agoutis. We found that seed removal, scatter hoarding, and seed predation increase with the abundance of agoutis. In contrast, the proportion of dispersed but non-cached seeds decreased with the abundance of agoutis. We did not find any effect of the abundance of agoutis on seed dispersal distance, but we did find a positive trend on the density of seedlings. We concluded that small rodents do not compensate the low abundance of agoutis on seed removal, scatter hoarding, and seed predation of this palm tree. Moreover, areas in which agoutis are already extinct did not present any seed removal or scatter hoarding, not even by small rodents. This study emphasizes both the importance of agoutis in dispersing seeds of A. aculeatissimum and the collapse in seed dispersal of this palm in areas where agoutis are already extinct.

Decomposition in tropical forests: a pan-tropical study of the effects of litter type, litter placement and mesofaunal exclusion across a precipitation gradient

1. Litter decomposition recycles nutrients and causes large fluxes of carbon dioxide into the atmosphere. It is typically assumed that climate, litter quality and decomposer communities determine litter decay rates, yet few comparative studies have examined their relative contributions in tropical forests. 2. We used a short-term litterbag experiment to quantify the effects of litter quality, placement and mesofaunal exclusion on decomposition in 23 tropical forests in 14 countries. Annual precipitation varied among sites (760-5797 mm). At each site, two standard substrates (Raphia farinifera and Laurus nobilis) were decomposed in fine- and coarse-mesh litterbags both above and below ground for approximately 1 year. 3. Decomposition was rapid, with >95% mass loss within a year at most sites. Litter quality, placement and mesofaunal exclusion all independently affected decomposition, but the magnitude depended upon site. Both the average decomposition rate at each site and the ratio of above- to below-ground decay increased linearly with annual precipitation, explaining 60-65% of among-site variation. Excluding mesofauna had the largest impact on decomposition, reducing decomposition rates by half on average, but the magnitude of decrease was largely independent of climate. This suggests that the decomposer community might play an important role in explaining patterns of decomposition among sites. Which litter type decomposed fastest varied by site, but was not related to climate. 4. Synthesis. A key goal of ecology is to identify general patterns across ecological communities, as well as relevant site-specific details to understand local dynamics. Our pan-tropical study shows that certain aspects of decomposition, including average decomposition rates and the ratio of above- to below-ground decomposition are highly correlated with a simple climatic index: mean annual precipitation. However, we found no relationship between precipitation and effects of mesofaunal exclusion or litter type, suggesting that site-specific details may also be required to understand how these factors affect decomposition at local scales.

The genetic effects of Late Quaternary climatic changes over a tropical latitudinal gradient: diversification of an Atlantic Forest passerine

The increase in biodiversity from high to low latitudes is a widely recognized biogeographical pattern. According to the latitudinal gradient hypothesis (LGH), this pattern was shaped by differential effects of Late Quaternary climatic changes across a latitudinal gradient. Here, we evaluate the effects of climatic changes across a tropical latitudinal gradient and its implications to diversification of an Atlantic Forest (AF) endemic passerine. We studied the intraspecific diversification and historical demography of Sclerurus scansor, based on mitochondrial (ND2, ND3 and cytb) and nuclear (FIB7) gene sequences. Phylogenetic analyses recovered three well-supported clades associated with distinct latitudinal zones. Coalescent-based methods were applied to estimate divergence times and changes in effective population sizes. Estimates of divergence times indicate that intraspecific diversification took place during Middle-Late Pleistocene. Distinct demographic scenarios were identified, with the southern lineage exhibiting a clear signature of demographic expansion, while the central one remained more stable. The northern lineage, contrasting with LGH predictions, exhibited a clear sign of a recent bottleneck. Our results suggest that different AF regions reacted distinctly, even in opposite ways, under the same climatic period, producing simultaneously favourable scenarios for isolation and contact among populations.

Evolution of Dendrocolaptes platyrostris (Aves: Furnariidae) between the South American open vegetation corridor and the Atlantic forest

The open vegetation corridor of South America is a region dominated by savanna biomes. It contains forests (i.e. riverine forests) that may act as corridors for rainforest specialists between the open vegetation corridor and its neighbouring biomes (i.e. the Amazonian and Atlantic forests). A prediction for this scenario is that populations of rainforest specialists in the open vegetation corridor and in the forested biomes show no significant genetic divergence. We addressed this hypothesis by studying plumage and genetic variation of the Planalto woodcreeper Dendrocolaptes platyrostris Spix (1824) (Aves: Furnariidae), a forest specialist that occurs in both open habitat and in the Atlantic forest. The study questions were: (1) is there any evidence of genetic continuity between populations of the open habitat and the Atlantic forest and (2) is plumage variation congruent with patterns of neutral genetic structure or with ecological factors related to habitat type? We used cytochrome b and mitochondrial DNA control region sequences to show that D. platyrostris is monophyletic and presents substantial intraspecific differentiation. We found two areas of plumage stability: one associated with Cerrado and the other associated with southern Atlantic Forest. Multiple Mantel tests showed that most of the plumage variation followed the transition of habitats but not phylogeographical gaps, suggesting that selection may be related to the evolution of the plumage of the species. The results were not compatible with the idea that forest specialists in the open vegetation corridor and in the Atlantic forest are linked at the population level because birds from each region were not part of the same genetic unit. Divergence in the presence of gene flow across the ecotone between both regions might explain our results. Also, our findings indicate that the southern Atlantic forest may have been significantly affected by Pleistocene climatic alteration, although such events did not cause local extinction of most taxa, as occurred in other regions of the globe where forests were significantly affected by global glaciations. Finally, our results neither support plumage stability areas, nor subspecies as full species. (C) 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103, 801-820.

Bayesian Estimation for Performance Measures of Two Diagnostic Tests in the Presence of Verification Bias

Sensitivity and specificity are measures that allow us to evaluate the performance of a diagnostic test. In practice, it is common to have situations where a proportion of selected individuals cannot have the real state of the disease verified, since the verification could be an invasive procedure, as occurs with biopsy. This happens, as a special case, in the diagnosis of prostate cancer, or in any other situation related to risks, that is, not practicable, nor ethical, or in situations with high cost. For this case, it is common to use diagnostic tests based only on the information of verified individuals. This procedure can lead to biased results or workup bias. In this paper, we introduce a Bayesian approach to estimate the sensitivity and the specificity for two diagnostic tests considering verified and unverified individuals, a result that generalizes the usual situation based on only one diagnostic test.

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.

Methods of Estimation in Multiple Linear Regression: Application to Clinical Data

Nesse artigo, tem-se o interesse em avaliar diferentes estratégias de estimação de parâmetros para um modelo de regressão linear múltipla. Para a estimação dos parâmetros do modelo foram utilizados dados de um ensaio clínico em que o interesse foi verificar se o ensaio mecânico da propriedade de força máxima (EM-FM) está associada com a massa femoral, com o diâmetro femoral e com o grupo experimental de ratas ovariectomizadas da raça Rattus norvegicus albinus, variedade Wistar. Para a estimação dos parâmetros do modelo serão comparadas três metodologias: a metodologia clássica, baseada no método dos mínimos quadrados; a metodologia Bayesiana, baseada no teorema de Bayes; e o método Bootstrap, baseado em processos de reamostragem.