Bayesian change-point analysis in linear regression model with scale mixtures of normal distributions

In this thesis, we consider Bayesian inference on the detection of variance change-point models with scale mixtures of normal (for short SMN) distributions. This class of distributions is symmetric and thick-tailed and includes as special cases: Gaussian, Student-t, contaminated normal, and slash distributions. The proposed models provide greater flexibility to analyze a lot of practical data, which often show heavy-tail and may not satisfy the normal assumption. As to the Bayesian analysis, we specify some prior distributions for the unknown parameters in the variance change-point models with the SMN distributions. Due to the complexity of the joint posterior distribution, we propose an efficient Gibbs-type with Metropolis- Hastings sampling algorithm for posterior Bayesian inference. Thereafter, following the idea of [1], we consider the problems of the single and multiple change-point detections. The performance of the proposed procedures is illustrated and analyzed by simulation studies. A real application to the closing price data of U.S. stock market has been analyzed for illustrative purposes.

Survival of bonded lingual retainers with chemical or photo polymerization over a 2-year period: a single-center, randomized controlled clinical trial

INTRODUCTION The objective of this trial was to compare the survival rates of mandibular lingual retainers bonded with either chemically cured or light-cured adhesive after orthodontic treatment. METHODS Patients having undergone orthodontic treatment at a private orthodontic office were randomly allocated to fixed retainers placed with chemically cured composite or light-cured composite. Eligibility criteria included no active caries, restorations, or fractures on the mandibular anterior teeth, and adequate oral hygiene. The main outcome was any type of first-time lingual retainer breakage; pattern of failure (adapted adhesive remnant index scores) was a secondary outcome. Randomization was accomplished with random permuted blocks of 20 patients with allocation concealed in sequentially numbered, opaque, sealed envelopes. Blinding was applicable for outcome assessment only. Patients were reviewed at 1, 3, and 6 months and then every 6 months after placement of the retainer until completion of the study. Data were analyzed using survival analysis including Cox regression; sensitivity analysis was carried out after data imputation for subjects lost to follow-up. RESULTS Two hundred twenty patients (median age, 16 years; interquartile range, 2; range, 12-47 years) were randomized in a 1:1 ratio to either chemical or light curing. Baseline characteristics were similar between groups, the median follow-up period was 2.19 years (range, 0.003-3.64 years), and 16 patients were lost to follow-up. At a minimum follow-up of 2 years, 47 of 110 (42.7%) and 55 of 110 (50.0%) retainers had some type of failure with chemically cured and light-cured adhesive, respectively (log-rank test, P = 0.35). Data were analyzed on an intention-to-treat basis, and the hazard ratio (HR) was 1.15 (95% confidence interval [CI], 0.88-1.70; P = 0.47). There was weak evidence that age is a significant predictor for lingual retainer failures (HR, 0.96; 95% CI, 0.93-1.00; P = 0.08). Adhesive remnant index scoring was possible for only 66 of the 102 (64.7%) failures and did not differ between composites (Fisher exact test, P = 0.16). No serious harm was observed other than gingivitis associated with plaque accumulation. CONCLUSIONS The results of this study indicated no evidence that survival of mandibular lingual retainers differs between chemically and light-cured adhesives. The overall failure rate was 46.4%; however, this included any type of failure, which may have exaggerated the overall failure rate.

Prognosis of Children with HIV-1 Infection Starting Antiretroviral Therapy in Southern Africa: A Collaborative Analysis of Treatment Programs

BACKGROUND: Prognostic models for children starting antiretroviral therapy (ART) in Africa are lacking. We developed models to estimate the probability of death during the first year receiving ART in Southern Africa. METHODS: We analyzed data from children ≤10 years old who started ART in Malawi, South Africa, Zambia or Zimbabwe from 2004-2010. Children lost to follow-up or transferred were excluded. The primary outcome was all-cause mortality in the first year of ART. We used Weibull survival models to construct two prognostic models: one with CD4%, age, WHO clinical stage, weight-for-age z-score (WAZ) and anemia and one without CD4%, because it is not routinely measured in many programs. We used multiple imputation to account for missing data. RESULTS: Among 12655 children, 877 (6.9%) died in the first year of ART. 1780 children were lost to follow-up/transferred and excluded from main analyses; 10875 children were included. With the CD4% model probability of death at 1 year ranged from 1.8% (95% CI: 1.5-2.3) in children 5-10 years with CD4% ≥10%, WHO stage I/II, WAZ ≥-2 and without severe anemia to 46.3% (95% CI: 38.2-55.2) in children <1 year with CD4% <5%, stage III/IV, WAZ< -3 and severe anemia. The corresponding range for the model without CD4% was 2.2% (95% CI: 1.8-2.7) to 33.4% (95% CI: 28.2-39.3). Agreement between predicted and observed mortality was good (C-statistics=0.753 and 0.745 for models with and without CD4% respectively). CONCLUSION: These models may be useful to counsel children/caregivers, for program planning and to assess program outcomes after allowing for differences in patient disease severity characteristics.

Heterotic T-folds with a small number of neutral moduli

We discuss non-geometric supersymmetric heterotic string models in D=4, in the framework of the free fermionic construction. We perform a systematic scan of models with four a priori left-right asymmetric Z2 projections and shifts. We analyze some 220 models, identifying 18 inequivalent classes and addressing variants generated by discrete torsions. They do not contain geometrical or trivial neutral moduli, apart from the dilaton. However, we show the existence of flat directions in the form of exactly marginal deformations and identify patterns of symmetry breaking where product gauge groups, realized at level one, are broken to their diagonal at higher level. We also describe an “inverse Gepner map” from Heterotic to Type II models that could be used, in certain non geometric settings, to define “effective” topological invariants.

THE MECHANISM OF TUMORIGENESIS IN THE IMMORTALIZED HUMAN PANCREATIC CELL LINES: CELL CULTURE MODELS OF HUMAN PANCREATIC CANCER

The mechanism of tumorigenesis in the immortalized human pancreatic cell lines: cell culture models of human pancreatic cancer Pancreatic ductal adenocarcinoma (PDAC) is the most lethal cancer in the world. The most common genetic lesions identified in PDAC include activation of K-ras (90%) and Her2 (70%), loss of p16 (95%) and p14 (40%), inactivation p53 (50-75%) and Smad4 (55%). However, the role of these signature gene alterations in PDAC is still not well understood, especially, how these genetic lesions individually or in combination contribute mechanistically to human pancreatic oncogenesis is still elusive. Moreover, a cell culture transformation model with sequential accumulation of signature genetic alterations in human pancreatic ductal cells that resembles the multiple-step human pancreatic carcinogenesis is still not established. In the present study, through the stepwise introduction of the signature genetic alterations in PDAC into the HPV16-E6E7 immortalized human pancreatic duct epithelial (HPDE) cell line and the hTERT immortalized human pancreatic ductal HPNE cell line, we developed the novel experimental cell culture transformation models with the most frequent gene alterations in PDAC and further dissected the molecular mechanism of transformation. We demonstrated that the combination of activation of K-ras and Her2, inactivation of p16/p14 and Smad4, or K-ras mutation plus p16 inactivation, was sufficient for the tumorigenic transformation of HPDE or HPNE cells respectively. We found that these transformed cells exhibited enhanced cell proliferation, anchorage-independent growth in soft agar, and grew tumors with PDAC histopathological features in orthotopic mouse model. Molecular analysis showed that the activation of K-ras and Her2 downstream effector pathways –MAPK, RalA, FAK, together with upregulation of cyclins and c-myc were involved in the malignant transformation. We discovered that MDM2, BMP7 and Bmi-1 were overexpressed in the tumorigenic HPDE cells, and that Smad4 played important roles in regulation of BMP7 and Bmi-1 gene expression and the tumorigenic transformation of HPDE cells. IPA signaling pathway analysis of microarray data revealed that abnormal signaling pathways are involved in transformation. This study is the first complete transformation model of human pancreatic ductal cells with the most common gene alterations in PDAC. Altogether, these novel transformation models more closely recapitulate the human pancreatic carcinogenesis from the cell origin, gene lesion, and activation of specific signaling pathway and histopathological features.

THE EFFECTS OF CURCUMIN ON BODY COMPOSITION OF PATIENTS WITH ADVANCED PANCREATIC CANCER

Cachexia is very common among patients with advanced pancreatic cancer and is a marker of poor prognosis. Weight loss in cachexia is due to both adipose and muscle compartments, and sarcopenia (severe muscle depletion) is associated with worse outcomes. Curcumin has shown a myriad of biological effects, including anti-cancer and anti-inflammatory. The ability of curcumin to attenuate cachexia and muscle loss has been tested in animal models, with conflicting results so far. The hypothesis of this study was that patients with advanced pancreatic cancer treated with curcumin for two months have less fat and muscle loss as compared to matched controls not treated with this compound. A matched 1:2 case-control retrospective study was conducted with 22 patients with pancreatic cancer who were treated with curcumin on a previous protocol and 44 untreated controls with the same diagnosis matched by age, gender, time from advanced cancer, body mass index, and number of prior therapies. Data was collected regarding oncologic treatment, medication use, weights, heights, and survival. Body composition was determined by computerized tomography analyses at two timepoints separated by 60±20 days. For treated patients, the first image was at the beginning of treatment and for controls it was determined by the matching time from advanced cancer. The evolution of body composition over time was quantitatively analyzed comparing both groups. All patients lost weight both due to fat and muscle losses, and patients treated with curcumin presented greater losses both in lean adipose body mass. Use of medications, chemotherapy, age, time from advanced cancer, baseline albumin, performance status, and number of prior therapies were not independently correlated with changes in body composition variables. Patients treated with curcumin had borderline shorter survival when compared with untreated patients. Sarcopenic treated patients had significantly shorter survival than non-sarcopenic counterparts, and sarcopenia status was not associated with survival among the controls. Treated patients with shorter survival showed a tendency to lose more lean and especially fat body mass as compared to untreated patients, maybe suggesting an effect of curcumin on shifting weight loss towards the end of life by impacting its mechanisms.

Sequential updating via dynamic hierarchical models: A Bayesian approach to longitudinal data analysis

Most statistical analysis, theory and practice, is concerned with static models; models with a proposed set of parameters whose values are fixed across observational units. Static models implicitly assume that the quantified relationships remain the same across the design space of the data. While this is reasonable under many circumstances this can be a dangerous assumption when dealing with sequentially ordered data. The mere passage of time always brings fresh considerations and the interrelationships among parameters, or subsets of parameters, may need to be continually revised. ^ When data are gathered sequentially dynamic interim monitoring may be useful as new subject-specific parameters are introduced with each new observational unit. Sequential imputation via dynamic hierarchical models is an efficient strategy for handling missing data and analyzing longitudinal studies. Dynamic conditional independence models offers a flexible framework that exploits the Bayesian updating scheme for capturing the evolution of both the population and individual effects over time. While static models often describe aggregate information well they often do not reflect conflicts in the information at the individual level. Dynamic models prove advantageous over static models in capturing both individual and aggregate trends. Computations for such models can be carried out via the Gibbs sampler. An application using a small sample repeated measures normally distributed growth curve data is presented. ^

Assessment of Skill and Portability in Regional Marine Biogeochemical Models: Role of Multiple Planktonic Groups

Application of biogeochemical models to the study of marine ecosystems is pervasive, yet objective quantification of these models' performance is rare. Here, 12 lower trophic level models of varying complexity are objectively assessed in two distinct regions (equatorial Pacific and Arabian Sea). Each model was run within an identical one-dimensional physical framework. A consistent variational adjoint implementation assimilating chlorophyll-a, nitrate, export, and primary productivity was applied and the same metrics were used to assess model skill. Experiments were performed in which data were assimilated from each site individually and from both sites simultaneously. A cross-validation experiment was also conducted whereby data were assimilated from one site and the resulting optimal parameters were used to generate a simulation for the second site. When a single pelagic regime is considered, the simplest models fit the data as well as those with multiple phytoplankton functional groups. However, those with multiple phytoplankton functional groups produced lower misfits when the models are required to simulate both regimes using identical parameter values. The cross-validation experiments revealed that as long as only a few key biogeochemical parameters were optimized, the models with greater phytoplankton complexity were generally more portable. Furthermore, models with multiple zooplankton compartments did not necessarily outperform models with single zooplankton compartments, even when zooplankton biomass data are assimilated. Finally, even when different models produced similar least squares model-data misfits, they often did so via very different element flow pathways, highlighting the need for more comprehensive data sets that uniquely constrain these pathways.

Dispersal Modeling of Fish Early Life Stages: Sensitivity with Application to Atlantic Cod in the Western Gulf of Maine

As an initial step in establishing mechanistic relationships between environmental variability and recruitment in Atlantic cod Gadhus morhua along the coast of the western Gulf of Maine, we assessed transport success of larvae from major spawning grounds to nursery areas with particle tracking using the unstructured grid model FVCOM (finite volume coastal ocean model). In coastal areas, dispersal of early planktonic life stages of fish and invertebrate species is highly dependent on the regional dynamics and its variability, which has to be captured by our models. With state-of-the-art forcing for the year 1995, we evaluate the sensitivity of particle dispersal to the timing and location of spawning, the spatial and temporal resolution of the model, and the vertical mixing scheme. A 3 d frequency for the release of particles is necessary to capture the effect of the circulation variability into an averaged dispersal pattern of the spawning season. The analysis of sensitivity to model setup showed that a higher resolution mesh, tidal forcing, and current variability do not change the general pattern of connectivity, but do tend to increase within-site retention. Our results indicate strong downstream connectivity among spawning grounds and higher chances for successful transport from spawning areas closer to the coast. The model run for January egg release indicates 1 to 19 % within-spawning ground retention of initial particles, which may be sufficient to sustain local populations. A systematic sensitivity analysis still needs to be conducted to determine the minimum mesh and forcing resolution that adequately resolves the complex dynamics of the western Gulf of Maine. Other sources of variability, i.e. large-scale upstream forcing and the biological environment, also need to be considered in future studies of the interannual variability in transport and survival of the early life stages of cod.

Recurrent stroke predictors differ in medically treated patients with pathogenic vs other PFOs

OBJECTIVE: To examine predictors of stroke recurrence in patients with a high vs a low likelihood of having an incidental patent foramen ovale (PFO) as defined by the Risk of Paradoxical Embolism (RoPE) score.METHODS: Patients in the RoPE database with cryptogenic stroke (CS) and PFO were classified as having a probable PFO-related stroke (RoPE score of >6, n = 647) and others (RoPE score of with stroke recurrence using Cox survival models with component database as a stratification factor. An interaction with RoPE score was checked for the variables that were significant.RESULTS: Follow-up was available for 92%, 79%, and 57% at 1, 2, and 3 years. Overall, a higher recurrence risk was associated with an index TIA. For all other predictors, effects were significantly different in the 2 RoPE score categories. For the low RoPE score group, but not the high RoPE score group, older age and antiplatelet (vs warfarin) treatment predicted recurrence. Conversely, echocardiographic features (septal hypermobility and a small shunt) and a prior (clinical) stroke/TIA were significant predictors in the high but not low RoPE score group.CONCLUSION: Predictors of recurrence differ when PFO relatedness is classified by the RoPE score, suggesting that patients with CS and PFO form a heterogeneous group with different stroke mechanisms. Echocardiographic features were only associated with recurrence in the high RoPE score group.

Multi-channel search for squarks and gluinos in sqrts=7 TeV pp collisions with the ATLAS detector

A search for supersymmetric particles in final states with zero, one, and two leptons, with and without jets identified as originating from b-quarks, in 4.7 fb(-1) of root s = 7 TeV pp collisions produced by the Large Hadron Collider and recorded by the ATLAS detector is presented. The search uses a set of variables carrying information on the event kinematics transverse and parallel to the beam line that are sensitive to several topologies expected in supersymmetry. Mutually exclusive final states are defined, allowing a combination of all channels to increase the search sensitivity. No deviation from the Standard Model expectation is observed. Upper limits at 95 % confidence level on visible cross-sections for the production of new particles are extracted. Results are interpreted in the context of the constrained minimal supersymmetric extension to the Standard Model and in supersymmetry-inspired models with diverse, high-multiplicity final states.

Search for supersymmetry in events with photons, bottom quarks, and missing transverse momentum in proton-proton collisions at a centre-of-mass energy of 7 TeV with the ATLAS detector

A search has been performed for the experimental signature of an isolated photon with high transverse momentum, at least one jet identified as originating from a bottom quark, and high missing transverse momentum. Such a final state may originate from supersymmetric models with gauge-mediated supersymmetry breaking in events in which one of a pair of higgsino-like neutralinos decays into a photon and a gravitino while the other decays into a Higgs boson and a gravitino. The search is performed using the full dataset of 7 TeV proton-proton collisions recorded with the ATLAS detector at the LHC in 2011, corresponding to an integrated luminosity of 4.7 fb(-1). A total of 7 candidate events are observed while 7.5 +/- 2.2 events are expected from the Standard Model background. The results of the search are interpreted in the context of general gauge mediation to exclude certain regions of a benchmark plane for higgsino-like neutralino production.

Search for dark matter candidates and large extra dimensions in events with a photon and missing transverse momentum in pp collision data at sqrts=7 TeV w the ATLAS detector

Results of a search for new phenomena in events with an energetic photon and large missing transverse momentum in proton-proton collisions at root s = 7 TeV are reported. Data collected by the ATLAS experiment at the LHC corresponding to an integrated luminosity of 4.6 fb(-1) are used. Good agreement is observed between the data and the standard model predictions. The results are translated into exclusion limits on models with large extra spatial dimensions and on pair production of weakly interacting dark matter candidates.

Search for squarks and gluinos with the ATLAS detector in final states with jets and missing transverse momentum using 4.7 fb⁻¹ of sqrts=7 TeV proton-pro collision data

A search for squarks and gluinos in final states containing jets, missing transverse momentum and no high-p(T) electrons or muons is presented. The data represent the complete sample recorded in 2011 by the ATLAS experiment in 7 TeV proton-proton collisions at the Large Hadron Collider, with a total integrated luminosity of 4.7 fb(-1). No excess above the Standard Model background expectation is observed. Gluino masses below 860 GeV and squark masses below 1320 GeV are excluded at the 95% confidence level in simplified models containing only squarks of the first two generations, a gluino octet and a massless neutralino, for squark or gluino masses below 2 TeV, respectively. Squarks and gluinos with equal masses below 1410 GeV are excluded. In minimal supergravity/constrained minimal supersymmetric Standard Model models with tan beta = 10, A(0) = 0 and mu > 0, squarks and gluinos of equal mass are excluded for masses below 1360 GeV. Constraints are also placed on the parameter space of supersymmetric models with compressed spectra. These limits considerably extend the region of supersymmetric parameter space excluded by previous measurements with the ATLAS detector.

Search for charginos nearly mass-degenerate with the lightest neutralino based on a disappearing-track signature in pp collisions at sqrts = 8 TeV with the ATLAS detector

A search is presented for direct chargino production based on a disappearing-track signature using 20.3  fb−1 of proton-proton collisions at s√=8  TeV collected with the ATLAS experiment at the LHC. In anomaly-mediated supersymmetry breaking (AMSB) models, the lightest chargino is nearly mass degenerate with the lightest neutralino and its lifetime is long enough to be detected in the tracking detectors by identifying decays that result in tracks with no associated hits in the outer region of the tracking system. Some models with supersymmetry also predict charginos with a significant lifetime. This analysis attains sensitivity for charginos with a lifetime between 0.1 and 10 ns, and significantly surpasses the reach of the LEP experiments. No significant excess above the background expectation is observed for candidate tracks with large transverse momentum, and constraints on chargino properties are obtained. In the AMSB scenarios, a chargino mass below 270 GeV is excluded at 95% confidence level.