Comparison of partial least squares regression and principal component regression for pelvic shape prediction.

This paper studied two different regression techniques for pelvic shape prediction, i.e., the partial least square regression (PLSR) and the principal component regression (PCR). Three different predictors such as surface landmarks, morphological parameters, or surface models of neighboring structures were used in a cross-validation study to predict the pelvic shape. Results obtained from applying these two different regression techniques were compared to the population mean model. In almost all the prediction experiments, both regression techniques unanimously generated better results than the population mean model, while the difference on prediction accuracy between these two regression methods is not statistically significant (α=0.01).

A BAYESIAN APPROACH TO DOSE-RESPONSE ASSESSMENT AND DRUG-DRUG INTERACTION ANALYSIS: APPLICATION TO IN VITRO STUDIES

The considerable search for synergistic agents in cancer research is motivated by the therapeutic benefits achieved by combining anti-cancer agents. Synergistic agents make it possible to reduce dosage while maintaining or enhancing a desired effect. Other favorable outcomes of synergistic agents include reduction in toxicity and minimizing or delaying drug resistance. Dose-response assessment and drug-drug interaction analysis play an important part in the drug discovery process, however analysis are often poorly done. This dissertation is an effort to notably improve dose-response assessment and drug-drug interaction analysis. The most commonly used method in published analysis is the Median-Effect Principle/Combination Index method (Chou and Talalay, 1984). The Median-Effect Principle/Combination Index method leads to inefficiency by ignoring important sources of variation inherent in dose-response data and discarding data points that do not fit the Median-Effect Principle. Previous work has shown that the conventional method yields a high rate of false positives (Boik, Boik, Newman, 2008; Hennessey, Rosner, Bast, Chen, 2010) and, in some cases, low power to detect synergy. There is a great need for improving the current methodology. We developed a Bayesian framework for dose-response modeling and drug-drug interaction analysis. First, we developed a hierarchical meta-regression dose-response model that accounts for various sources of variation and uncertainty and allows one to incorporate knowledge from prior studies into the current analysis, thus offering a more efficient and reliable inference. Second, in the case that parametric dose-response models do not fit the data, we developed a practical and flexible nonparametric regression method for meta-analysis of independently repeated dose-response experiments. Third, and lastly, we developed a method, based on Loewe additivity that allows one to quantitatively assess interaction between two agents combined at a fixed dose ratio. The proposed method makes a comprehensive and honest account of uncertainty within drug interaction assessment. Extensive simulation studies show that the novel methodology improves the screening process of effective/synergistic agents and reduces the incidence of type I error. We consider an ovarian cancer cell line study that investigates the combined effect of DNA methylation inhibitors and histone deacetylation inhibitors in human ovarian cancer cell lines. The hypothesis is that the combination of DNA methylation inhibitors and histone deacetylation inhibitors will enhance antiproliferative activity in human ovarian cancer cell lines compared to treatment with each inhibitor alone. By applying the proposed Bayesian methodology, in vitro synergy was declared for DNA methylation inhibitor, 5-AZA-2'-deoxycytidine combined with one histone deacetylation inhibitor, suberoylanilide hydroxamic acid or trichostatin A in the cell lines HEY and SKOV3. This suggests potential new epigenetic therapies in cell growth inhibition of ovarian cancer cells.

BAYESIAN PHASE I DOSE FINDING IN CANCER TRIALS

This dissertation explores phase I dose-finding designs in cancer trials from three perspectives: the alternative Bayesian dose-escalation rules, a design based on a time-to-dose-limiting toxicity (DLT) model, and a design based on a discrete-time multi-state (DTMS) model. We list alternative Bayesian dose-escalation rules and perform a simulation study for the intra-rule and inter-rule comparisons based on two statistical models to identify the most appropriate rule under certain scenarios. We provide evidence that all the Bayesian rules outperform the traditional ``3+3'' design in the allocation of patients and selection of the maximum tolerated dose. The design based on a time-to-DLT model uses patients' DLT information over multiple treatment cycles in estimating the probability of DLT at the end of treatment cycle 1. Dose-escalation decisions are made whenever a cycle-1 DLT occurs, or two months after the previous check point. Compared to the design based on a logistic regression model, the new design shows more safety benefits for trials in which more late-onset toxicities are expected. As a trade-off, the new design requires more patients on average. The design based on a discrete-time multi-state (DTMS) model has three important attributes: (1) Toxicities are categorized over a distribution of severity levels, (2) Early toxicity may inform dose escalation, and (3) No suspension is required between accrual cohorts. The proposed model accounts for the difference in the importance of the toxicity severity levels and for transitions between toxicity levels. We compare the operating characteristics of the proposed design with those from a similar design based on a fully-evaluated model that directly models the maximum observed toxicity level within the patients' entire assessment window. We describe settings in which, under comparable power, the proposed design shortens the trial. The proposed design offers more benefit compared to the alternative design as patient accrual becomes slower.

On the evolution of the oceanic component of the IPSL climate models from CMIP3 to CMIP5: A mean state comparison

This study analyses the impact on the oceanic mean state of the evolution of the oceanic component (NEMO) of the climate model developed at Institut Pierre Simon Laplace (IPSL-CM), from the version IPSL-CM4, used for third phase of the Coupled Model Intercomparison Project (CMIP3), to IPSL-CM5A, used for CMIP5. Several modifications have been implemented between these two versions, in particular an interactive coupling with a biogeochemical module, a 3-band model for the penetration of the solar radiation, partial steps at the bottom of the ocean and a set of physical parameterisations to improve the representation of the impact of turbulent and tidal mixing. A set of forced and coupled experiments is used to single out the effect of each of these modifications and more generally the evolution of the oceanic component on the IPSL coupled models family. Major improvements are located in the Southern Ocean, where physical parameterisations such as partial steps and tidal mixing reinforce the barotropic transport of water mass, in particular in the Antarctic Circumpolar Current) and ensure a better representation of Antarctic bottom water masses. However, our analysis highlights that modifications, which substantially improve ocean dynamics in forced configuration, can yield or amplify biases in coupled configuration. In particular, the activation of radiative biophysical coupling between biogeochemical cycle and ocean dynamics results in a cooling of the ocean mean state. This illustrates the difficulty to improve and tune coupled climate models, given the large number of degrees of freedom and the potential compensating effects masking some biases.

A convolution model for energy transport in a therapeutic fast neutron beam

A three-dimensional model has been proposed that uses Monte Carlo and fast Fourier transform convolution techniques to calculate the dose distribution from a fast neutron beam. This method transports scattered neutrons and photons in the forward, lateral, and backward directions and protons, electrons, and positrons in the forward and lateral directions by convolving energy spread kernels with initial interaction available energy distributions. The primary neutron and photon spectrums have been derived from narrow beam attenuation measurements. The positions and strengths of the effective primary neutron, scattered neutron, and photon sources have been derived from dual ion chamber measurements. The size of the effective primary neutron source has been measured using a copper activation technique. Heterogeneous tissue calculations require a weighted sum of two convolutions for each component since the kernels must be invariant for FFT convolution. Comparisons between calculations and measurements were performed for several water and heterogeneous phantom geometries. ^

Involvement of heparin-like glycosaminoglycans and expression of a novel heparan sulfate binding protein (p24) during human placentation and in a model for human implantation

Previous studies in our laboratory have indicated that heparan sulfate proteoglycans (HSPGs) play an important role in murine embryo implantation. To investigate the potential function of HSPGs in human implantation, two human cell lines (RL95 and JAR) were selected to model uterine epithelium and embryonal trophectoderm, respectively. A heterologous cell-cell adhesion assay showed that initial binding between JAR and RL95 cells is mediated by cell surface glycosaminoglycans (GAG) with heparin-like properties, i.e., heparan sulfate and dermatan sulfate. Furthermore, a single class of highly specific, protease-sensitive heparin/heparan sulfate binding sites exist on the surface of RL95 cells. Three heparin binding, tryptic peptide fragments were isolated from RL95 cell surfaces and their amino termini partially sequenced. Reverse transcription-polymerase chain reaction (RT-PCR) generated 1 to 4 PCR products per tryptic peptide. Northern blot analysis of RNA from RL95 cells using one of these RT-PCR products identified a 1.2 Kb mRNA species (p24). The amino acid sequence predicted from the cDNA sequence contains a putative heparin-binding domain. A synthetic peptide representing this putative heparin binding domain was used to generate a rabbit polyclonal antibody (anti-p24). Indirect immunofluorescence studies on RL95 and JAR cells as well as binding studies of anti-p24 to intact RL95 cells demonstrate that p24 is distributed on the cell surface. Western blots of RL95 membrane preparations identify a 24 kDa protein (p24) highly enriched in the 100,000 g pellet plasma membrane-enriched fraction. p24 eluted from membranes with 0.8 M NaCl, but not 0.6 M NaCl, suggesting that it is a peripheral membrane component. Solubilized p24 binds heparin by heparin affinity chromatography and $\sp{125}$I-heparin binding assays. Furthermore, indirect immunofluorescence studies indicate that cytotrophoblast of floating and attached villi of the human fetal-maternal interface are recognized by anti-p24. The study also indicates that the HSPG, perlecan, accumulates where chorionic villi are attached to uterine stroma and where p24-expressing cytotrophoblast penetrate the stroma. Collectively, these data indicate that p24 is a cell surface membrane-associated heparin/heparan sulfate binding protein found in cytotrophoblast, but not many other cell types of the fetal-maternal interface. Furthermore, p24 colocalizes with HSPGs in regions of cytotrophoblast invasion. These observations are consistent with a role for HSPGs and HSPG binding proteins in human trophoblast-uterine cell interactions. ^

Decadal Variability in the Northeast Pacific in a Physical-Ecosystem Model: Role of Mixed Layer Depth and Trophic Interactions

A basin-wide interdecadal change in both the physical state and the ecology of the North Pacific occurred near the end of 1976. Here we use a physical-ecosystem model to examine whether changes in the physical environment associated with the 1976-1977 transition influenced the lower trophic levels of the food web and if so by what means. The physical component is an ocean general circulation model, while the biological component contains 10 compartments: two phytoplankton, two zooplankton, two detritus pools, nitrate, ammonium, silicate, and carbon dioxide. The model is forced with observed atmospheric fields during 1960-1999. During spring, there is a similar to 40% reduction in plankton biomass in all four plankton groups during 1977-1988 relative to 1970-1976 in the central Gulf of Alaska (GOA). The epoch difference in plankton appears to be controlled by the mixed layer depth. Enhanced Ekman pumping after 1976 caused the halocline to shoal, and thus the mixed layer depth, which extends to the top of the halocline in late winter, did not penetrate as deep in the central GOA. As a result, more phytoplankton remained in the euphotic zone, and phytoplankton biomass began to increase earlier in the year after the 1976 transition. Zooplankton biomass also increased, but then grazing pressure led to a strong decrease in phytoplankton by April followed by a drop in zooplankton by May: Essentially, the mean seasonal cycle of plankton biomass was shifted earlier in the year. As the seasonal cycle progressed, the difference in plankton concentrations between epochs reversed sign again, leading to slightly greater zooplankton biomass during summer in the later epoch.

Social Support at the Workplace, Motivation to Transfer, and Training Transfer: A Multi-Level Indirect Effects Model

Supervisor support, peer support and transfer motivation have been identified as important predictors for training transfer. Transfer motivation is supposed to mediate the support–training transfer relationship. Especially after team training interventions that include all team members (i.e., intact-team training), individual perception of these factors might be shared among team members. However, an integration of the team level in the training transfer process is rare, yet still needed. Analyzing 194 employees from 34 teams in the context of intact-team training interventions, we found similar relationships and processes at both levels of analysis: Social support enhances transfer motivation at the individual and team levels. Furthermore, motivation to transfer increases training transfer and serves as a connecting mechanism in the social support–training transfer link. The results underline the importance of (1) considering multiple levels in theories and research about the training transfer process and (2) ensuring the practice of individual-directed support and a shared, supportive climate within teams.

General distress, hopelessness—suicidal ideation and worrying in adolescence: Concurrent and predictive validity of a symptom-level bifactor model for clinical diagnoses

BACKGROUND: Clinical disorders often share common symptoms and aetiological factors. Bifactor models acknowledge the role of an underlying general distress component and more specific sub-domains of psychopathology which specify the unique components of disorders over and above a general factor. METHODS: A bifactor model jointly calibrated data on subjective distress from The Mood and Feelings Questionnaire and the Revised Children's Manifest Anxiety Scale. The bifactor model encompassed a general distress factor, and specific factors for (a) hopelessness-suicidal ideation, (b) generalised worrying and (c) restlessness-fatigue at age 14 which were related to lifetime clinical diagnoses established by interviews at ages 14 (concurrent validity) and current diagnoses at 17 years (predictive validity) in a British population sample of 1159 adolescents. RESULTS: Diagnostic interviews confirmed the validity of a symptom-level bifactor model. The underlying general distress factor was a powerful but non-specific predictor of affective, anxiety and behaviour disorders. The specific factors for hopelessness-suicidal ideation and generalised worrying contributed to predictive specificity. Hopelessness-suicidal ideation predicted concurrent and future affective disorder; generalised worrying predicted concurrent and future anxiety, specifically concurrent generalised anxiety disorders. Generalised worrying was negatively associated with behaviour disorders. LIMITATIONS: The analyses of gender differences and the prediction of specific disorders was limited due to a low frequency of disorders other than depression. CONCLUSIONS: The bifactor model was able to differentiate concurrent and predict future clinical diagnoses. This can inform the development of targeted as well as non-specific interventions for prevention and treatment of different disorders.

Is the acetabular cup orientation after total hip arthroplasty on a two dimension or three dimension model accurate?

Purpose Malposition of the acetabular component in total hip arthroplasty (THA) is a common surgical problem that can lead to hip dislocation, reduced range of motion and may result in early loosening. The aim of this study is to validate the accuracy and reproducibility of a single x-ray image based 2D/3D reconstruction technique in determining cup inclination and anteversion against two different computer tomography (CT)-based measurement techniques. Methods Cup anteversion and inclination of 20 patients after cementless primary THA was measured on standard anteroposterior (AP) radiographs with the help of the single x-ray 2D/3D reconstruction program and compared with two different 3D CT-based analyses [Ground Truth (GT) and MeVis (MV) reconstruction model]. Results The measurements from the single x-ray 2D/3D reconstruction technique were strongly correlated with both types of CT image-processing protocols for both cup inclination [R²=0.69 (GT); R²=0.59 (MV)] and anteversion [R²=0.89 (GT); R²=0.80 (MV)]. Conclusions The single x-ray image based 2D/3D reconstruction technique is a feasible method to assess cup position on postoperative x-rays. CTscans remain the golden standard for a more complex biomechanical evaluation when a lower tolerance limit (+/-2 degrees) is required.

Transient Earth system responses to cumulative carbon dioxide emissions: linearities, uncertainties, and probabilities in an observation-constrained model ensemble

Information on the relationship between cumulative fossil CO2 emissions and multiple climate targets is essential to design emission mitigation and climate adaptation strategies. In this study, the transient response of a climate or environmental variable per trillion tonnes of CO2 emissions, termed TRE, is quantified for a set of impact-relevant climate variables and from a large set of multi-forcing scenarios extended to year 2300 towards stabilization. An  ∼ 1000-member ensemble of the Bern3D-LPJ carbon–climate model is applied and model outcomes are constrained by 26 physical and biogeochemical observational data sets in a Bayesian, Monte Carlo-type framework. Uncertainties in TRE estimates include both scenario uncertainty and model response uncertainty. Cumulative fossil emissions of 1000 Gt C result in a global mean surface air temperature change of 1.9 °C (68 % confidence interval (c.i.): 1.3 to 2.7 °C), a decrease in surface ocean pH of 0.19 (0.18 to 0.22), and a steric sea level rise of 20 cm (13 to 27 cm until 2300). Linearity between cumulative emissions and transient response is high for pH and reasonably high for surface air and sea surface temperatures, but less pronounced for changes in Atlantic meridional overturning, Southern Ocean and tropical surface water saturation with respect to biogenic structures of calcium carbonate, and carbon stocks in soils. The constrained model ensemble is also applied to determine the response to a pulse-like emission and in idealized CO2-only simulations. The transient climate response is constrained, primarily by long-term ocean heat observations, to 1.7 °C (68 % c.i.: 1.3 to 2.2 °C) and the equilibrium climate sensitivity to 2.9 °C (2.0 to 4.2 °C). This is consistent with results by CMIP5 models but inconsistent with recent studies that relied on short-term air temperature data affected by natural climate variability.

A vaccine formulation combining rhoptry proteins NcROP40 and NcROP2 improves pup survival in a pregnant mouse model of neosporosis.

Currently there are no effective vaccines for the control of bovine neosporosis. During the last years several subunit vaccines based on immunodominant antigens and other proteins involved in adhesion, invasion and intracellular proliferation of Neospora caninum have been evaluated as targets for vaccine development in experimental mouse infection models. Among them, the rhoptry antigen NcROP2 and the immunodominant NcGRA7 protein have been assessed with varying results. Recent studies have shown that another rhoptry component, NcROP40, and NcNTPase, a putative dense granule antigen, exhibit higher expression levels in tachyzoites of virulent N. caninum isolates, suggesting that these could be potential vaccine candidates to limit the effects of infection. In the present work, the safety and efficacy of these recombinant antigens formulated in Quil-A adjuvant as monovalent vaccines or pair-wise combinations (rNcROP40+rNcROP2 and rNcGRA7+rNcNTPase) were evaluated in a pregnant mouse model of neosporosis. All the vaccine formulations elicited a specific immune response against their respective native proteins after immunization. Mice vaccinated with rNcROP40 and rNcROP2 alone or in combination produced the highest levels of IFN-γ and exhibited low parasite burdens and low IgG antibody levels after the challenge. In addition, most of the vaccine formulations were able to increase the median survival time in the offspring. However, pup survival only ensued in the groups vaccinated with rNcROP40+rNcROP2 (16.2%) and rNcROP2 (6.3%). Interestingly, vertical transmission was not observed in those survivor pups immunized with rNcROP40+rNcROP2, as shown by PCR analyses. These results show a partial protection against N. caninum infection after vaccination with rNcROP40+rNcROP2, suggesting a synergistic effect of the two recombinant rhoptry antigens.

A Model Study of the Circulation in the Pearl River Estuary (PRE) and Its Adjacent Coastal Waters: 2. Sensitivity Experiments

The Princeton Ocean Model is used to study the circulation in the Pear River Estuary (PRE) and the adjacent coastal waters in the winter and summer seasons. Wong et al. [2003] compares the simulation results with the in situ measurements collected during the Pearl River Estuary Pollution Project (PREPP). In this paper, sensitivity experiments are carried out to examine the plume and the associated frontal dynamics in response to seasonal discharges and monsoon winds. During the winter, convergence between the seaward spreading plume water and the saline coastal water sets up a salinity front that aligns from the northeast to the southwest inside the PRE. During the summer the plume water fills the PRE at the surface and spreads eastward in the coastal waters in response to the prevailing southwesterly monsoon. The overall alignment of the plume is from the northwest to the southeast. The subsurface front is similar to that in the winter and summer except that the summer front is closer to the mouth and the winter front closer to the head of the estuary. Inside the PRE, bottom flows are always toward the head of the estuary, attributed to the density gradient associated with the plume front. In contrast, bottom flows in the shelf change from offshore in winter to onshore in summer, reflecting respectively the wintertime downwelling and summertime upwelling. Wind also plays an essential role in controlling the plume at the surface. An easterly wind drives the plume westward regardless winter or summer. The eastward spreading of the plume during the summer can be attributed to the southerly component of the wind. On the other hand, the surface area of the plume is positively proportional to the amount of discharge.

Using Leading Indicators to Forecast US Home Sales in a Bayesian VAR Framework

This paper uses Bayesian vector autoregressive models to examine the usefulness of leading indicators in predicting US home sales. The benchmark Bayesian model includes home sales, the price of homes, the mortgage rate, real personal disposable income, and the unemployment rate. We evaluate the forecasting performance of six alternative leading indicators by adding each, in turn, to the benchmark model. Out-of-sample forecast performance over three periods shows that the model that includes building permits authorized consistently produces the most accurate forecasts. Thus, the intention to build in the future provides good information with which to predict home sales. Another finding suggests that leading indicators with longer leads outperform the short-leading indicators.

Bayesian doubly optimal group sequential designs for randomized clinical trials

When conducting a randomized comparative clinical trial, ethical, scientific or economic considerations often motivate the use of interim decision rules after successive groups of patients have been treated. These decisions may pertain to the comparative efficacy or safety of the treatments under study, cost considerations, the desire to accelerate the drug evaluation process, or the likelihood of therapeutic benefit for future patients. At the time of each interim decision, an important question is whether patient enrollment should continue or be terminated; either due to a high probability that one treatment is superior to the other, or a low probability that the experimental treatment will ultimately prove to be superior. The use of frequentist group sequential decision rules has become routine in the conduct of phase III clinical trials. In this dissertation, we will present a new Bayesian decision-theoretic approach to the problem of designing a randomized group sequential clinical trial, focusing on two-arm trials with time-to-failure outcomes. Forward simulation is used to obtain optimal decision boundaries for each of a set of possible models. At each interim analysis, we use Bayesian model selection to adaptively choose the model having the largest posterior probability of being correct, and we then make the interim decision based on the boundaries that are optimal under the chosen model. We provide a simulation study to compare this method, which we call Bayesian Doubly Optimal Group Sequential (BDOGS), to corresponding frequentist designs using either O'Brien-Fleming (OF) or Pocock boundaries, as obtained from EaSt 2000. Our simulation results show that, over a wide variety of different cases, BDOGS either performs at least as well as both OF and Pocock, or on average provides a much smaller trial. ^