Bayesian Hierarchical Distributed Lag Models for Summer Ozone Exposure and Cardio-Respiratory Mortality

In this paper, we develop Bayesian hierarchical distributed lag models for estimating associations between daily variations in summer ozone levels and daily variations in cardiovascular and respiratory (CVDRESP) mortality counts for 19 U.S. large cities included in the National Morbidity Mortality Air Pollution Study (NMMAPS) for the period 1987 - 1994. At the first stage, we define a semi-parametric distributed lag Poisson regression model to estimate city-specific relative rates of CVDRESP associated with short-term exposure to summer ozone. At the second stage, we specify a class of distributions for the true city-specific relative rates to estimate an overall effect by taking into account the variability within and across cities. We perform the calculations with respect to several random effects distributions (normal, t-student, and mixture of normal), thus relaxing the common assumption of a two-stage normal-normal hierarchical model. We assess the sensitivity of the results to: 1) lag structure for ozone exposure; 2) degree of adjustment for long-term trends; 3) inclusion of other pollutants in the model;4) heat waves; 5) random effects distributions; and 6) prior hyperparameters. On average across cities, we found that a 10ppb increase in summer ozone level for every day in the previous week is associated with 1.25 percent increase in CVDRESP mortality (95% posterior regions: 0.47, 2.03). The relative rate estimates are also positive and statistically significant at lags 0, 1, and 2. We found that associations between summer ozone and CVDRESP mortality are sensitive to the confounding adjustment for PM_10, but are robust to: 1) the adjustment for long-term trends, other gaseous pollutants (NO_2, SO_2, and CO); 2) the distributional assumptions at the second stage of the hierarchical model; and 3) the prior distributions on all unknown parameters. Bayesian hierarchical distributed lag models and their application to the NMMAPS data allow us estimation of an acute health effect associated with exposure to ambient air pollution in the last few days on average across several locations. The application of these methods and the systematic assessment of the sensitivity of findings to model assumptions provide important epidemiological evidence for future air quality regulations.

Modeling multilevel sleep transitional data via Poisson log-linear multilevel models

This paper proposes Poisson log-linear multilevel models to investigate population variability in sleep state transition rates. We specifically propose a Bayesian Poisson regression model that is more flexible, scalable to larger studies, and easily fit than other attempts in the literature. We further use hierarchical random effects to account for pairings of individuals and repeated measures within those individuals, as comparing diseased to non-diseased subjects while minimizing bias is of epidemiologic importance. We estimate essentially non-parametric piecewise constant hazards and smooth them, and allow for time varying covariates and segment of the night comparisons. The Bayesian Poisson regression is justified through a re-derivation of a classical algebraic likelihood equivalence of Poisson regression with a log(time) offset and survival regression assuming piecewise constant hazards. This relationship allows us to synthesize two methods currently used to analyze sleep transition phenomena: stratified multi-state proportional hazards models and log-linear models with GEE for transition counts. An example data set from the Sleep Heart Health Study is analyzed.

CD4+ T-cell count increase in HIV-1-infected patients with suppressed viral load within 1 year after start of antiretroviral therapy

BACKGROUND: CD4+ T-cell recovery in patients with continuous suppression of plasma HIV-1 viral load (VL) is highly variable. This study aimed to identify predictive factors for long-term CD4+ T-cell increase in treatment-naive patients starting combination antiretroviral therapy (cART). METHODS: Treatment-naive patients in the Swiss HIV Cohort Study reaching two VL measurements <50 copies/ml >3 months apart during the 1st year of cART were included (n=1816 patients). We studied CD4+ T-cell dynamics until the end of suppression or up to 5 years, subdivided into three periods: 1st year, years 2-3 and years 4-5 of suppression. Multiple median regression adjusted for repeated CD4+ T-cell measurements was used to study the dependence of CD4+ T-cell slopes on clinical covariates and drug classes. RESULTS: Median CD4+ T-cell increases following VL suppression were 87, 52 and 19 cells/microl per year in the three periods. In the multiple regression model, median CD4+ T-cell increases over all three periods were significantly higher for female gender, lower age, higher VL at cART start, CD4+ T-cell <650 cells/microl at start of the period and low CD4+ T-cell increase in the previous period. Patients on tenofovir showed significantly lower CD4+ T-cell increases compared with stavudine. CONCLUSIONS: In our observational study, long-term CD4+ T-cell increase in drug-naive patients with suppressed VL was higher in regimens without tenofovir. The clinical relevance of these findings must be confirmed in, ideally, clinical trials or large, collaborative cohort projects but could influence treatment of older patients and those starting cART at low CD4+ T-cell levels.

Car ownership modeling and forecasts for China

A great increase of private car ownership took place in China from 1980 to 2009 with the development of the economy. To explain the relationship between car ownership and economic and social changes, an ordinary least squares linear regression model is developed using car ownership per capita as the dependent variable with GDP, savings deposits and highway mileages per capita as the independent variables. The model is tested and corrected for econometric problems such as spurious correlation and cointegration. Finally, the regression model is used to project oil consumption by the Chinese transportation sector through 2015. The result shows that about 2.0 million barrels of oil will be consumed by private cars in conservative scenario, and about 2.6 million barrels of oil per day in high case scenario in 2015. Both of them are much higher than the consumption level of 2009, which is 1.9 million barrels per day. It also shows that the annual growth rate of oil demand by transportation is 2.7% - 3.1% per year in the conservative scenario, and 6.9% - 7.3% per year in the high case forecast scenario from 2010 to 2015. As a result, actions like increasing oil efficiency need to be taken to deal with challenges of the increasing demand for oil.

Publication and non-publication of clinical trials: longitudinal study of applications submitted to a research ethics committee

BACKGROUND: Not all clinical trials are published, which may distort the evidence that is available in the literature. We studied the publication rate of a cohort of clinical trials and identified factors associated with publication and nonpublication of results. METHODS: We analysed the protocols of randomized clinical trials of drug interventions submitted to the research ethics committee of University Hospital (Inselspital) Bern, Switzerland from 1988 to 1998. We identified full articles published up to 2006 by searching the Cochrane CENTRAL database (issue 02/2006) and by contacting investigators. We analyzed factors associated with the publication of trials using descriptive statistics and logistic regression models. RESULTS: 451 study protocols and 375 corresponding articles were analyzed. 233 protocols resulted in at least one publication, a publication rate of 52%. A total of 366 (81%) trials were commercially funded, 47 (10%) had non-commercial funding. 346 trials (77%) were multi-centre studies and 272 of these (79%) were international collaborations. In the adjusted logistic regression model non-commercial funding (Odds Ratio [OR] 2.42, 95% CI 1.14-5.17), multi-centre status (OR 2.09, 95% CI 1.03-4.24), international collaboration (OR 1.87, 95% CI 0.99-3.55) and a sample size above the median of 236 participants (OR 2.04, 95% CI 1.23-3.39) were associated with full publication. CONCLUSIONS: In this cohort of applications to an ethics committee in Switzerland, only about half of clinical drug trials were published. Large multi-centre trials with non-commercial funding were more likely to be published than other trials, but most trials were funded by industry.

Androgen and gonadotropin patterns differ in HIV-1-infected men who develop lipoatrophy during antiretroviral therapy: a case-control study

OBJECTIVES: We compared androgen and gonadotropin values in HIV-infected men who did and did not develop lipoatrophy on combination antiretroviral therapy (cART). METHODS: From a population of 136 treatment-naïve male Caucasians under successful zidovudine/lamivudine-based cART, the 10 patients developing lipoatrophy (cases) were compared with 87 randomly chosen controls. Plasma levels of free testosterone (fT), dehydroepiandrosterone (DHEA), follicle-stimulating hormone and luteinizing hormone (LH) were measured at baseline and after 2 years of cART. RESULTS: At baseline, 60% of the cases and 71% of the controls showed abnormally low fT values. LH levels were normal or low in 67 and 94% of the patients, respectively, indicating a disturbance of the hypothalamic-pituitary-gonadal axis. fT levels did not significantly change after 2 years of cART. Cases showed a significant increase in LH levels, while controls showed a significant increase in DHEA levels. In a multivariate logistic regression model, lipoatrophy was associated with higher baseline DHEA levels (P=0.04), an increase in LH levels during cART (P=0.001), a lower body mass index and greater age. CONCLUSIONS: Hypogonadism is present in the majority of HIV-infected patients. The development of cART-related lipoatrophy is associated with an increase in LH and a lack of increase in DHEA levels.

Comparison between two caller groups of a medical call centre in Switzerland

OBJECTIVES: The incidence distribution of triage advice in the medical call centre Medi24 and the pattern of service utilisation were analysed with respect to two groups of callers with different insurance schemes. Individuals having contracted insurance of the Medi24 model could use the telephone consultation service of the medical call centre Medi24 (mainly part of the mandatory basic health insurance) voluntarily and free of charge whereas individuals holding an insurance policy of the Telmed model (special contract within the mandatory basic health insurance with a premium discount ranging from 8% to 12%) were obliged to have a telephone consultation before arranging an appointment with a medical doctor. METHODS: A cross-sectional study was carried out in the medical call centre Medi24 based on all triage datasets of the Medi24 and Telmed groups collected during the one year period from July 1st 2005 to June 30th 2006. The distribution of the six different urgency levels within the two groups and their respective pattern of service utilisation was determined. In a multivariable logistic regression model the Odds Ratio for every enquiry originating from the Telmed group versus those originating from the Medi24 group was calculated. RESULTS: During a one-year period 48 388 triage requests reached the medical call centre Medi24, 56% derived from the Telmed group and 44% from the Medi24 group. Within the Medi24 group more than 25% of the individuals received self-care advice, within the Telmed group, on the other hand, only about 18% received such advice. In contrast, 27% of the Telmed triage requests but only 18% of the Medi24 triage requests resulted in the advice to make a routine appointment with a medical doctor. The probability that an individual of the Telmed group obtained the advice to go to the accident and emergency department was lower than for an individual of the Medi24 group (OR 0.77, 95% CI 0.60-0.99). Likewise, the probability of self-care advice was decreased in regard to the Medi24 group (OR 0.80, 95% CI 0.75-0.85). However, regarding the advice to make a routine appointment with a medical doctor, the Telmed group was represented more frequently than the Medi24 group (OR 1.36, 95% CI 1.28-1.44). CONCLUSION: In respect of the triage advice, the Telmed group differed significantly from the Medi24 group within all urgency levels. The differences between the two groups in respect of the advice given were still less pronounced than expected against the background of their different contract conditions and the disparate temporal pattern of utilisation. We interprete this finding with the fact that appraising the urgency of health problems appropriately seems to be very difficult for the majority of people seeking advice.

Interrelation of peri-operative morbidity and ASA class assignment in patients undergoing gynaecological surgery

OBJECTIVE: The aim of this study was to estimate intra- and post-operative risk using the American Society of Anaesthesiologists (ASA) classification which is an important predictor of an intervention and of the entire operating programme. STUDY DESIGN: In this retrospective study, 4435 consecutive patients undergoing elective and emergency surgery at the Gynaecological Clinic of the University Hospital of Zurich were included. The ASA classification for pre-operative risk assessment was determined by an anaesthesiologist after a thorough physical examination. We observed several pre-, intra- and post-operative parameters, such as age, body-mass-index, duration of anaesthesia, duration of surgery, blood loss, duration of post-operative stay, complicated post-operative course, morbidity and mortality. The investigation of different risk factors was achieved by a multiple linear regression model for log-transformed duration of hospitalisation. RESULTS: Age and obesity were responsible for a higher ASA classification. ASA grade correlates with the duration of anaesthesia and the duration of the surgery itself. There was a significant difference in blood loss between ASA grades I (113+/-195 ml) and III (222+/-470 ml) and between classes II (176+/-432 ml) and III. The duration of post-operative hospitalisation could also be correlated with ASA class. ASA class I=1.7+/-3.0 days, ASA class II=3.6+/-4.3 days, ASA class III=6.8+/-8.2 days, and ASA class IV=6.2+/-3.9 days. The mean post-operative in-hospital stay was 2.5+/-4.0 days without complications, and 8.7+/-6.7 days with post-operative complications. Multiple linear regression model showed that not only the ASA classification contained an important information for the duration of hospitalisation. Parameters such as age, class of diagnosis, post-operative complications, etc. also have an influence on the duration of hospitalisation. CONCLUSION: This study shows that the ASA classification can be used as a good and early available predictor for the planning of an intervention in gynaecological surgery. The ASA classification helps the surgeon to assess the peri-operative risk profile of which important information can be derived for the planning of the operation programme.

Heart rate response determines long term exercise capacity after heart transplantation

BACKGROUND: Exercise capacity after heart transplantation (HTx) remains limited despite normal left ventricular systolic function of the allograft. Various clinical and haemodynamic parameters are predictive of exercise capacity following HTx. However, the predictive significance of chronotropic competence has not been demonstrated unequivocally despite its immediate relevance for cardiac output. AIMS: This study assesses the predictive value of various clinical and haemodynamic parameters for exercise capacity in HTx recipients with complete chronotropic competence evolving within the first 6 postoperative months. METHODS: 51 patients were enrolled in this exercise study. Patients were included when at least >6 months after HTx and without negative chronotropic medication or factors limiting exercise capacity such as significant transplant vasculopathy or allograft rejection. Clinical parameters were obtained by chart review, haemodynamic parameters from current cardiac catheterisation, and exercise capacity was assessed by treadmill stress testing. A stepwise multiple regression model analysed the proportion of the variance explained by the predictive parameters. RESULTS: The mean age of these 51 HTx recipients was 55.4 +/- 13.2 yrs on inclusion, 42 pts were male and the mean time interval after cardiac transplantation was 5.1 +/- 2.8 yrs. Five independent predictors explained 47.5% of the variance observed for peak exercise capacity (adjusted R2 = 0.475). In detail, heart rate response explained 31.6%, male gender 5.2%, age 4.1%, pulmonary vascular resistance 3.7%, and body-mass index 2.9%. CONCLUSION: Heart rate response is one of the most important predictors of exercise capacity in HTx recipients with complete chronotropic competence and without relevant transplant vasculopathy or acute allograft rejection.

MET Y1253D-activating point mutation and development of distant metastasis in advanced head and neck cancers

We investigated if the MET-activating point mutation Y1253D influences clinical outcomes in patients with advanced squamous cell carcinoma of the head and neck (HNSCC). The study population consisted of 152 HNSCC patients treated by hyperfractionated radiotherapy alone or concomitant with chemotherapy between September 1994 and July 2000. Tumors were screened for the presence of the MET-activating point mutation Y1253D. Seventy-eight patients (51%) received radiotherapy alone, 74 patients (49%) underwent radiotherapy concomitant with chemotherapy. Median patient age was 54 years and median follow-up was 5.5 years. Distant metastasis-free survival, local relapse-free survival and overall survival were compared with MET Y1253D status. During follow-up, 29 (19%) patients developed distant metastasis. MET Y1253D was detected in tumors of 21 out of 152 patients (14%). Distant metastasis-free survival (P = 0.008) was associated with MET Y1253D. In a multivariate Cox regression model, adjusted for T-category, only presence of MET Y1253D was associated with decreased distant metastasis-free survival: hazard ratio = 2.5 (95% confidence interval: 1.1, 5.8). The observed association between MET Y1253D-activating point mutation and decreased distant metastasis-free survival in advanced HNSCC suggests that MET may be a potential target for specific treatment interventions.

Infection with human herpesvirus 8 and transplant-associated gammopathy

BACKGROUND: The role of human herpesvirus (HHV)-8 in the pathogenesis of multiple myeloma and its pre-malignant state of monoclonal gammopathy is unclear. HHV-8 is transmitted by organ transplantation, representing a unique model with which to investigate primary HHV-8 infection. METHODS: The authors studied the incidence of clonal gammopathy in renal transplant recipients and correlated it with previous and recent HHV-8 infection. RESULTS: Clonal gammopathy was observed in 31 of 162 (19%) HHV-8-seronegative patients, in 5 of 17 (29%) HHV-8-seropositive patients, and in 9 of 24 (38%) HHV-8 seroconverters within 5 years after transplantation. Gammopathy was often transient, and no progression to myeloma was observed. Two patients with persistent gammopathy developed B-cell lymphoma. In a logistic regression model, HHV-8 serostatus of the graft recipient was significantly associated with subsequent development of gammopathy, with a relative risk (RR) of 1.9 and a 95% confidence interval (CI) of 0.5 to 6.4 for an HHV-8-seropositive recipient and an RR of 2.9 and a 95% CI of 1.01 to 8.0 for seroconverters as compared with baseline (HHV-8 seronegative). Other significant variables were cytomegalovirus (CMV) serostatus and the intensity of immunosuppression (RR of 10.4 and 95% CI of 2.6-41.7 for a CMV-negative recipient with a CMV-positive donor vs. a CMV-negative recipient with a CMV-negative donor and RR of 17.6 and 95% CI of 2.0-150.8 if OKT3 was used vs. no use of antilymphocytic substances). CONCLUSIONS: Transplant recipients with HHV-8 infection are more likely to develop clonal gammopathy. However, this risk is much lower than the risk conferred by CMV infection and antilymphocytic therapy, arguing against a major role of HHV-8 infection in the pathogenesis of clonal plasma cell proliferation.

Effects of Ownership on Hospital Efficiency in Germany

The objective of our study was to evaluate the efficiency of public, private for-profit, and private non-profit hospitals in Germany. First, bootstrapped data envelopment analysis (DEA) was used to evaluate the efficiency of a panel (n = 1,046) of public, private for-profit, and private non-profit hospitals between 2002 and 2006. This was followed by a second-step truncated linear regression model with bootstrapped DEA efficiency scores as dependent variable. The results show that public hospitals performed significantly better than their private for-profit and non-profit counterparts. In addition, we found a significant positive association between hospital size and efficiency, and that competitive pressure had a significant negative impact on hospital efficiency.

The influence of a continuous rate infusion of dexmedetomidine on the nociceptive withdrawal reflex and temporal summation during isoflurane anaesthesia in dogs

Objective To examine the influence of a low dose dexmedetomidine infusion on the nociceptive withdrawal reflex and temporal summation in dogs during isoflurane anaesthesia. Study design Prospective experimental blinded cross-over study. Animals Eight healthy mixed breed dogs, body weight Mean +/- SD 26.5 +/- 8.4 kg and age 25 +/- 16 months. Methods Anaesthesia was induced with propofol and maintained with isoflurane (Fe'ISO 1.3%) delivered in oxygen and air. After stabilization, baseline recordings (time 0) were obtained, then a dexmedetomidine bolus (1 mug kg(-1) IV) followed by a continuous rate infusion (1 mug kg(-1) hour(-1) ) or saline placebo were administered. At times 10, 30 and 60 minutes after the initial bolus, electrical stimulations of increasing intensity were applied over the lateral plantar digital nerve, and administered both as single and as repeated stimuli. The resulting reflex responses were recorded using electromyography. Data were analysed using a multivariable linear regression model and a Kruskal Wallis test for single stimulation data, and repeated measures anova and paired t-test for repeated stimulation data. Results The AUC for the stimulus-response curves after single stimulation were similar for both treatments at time 0. At times 10, 30 and 60 the AUCs for the stimulus-response curves were significantly lower with dexmedetomidine treatment than with placebo. Temporal summation was evident in both treatments at times 0, 10, 30 and 60 starting from a stimulation intensity of 10 mA. The magnitude of temporal summation was smaller in dexmedetomidine than in placebo treated dogs at time 10, 30 and 60, but not at time 0. Conclusions During isoflurane anaesthesia, low dose dexmedetomidine suppresses the nociceptive reflex responses after single and repeated electrical stimulation. Clinical relevance This experimental study confirms previous reports on its peri-operative efficacy under clinical conditions, and further indicates that dexmedetomidine might reduce the risk of post-operative chronic pain development.

Critical Control Points for Profitability in the Cow-Calf Enterprise

Financial, economic, and biological data collected from cow-calf producers who participated in the Illinois and Iowa Standardized Performance Analysis (SPA) programs were used in this study. Data used were collected for the 1996 through 1999 calendar years, with each herd within year representing one observation. This resulted in a final database of 225 observations (117 from Iowa and 108 from Illinois) from commercial herds with a range in size from 20 to 373 cows. Two analyses were conducted, one utilizing financial cost of production data, the other economic cost of production data. Each observation was analyzed as the difference from the mean for that given year. The independent variable utilized in both the financial and economic models as an indicator of profit was return to unpaid labor and management per cow (RLM). Used as dependent variables were the five factors that make up total annual cow cost: feed cost, operating cost, depreciation cost, capital charge, and hired labor, all on an annual cost per cow basis. In the economic analysis, family labor was also included. Production factors evaluated as dependent variables in both models were calf weight, calf price, cull weight, cull price, weaning percentage, and calving distribution. Herd size and investment were also analyzed. All financial factors analyzed were significantly correlated to RLM (P < .10) except cull weight, and cull price. All economic factors analyzed were significantly correlated to RLM (P < .10) except calf weight, cull weight and cull price. Results of the financial prediction equation indicate that there are eight measurements capable of explaining over 82 percent of the farm-to-farm variation in RLM. Feed cost is the overriding factor driving RLM in both the financial and economic stepwise regression analyses. In both analyses over 50 percent of the herd-to-herd variation in RLM could be explained by feed cost. Financial feed cost is correlated (P < .001) to operating cost, depreciation cost, and investment. Economic feed cost is correlated (P < .001) with investment and operating cost, as well as capital charge. Operating cost, depreciation, and capital charge were all negatively correlated (P < .10) to herd size, and positively correlated (P < .01) to feed cost in both analyses. Operating costs were positively correlated with capital charge and investment (P < .01) in both analyses. In the financial regression model, depreciation cost was the second critical factor explaining almost 9 percent of the herd-to-herd variation in RLM followed by operating cost (5 percent). Calf weight had a greater impact than calf price on RLM in both the financial and economic regression models. Calf weight was the fourth indicator of RLM in the financial model and was similar in magnitude to operating cost. Investment was not a significant variable in either regression model; however, it was highly correlated to a number of the significant cost variables including feed cost, depreciation cost, and operating cost (P < .001, financial; P < .10, economic). Cost factors were far more influential in driving RLM than production, reproduction, or producer controlled marketing factors. Of these cost factors, feed cost had by far the largest impact. As producers focus attention on factors that affect the profitability of the operation, feed cost is the most critical control point because it was responsible for over 50 percent of the herd-to-herd variation in profit.

On trend estimation under monotone Gaussian subordination with long-memory: application to fossil pollen series

Fossil pollen data from stratigraphic cores are irregularly spaced in time due to non-linear age-depth relations. Moreover, their marginal distributions may vary over time. We address these features in a nonparametric regression model with errors that are monotone transformations of a latent continuous-time Gaussian process Z(T). Although Z(T) is unobserved, due to monotonicity, under suitable regularity conditions, it can be recovered facilitating further computations such as estimation of the long-memory parameter and the Hermite coefficients. The estimation of Z(T) itself involves estimation of the marginal distribution function of the regression errors. These issues are considered in proposing a plug-in algorithm for optimal bandwidth selection and construction of confidence bands for the trend function. Some high-resolution time series of pollen records from Lago di Origlio in Switzerland, which go back ca. 20,000 years are used to illustrate the methods.