The equivalent-weights particle filter in a high-dimensional system

In general, particle filters need large numbers of model runs in order to avoid filter degeneracy in high-dimensional systems. The recently proposed, fully nonlinear equivalent-weights particle filter overcomes this requirement by replacing the standard model transition density with two different proposal transition densities. The first proposal density is used to relax all particles towards the high-probability regions of state space as defined by the observations. The crucial second proposal density is then used to ensure that the majority of particles have equivalent weights at observation time. Here, the performance of the scheme in a high, 65 500 dimensional, simplified ocean model is explored. The success of the equivalent-weights particle filter in matching the true model state is shown using the mean of just 32 particles in twin experiments. It is of particular significance that this remains true even as the number and spatial variability of the observations are changed. The results from rank histograms are less easy to interpret and can be influenced considerably by the parameter values used. This article also explores the sensitivity of the performance of the scheme to the chosen parameter values and the effect of using different model error parameters in the truth compared with the ensemble model runs.

Application of dried blood spots to determine vitamin D status in a large nutritional study with unsupervised sampling: the Food4Me project

An efficient and robust method to measure vitamin D (25-hydroxy vitamin D3 (25(OH)D3) and 25-hydroxy vitamin D2 in dried blood spots (DBS) has been developed and applied in the pan-European multi-centre, internet-based, personalised nutrition intervention study Food4Me. The method includes calibration with blood containing endogenous 25(OH)D3, spotted as DBS and corrected for haematocrit content. The methodology was validated following international standards. The performance characteristics did not reach those of the current gold standard liquid chromatography-MS/MS in plasma for all parameters, but were found to be very suitable for status-level determination under field conditions. DBS sample quality was very high, and 3778 measurements of 25(OH)D3 were obtained from 1465 participants. The study centre and the season within the study centre were very good predictors of 25(OH)D3 levels (P<0·001 for each case). Seasonal effects were modelled by fitting a sine function with a minimum 25(OH)D3 level on 20 January and a maximum on 21 July. The seasonal amplitude varied from centre to centre. The largest difference between winter and summer levels was found in Germany and the smallest in Poland. The model was cross-validated to determine the consistency of the predictions and the performance of the DBS method. The Pearson's correlation between the measured values and the predicted values was r 0·65, and the sd of their differences was 21·2 nmol/l. This includes the analytical variation and the biological variation within subjects. Overall, DBS obtained by unsupervised sampling of the participants at home was a viable methodology for obtaining vitamin D status information in a large nutritional study.

Mathematical relationships between scoring parameters in capped tendering

Mathematical relationships between Scoring Parameters can be used in Economic Scoring Formulas (ESF) in tendering to distribute the score among bidders in the economic part of a proposal. Each contracting authority must set an ESF when publishing tender specifications and the strategy of each bidder will differ depending on the ESF selected and the weight of the overall proposal scoring. This paper introduces the various mathematical relationships and density distributions that describe and inter-relate not only the main Scoring Parameters but the main Forecasting Parameters in any capped tender (those whose price is upper-limited). Forecasting Parameters, as variables that can be known in advance before the deadline of a tender is reached, together with Scoring Parameters constitute the basis of a future Bid Tender Forecasting Model.

Development of films based on blends of gelatin and poly(vinyl alcohol) cross linked with glutaraldehyde

Both gelatin and poly(vinyl alcohol) (PVA) can be cross linked with glutaraldehyde (GLU). In the case of gelatin, the GLU reacts with each e-NH2 functional group of adjacent lysine residues, while for PVA, the GLU reacts with two adjacent hydroxyl groups, forming acetal bridges. Thus it can be considered possible to cross link adjacent macromolecules of gelatin and PVA using GLU. In this context, the aims of this work were the development of biodegradable films based on blends of gelatin and poly(vinyl alcohol) cross linked with GLU, and the characterization of some of their main physical and functional properties. All the films were produced from film-forming solutions (FFS) containing 2 g macromolecules (PVA + gelatin)/100 g FFS, 25 g glycerol/100 g macromolecules, and 4 g GLU (25% solution)/100 g FFS. The FFS were prepared with two concentrations of PVA (20 or 50 g PVA/100 g macromolecules) and two reaction temperatures: 90 or 55 degrees C, applied for 30 min. The films were obtained after drying (30 degrees C/24 h) and conditioning at 25 degrees C and 58% of relative humidity for 7 days, and were then characterized. The results for the color parameters, mechanical properties, phase transitions and infrared spectra showed that some chemical modifications occurred, principally for the gelatin. However, in general, all the characteristics of the films were either typical of films based on blends of these macromolecules without cross linking, or slightly higher. A greater improvement in the properties of this material was probably not observed due to the crystallinity of the PVA, which has a melting point above 90 degrees C. The presence of microcrystals in the polymer chain probably reduced macromolecular mobility, hindering the reaction. Thus more research is necessary to produce biodegradable films with improved properties. (C) 2011 Elsevier Ltd. All rights reserved.

Further probing the nature of FSR 1767

With Two-Micron All-Sky Survey (2MASS) photometry and proper motions, Bonatto et al. suggested that FSR 1767 is a globular cluster (GC), while with J and K NTT/SOFI photometry Froebrich, Meusinger & Scholz concluded that it is not a star cluster. In this study, we combine previous and new evidence that are consistent with a GC. For instance, we show that the horizontal branch (HB) and red giant branch (RGB) stars, besides sharing a common proper motion, have radial density profiles that consistently follow the King`s law independently. Reddening maps around FSR 1767 are built using the bulge RGB as reference and also Schlegel`s extinction values to study local absorptions. Both approaches provide similar maps and show that FSR 1767 is not located in a dust window, which otherwise might have produced the stellar overdensity. Besides, neighbouring regions of similar reddening as FSR 1767 do not present the blue HB stars that are a conspicuous feature in the colour-magnitude diagram of FSR 1767. We report the presence of a compact group of stars located in the central parts of FSR 1767. It appears to be a detached post-collapse core, similar to those of other nearby low-luminosity GCs projected towards the bulge. We note that while the NTT/SOFI photometry of the star cluster FSR 1716 matches perfectly that from 2MASS, it shows a considerable offset for FSR 1767. We discuss the possible reasons why both photometries differ. We confirm our previous structural and photometric fundamental parameters for FSR 1767, which are consistent with a GC.

An analysis of the composite stellar population in M32

We obtained long-slit spectra of high signal-to-noise ratio of the galaxy M32 with the Gemini Multi-Object Spectrograph at the Gemini-North telescope. We analysed the integrated spectra by means of full spectral fitting in order to extract the mixture of stellar populations that best represents its composite nature. Three different galactic radii were analysed, from the nuclear region out to 2 arcmin from the centre. This allows us to compare, for the first time, the results of integrated light spectroscopy with those of resolved colour-magnitude diagrams from the literature. As a main result we propose that an ancient and an intermediate-age population co-exist in M32, and that the balance between these two populations change between the nucleus and outside one effective radius (1r(eff)) in the sense that the contribution from the intermediate population is larger at the nuclear region. We retrieve a smaller signal of a young population at all radii whose origin is unclear and may be a contamination from horizontal branch stars, such as the ones identified by Brown et al. in the nuclear region. We compare our metallicity distribution function for a region 1 to 2 arcmin from the centre to the one obtained with photometric data by Grillmair et al. Both distributions are broad, but our spectroscopically derived distribution has a significant component with [Z/Z(circle dot)] <= -1, which is not found by Grillmair et al.

A new species of Lauromacromia (Odonata: Corduliidae) from Southeastern Brazil, with a cladistic analysis of the genus and comments on Neotropical dragonfly biogeography

Lauromacromia melanica sp. nov. from Conceicao da Barra municipality, Espirito Santo State, Brazil, is described and illustrated based on two males (both in MNRJ n degrees 135). The new species is similar to L. picinguaba differing from it mainly by the absence of pale spots on S3-6 and by the ellipsoid shape of metepisternal pale stripe. A key for males of all species of the genus is provided. A cladistic analysis encompassing 43 external morphological male characters carried out in two distinct procedures, the first with all characters unordered and the second with two or three state characters ordered. The unordered analysis generated only one most-parsimonious tree (66 steps of length, CI = 0.69, RI = 0.62). The hypothesis of monophyly of Lauromacromia is supported and includes three groups, one formed by the Atlantic Forest species (L. melanica sp. nov. + L. picinguaba), and another by the Cerrado species (L. flaviae + (L. bedei + L. luismoojeni)), and L. dubitalis, positioned in polytomy with these two groups. The ordered analysis also generated only one most-parsimonious tree (68 steps of length, CI = 0.70, RI = 0.67), which maintained the monophyly of Lauromacromia but L. dubitalis positioned basally as sister-group to the Atlantic Forest + Cerrado species groups. The geographic distribution of Lauromacromia is updated with a new record of L. luismoojeni based on one adult male (Brazil: Mato Grosso do Sul State) and probable first Brazilian records for L. dubitalis (Amazonas and Para States) based on two larvae. A vicariance hypothesis is proposed to explain spatial evolution of Lauromacromia, and based on current biogeographical classifications we consider Gomphomacromia and Rialla apart from Neotropical biota. Some aspects of biology and ecology of Lauromacromia are also discussed.

Effects of low-level laser therapy and eccentric exercises in the treatment of recreational athletes with chronic achilles tendinopathy

Background: Eccentric exercises (EEs) are recommended for the treatment of Achilles tendinopathy, but the clinical effect from EE has a slow onset. Hypothesis: The addition of low-level laser therapy (LLLT) to EE may cause more rapid clinical improvement. Study Design: Randomized controlled trial; Level of evidence, 1. Methods: A total of 52 recreational athletes with chronic Achilles tendinopathy symptoms were randomized to groups receiving either EE + LLLT or EE + placebo LLLT over 8 weeks in a blinded manner. Low-level laser therapy (lambda = 820 nm) was administered in 12 sessions by irradiating 6 points along the Achilles tendon with a power density of 60 mW/cm(2) and a total dose of 5.4 J per session. Results: The results of the intention-to-treat analysis for the primary outcome, pain intensity during physical activity on the 100-mm visual analog scale, were significantly lower in the LLLT group than in the placebo LLLT group, with 53.6 mm versus 71.5 mm (P = .0003) at 4 weeks, 37.3 mm versus 62.8 mm (P = .0002) at 8 weeks, and 33.0 mm versus 53.0 mm (P =.007) at 12 weeks after randomization. Secondary outcomes of morning stiffness, active dorsiflexion, palpation tenderness, and crepitation showed the same pattern in favor of the LLLT group. Conclusion: Low-level laser therapy, with the parameters used in this study, accelerates clinical recovery from chronic Achilles tendinopathy when added to an EE regimen. For the LLLT group, the results at 4 weeks were similar to the placebo LLLT group results after 12 weeks.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Resilience of protein-protein interaction networks as determined by their large-scale topological features

The relationship between the structure and function of biological networks constitutes a fundamental issue in systems biology. Particularly, the structure of protein-protein interaction networks is related to important biological functions. In this work, we investigated how such a resilience is determined by the large scale features of the respective networks. Four species are taken into account, namely yeast Saccharomyces cerevisiae, worm Caenorhabditis elegans, fly Drosophila melanogaster and Homo sapiens. We adopted two entropy-related measurements (degree entropy and dynamic entropy) in order to quantify the overall degree of robustness of these networks. We verified that while they exhibit similar structural variations under random node removal, they differ significantly when subjected to intentional attacks (hub removal). As a matter of fact, more complex species tended to exhibit more robust networks. More specifically, we quantified how six important measurements of the networks topology (namely clustering coefficient, average degree of neighbors, average shortest path length, diameter, assortativity coefficient, and slope of the power law degree distribution) correlated with the two entropy measurements. Our results revealed that the fraction of hubs and the average neighbor degree contribute significantly for the resilience of networks. In addition, the topological analysis of the removed hubs indicated that the presence of alternative paths between the proteins connected to hubs tend to reinforce resilience. The performed analysis helps to understand how resilience is underlain in networks and can be applied to the development of protein network models.

A multivariate ultrastructural errors-in-variables model with equation error

This paper deals with asymptotic results on a multivariate ultrastructural errors-in-variables regression model with equation errors Sufficient conditions for attaining consistent estimators for model parameters are presented Asymptotic distributions for the line regression estimators are derived Applications to the elliptical class of distributions with two error assumptions are presented The model generalizes previous results aimed at univariate scenarios (C) 2010 Elsevier Inc All rights reserved

Inflated beta distributions

This paper considers the issue of modeling fractional data observed on [0,1), (0,1] or [0,1]. Mixed continuous-discrete distributions are proposed. The beta distribution is used to describe the continuous component of the model since its density can have quite different shapes depending on the values of the two parameters that index the distribution. Properties of the proposed distributions are examined. Also, estimation based on maximum likelihood and conditional moments is discussed. Finally, practical applications that employ real data are presented.

Testing hypotheses in the Birnbaum-Saunders distribution under type-II censored samples

The two-parameter Birnbaum-Saunders distribution has been used successfully to model fatigue failure times. Although censoring is typical in reliability and survival studies, little work has been published on the analysis of censored data for this distribution. In this paper, we address the issue of performing testing inference on the two parameters of the Birnbaum-Saunders distribution under type-II right censored samples. The likelihood ratio statistic and a recently proposed statistic, the gradient statistic, provide a convenient framework for statistical inference in such a case, since they do not require to obtain, estimate or invert an information matrix, which is an advantage in problems involving censored data. An extensive Monte Carlo simulation study is carried out in order to investigate and compare the finite sample performance of the likelihood ratio and the gradient tests. Our numerical results show evidence that the gradient test should be preferred. Further, we also consider the generalized Birnbaum-Saunders distribution under type-II right censored samples and present some Monte Carlo simulations for testing the parameters in this class of models using the likelihood ratio and gradient tests. Three empirical applications are presented. (C) 2011 Elsevier B.V. All rights reserved.

A DISCRETIZATION SCHEME FOR AN ONE-DIMENSIONAL REACTION-DIFFUSION EQUATION WITH DELAY AND ITS DYNAMICS

The goal of this paper is to present an approximation scheme for a reaction-diffusion equation with finite delay, which has been used as a model to study the evolution of a population with density distribution u, in such a way that the resulting finite dimensional ordinary differential system contains the same asymptotic dynamics as the reaction-diffusion equation.

Bias-corrected estimators for dispersion models with dispersion covariates

In this paper we discuss bias-corrected estimators for the regression and the dispersion parameters in an extended class of dispersion models (Jorgensen, 1997b). This class extends the regular dispersion models by letting the dispersion parameter vary throughout the observations, and contains the dispersion models as particular case. General formulae for the O(n(-1)) bias are obtained explicitly in dispersion models with dispersion covariates, which generalize previous results obtained by Botter and Cordeiro (1998), Cordeiro and McCullagh (1991), Cordeiro and Vasconcellos (1999), and Paula (1992). The practical use of the formulae is that we can derive closed-form expressions for the O(n(-1)) biases of the maximum likelihood estimators of the regression and dispersion parameters when the information matrix has a closed-form. Various expressions for the O(n(-1)) biases are given for special models. The formulae have advantages for numerical purposes because they require only a supplementary weighted linear regression. We also compare these bias-corrected estimators with two different estimators which are also bias-free to order O(n(-1)) that are based on bootstrap methods. These estimators are compared by simulation. (C) 2011 Elsevier B.V. All rights reserved.