Variance reduction methods for simulation of densities on Wiener space

We develop a general error analysis framework for the Monte Carlo simulationof densities for functionals in Wiener space. We also study variancereduction methods with the help of Malliavin derivatives. For this, wegive some general heuristic principles which are applied to diffusionprocesses. A comparison with kernel density estimates is made.

Regional differences in socio-economic health inequalities in Spain

This paper reports an analysis of income related health inequalities at the AutonomousCommunity level in Spain using the self assessed health measure in the 2001 edition of theEncuesta Nacional de Salud. We use recently developed methods in order to cardinalise andmodel self assessed health within a regression framework, decompose the sources ofinequality and explain the observed differences across regions. We find that the regions with the highest levels of mean health tend to enjoy the lowest degrees of income related health inequality and vice-versa. The main feature characterizing regions where income related health inequality is low is the absence of a positive gradient between income and health. In turn, the regions where income related health inequality is greater are characterized by a strong and significant positive gradient between health and income. These results suggest that policies aimed at eliminating the gradient between health and income can potentially lead to greate r reductions in socio-economic health inequalities than policies aimed at redistributing income.

Simple and fast methods based on solid-phase extraction coupled to liquid chromatography with UV detection for the monitoring of caffeine in natural and waste waters as marker of anthropogenic impact

Two concentration methods for fast and routine determination of caffeine (using HPLC-UV detection) in surface, and wastewater are evaluated. Both methods are based on solid-phase extraction (SPE) concentration with octadecyl silica sorbents. A common “offline” SPE procedure shows that quantitative recovery of caffeine is obtained with 2 mL of an elution mixture solvent methanol-water containing at least 60% methanol. The method detection limit is 0.1 μg L−1 when percolating 1 L samples through the cartridge. The development of an “online” SPE method based on a mini-SPE column, containing 100 mg of the same sorbent, directly connected to the HPLC system allows the method detection limit to be decreased to 10 ng L−1 with a sample volume of 100 mL. The “offline” SPE method is applied to the analysis of caffeine in wastewater samples, whereas the “on-line” method is used for analysis in natural waters from streams receiving significant water intakes from local wastewater treatment plants

A limited area model intercomparison on the "Montserrat-2000" flash-flood event using statistical and deterministic methods

In the scope of the European project Hydroptimet, INTERREG IIIB-MEDOCC programme, limited area model (LAM) intercomparison of intense events that produced many damages to people and territory is performed. As the comparison is limited to single case studies, the work is not meant to provide a measure of the different models' skill, but to identify the key model factors useful to give a good forecast on such a kind of meteorological phenomena. This work focuses on the Spanish flash-flood event, also known as "Montserrat-2000" event. The study is performed using forecast data from seven operational LAMs, placed at partners' disposal via the Hydroptimet ftp site, and observed data from Catalonia rain gauge network. To improve the event analysis, satellite rainfall estimates have been also considered. For statistical evaluation of quantitative precipitation forecasts (QPFs), several non-parametric skill scores based on contingency tables have been used. Furthermore, for each model run it has been possible to identify Catalonia regions affected by misses and false alarms using contingency table elements. Moreover, the standard "eyeball" analysis of forecast and observed precipitation fields has been supported by the use of a state-of-the-art diagnostic method, the contiguous rain area (CRA) analysis. This method allows to quantify the spatial shift forecast error and to identify the error sources that affected each model forecasts. High-resolution modelling and domain size seem to have a key role for providing a skillful forecast. Further work is needed to support this statement, including verification using a wider observational data set.

Statistical analysis and modelling of weather radar beam propagation conditions in the Po Valley (Italy)

Ground clutter caused by anomalous propagation (anaprop) can affect seriously radar rain rate estimates, particularly in fully automatic radar processing systems, and, if not filtered, can produce frequent false alarms. A statistical study of anomalous propagation detected from two operational C-band radars in the northern Italian region of Emilia Romagna is discussed, paying particular attention to its diurnal and seasonal variability. The analysis shows a high incidence of anaprop in summer, mainly in the morning and evening, due to the humid and hot summer climate of the Po Valley, particularly in the coastal zone. Thereafter, a comparison between different techniques and datasets to retrieve the vertical profile of the refractive index gradient in the boundary layer is also presented. In particular, their capability to detect anomalous propagation conditions is compared. Furthermore, beam path trajectories are simulated using a multilayer ray-tracing model and the influence of the propagation conditions on the beam trajectory and shape is examined. High resolution radiosounding data are identified as the best available dataset to reproduce accurately the local propagation conditions, while lower resolution standard TEMP data suffers from interpolation degradation and Numerical Weather Prediction model data (Lokal Model) are able to retrieve a tendency to superrefraction but not to detect ducting conditions. Observing the ray tracing of the centre, lower and upper limits of the radar antenna 3-dB half-power main beam lobe it is concluded that ducting layers produce a change in the measured volume and in the power distribution that can lead to an additional error in the reflectivity estimate and, subsequently, in the estimated rainfall rate.

Introduction of sensor spectral response into image fusion methods. Application to wavelet-based methods

Usual image fusion methods inject features from a high spatial resolution panchromatic sensor into every low spatial resolution multispectral band trying to preserve spectral signatures and improve spatial resolution to that of the panchromatic sensor. The objective is to obtain the image that would be observed by a sensor with the same spectral response (i.e., spectral sensitivity and quantum efficiency) as the multispectral sensors and the spatial resolution of the panchromatic sensor. But in these methods, features from electromagnetic spectrum regions not covered by multispectral sensors are injected into them, and physical spectral responses of the sensors are not considered during this process. This produces some undesirable effects, such as resolution overinjection images and slightly modified spectral signatures in some features. The authors present a technique which takes into account the physical electromagnetic spectrum responses of sensors during the fusion process, which produces images closer to the image obtained by the ideal sensor than those obtained by usual wavelet-based image fusion methods. This technique is used to define a new wavelet-based image fusion method.

Binding energy and momentum distribution of nuclear matter using Green's functions methods

The influence of hole-hole (h-h) propagation in addition to the conventional particle-particle (p-p) propagation, on the energy per particle and the momentum distribution is investigated for the v2 central interaction which is derived from Reid¿s soft-core potential. The results are compared to Brueckner-Hartree-Fock calculations with a continuous choice for the single-particle (SP) spectrum. Calculation of the energy from a self-consistently determined SP spectrum leads to a lower saturation density. This result is not corroborated by calculating the energy from the hole spectral function, which is, however, not self-consistent. A generalization of previous calculations of the momentum distribution, based on a Goldstone diagram expansion, is introduced that allows the inclusion of h-h contributions to all orders. From this result an alternative calculation of the kinetic energy is obtained. In addition, a direct calculation of the potential energy is presented which is obtained from a solution of the ladder equation containing p-p and h-h propagation to all orders. These results can be considered as the contributions of selected Goldstone diagrams (including p-p and h-h terms on the same footing) to the kinetic and potential energy in which the SP energy is given by the quasiparticle energy. The results for the summation of Goldstone diagrams leads to a different momentum distribution than the one obtained from integrating the hole spectral function which in general gives less depletion of the Fermi sea. Various arguments, based partly on the results that are obtained, are put forward that a self-consistent determination of the spectral functions including the p-p and h-h ladder contributions (using a realistic interaction) will shed light on the question of nuclear saturation at a nonrelativistic level that is consistent with the observed depletion of SP orbitals in finite nuclei.

Perturbing Hele-Shaw flow with a small gap gradient

A controlled perturbation is introduced into the Saffman-Taylor flow problem by adding a gradient to the gap of a Hele-Shaw cell. The stability of the single-finger steady state was found to be strongly affected by such a perturbation. Compared with patterns in a standard Hele-Shaw cell, the single Saffman-Taylor finger was stabilized or destabilized according to the sign of the gap gradient. While a linear stability analysis shows that this perturbation should have a negligible effect on the early-stage pattern formation, the experimental data indicate that the characteristic length for the initial breakup of a flat interface has been changed by the perturbation.

Application of fuzzy-rule-based postprocessing to correlation methods in pattern recognition

One of the most important problems in optical pattern recognition by correlation is the appearance of sidelobes in the correlation plane, which causes false alarms. We present a method that eliminate sidelobes of up to a given height if certain conditions are satisfied. The method can be applied to any generalized synthetic discriminant function filter and is capable of rejecting lateral peaks that are even higher than the central correlation. Satisfactory results were obtained in both computer simulations and optical implementation.

DnaSP, DNA polymorphism analyses by the coalescent and other methods

DnaSP is a software package for the analysis of DNA polymorphism data. Present version introduces several new modules and features which, among other options allow: (1) handling big data sets (~5 Mb per sequence); (2) conducting a large number of coalescent-based tests by Monte Carlo computer simulations; (3) extensive analyses of the genetic differentiation and gene flow among populations; (4) analysing the evolutionary pattern of preferred and unpreferred codons; (5) generating graphical outputs for an easy visualization of results. Availability: The software package, including complete documentation and examples, is freely available to academic users from: http://www.ub.es/dnasp

Rapid radiation in spiny lobsters (Palinurus spp) as revealed by classic and ABC methods using mtDNA and microsatellite data

Background: Molecular tools may help to uncover closely related and still diverging species from a wide variety of taxa and provide insight into the mechanisms, pace and geography of marine speciation. There is a certain controversy on the phylogeography and speciation modes of species-groups with an Eastern Atlantic-Western Indian Ocean distribution, with previous studies suggesting that older events (Miocene) and/or more recent (Pleistocene) oceanographic processes could have influenced the phylogeny of marine taxa. The spiny lobster genus Palinurus allows for testing among speciation hypotheses, since it has a particular distribution with two groups of three species each in the Northeastern Atlantic (P. elephas, P. mauritanicus and P. charlestoni) and Southeastern Atlantic and Southwestern Indian Oceans (P. gilchristi, P. delagoae and P. barbarae). In the present study, we obtain a more complete understanding of the phylogenetic relationships among these species through a combined dataset with both nuclear and mitochondrial markers, by testing alternative hypotheses on both the mutation rate and tree topology under the recently developed approximate Bayesian computation (ABC) methods. Results Our analyses support a North-to-South speciation pattern in Palinurus with all the South-African species forming a monophyletic clade nested within the Northern Hemisphere species. Coalescent-based ABC methods allowed us to reject the previously proposed hypothesis of a Middle Miocene speciation event related with the closure of the Tethyan Seaway. Instead, divergence times obtained for Palinurus species using the combined mtDNA-microsatellite dataset and standard mutation rates for mtDNA agree with known glaciation-related processes occurring during the last 2 my. Conclusion The Palinurus speciation pattern is a typical example of a series of rapid speciation events occurring within a group, with very short branches separating different species. Our results support the hypothesis that recent climate change-related oceanographic processes have influenced the phylogeny of marine taxa, with most Palinurus species originating during the last two million years. The present study highlights the value of new coalescent-based statistical methods such as ABC for testing different speciation hypotheses using molecular data.

Assessment of therapeutic benefit of antiviral therapy in chronic hepatitis C: is hepatic venous pressure gradient a better end point?

Chronic hepatitis C is a major healthcare problem. The response to antiviral therapy for patients with chronic hepatitis C has previously been defined biochemically and by PCR. However, changes in the hepatic venous pressure gradient (HVPG) may be considered as an adjunctive end point for the therapeutic evaluation of antiviral therapy in chronic hepatitis C. It is a validated technique which is safe, well tolerated, well established, and reproducible. Serial HVPG measurements may be the best way to evaluate response to therapy in chronic hepatitis C.

CO2 wedged hepatic venography in the evaluation of portal hypertension.

BACKGROUND/AIMS/METHODS During hepatic vein catheterisation, in addition to measurement of hepatic venous pressure gradient (HVPG), iodine wedged retrograde portography can be easily obtained. However, it rarely allows correct visualisation of the portal vein. Recently, CO2 has been suggested to allow better angiographic demonstration of the portal vein than iodine. In this study we investigated the efficacy of CO2 compared with iodinated contrast medium for portal vein imaging and its role in the evaluation of portal hypertension in a series of 100 patients undergoing hepatic vein catheterisation, 71 of whom had liver cirrhosis. RESULTS In the overall series, CO2 venography was markedly superior to iodine, allowing correct visualisation of the different segments of the portal venous system. In addition, CO2, but not iodine, visualised portal-systemic collaterals in 34 patients. In cirrhosis, non-visualisation of the portal vein on CO2 venography occurred in 11 cases; four had portal vein thrombosis and five had communications between different hepatic veins. Among non-cirrhotics, lack of portal vein visualisation had a 90% sensitivity, 88% specificity, 94% negative predictive value, and 83% positive predictive value in the diagnosis of pre-sinusoidal portal hypertension. CONCLUSIONS Visualisation of the venous portal system by CO2 venography is markedly superior to iodine. The use of CO2 wedged portography is a useful and safe complementary procedure during hepatic vein catheterisation which may help to detect portal thrombosis. Also, lack of demonstration of the portal vein in non-cirrhotic patients strongly suggests the presence of pre-sinusoidal portal hypertension.

Nadolol plus isosorbide mononitrate alone or associated with band ligation in the prevention of recurrent bleeding: A multicenter randomized controlled trial.

Background and aims: Previous clinical trials suggest that adding non-selective beta-blockers improves the efficacy of endoscopic band ligation (EBL) in the prevention of recurrent bleeding, but no study has evaluated whether EBL improves the efficacy of beta-blockers + isosorbide-5-mononitrate. The present study was aimed at evaluating this issue in a multicentre randomised controlled trial (RCT) and to correlate changes in hepatic venous pressure gradient (HVPG) during treatment with clinical outcomes. Methods: 158 patients with cirrhosis, admitted because of variceal bleeding, were randomised to receive nadolol+isosorbide-5-mononitrate alone (Drug: n=78) or combined with EBL (Drug+EBL; n=80). HVPG measurements were performed at randomisation and after 4¿6 weeks on medical therapy. Results: Median follow-up was 15 months. One-year probability of recurrent bleeding was similar in both groups (33% vs 26%: p=0.3). There were no significant differences in survival or need of rescue shunts. Overall adverse events or those requiring hospital admission were significantly more frequent in the Drug+EBL group. Recurrent bleeding was significantly more frequent in HVPG non-responders than in responders (HVPG reduction ¿20% or ¿12 mm Hg). Among non-responders recurrent bleeding was similar in patients treated with Drugs or Drugs+EBL. Conclusions: Adding EBL to pharmacological treatment did not reduce recurrent bleeding, the need for rescue therapy, or mortality, and was associated with more adverse events. Furthermore, associating EBL to drug therapy did not reduce the high rebleeding risk of HVPG non-responders.

The extent of the collateral circulation influences the postprandial increase in portal pressure in patients with cirrhosis.

Background: In cirrhosis, repeated flares of portal pressure and collateral blood flow provoked by postprandial hyperaemia may contribute to variceal dilation and rupture. Aim: To examine the effect of the extent of the collateral circulation on the postprandial increase in portal pressure observed in cirrhosis. Patients and methods: The hepatic venous pressure gradient (HVPG), hepatic blood flow and azygos blood flow were measured in 64 patients with cirrhosis before and after a standard liquid meal. Results: Peak increases in HVPG (median+14.9%), hepatic blood flow (median+25.4%), and azygos blood flow (median+32.2%) occurred at 30 min after the meal. Compared with patients with marked postprandial increase in HVPG (above the median, n¿=¿32), those showing mild (<15%, n¿=¿32) increase in HVPG had a higher baseline azygos flow (p<0.01) and underwent a greater postprandial increase in azygos flow (p<0.02). Hepatic blood flow increased similarly in both groups. Postprandial increases in HVPG were inversely correlated (p<0.001) with both baseline azygos flow (r¿=¿¿0.69) and its postprandial increase (r¿=¿¿0.72). Food intake increased nitric oxide products in the azygos (p<0.01), but not in the hepatic vein. Large varices (p<0.01) and previous variceal bleeding (p<0.001) were more frequent in patients with mild increase in HVPG. Conclusions: Postprandial hyperaemia simultaneously increases HVPG and collateral flow. The extent of the collateral circulation determines the HVPG response to food intake. Patients with extensive collateralisation show less pronounced postprandial increases in HVPG, but associated with marked flares in collateral flow. Collateral vessels preserve their ability to dilate in response to increased blood flow.