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

A full pedigree based method for the statistical assessment of genetic anticipation

Genetic anticipation is defined as a decrease in age of onset or increase in severity as the disorder is transmitted through subsequent generations. Anticipation has been noted in the literature for over a century. Recently, anticipation in several diseases including Huntington's Disease, Myotonic Dystrophy and Fragile X Syndrome were shown to be caused by expansion of triplet repeats. Anticipation effects have also been observed in numerous mental disorders (e.g. Schizophrenia, Bipolar Disorder), cancers (Li-Fraumeni Syndrome, Leukemia) and other complex diseases. ^ Several statistical methods have been applied to determine whether anticipation is a true phenomenon in a particular disorder, including standard statistical tests and newly developed affected parent/affected child pair methods. These methods have been shown to be inappropriate for assessing anticipation for a variety of reasons, including familial correlation and low power. Therefore, we have developed family-based likelihood modeling approaches to model the underlying transmission of the disease gene and penetrance function and hence detect anticipation. These methods can be applied in extended families, thus improving the power to detect anticipation compared with existing methods based only upon parents and children. The first method we have proposed is based on the regressive logistic hazard model. This approach models anticipation by a generational covariate. The second method allows alleles to mutate as they are transmitted from parents to offspring and is appropriate for modeling the known triplet repeat diseases in which the disease alleles can become more deleterious as they are transmitted across generations. ^ To evaluate the new methods, we performed extensive simulation studies for data simulated under different conditions to evaluate the effectiveness of the algorithms to detect genetic anticipation. Results from analysis by the first method yielded empirical power greater than 87% based on the 5% type I error critical value identified in each simulation depending on the method of data generation and current age criteria. Analysis by the second method was not possible due to the current formulation of the software. The application of this method to Huntington's Disease and Li-Fraumeni Syndrome data sets revealed evidence for a generation effect in both cases. ^

DNA nucleotide excision repair gene single nucleotide polymorphisms and hereditary nonpolyposis colorectal cancer

Hereditary nonpolyposis colorectal cancer (HNPCC) is an autosomal dominant disease caused by germline mutations in DNA mismatch repair(MMR) genes. The nucleotide excision repair(NER) pathway plays a very important role in cancer development. We systematically studied interactions between NER and MMR genes to identify NER gene single nucleotide polymorphism (SNP) risk factors that modify the effect of MMR mutations on risk for cancer in HNPCC. We analyzed data from polymorphisms in 10 NER genes that had been genotyped in HNPCC patients that carry MSH2 and MLH1 gene mutations. The influence of the NER gene SNPs on time to onset of colorectal cancer (CRC) was assessed using survival analysis and a semiparametric proportional hazard model. We found the median age of onset for CRC among MMR mutation carriers with the ERCC1 mutation was 3.9 years earlier than patients with wildtype ERCC1(median 47.7 vs 51.6, log-rank test p=0.035). The influence of Rad23B A249V SNP on age of onset of HNPCC is age dependent (likelihood ratio test p=0.0056). Interestingly, using the likelihood ratio test, we also found evidence of genetic interactions between the MMR gene mutations and SNPs in ERCC1 gene(C8092A) and XPG/ERCC5 gene(D1104H) with p-values of 0.004 and 0.042, respectively. An assessment using tree structured survival analysis (TSSA) showed distinct gene interactions in MLH1 mutation carriers and MSH2 mutation carriers. ERCC1 SNP genotypes greatly modified the age onset of HNPCC in MSH2 mutation carriers, while no effect was detected in MLH1 mutation carriers. Given the NER genes in this study play different roles in NER pathway, they may have distinct influences on the development of HNPCC. The findings of this study are very important for elucidation of the molecular mechanism of colon cancer development and for understanding why some mutation carriers of the MSH2 and MLH1 gene develop CRC early and others never develop CRC. Overall, the findings also have important implications for the development of early detection strategies and prevention as well as understanding the mechanism of colorectal carcinogenesis in HNPCC. ^

Statins and myasthenia gravis

Introduction. 3-hydroxy-3-methylglutaryl CoA reductase inhibitor ("statin") have been widely used for hypercholesteroremia and Statin induced myopathy is well known. Whether Statins contribute to exacerbation of Myasthenia Gravis (MG) requiring hospitalization is not well known. ^ Objectives. To determine the frequency of statin use in patients with MG seen at the neuromuscular division at University of Alabama in Birmingham (UAB) and to evaluate any association between use of statins and MG exacerbations requiring hospitalization in patients with an established diagnosis of Myasthenia Gravis. ^ Methods. We reviewed records of all current MG patients at the UAB neuromuscular department to obtain details on use of statins and any hospitalizations due to exacerbation of MG over the period from January 1, 2003 to December 31, 2006. ^ Results. Of the 113 MG patients on whom information was available for this period, 40 were on statins during at least one clinic visit. Statin users were more likely to be older (mean age 60.2 vs 53.8, p = 0.029), male (70.0% vs 43.8%, p = 0.008), and had a later onset of myasthenia gravis (mean age in years at onset 49.8 versus 42.9, p = 0.051). The total number of hospitalizations or the proportion of subjects who had at least one hospitalization during the study period did not differ in the statin versus no-statin group. However, when hospitalizations which occurred from a suspected precipitant were excluded ("event"), the proportion of subjects who had at least one such event during the study period was higher in the group using statins. In the final Cox proportional hazard model for cumulative time to event, statin use (OR = 6.44, p <0.01) and baseline immunosuppression (OR = 3.03, p = 0.07) were found to increase the odds of event. ^ Conclusions. Statin use may increase the rate of hospitalizations due to MG exacerbation, when excluding exacerbations precipitated by other suspected factors.^

Role of stress myocardial perfusion scan in pre-operative evaluation of cancer patients undergoing non-cardiac surgery

Background. Cardiac risk assessment in cancer patients has not extensively been studied. We evaluated the role of stress myocardial perfusion imaging (MPI) in predicting cardiovascular outcomes in cancer patients undergoing non-cardiac surgery. ^ Methods. A retrospective chart review was performed on 507 patients who had a MPI from 01/2002 - 03/2003 and underwent non-cardiac surgery. Median follow-up duration was 1.5 years. Cox proportional hazard model was used to determine the time-to-first event. End points included total cardiac events (cardiac death, myocardial infarction (MI) and coronary revascularization), cardiac death, and all cause mortality. ^ Results. Of all 507 MPI studies 146 (29%) were abnormal. There were significant differences in risk factors between normal and abnormal MPI groups. Mean age was 66±11 years, with 60% males and a median follow-up duration of 1.8 years (25th quartile=0.8 years, 75th quartile=2.2 years). The majority of patients had an adenosine stress study (53%), with fewer exercise (28%) and dobutamine stress (16%) studies. In the total group there were 39 total cardiac events, 31 cardiac deaths, and 223 all cause mortality events during the study. Univariate predictors of total cardiac events included CAD (p=0.005), previous MI (p=0.005), use of beta blockers (p=0.002), and not receiving chemotherapy (p=0.012). Similarly, the univariate predictors of cardiac death included previous MI (p=0.019) and use of beta blockers (p=0.003). In the multivariate model for total cardiac events, age at surgery (HR 1.04, p=0.030), use of beta blockers (HR 2.46; p=0.011), dobutamine MPI (HR 3.08; p=0.018) and low EF (HR 0.97; p=0.02) were significant predictors of worse outcomes. In the multivariate model for predictors of cardiac death, beta blocker use (HR=2.74; p=0.017) and low EF (HR=0.95; p<0.003) were predictors of cardiac death. The only univariate MPI predictor of total cardiac events was scar severity (p=0.005). While MPI predictors of cardiac death were scar severity (p= 0.001) and ischemia severity (p=0.02). ^ Conclusions. Stress MPI is a useful tool in predicting long term outcomes in cancer patients undergoing surgery. Ejection fraction and severity of myocardial scar are important factors determining long term outcomes in this group.^

Natural history of placebo treatment in adults acquiring travelers' diarrhea in Mexico

The natural history of placebo treated travelers' diarrhea and the prognostic factors of recovery from diarrhea were evaluated using 9 groups of placebo treated subjects from 9 clinical trial studies conducted since 1975, for use as a historical control in the future clinical trial of antidiarrheal agents. All of these studies were done by the same group of investigators in one site (Guadalajara, Mexico). The studies are similar in terms of population, measured parameters, microbiologic identification of enteropathogens and definitions of parameters. The studies had two different durations of followup. In some studies, subjects were followed for two days, and in some they were followed for five days.^ Using definitions established by the Infectious Diseases society of America and the Food and Drug Administration, the following efficacy parameters were evaluated: Time to last unformed stool (TLUS), number of unformed stools post-initiation of placebo treatment for five consecutive days of followup, microbiologic cure, and improvement of diarrhea. Among the groups that were followed for five days, the mean TLUS ranged from 59.1 to 83.5 hours. Fifty percent to 78% had diarrhea lasting more than 48 hours and 25% had diarrhea more than five days. The mean number of unformed stools passed on the first day post-initiation of therapy ranged from 3.6 to 5.8 and, for the fifth day ranged from 0.5 to 1.5. By the end of followup, diarrhea improved in 82.6% to 90% of the subjects. Subjects with enterotoxigenic E. coli had 21.6% to 90.0% microbiologic cure; and subjects with shigella species experienced 14.3% to 60.0% microbiologic cure.^ In evaluating the prognostic factors of recovery from diarrhea (primary efficacy parameter in evaluating the efficacy of antidiarrheal agents against travelers' diarrhea). The subjects from five studies were pooled and the Cox proportional hazard model was used to evaluate the predictors of prolonged diarrhea. After adjusting for design characteristics of each trial, fever with a rate ratio (RR) of 0.40, presence of invasive pathogens with a RR of 0.41, presence of severe abdominal pain and cramps with a RR of 0.50, number of watery stools more than five with a RR of 0.60, and presence of non-invasive pathogens with a RR of 0.84 predicted a longer duration of diarrhea. Severe vomiting with a RR of 2.53 predicted a shorter duration of diarrhea. The number of soft stools, presence of fecal leukocytes, presence of nausea, and duration of diarrhea before enrollment were not associated with duration of diarrhea. ^

Predictors of Contralateral Breast Cancer in BRCA Negative Women

Breast cancer is the most common cancer diagnosis and second leading cause of death in women. Risk factors associated with breast cancer include: increased age, alcohol consumption, cigarette smoking, white race, physical inactivity, benign breast conditions, reproductive and hormonal factors, dietary factors, and family history. Hereditary breast and ovarian cancer syndrome (HBOC) is caused by mutations in the BRCA1 and BRCA2 genes. Women carrying a mutation in these genes are at an increased risk to develop a second breast cancer. Contralateral breast cancer is the most common second primary cancer in patients treated for a first breast cancer. Other risk factors for developing contralateral breast cancer include a strong family history of breast cancer, age of onset of first primary breast cancer, and if the first primary was a lobular carcinoma, which has an increased risk of being bilateral. A retrospective chart review was performed on a select cohort of women in an IRB approved database at MD Anderson Cancer Center. The final cohort contained 572 women who tested negative for a BRCA1 or BRCA2 mutation, had their primary invasive breast cancer diagnosed under the age of 50, and had a BRCAPro risk assessment number over 10%. Of the 572 women, 97 women developed contralateral breast cancer. A number of predictors of contralateral breast cancer were looked at between the two groups. Using univariable Cox Proportional Hazard model, thirteen statistically interesting risk factors were found, defined as having a p-value under 0.2. Multivariable stepwise Cox Proportional Hazard model found four statistically significant variables out of the thirteen found in the univariable analysis. In our study population, the incidence of contralateral breast cancer was 17%. Four statistically significant variables were identified. Undergoing a prophylactic mastectomy was found to reduce the risk of developing contralateral breast cancer, while not having a prophylactic mastecomy, a young age at primary diagnosis, having a positive estrogen receptor status of the primary tumor, and having a family history of breast cancer increased a woman’s risk to develop contralateral breast cancer.

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

Comparison of clinical and economic outcomes between high-intensity care and hospice care for the end-of-life melanoma patients

The intensity of care for patients at the end-of-life is increasing in recent years. Publications have focused on intensity of care for many cancers, but none on melanoma patients. Substantial gaps exist in knowledge about intensive care and its alternative, hospice care, among the advanced melanoma patients at the end of life. End-of-life care may be used in quite different patterns and induce both intended and unintended clinical and economic consequences. We used the Surveillance, Epidemiology, and End Results (SEER)-Medicare linked databases to identify patients aged 65 years or older with metastatic melanoma who died between 2000 and 2007. We evaluated trends and associations between sociodemographic and health services characteristics and the use of hospice care, chemotherapy, surgery, and radiation therapy and costs. Survival, end-of-life costs, and incremental cost-effectiveness ratio were evaluated using propensity score methods. Costs were analyzed from the perspective of Medicare in 2009 dollars. In the first journal Article we found increasing use of surgery for patients with metastatic melanoma from 13% in 2000 to 30% in 2007 (P=0.03 for trend), no significant fluctuation in use of chemotherapy (P=0.43) or radiation therapy (P=0.46). Older patients were less likely to receive radiation therapy or chemotherapy. The use of hospice care increased from 61% in 2000 to 79% in 2007 (P =0.07 for trend). Enrollment in short-term (1-3 days) hospice care use increased, while long-term hospice care (≥ 4 days) remained stable. Patients living in the SEER Northeast and South regions were less likely to undergo surgery. Patients enrolled in long-term hospice care used significantly less chemotherapy, surgery and radiation therapy. In the second journal article, of 611 patients identified for this study, 358 (59%) received no hospice care after their diagnosis, 168 (27%) received 1 to 3 days of hospice care, and 85 (14%) received 4 or more days of hospice care. The median survival time was 181 days for patients with no hospice care, 196 days for patients enrolled in hospice for 1 to 3 days, and 300 days for patients enrolled for 4 or more days (log-rank test, P < 0.001). The estimated hazard ratios (HR) between 4 or more days hospice use and survival were similar within the original cohort Cox proportional hazard model (HR, 0.62; 95% CI, 0.49-0.78, P < 0.0001) and the propensity score-matched model (HR, 0.61; 95% CI, 0.47-0.78, P = 0.0001). Patients with ≥ 4 days of hospice care incurred lower end-of-life costs than the other two groups ($14,298 versus $19,380 for the 1- to 3-days hospice care, and $24,351 for patients with no hospice care; p < 0.0001). In conclusion, Surgery and hospice care use increased over the years of this study while the use of chemotherapy and radiation therapy remained consistent for patients diagnosed with metastatic melanoma. Patients diagnosed with advanced melanoma who enrolled in ≥ 4 days of hospice care experienced longer survival than those who had 1-3 days of hospice or no hospice care, and this longer overall survival was accompanied by lower end-of-life costs.^

Effects of informative censoring on the hazard function: Application to the proportional hazards regression model

The standard analyses of survival data involve the assumption that survival and censoring are independent. When censoring and survival are related, the phenomenon is known as informative censoring. This paper examines the effects of an informative censoring assumption on the hazard function and the estimated hazard ratio provided by the Cox model.^ The limiting factor in all analyses of informative censoring is the problem of non-identifiability. Non-identifiability implies that it is impossible to distinguish a situation in which censoring and death are independent from one in which there is dependence. However, it is possible that informative censoring occurs. Examination of the literature indicates how others have approached the problem and covers the relevant theoretical background.^ Three models are examined in detail. The first model uses conditionally independent marginal hazards to obtain the unconditional survival function and hazards. The second model is based on the Gumbel Type A method for combining independent marginal distributions into bivariate distributions using a dependency parameter. Finally, a formulation based on a compartmental model is presented and its results described. For the latter two approaches, the resulting hazard is used in the Cox model in a simulation study.^ The unconditional survival distribution formed from the first model involves dependency, but the crude hazard resulting from this unconditional distribution is identical to the marginal hazard, and inferences based on the hazard are valid. The hazard ratios formed from two distributions following the Gumbel Type A model are biased by a factor dependent on the amount of censoring in the two populations and the strength of the dependency of death and censoring in the two populations. The Cox model estimates this biased hazard ratio. In general, the hazard resulting from the compartmental model is not constant, even if the individual marginal hazards are constant, unless censoring is non-informative. The hazard ratio tends to a specific limit.^ Methods of evaluating situations in which informative censoring is present are described, and the relative utility of the three models examined is discussed. ^

The number and timing of time -dependent covariate measurements for the Cox regression model

The problem of analyzing data with updated measurements in the time-dependent proportional hazards model arises frequently in practice. One available option is to reduce the number of intervals (or updated measurements) to be included in the Cox regression model. We empirically investigated the bias of the estimator of the time-dependent covariate while varying the effect of failure rate, sample size, true values of the parameters and the number of intervals. We also evaluated how often a time-dependent covariate needs to be collected and assessed the effect of sample size and failure rate on the power of testing a time-dependent effect.^ A time-dependent proportional hazards model with two binary covariates was considered. The time axis was partitioned into k intervals. The baseline hazard was assumed to be 1 so that the failure times were exponentially distributed in the ith interval. A type II censoring model was adopted to characterize the failure rate. The factors of interest were sample size (500, 1000), type II censoring with failure rates of 0.05, 0.10, and 0.20, and three values for each of the non-time-dependent and time-dependent covariates (1/4,1/2,3/4).^ The mean of the bias of the estimator of the coefficient of the time-dependent covariate decreased as sample size and number of intervals increased whereas the mean of the bias increased as failure rate and true values of the covariates increased. The mean of the bias of the estimator of the coefficient was smallest when all of the updated measurements were used in the model compared with two models that used selected measurements of the time-dependent covariate. For the model that included all the measurements, the coverage rates of the estimator of the coefficient of the time-dependent covariate was in most cases 90% or more except when the failure rate was high (0.20). The power associated with testing a time-dependent effect was highest when all of the measurements of the time-dependent covariate were used. An example from the Systolic Hypertension in the Elderly Program Cooperative Research Group is presented. ^