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. ^

The role of smoking status on postoperative complications and survival in lung cancer patients presenting to M.D. Anderson Cancer Center

It is estimated that 50% of all lung cancer patients continue to smoke after diagnosis. Many of these lung cancer patients who are current smokers often experience tremendous guilt and responsibility for their disease, and feel it might be too late for them to quit smoking. In addition, many oncologists may be heard to say that it is 'too late', 'it doesn't matter', 'it is too difficult', 'it is too stressful' for their patients to stop smoking, or they never identify the smoking status of the patient. Many oncologists feel unprepared to address smoking cessation as part of their clinical practice. In reality, physicians can have tremendous effects on motivating patients, particularly when patients are initially being diagnosed with cancer. More information is needed to convince patients to quit smoking and to encourage clinicians to assist patients with their smoking cessation. ^ In this current study, smoking status at time of lung cancer diagnosis was assessed to examine its impact on complications and survival, after exploring the reliability of smoking data that is self-reported. Logistic Regression was used to determine the risks of smoking prior to lung resection. In addition, survival analysis was performed to examine the impact of smoking on survival. ^ The reliability of how patients report their smoking status was high, but there was some discordance between current smokers and recent quitters. In addition, we found that cigarette pack-year history and duration of smoking cessation were directly related to the rate of a pulmonary complication. In regards to survival, we found that current smoking at time of lung cancer diagnosis was an independent predictor of early stage lung cancer. This evidence supports the idea that it is "never too late" for patients to quit smoking and health care providers should incorporate smoking status regularly into their clinical practice.^

Modeling survival and the development of second primary tumors among Hodgkin's disease patients

Hodgkin's disease (HD) is a cancer of the lymphatic system. Survivors of HD face varieties of consequent adverse effects, in which secondary primary tumors (SPT) is one of the most serious consequences. This dissertation is aimed to model time-to-SPT in the presence of death and HD relapses during follow-up.^ The model is designed to handle a mixture phenomenon of SPT and the influence of death. Relapses of HD are adjusted as a covariate. Proportional hazards framework is used to define SPT intensity function, which includes an exponential term to estimate explanatory variables. Death as a competing risk is considered according to different scenarios, depending on which terminal event comes first. Newton-Raphson method is used to estimate the parameter estimates in the end.^ The proposed method is applied to a real data set containing a group of HD patients. Several risk factors for the development of SPT are identified and the findings are noteworthy in the development of healthcare guidelines that may lead to the early detection or prevention of SPT.^

Racial disparities in patient survival after liver transplantation in the United States

Background. Racial disparities in healthcare span such areas as access, outcomes after procedures, and patient satisfaction. Previous work suggested that minorities experience less healthcare and worse survival rates. In adult orthotopic liver transplantation (OLT) mixed results have been reported, with some showing African-American recipients having poor survival compared to Caucasians, and others finding no such discrepancy. ^ Purpose. This study’s purpose was to analyze the most recent United Network for Organ Sharing (UNOS) data, both before and after the implementation of the Model for End-Stage Liver Disease (MELD)/Pediatric End-Stage Liver Disease (PELD) scoring system, to determine if minority racial groups still experience poor outcomes after OLT. ^ Methods. The UNOS dataset for 1992-2001 (Era I) and 2002-2007 (Era II) was used. Patient survival rates for each Era and for adult and pediatric recipients were analyzed with adjustment. A separate multivariate analysis was performed on African-American adult patients in Era II in order to identify unique predictors for poor patient survival. ^ Results. The overall study included 66,118 OLT recipients. The majority were Caucasian (78%), followed by Hispanics (13%) and African-Americans (9%). Hispanic and African-American adults were more likely to be female, have Hepatitis C, to be in the intensive care unit (ICU) or ventilated at time of OLT, to have a MELD score ≥23, to have a lower education level, and to have public insurance when compared to Caucasian adults (all p-values < 0.05). Hispanic and African-American pediatric recipients were more likely have public insurance and less likely to receive a living donor OLT than were Caucasian pediatric OLT recipients (p <0.05). There was no difference in the likelihood of having a PELD score ≥21 among racial groups (p >0.40). African-American adults in Era I and Era II had worse patient survival rates than both Caucasians and Hispanic (pair-wise p-values <0.05). This same disparity was seen for pediatric recipients in Era I, but not in Era II. Multivariate analysis of African-American recipients revealed no unique predictors of patient death. ^ Conclusions. African-American race is still a predictor of poor outcome after adult OLT, even after adjustment for multiple clinical, demographic, and liver disease severity variables. Although African-American and Hispanic subgroups share many characteristics previously thought to increase risk of post-OLT death, only African-American patients have poor survival rates when compared to Caucasians. ^

Ethnic differences in tumor characteristics, treatment and survival in glioma patients

Purpose. A descriptive analysis of glioma patients by race was carried out in order to better elucidate potential differences between races in demographics, treatment, characteristics, prognosis and survival. ^ Patients and Methods. Among 1,967 patients ≥ 18 years diagnosed with glioma seen between July 2000 and September 2006 at The University of Texas M.D. Anderson Cancer Center (UTMDACC). Data were collated from the UTMDACC Patient History Database (PHDB) and the UTMDACC Tumor Registry Database (TRDB). Chi-square analysis, uni- /multivariate Cox proportional hazards modeling and survival analysis were used to analyze differences by race. ^ Results. Demographic, treatment and histologic differences exist between races. Though risk differences were seen between races, race was not found to be a significant predictor in multivariate regression analysis after accounting for age, surgery, chemotherapy, radiation, tumor type as stratified by WHO tumor grade. Age was the most consistent predictor in risk for death. Overall survival by race was significantly different (p=0.0049) only in low-grade gliomas after adjustment for age although survival differences were very slight. ^ Conclusion. Among this cohort of glioma patients, age was the strongest predictor for survival. It is likely that survival is more influenced by age, time to treatment, tumor grade and surgical expertise rather than racial differences. However, age at diagnosis, gender ratios, histology and history of cancer differed significantly between race and genetic differences to this effect cannot be excluded. ^

The role of cardiovascular comorbidities in ovarian cancer survival

Objective. One facet of cancer care that often goes ignored is comorbidities, or diseases that exist in concert with cancer. Comorbid conditions may affect survival by influencing treatment decisions and prognosis. The purpose of this secondary data analysis was to identify whether a history of cardiovascular comorbidities among ovarian cancer patients influenced survival time at the University of Texas M. D. Anderson Cancer Center. The parent study, Project Peace, has a longitudinal design with an embedded randomized efficacy study which seeks to improve detection of depressive disorders in ovarian, peritoneal, and fallopian tube cancers. ^ Methods. Survival time was calculated for the 249 ovarian cancer patients abstracted by Project Peace staff. Cardiovascular comorbidities were documented as present, based upon information from medical records in addition to self reported comorbidities in a baseline study questionnaire. Kaplan-Meier survival curves were used to compare survival time among patients with a presence or absence of particular cardiovascular comorbidities. Cox Regression proportional models accounted for multivariable factors such as age, staging, family history of cardiovascular comorbidities, and treatment. ^ Results. Among our patient population, there was a statistically significant relationship between shorter survival time and a history of thrombosis, pericardial disease/tamponade, or COPD/pulmonary hypertension. Ovarian cancer patients with a history of thrombosis lived approximately half as long as patients without thrombosis (58.06 months vs. 121.55 months; p=.001). In addition, patients who suffered from pericardial disease/tamponade had poorer survival than those without a history of pericardial disease/tamponade (48 months vs. 80.07 months; p=.002). Ovarian cancer patients with a history of COPD or pulmonary hypertension had a median survival of 60.2 months, while the median survival for patients without these comorbidities was 80.2 months (p=.014). ^ Conclusion. Especially because of its relatively lower survival rate, greater emphasis needs to be placed on the potential influence of cardiovascular comorbid conditions in ovarian cancer.^

Hormone therapy and physical activity as predictors of lung cancer risk and survival: The California Teachers Study

Cigarette smoking is responsible for the majority of lung cancer cases worldwide; however, a proportion of never smokers still develop lung cancer over their lifetime, prompting investigation into additional factors that may modify lung cancer incidence, as well as mortality. Although hormone therapy (HT), physical activity (PA), and lung cancer have been previously examined, the associations remain unclear. This study investigated exposure to HT and PA that may modulate underlying mechanisms of lung cancer etiology and progression among women by using existing, de-identified data from the California Teachers Study (CTS).^ The CTS cohort, established in 1995–1996, has 133,479 active and retired female teachers and administrators, recruited through the California State Teachers Retirement System, and followed annually for cancer diagnosis, death, and change of address. Each woman enrolled in the CTS returned a questionnaire covering a wide variety of issues related to cancer risk and women's health, including recent and past HT use and physical activity, as well as active and environmental cigarette smoke exposure. Complete data to assess the associations between HT and lung cancer risk and survival were available for 60,592 postmenopausal women. Between 1995 and 2007, 727 of these women were diagnosed with invasive lung cancer; 441 of these died. Complete data to assess the associations between PA and lung cancer risk and survival were available for 118,513 women. Between 1995 and 2007, 853 of these women were diagnosed with invasive lung cancer; 516 of these died.^ After careful adjustment for smoking habits and other potential confounders, no measure of HT use was associated with lung cancer risk; however, any HT use (vs. no use) was associated with a decrease in lung-cancer-specific mortality. Specifically, among women who only used estrogen (E-only), decreases in lung cancer mortality were seen for recent use, but not for former use; no association was observed for estrogen plus progestin (E+P). Furthermore, among former users of HT, a statistically significant decrease in lung cancer mortality was observed for E-only use within 5 years prior to baseline, but not for E-only use >5 years prior to baseline. Neither long-term recreational PA nor recent recreational PA alone were associated with lung cancer risk; however, among women with a BMI<25 and ever smokers, high long-term moderate+strenuous PA was associated with a decrease in lung cancer risk. Women with non-local disease showed a decrease in lung cancer mortality associated with increasing duration of strenuous long-term activity, and 1.50-3.00 h/wk/y of recent moderate or recent strenuous PA. Long-term moderate PA was associated with decreased lung cancer mortality in never smokers, whereas recent moderate PA was associated with increased lung cancer mortality in current smokers. ^ Placing our findings in the context of the current literature, HT does not appear to be associated with lung cancer risk and previous studies reporting a protective effect of HT use on lung cancer risk may be subject to residual confounding by smoking. Looking at our findings regarding PA overall, the evidence still remains inconclusive regarding whether or not physical activity influence lung cancer risk or mortality. Our results suggest that recreational PA may associated with decreased lung cancer risk among women with BMI<25 and ever smoking-women; however, residual confounding by smoking should be strongly considered. To our knowledge, this is the first study to investigate lifetime recreational PA and lung cancer mortality among women. Our results contribute to the growing body of knowledge regarding non-smoking-related risk factors for lung cancer incidence and mortality among women. Given the potential clinical and interventional significance, further study and validation of these findings is warranted.^

Demographic analysis of Certified Industrial Hygienists (CIH) with the use of survival analysis techniques

Objective. The goal of this study is to characterize the current workforce of CIHs, the lengths of professional practice careers of the past and current CIHs.^ Methods. This is a secondary data analysis of data compiled from all of the nearly 50 annual roster listings of the American Board of Industrial Hygiene (ABIH) for Certified Industrial Hygienists active in each year since 1960. Survival analysis was performed as a technique to measure the primary outcome of interest. The technique which was involved in this study was the Kaplan-Meier method for estimating the survival function.^ Study subjects: The population to be studied is all Certified Industrial Hygienists (CIHs). A CIH is defined by the ABIH as an individual who has achieved the minimum requirements for education, working experience and through examination, has demonstrated a minimum level of knowledge and competency in the prevention of occupational illnesses. ^ Results. A Cox-proportional hazards model analysis was performed by different start-time cohorts of CIHs. In this model we chose cohort 1 as the reference cohort. The estimated relative risk of the event (defined as retirement, or absent from 5 consecutive years of listing) occurred for CIHs for cohorts 2,3,4,5 relative to cohort 1 is 0.385, 0.214, 0.234, 0.299 relatively. The result show that cohort 2 (CIHs issued from 1970-1980) has the lowest hazard ratio which indicates the lowest retirement rate.^ Conclusion. The manpower of CIHs (still actively practicing up to the end of 2009) increased tremendously starting in 1980 and grew into a plateau in recent decades. This indicates that the supply and demand of the profession may have reached equilibrium. More demographic information and variables are needed to actually predict the future number of CIHs needed. ^

Modeling the survival of T-cell lymphocytes using compound branching processes

Radiotherapy has been a method of choice in cancer treatment for a number of years. Mathematical modeling is an important tool in studying the survival behavior of any cell as well as its radiosensitivity. One particular cell under investigation is the normal T-cell, the radiosensitivity of which may be indicative to the patient's tolerance to radiation doses.^ The model derived is a compound branching process with a random initial population of T-cells that is assumed to have compound distribution. T-cells in any generation are assumed to double or die at random lengths of time. This population is assumed to undergo a random number of generations within a period of time. The model is then used to obtain an estimate for the survival probability of T-cells for the data under investigation. This estimate is derived iteratively by applying the likelihood principle. Further assessment of the validity of the model is performed by simulating a number of subjects under this model.^ This study shows that there is a great deal of variation in T-cells survival from one individual to another. These variations can be observed under normal conditions as well as under radiotherapy. The findings are in agreement with a recent study and show that genetic diversity plays a role in determining the survival of T-cells. ^

Geographic variation in the utilization and survival associated with oxaliplatin chemotherapy in patients with stage III colon cancer

Background: Overall objectives of this dissertation are to examine the geographic variation and socio-demographic disparities (by age, race and gender) in the utilization and survival of newly FDA-approved chemotherapy agents (Oxaliplatin-containing regimens) as well as to determine the cost-effectiveness of Oxaliplatin in a large nationwide and population-based cohort of Medicare patients with resected stage-III colon cancer. Methods: A retrospective cohort of 7,654 Medicare patients was identified from the Surveillance, Epidemiology and End Results – Medicare linked database. Multiple logistic regression was performed to examine the relationship between receipt of Oxaliplatin-containing chemotherapy and geographic regions while adjusting for other patient characteristics. Cox proportional hazard model was used to estimate the effect of Oxaliplatin-containing chemotherapy on the survival variation across regions using 2004-2005 data. Propensity score adjustments were also made to control for potential bias related to non-random allocation of the treatment group. We used Kaplan-Meier sample average estimator to calculate the cost of disease after cancer-specific surgery to death, loss-to follow-up or censorship. Results: Only 51% of the stage-III patients received adjuvant chemotherapy within three to six months of colon-cancer specific surgery. Patients in the rural regions were approximately 30% less likely to receive Oxaliplatin chemotherapy than those residing in a big metro region (OR=0.69, p=0.033). The hazard ratio for patients residing in metro region was comparable to those residing in big metro region (HR: 1.05, 95% CI: 0.49-2.28). Patients who received Oxalipaltin chemotherapy were 33% less likely to die than those received 5-FU only chemotherapy (adjusted HR=0.67, 95% CI: 0.41-1.11). KMSA-adjusted mean payments were almost 2.5 times higher in the Oxaliplatin-containing group compared to 5-FU only group ($45,378 versus $17,856). When compared to no chemotherapy group, ICER of 5-FU based regimen was $12,767 per LYG, and ICER of Oxaliplatin-chemotherapy was $60,863 per LYG. Oxaliplatin was found economically dominated by 5-FU only chemotherapy in this study population. Conclusion: Chemotherapy use varies across geographic regions. We also observed considerable survival differences across geographic regions; the difference remained even after adjusting for socio-demographic characteristics. The cost-effectiveness of Oxaliplatin in Medicare patients may be over-estimated in the clinical trials. Our study found 5-FU only chemotherapy cost-effective in adjuvant settings in patients with stage-III colon cancer.^

Data quality in trauma transfusion studies and the impact of missing data on predicting massive transfusion

Maximizing data quality may be especially difficult in trauma-related clinical research. Strategies are needed to improve data quality and assess the impact of data quality on clinical predictive models. This study had two objectives. The first was to compare missing data between two multi-center trauma transfusion studies: a retrospective study (RS) using medical chart data with minimal data quality review and the PRospective Observational Multi-center Major Trauma Transfusion (PROMMTT) study with standardized quality assurance. The second objective was to assess the impact of missing data on clinical prediction algorithms by evaluating blood transfusion prediction models using PROMMTT data. RS (2005-06) and PROMMTT (2009-10) investigated trauma patients receiving ≥ 1 unit of red blood cells (RBC) from ten Level I trauma centers. Missing data were compared for 33 variables collected in both studies using mixed effects logistic regression (including random intercepts for study site). Massive transfusion (MT) patients received ≥ 10 RBC units within 24h of admission. Correct classification percentages for three MT prediction models were evaluated using complete case analysis and multiple imputation based on the multivariate normal distribution. A sensitivity analysis for missing data was conducted to estimate the upper and lower bounds of correct classification using assumptions about missing data under best and worst case scenarios. Most variables (17/33=52%) had <1% missing data in RS and PROMMTT. Of the remaining variables, 50% demonstrated less missingness in PROMMTT, 25% had less missingness in RS, and 25% were similar between studies. Missing percentages for MT prediction variables in PROMMTT ranged from 2.2% (heart rate) to 45% (respiratory rate). For variables missing >1%, study site was associated with missingness (all p≤0.021). Survival time predicted missingness for 50% of RS and 60% of PROMMTT variables. MT models complete case proportions ranged from 41% to 88%. Complete case analysis and multiple imputation demonstrated similar correct classification results. Sensitivity analysis upper-lower bound ranges for the three MT models were 59-63%, 36-46%, and 46-58%. Prospective collection of ten-fold more variables with data quality assurance reduced overall missing data. Study site and patient survival were associated with missingness, suggesting that data were not missing completely at random, and complete case analysis may lead to biased results. Evaluating clinical prediction model accuracy may be misleading in the presence of missing data, especially with many predictor variables. The proposed sensitivity analysis estimating correct classification under upper (best case scenario)/lower (worst case scenario) bounds may be more informative than multiple imputation, which provided results similar to complete case analysis.^

Analyzing left-truncated right-censored data with uncertain onset time with parametric models

Prevalent sampling is an efficient and focused approach to the study of the natural history of disease. Right-censored time-to-event data observed from prospective prevalent cohort studies are often subject to left-truncated sampling. Left-truncated samples are not randomly selected from the population of interest and have a selection bias. Extensive studies have focused on estimating the unbiased distribution given left-truncated samples. However, in many applications, the exact date of disease onset was not observed. For example, in an HIV infection study, the exact HIV infection time is not observable. However, it is known that the HIV infection date occurred between two observable dates. Meeting these challenges motivated our study. We propose parametric models to estimate the unbiased distribution of left-truncated, right-censored time-to-event data with uncertain onset times. We first consider data from a length-biased sampling, a specific case in left-truncated samplings. Then we extend the proposed method to general left-truncated sampling. With a parametric model, we construct the full likelihood, given a biased sample with unobservable onset of disease. The parameters are estimated through the maximization of the constructed likelihood by adjusting the selection bias and unobservable exact onset. Simulations are conducted to evaluate the finite sample performance of the proposed methods. We apply the proposed method to an HIV infection study, estimating the unbiased survival function and covariance coefficients. ^

The risk prediction for SPM and its impact on the survival in patients with head and neck squamous cell carcinoma

Head and Neck Squamous Cell Carcinoma (HNSCC) is the sixth common malignancy in the world, with high rates of developing second primary malignancy (SPM) and moderately low survival rates. This disease has become an enormous challenge in the cancer research and treatments. For HNSCC patients, a highly significant cause of post-treatment mortality and morbidity is the development of SPM. Hence, assessment of predicting the risk for the development of SPM would be very helpful for patients, clinicians and policy makers to estimate the survival of patients with HNSCC. In this study, we built a prognostic model to predict the risk of developing SPM in patients with newly diagnosed HNSCC. The dataset used in this research was obtained from The University of Texas MD Anderson Cancer Center. For the first aim, we used stepwise logistic regression to identify the prognostic factors for the development of SPM. Our final model contained cancer site and overall cancer stage as our risk factors for SPM. The Hosmer-Lemeshow test (p-value= 0.15>0.05) showed the final prognostic model fit the data well. The area under the ROC curve was 0.72 that suggested the discrimination ability of our model was acceptable. The internal validation confirmed the prognostic model was a good fit and the final prognostic model would not over optimistically predict the risk of SPM. This model needs external validation by using large data sample size before it can be generalized to predict SPM risk for other HNSCC patients. For the second aim, we utilized a multistate survival analysis approach to estimate the probability of death for HNSCC patients taking into consideration of the possibility of SPM. Patients without SPM were associated with longer survival. These findings suggest that the development of SPM could be a predictor of survival rates among the patients with HNSCC.^

Assessing the patterns of colorectal cancer treatment, adherence to treatment guidelines and survival in Texas

Colorectal cancer (CRC) is the third leading cancer in both incidence and mortality in Texas. This study investigated the adherence of CRC treatment to standard treatment guidelines and the association between standard treatment and CRC survival in Texas. The author used Texas Cancer Registry (TCR) and Medicare linked data to study the CRC treatment patterns and factors associated with standard treatment in patients who were more than 65 years old and were diagnosed in 2001 through 2007. We also determined whether adherence to standard treatment affect patients' survival. Multiple logistic regression and Cox regression analysis were used to analyze our data. Both regression models are adjusted for demographic characteristics and tumor characteristics. We found that for the 3977 regional colon cancer patients 80 years old or younger, 60.2% of them received chemotherapy, in adherence to the recommended treatment guidelines. People with younger age, female gender, higher education and lower comorbidity score are more likely adherent to this surgery guideline. Patients' adherence to chemotherapy in this cohort have better survival compared to those who are not (HR: 0.76, 95% CI: 0.68-0.84). For the 12709 colon cancer patients treated with surgery, 49.3% have more than 12 lymph nodes removed, in adherence to the treatment guidelines. People with younger age, female gender, higher education, regional stage, lager tumor size and lower comorbidity score are more likely to adherent to this surgery guideline. Patients with more than 12 lymph nodes removed in this cohort have better survival (HR: 0.86, 95% CI: 0.82-0.91). For the 1211 regional rectal cancer patients 80 years old or younger, 63.2% of them were adherent to radiation treatment. People with smaller tumor size and lower comorbidity score are more likely to adherent to this radiation guideline. There is no significant survival difference between radiation adherent patients and non-adherent patients (HR: 1.03, 95% CI: 0.82-1.29). For the 1122 regional rectal cancer patients 80 years old or younger who were treated with surgery, 76.0% of them received postoperative chemotherapy, in adherence to the treatment guidelines. People with younger age and smaller comorbidity score are related with higher adherence rate. Patients adherent with adjuvant chemotherapy in this cohort have better survival than those were not adherent (HR: 0.60, 95% CI: 0.45-0.79).^