MODELING DIFFERENTIATED TREATMENT EFFECTS FOR MULTIPLE OUTCOMES DATA

Multiple outcomes data are commonly used to characterize treatment effects in medical research, for instance, multiple symptoms to characterize potential remission of a psychiatric disorder. Often either a global, i.e. symptom-invariant, treatment effect is evaluated. Such a treatment effect may over generalize the effect across the outcomes. On the other hand individual treatment effects, varying across all outcomes, are complicated to interpret, and their estimation may lose precision relative to a global summary. An effective compromise to summarize the treatment effect may be through patterns of the treatment effects, i.e. "differentiated effects." In this paper we propose a two-category model to differentiate treatment effects into two groups. A model fitting algorithm and simulation study are presented, and several methods are developed to analyze heterogeneity presenting in the treatment effects. The method is illustrated using an analysis of schizophrenia symptom data.

ON MARGINALIZED MULTILEVEL MODELS AND THEIR COMPUTATION

Clustered data analysis is characterized by the need to describe both systematic variation in a mean model and cluster-dependent random variation in an association model. Marginalized multilevel models embrace the robustness and interpretations of a marginal mean model, while retaining the likelihood inference capabilities and flexible dependence structures of a conditional association model. Although there has been increasing recognition of the attractiveness of marginalized multilevel models, there has been a gap in their practical application arising from a lack of readily available estimation procedures. We extend the marginalized multilevel model to allow for nonlinear functions in both the mean and association aspects. We then formulate marginal models through conditional specifications to facilitate estimation with mixed model computational solutions already in place. We illustrate this approach on a cerebrovascular deficiency crossover trial.

RANDOM EFFECTS MODELS IN A META-ANALYSIS OF THE ACCURACY OF DIAGNOSTIC TESTS WITHIN A GOLD STANDARD IN THE PRESENCE OF MISSING DATA

In evaluating the accuracy of diagnosis tests, it is common to apply two imperfect tests jointly or sequentially to a study population. In a recent meta-analysis of the accuracy of microsatellite instability testing (MSI) and traditional mutation analysis (MUT) in predicting germline mutations of the mismatch repair (MMR) genes, a Bayesian approach (Chen, Watson, and Parmigiani 2005) was proposed to handle missing data resulting from partial testing and the lack of a gold standard. In this paper, we demonstrate an improved estimation of the sensitivities and specificities of MSI and MUT by using a nonlinear mixed model and a Bayesian hierarchical model, both of which account for the heterogeneity across studies through study-specific random effects. The methods can be used to estimate the accuracy of two imperfect diagnostic tests in other meta-analyses when the prevalence of disease, the sensitivities and/or the specificities of diagnostic tests are heterogeneous among studies. Furthermore, simulation studies have demonstrated the importance of carefully selecting appropriate random effects on the estimation of diagnostic accuracy measurements in this scenario.

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.

Farmer land allocation for maize, groundnut and cotton production in Chipata District, Eastern Province, Zambia

Small-scale farmers in the Chipata District of Zambia rely on their farm fields to grow maize and groundnuts for food security. Cotton production and surplus food security crops are used to generate income to provide for their families. With increasing population pressure, available land has decreased and farmers struggle to provide the necessary food requirements and income to meet their family’s needs. The purpose of the study was to determine how a farmer can best allocate his land to produce maize, groundnuts and cotton when constrained by labor and capital resources to generate the highest potential for food security and financial gains. Data from the 2008-2009 growing season was compiled and analyzed using a linear programming model. The study determined that farmers make the most profit by allocating all additional land and resources to cotton after meeting their minimum food security requirements. The study suggests growing cotton is a beneficial practice for small-scale subsistence farmers to generate income when restricted by limited resources.

Using water balance models to approximate the effects of climate change on spring catchment discharge : Mt. Hanang, Tanzania

This project addresses the potential impacts of changing climate on dry-season water storage and discharge from a small, mountain catchment in Tanzania. Villagers and water managers around the catchment have experienced worsening water scarcity and attribute it to increasing population and demand, but very little has been done to understand the physical characteristics and hydrological behavior of the spring catchment. The physical nature of the aquifer was characterized and water balance models were calibrated to discharge observations so as to be able to explore relative changes in aquifer storage resulting from climate changes. To characterize the shallow aquifer supplying water to the Jandu spring, water quality and geochemistry data were analyzed, discharge recession analysis was performed, and two water balance models were developed and tested. Jandu geochemistry suggests a shallow, meteorically-recharged aquifer system with short circulation times. Baseflow recession analysis showed that the catchment behavior could be represented by a linear storage model with an average recession constant of 0.151/month from 2004-2010. Two modified Thornthwaite-Mather Water Balance (TMWB) models were calibrated using historic rainfall and discharge data and shown to reproduce dry-season flows with Nash-Sutcliffe efficiencies between 0.86 and 0.91. The modified TMWB models were then used to examine the impacts of nineteen, perturbed climate scenarios to test the potential impacts of regional climate change on catchment storage during the dry season. Forcing the models with realistic scenarios for average monthly temperature, annual precipitation, and seasonal rainfall distribution demonstrated that even small climate changes might adversely impact aquifer storage conditions at the onset of the dry season. The scale of the change was dependent on the direction (increasing vs. decreasing) and magnitude of climate change (temperature and precipitation). This study demonstrates that small, mountain aquifer characterization is possible using simple water quality parameters, recession analysis can be integrated into modeling aquifer storage parameters, and water balance models can accurately reproduce dry-season discharges and might be useful tools to assess climate change impacts. However, uncertainty in current climate projections and lack of data for testing the predictive capabilities of the model beyond the present data set, make the forecasts of changes in discharge also uncertain. The hydrologic tools used herein offer promise for future research in understanding small, shallow, mountainous aquifers and could potentially be developed and used by water resource professionals to assess climatic influences on local hydrologic systems.

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.

A STUDY ON REGIONAL CIRCULAR ECONOMY SYSTEM AND ITS CONSTRUCTION, OPERATION AND SUGGESTION FOR SHANGHAI

The characteristics of the traditional linear economic model are high consumption, high emission and low efficiency. Economic development is still largely at the expense of the environment and requires a natural resource investment. This can realize rapid economic development but resource depletion and environmental pollution become increasingly serious. In the 1990's a new economic model, circular economics, began to enter our vision. The circular economy maximizes production and minimizes the impact of economic activities on the ecological environment through organizing the activities through the closed-loop feedback cycle of "resources - production - renewable resource". Circular economy is a better way to solve the contradictions between the economic development and resource shortages. Developing circular economy has become the major strategic initiatives to achieving sustainable development in countries all over the world. The evaluation of the development of circular economics is a necessary step for regional circular economy development. Having a quantitative evaluation of circular economy can better monitor and reveal the contradictions and problems in the process of the development of recycling economy. This thesis will: 1) Create an evaluation model framework and new types of industries and 2) Make an evaluation of the Shanghai circular economy currently to analyze the situation of Shanghai in the development of circular economy. I will then propose suggestions about the structure and development of Shanghai circular economy.

Long-term gas exchange characteristics as marker of deterioration in patients with cystic fibrosis

Background and Aim In patients with cystic fibrosis (CF) the architecture of the developing lungs and the ventilation of lung units are progressively affected, influencing intrapulmonary gas mixing and gas exchange. We examined the long-term course of blood gas measurements in relation to characteristics of lung function and the influence of different CFTR genotype upon this process. Methods Serial annual measurements of PaO2 and PaCO2 assessed in relation to lung function, providing functional residual capacity (FRCpleth), lung clearance index (LCI), trapped gas (VTG), airway resistance (sReff), and forced expiratory indices (FEV1, FEF50), were collected in 178 children (88 males; 90 females) with CF, over an age range of 5 to 18 years. Linear mixed model analysis and binary logistic regression analysis were used to define predominant lung function parameters influencing oxygenation and carbon dioxide elimination. Results PaO2 decreased linearly from age 5 to 18 years, and was mainly associated with FRCpleth, (p < 0.0001), FEV1 (p < 0.001), FEF50 (p < 0.002), and LCI (p < 0.002), indicating that oxygenation was associated with the degree of pulmonary hyperinflation, ventilation inhomogeneities and impeded airway function. PaCO2 showed a transitory phase of low PaCO2 values, mainly during the age range of 5 to 12 years. Both PaO2 and PaCO2 presented with different progression slopes within specific CFTR genotypes. Conclusion In the long-term evaluation of gas exchange characteristics, an association with different lung function patterns was found and was closely related to specific genotypes. Early examination of blood gases may reveal hypocarbia, presumably reflecting compensatory mechanisms to improve oxygenation.

Mechanical characterization of brushite and hydroxyapatite cements

Compression, tension and torsion tests were designed and completed successfully on a brushite and a precipitated hydroxyapatite cement in moist condition. Elastic and strength properties were measured for these three loading cases. For each cement, the full set of strength data was fitted to an isotropic Tsai-Wu criterion and the associated coefficients identified. Since the compressive Young's moduli were about 10% larger than the tensile moduli, the full set of elastic data of each cement was fitted to a conewise linear elastic model. Hysteresis of the stress-strain curves was also observed, indicating dissipation mechanisms within these cement microstructures. A comparison of the measured mechanical properties with human cancellous bone confirmed the indication of brushite as a bone filling material and the potential of the hydroxyapatite cement as a structural biomaterial.

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.

Brain metabolite composition during early human brain development as measured by quantitative in vivo 1H magnetic resonance spectroscopy

Biochemical maturation of the brain can be studied noninvasively by (1)H magnetic resonance spectroscopy (MRS) in human infants. Detailed time courses of cerebral tissue contents are known for the most abundant metabolites only, and whether or not premature birth affects biochemical maturation of the brain is disputed. Hence, the last trimester of gestation was observed in infants born prematurely, and their cerebral metabolite contents at birth and at expected term were compared with those of fullterm infants. Successful quantitative short-TE (1)H MRS was performed in three cerebral locations in 21 infants in 28 sessions (gestational age 32-43 weeks). The spectra were analyzed with linear combination model fitting, considerably extending the range of observable metabolites to include acetate, alanine, aspartate, cholines, creatines, gamma-aminobutyrate, glucose, glutamine, glutamate, glutathione, glycine, lactate, myo-inositol, macromolecular contributions, N-acetylaspartate, N-acetylaspartylglutamate, o-phosphoethanolamine, scyllo-inositol, taurine, and threonine. Significant effects of age and location were found for many metabolites, including the previously observed neuronal maturation reflected by an increase in N-acetylaspartate. Absolute brain metabolite content in premature infants at term was not considerably different from that in fullterm infants, indicating that prematurity did not affect biochemical brain maturation substantially in the studied population, which did not include infants of extremely low birthweight.

Patient classification system: An integrated method for measuring nursing intensity and optimizing resource allocation

A patient classification system was developed integrating a patient acuity instrument with a computerized nursing distribution method based on a linear programming model. The system was designed for real-time measurement of patient acuity (workload) and allocation of nursing personnel to optimize the utilization of resources.^ The acuity instrument was a prototype tool with eight categories of patients defined by patient severity and nursing intensity parameters. From this tool, the demand for nursing care was defined in patient points with one point equal to one hour of RN time. Validity and reliability of the instrument was determined as follows: (1) Content validity by a panel of expert nurses; (2) predictive validity through a paired t-test analysis of preshift and postshift categorization of patients; (3) initial reliability by a one month pilot of the instrument in a practice setting; and (4) interrater reliability by the Kappa statistic.^ The nursing distribution system was a linear programming model using a branch and bound technique for obtaining integer solutions. The objective function was to minimize the total number of nursing personnel used by optimally assigning the staff to meet the acuity needs of the units. A penalty weight was used as a coefficient of the objective function variables to define priorities for allocation of staff.^ The demand constraints were requirements to meet the total acuity points needed for each unit and to have a minimum number of RNs on each unit. Supply constraints were: (1) total availability of each type of staff and the value of that staff member (value was determined relative to that type of staff's ability to perform the job function of an RN (i.e., value for eight hours RN = 8 points, LVN = 6 points); (2) number of personnel available for floating between units.^ The capability of the model to assign staff quantitatively and qualitatively equal to the manual method was established by a thirty day comparison. Sensitivity testing demonstrated appropriate adjustment of the optimal solution to changes in penalty coefficients in the objective function and to acuity totals in the demand constraints.^ Further investigation of the model documented: correct adjustment of assignments in response to staff value changes; and cost minimization by an addition of a dollar coefficient to the objective function. ^

PREDICTIVE VARIABLES FOR WEIGHT CHANGE 2 YEARS FOLLOWING PARTICIPATION IN A WEIGHT LOSS INTERVENTION PROGRAM

A historical prospective study was designed to assess the man weight status of subjects who participated in a behavioral weight reduction program in 1983 and to determine whether there was an association between the dependent variable weight change and any of 31 independent variables after a 2 year follow-up period. Data was obtained by abstracting the subjects records and from a follow-up questionnaire administered 2 years following program participation. Five hundred nine subjects (386 females and 123 males) of 1460 subjects who participated in the program, completed and returned the questionnaire. Results showed that mean weight was significantly different (p < 0.001) between the measurement at baseline and after a 2 year follow-up period. The mean weight loss of the group was 5.8 pounds, 10.7 pounds for males and 4.2 pounds for females after a 2 year follow-up period. A total of 63.9% of the group, 69.9% of males and 61.9% of females were still below their initial weight after the 2 year follow-up period. Sixteen of the 31 variables assessed utilizing bivariate analyses were found to be significantly (p (LESSTHEQ) 0.05) associated with weight change after a 2 year follow-up period. These variables were then entered into a multivariate linear regression model. A total of 37.9% of the variance of the dependent variable, weight change, was accounted for by all 16 variables. Eight of these variables were found to be significantly (p (LESSTHEQ) 0.05) predictive of weight change in the stepwise multivariate process accounting for 37.1% of the variance. These variables included: Two baseline variables (percent over ideal body weight at enrollment and occupation) and six follow-up variables (feeling in control of eating habits, percent of body weight lost during treatment, frequency of weight measurement, physical activity, eating in response to emotions, and number of pounds of weight gain needed to resume a diet). It was concluded that a greater amount of emphasis should be placed on the six follow-up variables by clinicians involved in the treatment of obesity, and by the subjects themselves to enhance their chances of success at long-term weight loss. ^