The Effect of Four-Lane to Three-Lane Conversion on the Number of Crashes and Crash Rates in Iowa Roads, June 2005

We analyze crash data collected by the Iowa Department of Transportation using Bayesian methods. The data set includes monthly crash numbers, estimated monthly traffic volumes, site length and other information collected at 30 paired sites in Iowa over more than 20 years during which an intervention experiment was set up. The intervention consisted in transforming 15 undivided road segments from four-lane to three lanes, while an additional 15 segments, thought to be comparable in terms of traffic safety-related characteristics were not converted. The main objective of this work is to find out whether the intervention reduces the number of crashes and the crash rates at the treated sites. We fitted a hierarchical Poisson regression model with a change-point to the number of monthly crashes per mile at each of the sites. Explanatory variables in the model included estimated monthly traffic volume, time, an indicator for intervention reflecting whether the site was a “treatment” or a “control” site, and various interactions. We accounted for seasonal effects in the number of crashes at a site by including smooth trigonometric functions with three different periods to reflect the four seasons of the year. A change-point at the month and year in which the intervention was completed for treated sites was also included. The number of crashes at a site can be thought to follow a Poisson distribution. To estimate the association between crashes and the explanatory variables, we used a log link function and added a random effect to account for overdispersion and for autocorrelation among observations obtained at the same site. We used proper but non-informative priors for all parameters in the model, and carried out all calculations using Markov chain Monte Carlo methods implemented in WinBUGS. We evaluated the effect of the four to three-lane conversion by comparing the expected number of crashes per year per mile during the years preceding the conversion and following the conversion for treatment and control sites. We estimated this difference using the observed traffic volumes at each site and also on a per 100,000,000 vehicles. We also conducted a prospective analysis to forecast the expected number of crashes per mile at each site in the study one year, three years and five years following the four to three-lane conversion. Posterior predictive distributions of the number of crashes, the crash rate and the percent reduction in crashes per mile were obtained for each site for the months of January and June one, three and five years after completion of the intervention. The model appears to fit the data well. We found that in most sites, the intervention was effective and reduced the number of crashes. Overall, and for the observed traffic volumes, the reduction in the expected number of crashes per year and mile at converted sites was 32.3% (31.4% to 33.5% with 95% probability) while at the control sites, the reduction was estimated to be 7.1% (5.7% to 8.2% with 95% probability). When the reduction in the expected number of crashes per year, mile and 100,000,000 AADT was computed, the estimates were 44.3% (43.9% to 44.6%) and 25.5% (24.6% to 26.0%) for converted and control sites, respectively. In both cases, the difference in the percent reduction in the expected number of crashes during the years following the conversion was significantly larger at converted sites than at control sites, even though the number of crashes appears to decline over time at all sites. Results indicate that the reduction in the expected number of sites per mile has a steeper negative slope at converted than at control sites. Consistent with this, the forecasted reduction in the number of crashes per year and mile during the years after completion of the conversion at converted sites is more pronounced than at control sites. Seasonal effects on the number of crashes have been well-documented. In this dataset, we found that, as expected, the expected number of monthly crashes per mile tends to be higher during winter months than during the rest of the year. Perhaps more interestingly, we found that there is an interaction between the four to three-lane conversion and season; the reduction in the number of crashes appears to be more pronounced during months, when the weather is nice than during other times of the year, even though a reduction was estimated for the entire year. Thus, it appears that the four to three-lane conversion, while effective year-round, is particularly effective in reducing the expected number of crashes in nice weather.

Intravenous streptomycin dosing regimen in a patient undergoing hemodialysis: Plasma level monitoring and pharmacokinetic simulation

Introduction: Streptomycin, as other aminoglycosides, exhibits concentration-dependent bacterial killing but has a narrow therapeutic window. It is primarily eliminated unchanged by the kidneys. Data and dosing information to achieve a safe regimen in patients with chronic renal failure undergoing hemodialysis (HD) are scarce. Although main adverse reactions are related to prolonged, elevated serum concentrations, literature recommendation is to administer streptomycin after each HD. Patients (or Materials) and Methods: We report the case of a patient with end-stage renal failure, undergoing HD, who was successfully treated with streptomycin for gentamicin-resistant Enterococcus faecalis bacteremia with prosthetic arteriovenous fistula infection. Streptomycin was administered intravenously 7.5 mg/kg, 3 hours before each dialysis (3 times a week) during 6 weeks in combination with amoxicillin. Streptomycin plasma levels were monitored with repeated blood sampling before, after, and between HD sessions. A 2-compartment model was used to reconstruct the concentration time profile over days on and off HD. Results: Streptomycin trough plasma-concentration was 2.8 mg/L. It peaked to 21.4 mg/L 30 minutes after intravenous administration, decreased to 18.2 mg/L immediately before HD, and dropped to 4.5 mg/L at the end of a 4-hour HD session. Plasma level increased again to 5.7 mg/L 2 hours after the end of HD and was 2.8 mg/L 48 hours later, before the next administration and HD. The pharmacokinetics of streptomycin was best described with a 2-compartment model. The computer simulation fitted fairly well to the observed concentrations during or between HD sessions. Redistribution between the 2 compartments after the end of HD reproduced the rebound of plasma concentrations after HD. No significant toxicity was observed during treatment. The outcome of the infection was favorable, and no sign of relapse was observed after a follow-up of 3 months. Conclusion: Streptomycin administration of 7.5 mg/kg 3 hours before HD sessions in a patient with end-stage renal failure resulted in an effective and safe dosing regimen. Monitoring plasma levels along with pharmacokinetic simulation document the suitability of this dosing scheme, which should replace current dosage recommendations for streptomycin in HD.

Resurgence of HIV infection among men who have sex with men in Switzerland: mathematical modelling study.

BACKGROUND: New HIV infections in men who have sex with men (MSM) have increased in Switzerland since 2000 despite combination antiretroviral therapy (cART). The objectives of this mathematical modelling study were: to describe the dynamics of the HIV epidemic in MSM in Switzerland using national data; to explore the effects of hypothetical prevention scenarios; and to conduct a multivariate sensitivity analysis. METHODOLOGY/PRINCIPAL FINDINGS: The model describes HIV transmission, progression and the effects of cART using differential equations. The model was fitted to Swiss HIV and AIDS surveillance data and twelve unknown parameters were estimated. Predicted numbers of diagnosed HIV infections and AIDS cases fitted the observed data well. By the end of 2010, an estimated 13.5% (95% CI 12.5, 14.6%) of all HIV-infected MSM were undiagnosed and accounted for 81.8% (95% CI 81.1, 82.4%) of new HIV infections. The transmission rate was at its lowest from 1995-1999, with a nadir of 46 incident HIV infections in 1999, but increased from 2000. The estimated number of new infections continued to increase to more than 250 in 2010, although the reproduction number was still below the epidemic threshold. Prevention scenarios included temporary reductions in risk behaviour, annual test and treat, and reduction in risk behaviour to levels observed earlier in the epidemic. These led to predicted reductions in new infections from 2 to 26% by 2020. Parameters related to disease progression and relative infectiousness at different HIV stages had the greatest influence on estimates of the net transmission rate. CONCLUSIONS/SIGNIFICANCE: The model outputs suggest that the increase in HIV transmission amongst MSM in Switzerland is the result of continuing risky sexual behaviour, particularly by those unaware of their infection status. Long term reductions in the incidence of HIV infection in MSM in Switzerland will require increased and sustained uptake of effective interventions.

Improving generalized regression analysis for the spatial prediction of forest communities

Aim This study used data from temperate forest communities to assess: (1) five different stepwise selection methods with generalized additive models, (2) the effect of weighting absences to ensure a prevalence of 0.5, (3) the effect of limiting absences beyond the environmental envelope defined by presences, (4) four different methods for incorporating spatial autocorrelation, and (5) the effect of integrating an interaction factor defined by a regression tree on the residuals of an initial environmental model. Location State of Vaud, western Switzerland. Methods Generalized additive models (GAMs) were fitted using the grasp package (generalized regression analysis and spatial predictions, http://www.cscf.ch/grasp). Results Model selection based on cross-validation appeared to be the best compromise between model stability and performance (parsimony) among the five methods tested. Weighting absences returned models that perform better than models fitted with the original sample prevalence. This appeared to be mainly due to the impact of very low prevalence values on evaluation statistics. Removing zeroes beyond the range of presences on main environmental gradients changed the set of selected predictors, and potentially their response curve shape. Moreover, removing zeroes slightly improved model performance and stability when compared with the baseline model on the same data set. Incorporating a spatial trend predictor improved model performance and stability significantly. Even better models were obtained when including local spatial autocorrelation. A novel approach to include interactions proved to be an efficient way to account for interactions between all predictors at once. Main conclusions Models and spatial predictions of 18 forest communities were significantly improved by using either: (1) cross-validation as a model selection method, (2) weighted absences, (3) limited absences, (4) predictors accounting for spatial autocorrelation, or (5) a factor variable accounting for interactions between all predictors. The final choice of model strategy should depend on the nature of the available data and the specific study aims. Statistical evaluation is useful in searching for the best modelling practice. However, one should not neglect to consider the shapes and interpretability of response curves, as well as the resulting spatial predictions in the final assessment.

Increasing prevalence of gastroschisis in Europe 1980-2002: a phenomenon restricted to younger mothers?

Gastroschisis is an abdominal wall defect more prevalent in offspring of young mothers. It is known to be increasing in prevalence despite the general decrease in the proportion of births to young European women. We investigated whether the increase in prevalence was restricted to the high-risk younger mothers. We analysed 936 cases of gastroschisis from 25 population-based registries in 15 European countries, 1980-2002. We fitted a Bayesian Hierarchical Model which allowed us to estimate trend, to estimate which registries were significantly different from the common distribution, and to adjust simultaneously for maternal age, time (in grouped years) and the random variation between registries. The maternal age-standardised prevalence (standardised to the year 2000 European maternal age structure) increased almost fourfold from 0.54 [95% Credible Interval (CrI) 0.37, 0.75] per 10,000 births in 1980-84 to 2.12 [95% CrI 1.85, 2.40] per 10,000 births in 2000-02. The relative risk of gastroschisis for mothers <20 years of age in 1995-2002 was 7.0 [95% CrI 5.6, 8.7]. There were geographical differences within Europe, with higher rates of gastroschisis in the UK, and lower rates in Italy after adjusting for maternal age. After standardising for regional variation, our results showed that the increase in risk over time was the same for mothers of all ages--the increase for mothers <20 years was 3.96-fold compared with an increase of 3.95-fold for mothers in the other age groups. These findings indicate that the phenomenon of increasing gastroschisis prevalence is not restricted to younger mothers only.

Forecasting with missing data: Application to a real case

This paper presents a comparative analysis of linear and mixed modelsfor short term forecasting of a real data series with a high percentage of missing data. Data are the series of significant wave heights registered at regular periods of three hours by a buoy placed in the Bay of Biscay.The series is interpolated with a linear predictor which minimizes theforecast mean square error. The linear models are seasonal ARIMA models and themixed models have a linear component and a non linear seasonal component.The non linear component is estimated by a non parametric regression of dataversus time. Short term forecasts, no more than two days ahead, are of interestbecause they can be used by the port authorities to notice the fleet.Several models are fitted and compared by their forecasting behavior.

Rate of development of forensically-important Dipterain southern Brazil

Dipteran larvae were collected from rabbit (Oryctolagus cunniculus L.) carcasses during the four seasons in 2005 in the southernmost state of Rio Grande do Sul, Brazil. The larvae were fed ground beef at ambient temperatures following collection from carcasses. The development of each species under these conditions was estimated. The most abundant species in the carcasses were Lucilia eximia (Wiedemann) and Chrysomya albiceps (Wiedemann) (Calliphoridae), and they were found in all seasons. The data were fitted to a linear model that describes the relationship between temperature and linear developmental rating. These two species are primary forensic indicators in southern Brazil. Other species such as Hemilucilia semidiaphana (Rondani) (Calliphoridae), Synthesiomyia nudiseta (Wulp), Muscina stabulans (Fallen) (Muscidae), and Fannia pusio (Wiedemann) (Fanniidae) were forensically less important because they only occurred in high frequency in certain seasons and during the first days of carcass decomposition.

Estimating learning models from experimental data

We study the statistical properties of three estimation methods for a model of learning that is often fitted to experimental data: quadratic deviation measures without unobserved heterogeneity, and maximum likelihood withand without unobserved heterogeneity. After discussing identification issues, we show that the estimators are consistent and provide their asymptotic distribution. Using Monte Carlo simulations, we show that ignoring unobserved heterogeneity can lead to seriously biased estimations in samples which have the typical length of actual experiments. Better small sample properties areobtained if unobserved heterogeneity is introduced. That is, rather than estimating the parameters for each individual, the individual parameters are considered random variables, and the distribution of those random variables is estimated.

CABO VERDE: DE UM “ESTADO INVIÁVEL” AO PRAGMATISMO NA POLÍTICA EXTERNA

O artigo procura compreender de forma sincrónica as linhas fundamentais da política externa de Cabo Verde e seus reflexos nas relações internacionais, abarcando o período, desde a independência à actualidade. A perspectiva central é a de que a independência em 1975 representou uma nova era para o Estado de Cabo Verde, embora se tenha propalado a ideia de um Estado inviável. A estratégia para ultrapassar esta situação passava por uma aposta activa nas relações internacionais, nomeadamente numa política externa pragmática que se ajustasse num quadro de concertação e interacção com outros Estados. Procuramos, através de uma abordagem interdisciplinar e interpretativa, analisar os fenómenos sociais-internacionais que se configuram na política externa cabo-verdiana, sabendo que estes se relacionam tanto no campo da história e da sociologia como no campo da ciência política e das relações internacionais. Propõe-se novos desafios para a política externa cabo-verdiana, numa conjuntura internacional abalada pela crise financeira, uma aposta para uma maior integração regional no continente africano, aproveitando de forma estratégica o continente como nova fronteira de desenvolvimento.

Evaluation of the Adolescent Drug Abuse Diagnosis instrument in a Swiss sample of drug abusers.

OBJECTIVE: The aim of this study was to evaluate a French language version of the Adolescent Drug Abuse Diagnosis (ADAD) instrument in a Swiss sample of adolescent illicit drug and/or alcohol users. PARTICIPANTS AND SETTING: The participants in the study were 102 French-speaking adolescents aged 13-19 years who fitted the criteria of illicit drug or alcohol use (at least one substance--except tobacco--once a week during the last 3 months). They were recruited in hospitals, institutions and leisure places. Procedure. The ADAD was administered individually by trained psychologists. It was integrated into a broader protocol including alcohol and drug abuse DSM-IV diagnoses, the BDI-13 (Beck Depression Inventory), life events and treatment trajectories. RESULTS: The ADAD appears to show good inter-rater reliability; the subscales showed good internal coherence and the correlations between the composite scores and the severity ratings were moderate to high. Finally, the results confirmed good concurrent validity for three out of eight ADAD dimensions. CONCLUSIONS: The French language version of the ADAD appears to be an adequate instrument for assessing drug use and associated problems in adolescents. Despite its complexity, the instrument has acceptable validity, reliability and usefulness criteria, enabling international and transcultural comparisons.

Crecimiento de Platyxanthus orbignyi "Cangrejo violaceo" procedente de la caleta de Huanchaco - La Libertad durante el 2015

Esta investigación se realizó para determinar el crecimiento de Platyxanthus orbignyi “cangrejo violáceo” capturados en la caleta de Huanchaco, ubicado en el distrito de Huanchaco, Provincia de Trujillo, Departamento La Libertad, durante el 2015. El objetivo fue determinar el crecimiento de P. orbignyi, la muestra estuvo constituida por 1477 ejemplares procedentes de las capturas que realizan los pescadores artesanales empleando caballitos de totora y trampas denominadas nasas. Estos muestreos fueron mensuales durante el año 2015. La relación entre longitud cefalotorácica (mm) y ancho cefalotorácico (mm) presentó diferencia estadística entre las líneas de regresión de ambos sexos, la relación entre el ancho cefalotorácico (mm) y peso total (g) mostró alometría en el crecimiento de hembras y machos, estadísticamente no se encontró diferencia entre los sexos. Mediante el método de Bhattacharya (1967) modificado por Pauly y Caddy (1985) se establecieron siete estados de muda para las hembras y nueve estados de muda para los machos. Según el criterio de Kurata (1962) se estableció que el crecimiento en ambos sexos es de tipo geométrico retrogresivo. Los parámetros calculados se ajustaron a la ecuación de von Bertalanffy.

Comparison of the glomerular filtration rate in children by the new revised Schwartz formula and a new generalized formula.

The most widely used formula for estimating glomerular filtration rate (eGFR) in children is the Schwartz formula. It was revised in 2009 using iohexol clearances with measured GFR (mGFR) ranging between 15 and 75 ml/min × 1.73 m(2). Here we assessed the accuracy of the Schwartz formula using the inulin clearance (iGFR) method to evaluate its accuracy for children with less renal impairment comparing 551 iGFRs of 392 children with their Schwartz eGFRs. Serum creatinine was measured using the compensated Jaffe method. In order to find the best relationship between iGFR and eGFR, a linear quadratic regression model was fitted and a more accurate formula was derived. This quadratic formula was: 0.68 × (Height (cm)/serum creatinine (mg/dl))-0.0008 × (height (cm)/serum creatinine (mg/dl))(2)+0.48 × age (years)-(21.53 in males or 25.68 in females). This formula was validated using a split-half cross-validation technique and also externally validated with a new cohort of 127 children. Results show that the Schwartz formula is accurate until a height (Ht)/serum creatinine value of 251, corresponding to an iGFR of 103 ml/min × 1.73 m(2), but significantly unreliable for higher values. For an accuracy of 20 percent, the quadratic formula was significantly better than the Schwartz formula for all patients and for patients with a Ht/serum creatinine of 251 or greater. Thus, the new quadratic formula could replace the revised Schwartz formula, which is accurate for children with moderate renal failure but not for those with less renal impairment or hyperfiltration.

Numerical simulation of ambient seismic wavefield modification caused by pore-fluid effects in an oil reservoir

We have modeled numerically the seismic response of a poroelastic inclusion with properties applicable to an oil reservoir that interacts with an ambient wavefield. The model includes wave-induced fluid flow caused by pressure differences between mesoscopic-scale (i.e., in the order of centimeters to meters) heterogeneities. We used a viscoelastic approximation on the macroscopic scale to implement the attenuation and dispersion resulting from this mesoscopic-scale theory in numerical simulations of wave propagation on the kilometer scale. This upscaling method includes finite-element modeling of wave-induced fluid flow to determine effective seismic properties of the poroelastic media, such as attenuation of P- and S-waves. The fitted, equivalent, viscoelastic behavior is implemented in finite-difference wave propagation simulations. With this two-stage process, we model numerically the quasi-poroelastic wave-propagation on the kilometer scale and study the impact of fluid properties and fluid saturation on the modeled seismic amplitudes. In particular, we addressed the question of whether poroelastic effects within an oil reservoir may be a plausible explanation for low-frequency ambient wavefield modifications observed at oil fields in recent years. Our results indicate that ambient wavefield modification is expected to occur for oil reservoirs exhibiting high attenuation. Whether or not such modifications can be detected in surface recordings, however, will depend on acquisition design and noise mitigation processing as well as site-specific conditions, such as the geologic complexity of the subsurface, the nature of the ambient wavefield, and the amount of surface noise.

Uncertainty Quantification with Support Vector Regression Prediction Models

Uncertainty quantification of petroleum reservoir models is one of the present challenges, which is usually approached with a wide range of geostatistical tools linked with statistical optimisation or/and inference algorithms. The paper considers a data driven approach in modelling uncertainty in spatial predictions. Proposed semi-supervised Support Vector Regression (SVR) model has demonstrated its capability to represent realistic features and describe stochastic variability and non-uniqueness of spatial properties. It is able to capture and preserve key spatial dependencies such as connectivity, which is often difficult to achieve with two-point geostatistical models. Semi-supervised SVR is designed to integrate various kinds of conditioning data and learn dependences from them. A stochastic semi-supervised SVR model is integrated into a Bayesian framework to quantify uncertainty with multiple models fitted to dynamic observations. The developed approach is illustrated with a reservoir case study. The resulting probabilistic production forecasts are described by uncertainty envelopes.

Biokinetics and dosimetry of (111)In-DOTA-NOC-ATE compared with (111)In-DTPA-octreotide.

PURPOSE: The biokinetics and dosimetry of (111)In-DOTA-NOC-ATE (NOCATE), a high-affinity ligand of SSTR-2 and SSTR-5, and (111)In-DTPA-octreotide (Octreoscan?, OCTREO) were compared in the same patients. METHODS: Seventeen patients (10 men, 7 women; mean age 60 years), referred for an OCTREO scan for imaging of a neuroendocrine tumour (15), thymoma (1) or medullary thyroid carcinoma (1), agreed to undergo a second study with NOCATE. Whole-body anterior-posterior scans were recorded 0.5 (100 % reference scan), 4, 24 and 48 h (17 patients) and 120 h (5 patients) after injection. In 16 patients the OCTREO scan (178 ± 15 MBq) was performed 16 ± 5 days before the NOCATE scan (108 ± 14 MBq) with identical timing; 1 patient had the NOCATE scan before the OCTREO scan. Blood samples were obtained from 14 patients 5 min to 48 h after injection. Activities expressed as percent of the initial (reference) activity in the whole body, lung, kidney, liver, spleen and blood were fitted to biexponential or single exponential functions. Dosimetry was performed using OLINDA/EXM. RESULTS: Initial whole-body, lung and kidney activities were similar, but retention of NOCATE was higher than that of OCTREO. Liver and spleen uptakes of NOCATE were higher from the start (p < 0.001) and remained so over time. Whole-body activity showed similar α and β half-lives, but the β fraction of NOCATE was double that of OCTREO. Blood T (1/2)β for NOCATE was longer (19 vs. 6 h). As a result, the effective dose of NOCATE (105 μSv/MBq) exceeded that of OCTREO (52 μSv/MBq), and the latter result was similar to the ICRP 106 value of 54 μSv/MBq. Differential activity measurement in blood cells and plasma showed an average of <5 % of NOCATE and OCTREO attached to globular blood components. CONCLUSION: NOCATE showed a slower clearance from normal tissues and its effective dose was roughly double that of OCTREO.