Comparison of time-lapse GPR data collected under natural and forced infiltration conditions to estimate vadose zone hydraulic parameters

Time-lapse crosshole ground-penetrating radar (GPR) data, collected while infiltration occurs, can provide valuable information regarding the hydraulic properties of the unsaturated zone. In particular, the stochastic inversion of such data provides estimates of parameter uncertainties, which are necessary for hydrological prediction and decision making. Here, we investigate the effect of different infiltration conditions on the stochastic inversion of time-lapse, zero-offset-profile, GPR data. Inversions are performed using a Bayesian Markov-chain-Monte-Carlo methodology. Our results clearly indicate that considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions

Estimation of dynamic latent variable models using simulated nonparametric moments

Abstract. Given a model that can be simulated, conditional moments at a trial parameter value can be calculated with high accuracy by applying kernel smoothing methods to a long simulation. With such conditional moments in hand, standard method of moments techniques can be used to estimate the parameter. Because conditional moments are calculated using kernel smoothing rather than simple averaging, it is not necessary that the model be simulable subject to the conditioning information that is used to define the moment conditions. For this reason, the proposed estimator is applicable to general dynamic latent variable models. It is shown that as the number of simulations diverges, the estimator is consistent and a higher-order expansion reveals the stochastic difference between the infeasible GMM estimator based on the same moment conditions and the simulated version. In particular, we show how to adjust standard errors to account for the simulations. Monte Carlo results show how the estimator may be applied to a range of dynamic latent variable (DLV) models, and that it performs well in comparison to several other estimators that have been proposed for DLV models.

Beurling-Landau densities of weighted Fekete sets and correlation kernel estimates

Vegeu el resum a l'inici del document del fitxer adjunt.

Error estimate and control in the MSFV method for multiphase flow in porous media

n this paper the iterative MSFV method is extended to include the sequential implicit simulation of time dependent problems involving the solution of a system of pressure-saturation equations. To control numerical errors in simulation results, an error estimate, based on the residual of the MSFV approximate pressure field, is introduced. In the initial time steps in simulation iterations are employed until a specified accuracy in pressure is achieved. This initial solution is then used to improve the localization assumption at later time steps. Additional iterations in pressure solution are employed only when the pressure residual becomes larger than a specified threshold value. Efficiency of the strategy and the error control criteria are numerically investigated. This paper also shows that it is possible to derive an a-priori estimate and control based on the allowed pressure-equation residual to guarantee the desired accuracy in saturation calculation.

2003 Revised Estimates for Public Services - Estimate for the Department of Health and Children

2003 Revised Estimates for Public Services – Estimate for the Department of Health and Children Estimate of the amount required in the year ending 31 December, 2003 for the salaries and expenses of the Office of the Minister for Health and Children (including Oifig an Ard- Chláraitheora), and certain other services administered by that Office, including grants to Health Boards and miscellaneous grants. Click here to download PDF 52kb

2003 Revised Estimates for Public Services Estimate for the Department of Health and Children

Click here to download PDF

2006 Drugs estimate refuels debate on governance, Issue 17, Spring 2006

In November 2005 €34.027 million was voted for the Drugs Initiative, which funds the local and regional drugs task forces, and the Young People’s Facilities and Services Fund (YPFSF), in 2006 – an 8 per cent increase on the allocation for 2005 and ‘well in excess of inflation’ according to the government. The community sector strongly criticised this estimate, calling for an additional €8 million to €15 million, to fund the projects identified following the creation of the Emerging Needs Fund in January 2005. In February 2006, the government revised its drugs estimate upwards by a further €8.979 million. As well as criticising the 2006 drugs estimate, the community sector raised a series of concerns with regard to governance – the rules, processes and behaviour that affect the way in which powers are exercised4 – or, in short, decision-making with regard to drugs policy in Ireland.This resource was contributed by The National Documentation Centre on Drug Use.

Alternative ways to estimate change points in multinomial sequences. An application to an authorship attribution problem

The statistical analysis of literary style is the part of stylometry that compares measurable characteristicsin a text that are rarely controlled by the author, with those in other texts. When thegoal is to settle authorship questions, these characteristics should relate to the author’s style andnot to the genre, epoch or editor, and they should be such that their variation between authors islarger than the variation within comparable texts from the same author.For an overview of the literature on stylometry and some of the techniques involved, see for exampleMosteller and Wallace (1964, 82), Herdan (1964), Morton (1978), Holmes (1985), Oakes (1998) orLebart, Salem and Berry (1998).Tirant lo Blanc, a chivalry book, is the main work in catalan literature and it was hailed to be“the best book of its kind in the world” by Cervantes in Don Quixote. Considered by writterslike Vargas Llosa or Damaso Alonso to be the first modern novel in Europe, it has been translatedseveral times into Spanish, Italian and French, with modern English translations by Rosenthal(1996) and La Fontaine (1993). The main body of this book was written between 1460 and 1465,but it was not printed until 1490.There is an intense and long lasting debate around its authorship sprouting from its first edition,where its introduction states that the whole book is the work of Martorell (1413?-1468), while atthe end it is stated that the last one fourth of the book is by Galba (?-1490), after the death ofMartorell. Some of the authors that support the theory of single authorship are Riquer (1990),Chiner (1993) and Badia (1993), while some of those supporting the double authorship are Riquer(1947), Coromines (1956) and Ferrando (1995). For an overview of this debate, see Riquer (1990).Neither of the two candidate authors left any text comparable to the one under study, and thereforediscriminant analysis can not be used to help classify chapters by author. By using sample textsencompassing about ten percent of the book, and looking at word length and at the use of 44conjunctions, prepositions and articles, Ginebra and Cabos (1998) detect heterogeneities that mightindicate the existence of two authors. By analyzing the diversity of the vocabulary, Riba andGinebra (2000) estimates that stylistic boundary to be near chapter 383.Following the lead of the extensive literature, this paper looks into word length, the use of the mostfrequent words and into the use of vowels in each chapter of the book. Given that the featuresselected are categorical, that leads to three contingency tables of ordered rows and therefore tothree sequences of multinomial observations.Section 2 explores these sequences graphically, observing a clear shift in their distribution. Section 3describes the problem of the estimation of a suden change-point in those sequences, in the followingsections we propose various ways to estimate change-points in multinomial sequences; the methodin section 4 involves fitting models for polytomous data, the one in Section 5 fits gamma modelsonto the sequence of Chi-square distances between each row profiles and the average profile, theone in Section 6 fits models onto the sequence of values taken by the first component of thecorrespondence analysis as well as onto sequences of other summary measures like the averageword length. In Section 7 we fit models onto the marginal binomial sequences to identify thefeatures that distinguish the chapters before and after that boundary. Most methods rely heavilyon the use of generalized linear models

Official Positions for FRAX(®) clinical regarding glucocorticoids: the impact of the use of glucocorticoids on the estimate by FRAX(®) of the 10 year risk of fracture from Joint Official Positions Development Conference of the International Society for Clinical Densitometry and International Osteoporosis Foundation on FRAX(®).

Given the significant impact the use of glucocorticoids can have on fracture risk independent of bone density, their use has been incorporated as one of the clinical risk factors for calculating the 10-year fracture risk in the World Health Organization's Fracture Risk Assessment Tool (FRAX(®)). Like the other clinical risk factors, the use of glucocorticoids is included as a dichotomous variable with use of steroids defined as past or present exposure of 3 months or more of use of a daily dose of 5 mg or more of prednisolone or equivalent. The purpose of this report is to give clinicians guidance on adjustments which should be made to the 10-year risk based on the dose, duration of use and mode of delivery of glucocorticoids preparations. A subcommittee of the International Society for Clinical Densitometry and International Osteoporosis Foundation joint Position Development Conference presented its findings to an expert panel and the following recommendations were selected. 1) There is a dose relationship between glucocorticoid use of greater than 3 months and fracture risk. The average dose exposure captured within FRAX(®) is likely to be a prednisone dose of 2.5-7.5 mg/day or its equivalent. Fracture probability is under-estimated when prednisone dose is greater than 7.5 mg/day and is over-estimated when the prednisone dose is less than 2.5 mg/day. 2) Frequent intermittent use of higher doses of glucocorticoids increases fracture risk. Because of the variability in dose and dosing schedule, quantification of this risk is not possible. 3) High dose inhaled glucocorticoids may be a risk factor for fracture. FRAX(®) may underestimate fracture probability in users of high dose inhaled glucocorticoids. 4) Appropriate glucocorticoid replacement in individuals with adrenal insufficiency has not been found to increase fracture risk. In such patients, use of glucocorticoids should not be included in FRAX(®) calculations.

Plasma and urine profiles of Delta9-tetrahydrocannabinol and its metabolites 11-hydroxy-Delta9-tetrahydrocannabinol and 11-nor-9-carboxy-Delta9-tetrahydrocannabinol after cannabis smoking by male volunteers to estimate recent consumption by athletes.

Since 2004, cannabis has been prohibited by the World Anti-Doping Agency for all sports competitions. In the years since then, about half of all positive doping cases in Switzerland have been related to cannabis consumption. In doping urine analysis, the target analyte is 11-nor-9-carboxy-Delta(9)-tetrahydrocannabinol (THC-COOH), the cutoff being 15 ng/mL. However, the wide urinary detection window of the long-term metabolite of Delta(9)-tetrahydrocannabinol (THC) does not allow a conclusion to be drawn regarding the time of consumption or the impact on the physical performance. The purpose of the present study on light cannabis smokers was to evaluate target analytes with shorter urinary excretion times. Twelve male volunteers smoked a cannabis cigarette standardized to 70 mg THC per cigarette. Plasma and urine were collected up to 8 h and 11 days, respectively. Total THC, 11-hydroxy-Delta(9)-tetrahydrocannabinol (THC-OH), and THC-COOH were determined after hydrolysis followed by solid-phase extraction and gas chromatography/mass spectrometry. The limits of quantitation were 0.1-1.0 ng/mL. Eight puffs delivered a mean THC dose of 45 mg. Plasma levels of total THC, THC-OH, and THC-COOH were measured in the ranges 0.2-59.1, 0.1-3.9, and 0.4-16.4 ng/mL, respectively. Peak concentrations were observed at 5, 5-20, and 20-180 min. Urine levels were measured in the ranges 0.1-1.3, 0.1-14.4, and 0.5-38.2 ng/mL, peaking at 2, 2, and 6-24 h, respectively. The times of the last detectable levels were 2-8, 6-96, and 48-120 h. Besides high to very high THC-COOH levels (245 +/- 1,111 ng/mL), THC (3 +/- 8 ng/mL) and THC-OH (51 +/- 246 ng/mL) were found in 65 and 98% of cannabis-positive athletes' urine samples, respectively. In conclusion, in addition to THC-COOH, the pharmacologically active THC and THC-OH should be used as target analytes for doping urine analysis. In the case of light cannabis use, this may allow the estimation of more recent consumption, probably influencing performance during competitions. However, it is not possible to discriminate the intention of cannabis use, i.e., for recreational or doping purposes. Additionally, pharmacokinetic data of female volunteers are needed to interpret cannabis-positive doping cases of female athletes.

Nonparametric estimation of Value-at-Risk

A method to estimate an extreme quantile that requires no distributional assumptions is presented. The approach is based on transformed kernel estimation of the cumulative distribution function (cdf). The proposed method consists of a double transformation kernel estimation. We derive optimal bandwidth selection methods that have a direct expression for the smoothing parameter. The bandwidth can accommodate to the given quantile level. The procedure is useful for large data sets and improves quantile estimation compared to other methods in heavy tailed distributions. Implementation is straightforward and R programs are available.

Use of variability in national and regional data to estimate the prevalence of lymphangioleiomyomatosis.

Background: Understanding the true prevalence of lymphangioleiomyomatosis (LAM) is important in estimating disease burden and targeting specific interventions. As with all rare diseases, obtaining reliable epidemiological data is difficult and requires innovative approaches.Aim: To determine the prevalence and incidence of LAM using data from patient organizations in seven countries, and to use the extent to which the prevalence of LAM varies regionally and nationally to determine whether prevalence estimates are related to health-care provision.Methods: Numbers of women with LAM were obtained from patient groups and national databases from seven countries (n = 1001). Prevalence was calculated for regions within countries using female population figures from census data. Incidence estimates were calculated for the USA, UK and Switzerland. Regional variation in prevalence and changes in incidence over time were analysed using Poisson regression and linear regression.Results: Prevalence of LAM in the seven countries ranged from 3.4 to 7.8/million women with significant variation, both between countries and between states in the USA. This variation did not relate to the number of pulmonary specialists in the region nor the percentage of population with health insurance, but suggests a large number of patients remain undiagnosed. The incidence of LAM from 2004 to 2008 ranged from 0.23 to 0.31/million women/per year in the USA, UK and Switzerland.Conclusions: Using this method, we have found that the prevalence of LAM is higher than that previously recorded and that many patients with LAM are undiagnosed.

A comparison of the alr and ilr transformations for kernel density estimation of compositional data

In a seminal paper, Aitchison and Lauder (1985) introduced classical kernel densityestimation techniques in the context of compositional data analysis. Indeed, they gavetwo options for the choice of the kernel to be used in the kernel estimator. One ofthese kernels is based on the use the alr transformation on the simplex SD jointly withthe normal distribution on RD-1. However, these authors themselves recognized thatthis method has some deficiencies. A method for overcoming these dificulties based onrecent developments for compositional data analysis and multivariate kernel estimationtheory, combining the ilr transformation with the use of the normal density with a fullbandwidth matrix, was recently proposed in Martín-Fernández, Chacón and Mateu-Figueras (2006). Here we present an extensive simulation study that compares bothmethods in practice, thus exploring the finite-sample behaviour of both estimators

Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten-Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.