Simulation study of joint transition and Weibull survival model with shared parameters

This study investigates a theoretical model where a longitudinal process, that is a stationary Markov-Chain, and a Weibull survival process share a bivariate random effect. Furthermore, a Quality-of-Life adjusted survival is calculated as the weighted sum of survival time. Theoretical values of population mean adjusted survival of the described model are computed numerically. The parameters of the bivariate random effect do significantly affect theoretical values of population mean. Maximum-Likelihood and Bayesian methods are applied on simulated data to estimate the model parameters. Based on the parameter estimates, predicated population mean adjusted survival can then be calculated numerically and compared with the theoretical values. Bayesian method and Maximum-Likelihood method provide parameter estimations and population mean prediction with comparable accuracy; however Bayesian method suffers from poor convergence due to autocorrelation and inter-variable correlation. ^

Bayesian and survival model for tumor relapse status and disease-specific survival, with applications to breast cancer

Breast cancer is the most common non-skin cancer and the second leading cause of cancer-related death in women in the United States. Studies on ipsilateral breast tumor relapse (IBTR) status and disease-specific survival will help guide clinic treatment and predict patient prognosis.^ After breast conservation therapy, patients with breast cancer may experience breast tumor relapse. This relapse is classified into two distinct types: true local recurrence (TR) and new ipsilateral primary tumor (NP). However, the methods used to classify the relapse types are imperfect and are prone to misclassification. In addition, some observed survival data (e.g., time to relapse and time from relapse to death)are strongly correlated with relapse types. The first part of this dissertation presents a Bayesian approach to (1) modeling the potentially misclassified relapse status and the correlated survival information, (2) estimating the sensitivity and specificity of the diagnostic methods, and (3) quantify the covariate effects on event probabilities. A shared frailty was used to account for the within-subject correlation between survival times. The inference was conducted using a Bayesian framework via Markov Chain Monte Carlo simulation implemented in softwareWinBUGS. Simulation was used to validate the Bayesian method and assess its frequentist properties. The new model has two important innovations: (1) it utilizes the additional survival times correlated with the relapse status to improve the parameter estimation, and (2) it provides tools to address the correlation between the two diagnostic methods conditional to the true relapse types.^ Prediction of patients at highest risk for IBTR after local excision of ductal carcinoma in situ (DCIS) remains a clinical concern. The goals of the second part of this dissertation were to evaluate a published nomogram from Memorial Sloan-Kettering Cancer Center, to determine the risk of IBTR in patients with DCIS treated with local excision, and to determine whether there is a subset of patients at low risk of IBTR. Patients who had undergone local excision from 1990 through 2007 at MD Anderson Cancer Center with a final diagnosis of DCIS (n=794) were included in this part. Clinicopathologic factors and the performance of the Memorial Sloan-Kettering Cancer Center nomogram for prediction of IBTR were assessed for 734 patients with complete data. Nomogram for prediction of 5- and 10-year IBTR probabilities were found to demonstrate imperfect calibration and discrimination, with an area under the receiver operating characteristic curve of .63 and a concordance index of .63. In conclusion, predictive models for IBTR in DCIS patients treated with local excision are imperfect. Our current ability to accurately predict recurrence based on clinical parameters is limited.^ The American Joint Committee on Cancer (AJCC) staging of breast cancer is widely used to determine prognosis, yet survival within each AJCC stage shows wide variation and remains unpredictable. For the third part of this dissertation, biologic markers were hypothesized to be responsible for some of this variation, and the addition of biologic markers to current AJCC staging were examined for possibly provide improved prognostication. The initial cohort included patients treated with surgery as first intervention at MDACC from 1997 to 2006. Cox proportional hazards models were used to create prognostic scoring systems. AJCC pathologic staging parameters and biologic tumor markers were investigated to devise the scoring systems. Surveillance Epidemiology and End Results (SEER) data was used as the external cohort to validate the scoring systems. Binary indicators for pathologic stage (PS), estrogen receptor status (E), and tumor grade (G) were summed to create PS+EG scoring systems devised to predict 5-year patient outcomes. These scoring systems facilitated separation of the study population into more refined subgroups than the current AJCC staging system. The ability of the PS+EG score to stratify outcomes was confirmed in both internal and external validation cohorts. The current study proposes and validates a new staging system by incorporating tumor grade and ER status into current AJCC staging. We recommend that biologic markers be incorporating into revised versions of the AJCC staging system for patients receiving surgery as the first intervention.^ Chapter 1 focuses on developing a Bayesian method to solve misclassified relapse status and application to breast cancer data. Chapter 2 focuses on evaluation of a breast cancer nomogram for predicting risk of IBTR in patients with DCIS after local excision gives the statement of the problem in the clinical research. Chapter 3 focuses on validation of a novel staging system for disease-specific survival in patients with breast cancer treated with surgery as the first intervention. ^

Mixture modeling for joint analysis of survival, discrete, and continuous data

Mixture modeling is commonly used to model categorical latent variables that represent subpopulations in which population membership is unknown but can be inferred from the data. In relatively recent years, the potential of finite mixture models has been applied in time-to-event data. However, the commonly used survival mixture model assumes that the effects of the covariates involved in failure times differ across latent classes, but the covariate distribution is homogeneous. The aim of this dissertation is to develop a method to examine time-to-event data in the presence of unobserved heterogeneity under a framework of mixture modeling. A joint model is developed to incorporate the latent survival trajectory along with the observed information for the joint analysis of a time-to-event variable, its discrete and continuous covariates, and a latent class variable. It is assumed that the effects of covariates on survival times and the distribution of covariates vary across different latent classes. The unobservable survival trajectories are identified through estimating the probability that a subject belongs to a particular class based on observed information. We applied this method to a Hodgkin lymphoma study with long-term follow-up and observed four distinct latent classes in terms of long-term survival and distributions of prognostic factors. Our results from simulation studies and from the Hodgkin lymphoma study demonstrated the superiority of our joint model compared with the conventional survival model. This flexible inference method provides more accurate estimation and accommodates unobservable heterogeneity among individuals while taking involved interactions between covariates into consideration.^

A multivariate frailty model for disease recurrences and survival

A multivariate frailty hazard model is developed for joint-modeling of three correlated time-to-event outcomes: (1) local recurrence, (2) distant recurrence, and (3) overall survival. The term frailty is introduced to model population heterogeneity. The dependence is modeled by conditioning on a shared frailty that is included in the three hazard functions. Independent variables can be included in the model as covariates. The Markov chain Monte Carlo methods are used to estimate the posterior distributions of model parameters. The algorithm used in present application is the hybrid Metropolis-Hastings algorithm, which simultaneously updates all parameters with evaluations of gradient of log posterior density. The performance of this approach is examined based on simulation studies using Exponential and Weibull distributions. We apply the proposed methods to a study of patients with soft tissue sarcoma, which motivated this research. Our results indicate that patients with chemotherapy had better overall survival with hazard ratio of 0.242 (95% CI: 0.094 - 0.564) and lower risk of distant recurrence with hazard ratio of 0.636 (95% CI: 0.487 - 0.860), but not significantly better in local recurrence with hazard ratio of 0.799 (95% CI: 0.575 - 1.054). The advantages and limitations of the proposed models, and future research directions are discussed. ^

Evaluation and Model of a Clinical Librarian Program

Clinical medical librarianship is entering its second decade, but little evaluative data has accrued in the literature. Variations from the original programs and novel new approaches have insured the survival of the program so far. The clinical librarian (CL) forms a vital link between the library and the health care professional, operating as an important information transfer agent. However, to further insure the survival of these vital programs, hard evaluative evidence is needed. The University of Texas Medical Branch (UTMB) at Galveston began a CL Program in 1978/79. An extensive three-year pre/post evaluation study was conducted using a specifically developed evaluation model, which, if adopted by others, will provide the needed comparative data. Both a pilot study, or formative evaluation, and a summative evaluation were conducted. The results of this evaluation confirmed many of the conclusions reported by other CL programs. Eight hypotheses were proposed at the beginning of this study. Data were collected and used to support acceptance or rejection of the null hypotheses, and conclusions were drawn according to the results. Implications relevant to the study conclusions and future trends in medical librarianship are also discussed in the closing chapter.

SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

Detailed analysis of sequence changes occurring during vlsE antigenic variation in the mouse model of Borrelia burgdorferi infection.

Lyme disease Borrelia can infect humans and animals for months to years, despite the presence of an active host immune response. The vls antigenic variation system, which expresses the surface-exposed lipoprotein VlsE, plays a major role in B. burgdorferi immune evasion. Gene conversion between vls silent cassettes and the vlsE expression site occurs at high frequency during mammalian infection, resulting in sequence variation in the VlsE product. In this study, we examined vlsE sequence variation in B. burgdorferi B31 during mouse infection by analyzing 1,399 clones isolated from bladder, heart, joint, ear, and skin tissues of mice infected for 4 to 365 days. The median number of codon changes increased progressively in C3H/HeN mice from 4 to 28 days post infection, and no clones retained the parental vlsE sequence at 28 days. In contrast, the decrease in the number of clones with the parental vlsE sequence and the increase in the number of sequence changes occurred more gradually in severe combined immunodeficiency (SCID) mice. Clones containing a stop codon were isolated, indicating that continuous expression of full-length VlsE is not required for survival in vivo; also, these clones continued to undergo vlsE recombination. Analysis of clones with apparent single recombination events indicated that recombinations into vlsE are nonselective with regard to the silent cassette utilized, as well as the length and location of the recombination event. Sequence changes as small as one base pair were common. Fifteen percent of recovered vlsE variants contained "template-independent" sequence changes, which clustered in the variable regions of vlsE. We hypothesize that the increased frequency and complexity of vlsE sequence changes observed in clones recovered from immunocompetent mice (as compared with SCID mice) is due to rapid clearance of relatively invariant clones by variable region-specific anti-VlsE antibody responses.

Long-term survival after multidisciplinary management of resected pancreatic adenocarcinoma.

INTRODUCTION: Actual 5-year survival rates of 10-18% have been reported for patients with resected pancreatic adenocarcinoma (PC), but the use of multimodality therapy was uncommon in these series. We evaluated long-term survival and patterns of recurrence in patients treated for PC with contemporary staging and multimodality therapy. METHODS: We analyzed 329 consecutive patients with PC evaluated between 1990 and 2002 who underwent resection. Each received a multidisciplinary evaluation and a standard operative approach. Pre- or postoperative chemotherapy and/or chemoradiation were routine. Surgical specimens of 5-year survivors were re-reviewed. A multivariate model of factors associated with long-term survival was constructed. RESULTS: Patients underwent pancreaticoduodenectomy (n = 302; 92%), distal (n = 20; 6%), or total pancreatectomy (n = 7; 2%). A total of 108 patients (33%) underwent vascular reconstruction, 301 patients (91%) received neoadjuvant or adjuvant therapy, 157 specimens (48%) were node positive, and margins were microscopically positive in 52 patients (16%). Median overall survival and disease-specific survival was 23.9 and 26.5 months. Eighty-eight patients (27%) survived a minimum of 5 years and had a median overall survival of 11 years. Of these, 21 (24%) experienced recurrence, 7 (8%) after 5 years. Late recurrences occurred most frequently in the lungs, the latest at 6.7 years. Multivariate analysis identified disease-negative lymph nodes (P = .02) and no prior attempt at resection (P = 0.01) as associated with 5-year survival. CONCLUSIONS: Our 27% actual 5-year survival rate for patients with resected PC is superior to that previously reported, and it is influenced by our emphasis on detailed staging and patient selection, a standardized operative approach, and routine use of multimodality therapy.

Detailed analysis of sequence changes occurring during vlsE antigenic variation in the mouse model of Borrelia burgdorferi infection.

Lyme disease Borrelia can infect humans and animals for months to years, despite the presence of an active host immune response. The vls antigenic variation system, which expresses the surface-exposed lipoprotein VlsE, plays a major role in B. burgdorferi immune evasion. Gene conversion between vls silent cassettes and the vlsE expression site occurs at high frequency during mammalian infection, resulting in sequence variation in the VlsE product. In this study, we examined vlsE sequence variation in B. burgdorferi B31 during mouse infection by analyzing 1,399 clones isolated from bladder, heart, joint, ear, and skin tissues of mice infected for 4 to 365 days. The median number of codon changes increased progressively in C3H/HeN mice from 4 to 28 days post infection, and no clones retained the parental vlsE sequence at 28 days. In contrast, the decrease in the number of clones with the parental vlsE sequence and the increase in the number of sequence changes occurred more gradually in severe combined immunodeficiency (SCID) mice. Clones containing a stop codon were isolated, indicating that continuous expression of full-length VlsE is not required for survival in vivo; also, these clones continued to undergo vlsE recombination. Analysis of clones with apparent single recombination events indicated that recombinations into vlsE are nonselective with regard to the silent cassette utilized, as well as the length and location of the recombination event. Sequence changes as small as one base pair were common. Fifteen percent of recovered vlsE variants contained "template-independent" sequence changes, which clustered in the variable regions of vlsE. We hypothesize that the increased frequency and complexity of vlsE sequence changes observed in clones recovered from immunocompetent mice (as compared with SCID mice) is due to rapid clearance of relatively invariant clones by variable region-specific anti-VlsE antibody responses.

Hepatocellular carcinoma survival in uninsured and underinsured patients.

BACKGROUND: The incidence of hepatitis C virus (HCV) and hepatocellular carcinoma (HCC) is increasing. The purpose of this study is to establish baseline survival in a medically-underserved population and to evaluate the effect of HCV seropositivity on our patient population. MATERIALS AND METHODS: We reviewed clinicopathologic parameters from a prospective tumor registry and medical records from the Harris County Hospital District (HCHD). Outcomes were compared using Kaplan-Meier survival analysis and log-rank tests. RESULTS: A total of 298 HCC patients were identified. The median survival for the entire cohort was 3.4 mo. There was no difference in survival between the HCV seropositive and the HCV seronegative groups (3.6 mo versus 2.6 mo, P = 0.7). Patients with a survival <1 mo had a significant increase in>αfetoprotein (AFP), international normalized ratio (INR), model for end-stage liver disease (MELD) score, and total bilirubin and decrease in albumin compared with patients with a survival ≥ 1 mo. CONCLUSIONS: Survival for HCC patients in the HCHD is extremely poor compared with an anticipated median survival of 7 mo reported in other studies. HCV seropositive patients have no survival advantage over HCV seronegative patients. Poorer liver function at diagnosis appears to be related to shorter survival. Further analysis into variables contributing to decreased survival is needed.

Health-related quality of life and survival in liver transplant candidates.

Health-related quality of life (HRQOL) is an important measure of the effects of chronic liver disease in affected patients that helps guide interventions to improve well-being. However, the relationship between HRQOL and survival in liver transplant candidates remains unclear. We examined whether the Physical Component Summary (PCS) and Mental Component Summary (MCS) scores from the Short Form 36 (SF-36) Health Survey were associated with survival in liver transplant candidates. We administered the SF-36 questionnaire (version 2.0) to patients in the Pulmonary Vascular Complications of Liver Disease study, a multicenter prospective cohort of patients evaluated for liver transplantation in 7 academic centers in the United States between 2003 and 2006. Cox proportional hazards models were used with death as the primary outcome and adjustment for liver transplantation as a time-varying covariate. The mean age of the 252 participants was 54 +/- 10 years, 64% were male, and 94% were white. During the 422 person years of follow-up, 147 patients (58%) were listed, 75 patients (30%) underwent transplantation, 49 patients (19%) died, and 3 patients were lost to follow-up. Lower baseline PCS scores were associated with an increased mortality rate despite adjustments for age, gender, Model for End-Stage Liver Disease score, and liver transplantation (P for the trend = 0.0001). The MCS score was not associated with mortality (P for the trend = 0.53). In conclusion, PCS significantly predicts survival in liver transplant candidates, and interventions directed toward improving the physical status may be helpful in improving outcomes in liver transplant candidates.

Activity of anidulafungin in a murine model of Candida krusei infection: evaluation of mortality and disease burden by quantitative tissue cultures and measurement of serum (1,3)-beta-D-glucan levels.

Experience with anidulafungin against Candida krusei is limited. Immunosuppressed mice were injected with 1.3 x 10(7) to 1.5 x 10(7) CFU of C. krusei. Animals were treated with saline, 40 mg/kg fluconazole, 1 mg/kg amphotericin B, or 10 and 20 mg/kg anidulafungin for 5 days. Anidulafungin improved survival and significantly reduced the number of CFU/g in kidneys and serum beta-glucan levels.

Characterization of recombinant class I MHC alloantigen action upon the survival of disparate heart allografts

Previous studies have led to the development of allochimeric class I MHC proteins as agents that effectively induce donor-specific transplantation tolerance in a rat system with or without additional immunosuppression. Within the α1-helical region of RT1.Au, an epitope that conferred immunologic tolerance was discovered. Studies presented herein were designed to test our central hypothesis that allochimeric proteins onfer tolerance in a mouse model. To test this hypothesis, portal vein (PV) injection of wild-type H2Kd and H2Dd proteins were produced in a bacterial expression system and found to specifically prolong the survival of BALB/c (H2d) heart allografts in C57BL/10 (H2b) recipients. Although a single PV injection of 50 μg α1–α 3 H2Kd alone was ineffective, 50 μg α1 –α3 alone slightly prolonged BALB/c heart allograft survivals. In contrast, the combination of 25 μg α1–α 3 H2Kd and 25 μg α1–α 3 H2Dd proteins prolonged BALB/c graft survivals to 20.2 ± 6.4 days (p < 0.004). The effect was donor-specific, since a combination of 25 μg α1–α3 H2Kd and 25 μg α1–α3 H2Dd proteins failed to affect survivals of third-party C3H (H2k k) heart allografts, namely 9.0 ± 0.0 days in treated versus 7.8 ± 0.5 days in untreated hosts. Thus, the combination of two H2K d and H2Dd proteins is more effective in prolonging allograft survival than a single protein produced in a bacterial expression system. A single PV injection (day 0) of 25 μg α1–α 2 H2Kd and 25 μg α1–α 2 H2Dd proteins to C57BL/10 mice prolonged the survival of BALB/c heart allografts to 22.4 ± 4.5 days. Within a WF to ACI rat heart allograft system, a single PV injection of 20 μg 70–77 u-RT1.Aa induced specific tolerance of allografts. This therapy could be combined with CsA to induce transplantation tolerance. However, combination of 70–77u-RT1.Aa with CTLA4Ig, rapamycin, or AG-490 effectively blocked the induction of transplantation tolerance. Tolerance generated by allochimeric protein could be adoptively transferred to naive recipients. Intragraft cytokine mRNA levels showed a bias towards a Th2-type phenotype. Additionally, studies of cytokine signaling and activation of transcription factors revealed a requirement that these pathways remain available for signaling in order for transplantation tolerance to occur. These studies suggest that the generation of regulatory cells are required for the induction of transplantation tolerance through the use of allochimeric proteins. ^

Comparison of a traditional and a multilevel Cox Proportional Hazards model

Hierarchically clustered populations are often encountered in public health research, but the traditional methods used in analyzing this type of data are not always adequate. In the case of survival time data, more appropriate methods have only begun to surface in the last couple of decades. Such methods include multilevel statistical techniques which, although more complicated to implement than traditional methods, are more appropriate. ^ One population that is known to exhibit a hierarchical structure is that of patients who utilize the health care system of the Department of Veterans Affairs where patients are grouped not only by hospital, but also by geographic network (VISN). This project analyzes survival time data sets housed at the Houston Veterans Affairs Medical Center Research Department using two different Cox Proportional Hazards regression models, a traditional model and a multilevel model. VISNs that exhibit significantly higher or lower survival rates than the rest are identified separately for each model. ^ In this particular case, although there are differences in the results of the two models, it is not enough to warrant using the more complex multilevel technique. This is shown by the small estimates of variance associated with levels two and three in the multilevel Cox analysis. Much of the differences that are exhibited in identification of VISNs with high or low survival rates is attributable to computer hardware difficulties rather than to any significant improvements in the model. ^