1
|
Misetic H, Keddar MR, Jeannon JP, Ciccarelli FD. Mechanistic insights into the interactions between cancer drivers and the tumour immune microenvironment. Genome Med 2023; 15:40. [PMID: 37277866 DOI: 10.1186/s13073-023-01197-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Accepted: 05/25/2023] [Indexed: 06/07/2023] Open
Abstract
BACKGROUND The crosstalk between cancer and the tumour immune microenvironment (TIME) has attracted significant interest in the latest years because of its impact on cancer evolution and response to treatment. Despite this, cancer-specific tumour-TIME interactions and their mechanistic insights are still poorly understood. METHODS Here, we compute the significant interactions occurring between cancer-specific genetic drivers and five anti- and pro-tumour TIME features in 32 cancer types using Lasso regularised ordinal regression. Focusing on head and neck squamous cancer (HNSC), we rebuild the functional networks linking specific TIME driver alterations to the TIME state they associate with. RESULTS The 477 TIME drivers that we identify are multifunctional genes whose alterations are selected early in cancer evolution and recur across and within cancer types. Tumour suppressors and oncogenes have an opposite effect on the TIME and the overall anti-tumour TIME driver burden is predictive of response to immunotherapy. TIME driver alterations predict the immune profiles of HNSC molecular subtypes, and perturbations in keratinization, apoptosis and interferon signalling underpin specific driver-TIME interactions. CONCLUSIONS Overall, our study delivers a comprehensive resource of TIME drivers, gives mechanistic insights into their immune-regulatory role, and provides an additional framework for patient prioritisation to immunotherapy. The full list of TIME drivers and associated properties are available at http://www.network-cancer-genes.org .
Collapse
Affiliation(s)
- Hrvoje Misetic
- Cancer Systems Biology Laboratory, The Francis Crick Institute, London, NW1 1AT, UK
- School of Cancer and Pharmaceutical Sciences, King's College London, London, SE11UL, UK
| | - Mohamed Reda Keddar
- Cancer Systems Biology Laboratory, The Francis Crick Institute, London, NW1 1AT, UK
- School of Cancer and Pharmaceutical Sciences, King's College London, London, SE11UL, UK
| | - Jean-Pierre Jeannon
- Department of Head & Neck Surgery, Great Maze Pond, Guy's Hospital, London, SE1 9RT, UK
| | - Francesca D Ciccarelli
- Cancer Systems Biology Laboratory, The Francis Crick Institute, London, NW1 1AT, UK.
- School of Cancer and Pharmaceutical Sciences, King's College London, London, SE11UL, UK.
| |
Collapse
|
2
|
Zhang H, Zhang N, Wu W, Zhou R, Li S, Wang Z, Dai Z, Zhang L, Liu Z, Zhang J, Luo P, Liu Z, Cheng Q. Machine learning-based tumor-infiltrating immune cell-associated lncRNAs for predicting prognosis and immunotherapy response in patients with glioblastoma. Brief Bioinform 2022; 23:6711411. [PMID: 36136350 DOI: 10.1093/bib/bbac386] [Citation(s) in RCA: 46] [Impact Index Per Article: 23.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Revised: 07/29/2022] [Accepted: 08/10/2022] [Indexed: 12/14/2022] Open
Abstract
Long noncoding ribonucleic acids (RNAs; lncRNAs) have been associated with cancer immunity regulation. However, the roles of immune cell-specific lncRNAs in glioblastoma (GBM) remain largely unknown. In this study, a novel computational framework was constructed to screen the tumor-infiltrating immune cell-associated lncRNAs (TIIClnc) for developing TIIClnc signature by integratively analyzing the transcriptome data of purified immune cells, GBM cell lines and bulk GBM tissues using six machine learning algorithms. As a result, TIIClnc signature could distinguish survival outcomes of GBM patients across four independent datasets, including the Xiangya in-house dataset, and more importantly, showed superior performance than 95 previously established signatures in gliomas. TIIClnc signature was revealed to be an indicator of the infiltration level of immune cells and predicted the response outcomes of immunotherapy. The positive correlation between TIIClnc signature and CD8, PD-1 and PD-L1 was verified in the Xiangya in-house dataset. As a newly demonstrated predictive biomarker, the TIIClnc signature enabled a more precise selection of the GBM population who would benefit from immunotherapy and should be validated and applied in the near future.
Collapse
Affiliation(s)
- Hao Zhang
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China.,Department of Neurosurgery, The Second Affiliated Hospital, Chongqing Medical University, China
| | - Nan Zhang
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,One-third Lab, College of Bioinformatics Science and Technology, Harbin Medical University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China
| | - Wantao Wu
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,Department of Oncology, Xiangya Hospital, Central South University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China
| | - Ran Zhou
- Division of Neuroscience and Experimental Psychology, Faculty of Biology, Medicine and Health, University of Manchester, UK
| | - Shuyu Li
- Department of Thyroid and Breast Surgery, Tongji Hospital, Tongji Medical College of Huazhong University of Science and Technology, China
| | - Zeyu Wang
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China
| | - Ziyu Dai
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China
| | - Liyang Zhang
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China
| | - Zaoqu Liu
- Department of Interventional Radiology, The First Affiliated Hospital of Zhengzhou, China
| | - Jian Zhang
- Department of Oncology, Zhujiang Hospital, Southern Medical University, China
| | - Peng Luo
- Department of Oncology, Zhujiang Hospital, Southern Medical University, China
| | - Zhixiong Liu
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China
| | - Quan Cheng
- Department of Neurosurgery, Xiangya Hospital, Central South University, China.,National Clinical Research Center for Geriatric Disorders, Xiangya Hospital, Central South University, China
| |
Collapse
|
3
|
Iannario M, Monti AC, Scalera P. The number of response categories in ordered response models. Int J Biostat 2022; 18:593-611. [PMID: 34547198 DOI: 10.1515/ijb-2021-0013] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Accepted: 08/23/2021] [Indexed: 01/10/2023]
Abstract
The choice of the number m of response categories is a crucial issue in categorization of a continuous response. The paper exploits the Proportional Odds Models' property which allows to generate ordinal responses with a different number of categories from the same underlying variable. It investigates the asymptotic efficiency of the estimators of the regression coefficients and the accuracy of the derived inferential procedures when m varies. The analysis is based on models with closed-form information matrices so that the asymptotic efficiency can be analytically evaluated without need of simulations. The paper proves that a finer categorization augments the information content of the data and consequently shows that the asymptotic efficiency and the power of the tests on the regression coefficients increase with m. The impact of the loss of information produced by merging categories on the efficiency of the estimators is also considered, highlighting its risks especially when performed in its extreme form of dichotomization. Furthermore, the appropriate value of m for various sample sizes is explored, pointing out that a large number of categories can offset the limited amount of information of a small sample by a better quality of the data. Finally, two case studies on the quality of life of chemotherapy patients and on the perception of pain, based on discretized continuous scales, illustrate the main findings of the paper.
Collapse
Affiliation(s)
- Maria Iannario
- Department of Political Sciences, University of Naples Federico II, Napoli, Italy
| | - Anna Clara Monti
- Department of Law, Economics, Management and Quantitative Methods, University of Sannio, Benevento, Italy
| | - Pietro Scalera
- Department of Political Sciences, University of Naples Federico II, Napoli, Italy
| |
Collapse
|
4
|
Macnair W, Gupta R, Claassen M. psupertime: supervised pseudotime analysis for time-series single-cell RNA-seq data. Bioinformatics 2022; 38:i290-i298. [PMID: 35758781 PMCID: PMC9235474 DOI: 10.1093/bioinformatics/btac227] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Motivation Improvements in single-cell RNA-seq technologies mean that studies measuring multiple experimental conditions, such as time series, have become more common. At present, few computational methods exist to infer time series-specific transcriptome changes, and such studies have therefore typically used unsupervised pseudotime methods. While these methods identify cell subpopulations and the transitions between them, they are not appropriate for identifying the genes that vary coherently along the time series. In addition, the orderings they estimate are based only on the major sources of variation in the data, which may not correspond to the processes related to the time labels. Results We introduce psupertime, a supervised pseudotime approach based on a regression model, which explicitly uses time-series labels as input. It identifies genes that vary coherently along a time series, in addition to pseudotime values for individual cells, and a classifier that can be used to estimate labels for new data with unknown or differing labels. We show that psupertime outperforms benchmark classifiers in terms of identifying time-varying genes and provides better individual cell orderings than popular unsupervised pseudotime techniques. psupertime is applicable to any single-cell RNA-seq dataset with sequential labels (e.g. principally time series but also drug dosage and disease progression), derived from either experimental design and provides a fast, interpretable tool for targeted identification of genes varying along with specific biological processes. Availability and implementation R package available at github.com/wmacnair/psupertime and code for results reproduction at github.com/wmacnair/psupplementary. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Will Macnair
- Institute of Molecular Systems Biology, Department of Biology, ETH Zurich, Zurich 8093, Switzerland
| | - Revant Gupta
- Inner Medicine I, Faculty of Medicine, University of Tübingen, University Hospital Tübingen, 72074, Germany
| | - Manfred Claassen
- Inner Medicine I, Faculty of Medicine, University of Tübingen, University Hospital Tübingen, 72074, Germany.,Department of Computer Science, University of Tübingen, Tübingen 72074, Germany
| |
Collapse
|
5
|
ordinalbayes: Fitting Ordinal Bayesian Regression Models to High-Dimensional Data Using R. STATS 2022; 5:371-384. [PMID: 35574500 PMCID: PMC9097970 DOI: 10.3390/stats5020021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
Abstract
The stage of cancer is a discrete ordinal response that indicates the aggressiveness of disease and is often used by physicians to determine the type and intensity of treatment to be administered. For example, the FIGO stage in cervical cancer is based on the size and depth of the tumor as well as the level of spread. It may be of clinical relevance to identify molecular features from high-throughput genomic assays that are associated with the stage of cervical cancer to elucidate pathways related to tumor aggressiveness, identify improved molecular features that may be useful for staging, and identify therapeutic targets. High-throughput RNA-Seq data and corresponding clinical data (including stage) for cervical cancer patients have been made available through The Cancer Genome Atlas Project (TCGA). We recently described penalized Bayesian ordinal response models that can be used for variable selection for over-parameterized datasets, such as the TCGA-CESC dataset. Herein, we describe our ordinalbayes R package, available from the Comprehensive R Archive Network (CRAN), which enhances the runjags R package by enabling users to easily fit cumulative logit models when the outcome is ordinal and the number of predictors exceeds the sample size, P > N, such as for TCGA and other high-throughput genomic data. We demonstrate the use of this package by applying it to the TCGA cervical cancer dataset. Our ordinalbayes package can be used to fit models to high-dimensional datasets, and it effectively performs variable selection.
Collapse
|
6
|
Titz B, Sewer A, Luettich K, Wong ET, Guedj E, Nury C, Schneider T, Xiang Y, Trivedi K, Vuillaume G, Leroy P, Büttner A, Martin F, Ivanov NV, Vanscheeuwijck P, Hoeng J, Peitsch MC. Respiratory Effects of Exposure to Aerosol From the Candidate Modified-Risk Tobacco Product THS 2.2 in an 18-Month Systems Toxicology Study With A/J Mice. Toxicol Sci 2021; 178:138-158. [PMID: 32780831 PMCID: PMC7657339 DOI: 10.1093/toxsci/kfaa132] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Smoking cessation is the most effective measure for reducing the risk of smoking-related diseases. However, switching to less harmful products (modified-risk tobacco products [MRTP]) can be an alternative to help reduce the risk for adult smokers who would otherwise continue to smoke. In an 18-month chronic carcinogenicity/toxicity study in A/J mice (OECD Test Guideline 453), we assessed the aerosol of Tobacco Heating System 2.2 (THS 2.2), a candidate MRTP based on the heat-not-burn principle, compared with 3R4F cigarette smoke (CS). To capture toxicity- and disease-relevant mechanisms, we complemented standard toxicology endpoints with in-depth systems toxicology analyses. In this part of our publication series, we report on integrative assessment of the apical and molecular exposure effects on the respiratory tract (nose, larynx, and lungs). Across the respiratory tract, we found changes in inflammatory response following 3R4F CS exposure (eg, antimicrobial peptide response in the nose), with both shared and distinct oxidative and xenobiotic responses. Compared with 3R4F CS, THS 2.2 aerosol exerted far fewer effects on respiratory tract histology, including adaptive tissue changes in nasal and laryngeal epithelium and inflammation and emphysematous changes in the lungs. Integrative analysis of molecular changes confirmed the substantially lower impact of THS 2.2 aerosol than 3R4F CS on toxicologically and disease-relevant molecular processes such as inflammation, oxidative stress responses, and xenobiotic metabolism. In summary, this work exemplifies how apical and molecular endpoints can be combined effectively for toxicology assessment and further supports findings on the reduced respiratory health risks of THS 2.2 aerosol.
Collapse
Affiliation(s)
- Bjoern Titz
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | - Alain Sewer
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | - Karsta Luettich
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | - Ee Tsin Wong
- Philip Morris International Research Laboratories Pte. Ltd, Singapore 117406
| | - Emmanuel Guedj
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | - Catherine Nury
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | | | - Yang Xiang
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | - Keyur Trivedi
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | | | - Patrice Leroy
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | | | - Florian Martin
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | - Nikolai V Ivanov
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | | | - Julia Hoeng
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| | - Manuel C Peitsch
- PMI R&D, Philip Morris Products S.A, CH-2000 Neuchâtel, Switzerland
| |
Collapse
|
7
|
Abstract
The so-called proportional odds assumption is popular in cumulative, ordinal regression. In practice, however, such an assumption is sometimes too restrictive. For instance, when modeling the perception of boar taint on an individual level, it turns out that, at least for some subjects, the effects of predictors (androstenone and skatole) vary between response categories. For more flexible modeling, we consider the use of a ‘smooth-effects-on-response penalty’ (SERP) as a connecting link between proportional and fully non-proportional odds models, assuming that parameters of the latter vary smoothly over response categories. The usefulness of SERP is further demonstrated through a simulation study. Besides flexible and accurate modeling, SERP also enables fitting of parameters in cases where the pure, unpenalized non-proportional odds model fails to converge.
Collapse
|
8
|
Liu H, Wang X, Tang K, Peng E, Xia D, Chen Z. Machine learning-assisted decision-support models to better predict patients with calculous pyonephrosis. Transl Androl Urol 2021; 10:710-723. [PMID: 33718073 PMCID: PMC7947454 DOI: 10.21037/tau-20-1208] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
Background To develop a machine learning (ML)-assisted model capable of accurately identifying patients with calculous pyonephrosis before making treatment decisions by integrating multiple clinical characteristics. Methods We retrospectively collected data from patients with obstructed hydronephrosis who underwent retrograde ureteral stent insertion, percutaneous nephrostomy (PCN), or percutaneous nephrolithotomy (PCNL). The study cohort was divided into training and testing datasets in a 70:30 ratio for further analysis. We developed 5 ML-assisted models from 22 clinical features using logistic regression (LR), LR optimized by least absolute shrinkage and selection operator (Lasso) regularization (Lasso-LR), support vector machine (SVM), extreme gradient boosting (XGBoost), and random forest (RF). The area under the curve (AUC) was applied to determine the model with the highest discrimination. Decision curve analysis (DCA) was used to investigate the clinical net benefit associated with using the predictive models. Results A total of 322 patients were included, with 225 patients in the training dataset, and 97 patients in the testing dataset. The XGBoost model showed good discrimination with the AUC, accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) of 0.981, 0.991, 0.962, 1.000, 1.000, and 0.989, respectively, followed by SVM [AUC =0.985, 95% confidence interval (CI): 0.970–1.000], Lasso-LR (AUC =0.977, 95% CI: 0.958–0.996), LR (AUC =0.936, 95% CI: 0.905–0.968), and RF (AUC =0.920, 95% CI: 0.870–0.970). Validation of the model showed that SVM yielded the highest AUC (0.977, 95% CI: 0.952–1.000), followed by Lasso-LR (AUC =0.959, 95% CI: 0.921–0.997), XGBoost (AUC =0.958, 95% CI: 0.902–1.000), LR (AUC =0.932, 95% CI: 0.878–0.987), and RF (AUC =0.868, 95% CI: 0.779–0.958) in the testing dataset. Conclusions Our ML-based models had good discrimination in predicting patients with obstructed hydronephrosis at high risk of harboring pyonephrosis, and the use of these models may be greatly beneficial to urologists in treatment planning, patient selection, and decision-making.
Collapse
Affiliation(s)
- Hailang Liu
- Department of Urology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Xinguang Wang
- Department of Urology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Kun Tang
- Department of Urology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Ejun Peng
- Department of Urology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Ding Xia
- Department of Urology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China
| | - Zhiqiang Chen
- Department of Urology, Tongji Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, Hubei, China
| |
Collapse
|
9
|
Zhang Y, Archer KJ. Bayesian penalized cumulative logit model for high-dimensional data with an ordinal response. Stat Med 2020; 40:1453-1481. [PMID: 33336826 DOI: 10.1002/sim.8851] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2019] [Revised: 11/23/2020] [Accepted: 11/23/2020] [Indexed: 01/15/2023]
Abstract
Many previous studies have identified associations between gene expression, measured using high-throughput genomic platforms, and quantitative or dichotomous traits. However, we note that health outcome and disease status measurements frequently appear on an ordinal scale, that is, the outcome is categorical but has inherent ordering. Identification of important genes may be useful for developing novel diagnostic and prognostic tools to predict or classify stage of disease. Gene expression data are usually high-dimensional, meaning that the number of genes is much larger than the sample size or number of patients. Herein we describe some existing frequentist methods for modeling an ordinal response in a high-dimensional predictor space. Following Tibshirani (1996), who described the LASSO estimate as the Bayesian posterior mode when the regression coefficients have independent Laplace priors, we propose a new approach for high-dimensional data with an ordinal response that is rooted in the Bayesian paradigm. We show that our proposed Bayesian approach outperforms existing frequentist methods through simulation studies. We then compare the performance of frequentist and Bayesian approaches using a study evaluating progression to hepatocellular carcinoma in hepatitis C infected patients.
Collapse
Affiliation(s)
- Yiran Zhang
- College of Public Health, The Ohio State University, Columbus, Ohio, USA
| | - Kellie J Archer
- College of Public Health, The Ohio State University, Columbus, Ohio, USA
| |
Collapse
|
10
|
Developing an instrument for an early prediction model of long-term functional outcomes in people with acquired injuries of the central nervous system: protocol and methodological aspects. Neurol Sci 2020; 42:2441-2446. [PMID: 33078248 PMCID: PMC8159777 DOI: 10.1007/s10072-020-04821-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2020] [Accepted: 10/10/2020] [Indexed: 11/24/2022]
Abstract
Severe acquired brain injury (ABI) is a major cause of long-term disability and is the main determinant of health and societal costs. Early identification of favourable long-term recovery would allow personalized rehabilitative programs and better health care resources allocation. In light of the higher survival rate from intensive care units (ICU) in recent years, there is a growing need for early prognostication markers of functional recovery; to date, these data have been mainly collected at rehabilitation unit admission and not during the acute phase. We present the protocol and methodology to develop prediction models in people with severe acquired brain injury (GCS at admission to ICU < 8) for the functional and cognitive outcome at 12 months from the event. Predictors will be collected during the acute stage. Participants will be recruited within the first 72 h from the event in the ICUs of two teaching hospitals (Padova and Treviso). Participants will be followed up at discharge from ICU, admission and discharge from Neurorehabilitation and after 12 months from the event. Clinical and functional scales, electroencephalography, evoked potentials, magnetic resonance imaging and serological markers will be entered into a digital registry. Survival will be estimated using the Cox proportional hazard model. A multivariate prediction model will be developed for each of the functional and cognitive outcomes at 12 months from the event.
Collapse
|
11
|
Predicting Ewe Body Condition Score Using Lifetime Liveweight and Liveweight Change, and Previous Body Condition Score Record. Animals (Basel) 2020; 10:ani10071182. [PMID: 32668688 PMCID: PMC7401657 DOI: 10.3390/ani10071182] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 06/25/2020] [Accepted: 07/09/2020] [Indexed: 11/17/2022] Open
Abstract
The body condition score (BCS) in sheep (Ovis aries) is a widely used subjective measure of body condition. Body condition score and liveweight have been reported to be statistically and often linearly related in ewes. Therefore, it was hypothesized that current BCS could be accurately and indirectly predicted using a ewe's lifetime liveweight, liveweight change, and previous BCS record. Ewes born between 2011 and 2012 (n = 11,798) were followed from 8 months to approximately 67 months of age in New Zealand. Individual ewe data was collected on liveweight and body condition scores at each stage of the annual cycle (pre-breeding, pregnancy diagnosis, pre-lambing, and weaning). Linear regression models were fitted to predict BCS at a given ewe age and stage of the annual cycle using a ewe's lifetime liveweight records (liveweight alone models). Further, linear models were then fitted using previous BCS and changes in liveweight, in addition to the lifetime liveweight records (combined models). Using the combined models improved (p < 0.01) the R2 value by 39.8% (from 0.32 to 0.45) and lowered the average prediction error by 10% to 12% (from 0.29 to 0.26 body condition scores). However, a significant portion of the variability in BCS remained unaccounted for (39% to 89%) even in the combined models. The procedures found in this study, therefore, may overestimate or underestimate measures by 0.23 to 0.32 BCS, which could substantially change the status of the ewe, leading to incorrect management decisions. However, the findings do still suggest that there is potential for predicting ewe BCS from liveweight using linear regression if the key variables affecting the relationship between BCS and liveweight are accounted for.
Collapse
|
12
|
Fu H, Archer KJ. High-dimensional variable selection for ordinal outcomes with error control. Brief Bioinform 2020; 22:334-345. [PMID: 32031572 DOI: 10.1093/bib/bbaa007] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Revised: 01/06/2020] [Indexed: 12/24/2022] Open
Abstract
Many high-throughput genomic applications involve a large set of potential covariates and a response which is frequently measured on an ordinal scale, and it is crucial to identify which variables are truly associated with the response. Effectively controlling the false discovery rate (FDR) without sacrificing power has been a major challenge in variable selection research. This study reviews two existing variable selection frameworks, model-X knockoffs and a modified version of reference distribution variable selection (RDVS), both of which utilize artificial variables as benchmarks for decision making. Model-X knockoffs constructs a 'knockoff' variable for each covariate to mimic the covariance structure, while RDVS generates only one null variable and forms a reference distribution by performing multiple runs of model fitting. Herein, we describe how different importance measures for ordinal responses can be constructed that fit into these two selection frameworks, using either penalized regression or machine learning techniques. We compared these measures in terms of the FDR and power using simulated data. Moreover, we applied these two frameworks to high-throughput methylation data for identifying features associated with the progression from normal liver tissue to hepatocellular carcinoma to further compare and contrast their performances.
Collapse
|
13
|
Connell SP, Yazbek‐Hanna M, McCarthy F, Hurst R, Webb M, Curley H, Walker H, Mills R, Ball RY, Sanda MG, Pellegrini KL, Patil D, Perry AS, Schalken J, Pandha H, Whitaker H, Dennis N, Stuttle C, Mills IG, Guldvik I, Parker C, Brewer DS, Cooper CS, Clark J. A four-group urine risk classifier for predicting outcomes in patients with prostate cancer. BJU Int 2019; 124:609-620. [PMID: 31106513 PMCID: PMC6851983 DOI: 10.1111/bju.14811] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
OBJECTIVES To develop a risk classifier using urine-derived extracellular vesicle (EV)-RNA capable of providing diagnostic information on disease status prior to biopsy, and prognostic information for men on active surveillance (AS). PATIENTS AND METHODS Post-digital rectal examination urine-derived EV-RNA expression profiles (n = 535, multiple centres) were interrogated with a curated NanoString panel. A LASSO-based continuation ratio model was built to generate four prostate urine risk (PUR) signatures for predicting the probability of normal tissue (PUR-1), D'Amico low-risk (PUR-2), intermediate-risk (PUR-3), and high-risk (PUR-4) prostate cancer. This model was applied to a test cohort (n = 177) for diagnostic evaluation, and to an AS sub-cohort (n = 87) for prognostic evaluation. RESULTS Each PUR signature was significantly associated with its corresponding clinical category (P < 0.001). PUR-4 status predicted the presence of clinically significant intermediate- or high-risk disease (area under the curve = 0.77, 95% confidence interval [CI] 0.70-0.84). Application of PUR provided a net benefit over current clinical practice. In an AS sub-cohort (n = 87), groups defined by PUR status and proportion of PUR-4 had a significant association with time to progression (interquartile range hazard ratio [HR] 2.86, 95% CI 1.83-4.47; P < 0.001). PUR-4, when used continuously, dichotomized patient groups with differential progression rates of 10% and 60% 5 years after urine collection (HR 8.23, 95% CI 3.26-20.81; P < 0.001). CONCLUSION Urine-derived EV-RNA can provide diagnostic information on aggressive prostate cancer prior to biopsy, and prognostic information for men on AS. PUR represents a new and versatile biomarker that could result in substantial alterations to current treatment of patients with prostate cancer.
Collapse
Affiliation(s)
| | | | | | - Rachel Hurst
- Norwich Medical SchoolUniversity of East AngliaNorwichUK
| | - Martyn Webb
- Norwich Medical SchoolUniversity of East AngliaNorwichUK
| | - Helen Curley
- Norwich Medical SchoolUniversity of East AngliaNorwichUK
| | - Helen Walker
- Norfolk and Norwich University Hospitals NHS Foundation TrustNorwichUK
| | - Rob Mills
- Norfolk and Norwich University Hospitals NHS Foundation TrustNorwichUK
| | - Richard Y. Ball
- Norfolk and Norwich University Hospitals NHS Foundation TrustNorwichUK
| | - Martin G. Sanda
- Department of UrologyWinship Cancer InstituteEmory University School of MedicineAtlantaGAUSA
| | - Kathryn L. Pellegrini
- Department of UrologyWinship Cancer InstituteEmory University School of MedicineAtlantaGAUSA
| | - Dattatraya Patil
- Department of UrologyWinship Cancer InstituteEmory University School of MedicineAtlantaGAUSA
| | - Antoinette S. Perry
- School of Biology and Environmental ScienceScience WestUniversity College DublinDublin 4Ireland
| | - Jack Schalken
- Nijmegen Medical CentreRadboud University Medical CentreNijmegenThe Netherlands
| | - Hardev Pandha
- Faculty of Health and Medical SciencesThe University of SurreyGuildfordUK
| | - Hayley Whitaker
- Molecular Diagnostics and Therapeutics GroupUniversity College LondonLondonUK
| | | | | | - Ian G. Mills
- School of MedicineDentistry and Biomedical SciencesInstitute for Health SciencesCentre for Cancer Research and Cell BiologyQueen's University BelfastBelfastUK,Centre for Molecular MedicineUniversity of OsloOsloNorway,Nuffield Department of Surgical SciencesUniversity of OxfordOxfordUK
| | - Ingrid Guldvik
- Centre for Molecular MedicineUniversity of OsloOsloNorway
| | | | | | - Daniel S. Brewer
- Norwich Medical SchoolUniversity of East AngliaNorwichUK,Earlham InstituteNorwichUK
| | | | - Jeremy Clark
- Norwich Medical SchoolUniversity of East AngliaNorwichUK
| |
Collapse
|
14
|
Ma J, Stingo FC, Hobbs BP. Bayesian personalized treatment selection strategies that integrate predictive with prognostic determinants. Biom J 2019; 61:902-917. [PMID: 30786040 PMCID: PMC7341533 DOI: 10.1002/bimj.201700323] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2017] [Revised: 09/28/2018] [Accepted: 12/04/2018] [Indexed: 01/13/2023]
Abstract
The evolution of "informatics" technologies has the potential to generate massive databases, but the extent to which personalized medicine may be effectuated depends on the extent to which these rich databases may be utilized to advance understanding of the disease molecular profiles and ultimately integrated for treatment selection, necessitating robust methodology for dimension reduction. Yet, statistical methods proposed to address challenges arising with the high-dimensionality of omics-type data predominately rely on linear models and emphasize associations deriving from prognostic biomarkers. Existing methods are often limited for discovering predictive biomarkers that interact with treatment and fail to elucidate the predictive power of their resultant selection rules. In this article, we present a Bayesian predictive method for personalized treatment selection that is devised to integrate both the treatment predictive and disease prognostic characteristics of a particular patient's disease. The method appropriately characterizes the structural constraints inherent to prognostic and predictive biomarkers, and hence properly utilizes these complementary sources of information for treatment selection. The methodology is illustrated through a case study of lower grade glioma. Theoretical considerations are explored to demonstrate the manner in which treatment selection is impacted by prognostic features. Additionally, simulations based on an actual leukemia study are provided to ascertain the method's performance with respect to selection rules derived from competing methods.
Collapse
Affiliation(s)
- Junsheng Ma
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center Houston, Texas 77030
| | - Francesco C. Stingo
- Department of Statistica, Informatica, Applicazioni “G.Parenti”, University of Florence, Florence, 50134, Italy
| | - Brian P. Hobbs
- Quantitative Health Sciences and The Taussig Cancer Institute, Cleveland Clinic, Cleveland, Ohio 44195
| |
Collapse
|
15
|
Przybyła P, Brockmeier AJ, Ananiadou S. Quantifying risk factors in medical reports with a context-aware linear model. J Am Med Inform Assoc 2019; 26:537-546. [PMID: 30840055 PMCID: PMC6515525 DOI: 10.1093/jamia/ocz004] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2018] [Revised: 12/14/2018] [Accepted: 01/03/2019] [Indexed: 12/03/2022] Open
Abstract
OBJECTIVE We seek to quantify the mortality risk associated with mentions of medical concepts in textual electronic health records (EHRs). Recognizing mentions of named entities of relevant types (eg, conditions, symptoms, laboratory tests or behaviors) in text is a well-researched task. However, determining the level of risk associated with them is partly dependent on the textual context in which they appear, which may describe severity, temporal aspects, quantity, etc. METHODS To take into account that a given word appearing in the context of different risk factors (medical concepts) can make different contributions toward risk level, we propose a multitask approach, called context-aware linear modeling, which can be applied using appropriately regularized linear regression. To improve the performance for risk factors unseen in training data (eg, rare diseases), we take into account their distributional similarity to other concepts. RESULTS The evaluation is based on a corpus of 531 reports from EHRs with 99 376 risk factors rated manually by experts. While context-aware linear modeling significantly outperforms single-task models, taking into account concept similarity further improves performance, reaching the level of human annotators' agreements. CONCLUSION Our results show that automatic quantification of risk factors in EHRs can achieve performance comparable to human assessment, and taking into account the multitask structure of the problem and the ability to handle rare concepts is crucial for its accuracy.
Collapse
Affiliation(s)
- Piotr Przybyła
- National Centre for Text Mining, School of Computer Science, University of Manchester, Manchester, United Kingdom
| | - Austin J Brockmeier
- National Centre for Text Mining, School of Computer Science, University of Manchester, Manchester, United Kingdom
| | - Sophia Ananiadou
- National Centre for Text Mining, School of Computer Science, University of Manchester, Manchester, United Kingdom
| |
Collapse
|
16
|
Lehman RR, Archer KJ. Penalized negative binomial models for modeling an overdispersed count outcome with a high-dimensional predictor space: Application predicting micronuclei frequency. PLoS One 2019; 14:e0209923. [PMID: 30620740 PMCID: PMC6324811 DOI: 10.1371/journal.pone.0209923] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2018] [Accepted: 12/13/2018] [Indexed: 11/19/2022] Open
Abstract
Chromosomal aberrations, such as micronuclei (MN), have served as biomarkers of genotoxic exposure and cancer risk. Guidelines for the process of scoring MN have been presented by the HUman MicroNucleus (HUMN) project. However, these guidelines were developed for assay performance but do not address how to statistically analyze the data generated by the assay. This has led to the application of various statistical methods that may render different interpretations and conclusions. By combining MN with data from other high-throughput genomic technologies such as gene expression microarray data, we may elucidate molecular features involved in micronucleation. Traditional methods that can model discrete (synonymously, count) data, such as MN frequency, require that the number of explanatory variables (P) is less than the number of samples (N). Due to this limitation, penalized models have been developed to enable model fitting for such over-parameterized datasets. Because penalized models in the discrete response setting are lacking, particularly when the count outcome is over-dispersed, herein we present our penalized negative binomial regression model that can be fit when P > N. Using simulation studies we demonstrate the performance of our method in comparison to commonly used penalized Poisson models when the outcome is over-dispersed and applied it to MN frequency and gene expression data collected as part of the Norwegian Mother and Child Cohort Study. Our countgmifs R package is available for download from the Comprehensive R Archive Network and can be applied to datasets having a discrete outcome that is either Poisson or negative binomial distributed and a high-dimensional covariate space.
Collapse
Affiliation(s)
- Rebecca R. Lehman
- United Network for Organ Sharing, Richmond, VA, United States of America
| | - Kellie J. Archer
- Division of Biostatistics, College of Public Health, The Ohio State University, Columbus, OH, United States of America
- * E-mail:
| |
Collapse
|
17
|
Mankad S, Hu S, Gopal A. Single stage prediction with embedded topic modeling of online reviews for mobile app management. Ann Appl Stat 2018. [DOI: 10.1214/18-aoas1152] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
18
|
|
19
|
Starr TN, Picton LK, Thornton JW. Alternative evolutionary histories in the sequence space of an ancient protein. Nature 2017; 549:409-413. [PMID: 28902834 PMCID: PMC6214350 DOI: 10.1038/nature23902] [Citation(s) in RCA: 112] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2017] [Accepted: 08/08/2017] [Indexed: 12/28/2022]
Abstract
To understand why molecular evolution turned out as it did, we must characterize not only the path that evolution followed across the space of possible molecular sequences but also the many alternative trajectories that could have been taken but were not. A large-scale comparison of real and possible histories would establish whether the outcome of evolution represents an optimal state driven by natural selection or the contingent product of historical chance events; it would also reveal how the underlying distribution of functions across sequence space shaped historical evolution. Here we combine ancestral protein reconstruction with deep mutational scanning to characterize alternative histories in the sequence space around an ancient transcription factor, which evolved a novel biological function through well-characterized mechanisms. We find hundreds of alternative protein sequences that use diverse biochemical mechanisms to perform the derived function at least as well as the historical outcome. These alternatives all require prior permissive substitutions that do not enhance the derived function, but not all require the same permissive changes that occurred during history. We find that if evolution had begun from a different starting point within the network of sequences encoding the ancestral function, outcomes with different genetic and biochemical forms would probably have resulted; this contingency arises from the distribution of functional variants in sequence space and epistasis between residues. Our results illuminate the topology of the vast space of possibilities from which history sampled one path, highlighting how the outcome of evolution depends on a serial chain of compounding chance events.
Collapse
Affiliation(s)
- Tyler N Starr
- Department of Biochemistry and Molecular Biology, University of Chicago, Chicago, Illinois 60637, USA
| | - Lora K Picton
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois 60637, USA
| | - Joseph W Thornton
- Department of Ecology and Evolution, University of Chicago, Chicago, Illinois 60637, USA
- Department of Human Genetics, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
20
|
Schauberger G, Tutz G. Subject-specific modelling of paired comparison data: A lasso-type penalty approach. STAT MODEL 2017. [DOI: 10.1177/1471082x17693086] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
In traditional paired comparison models heterogeneity in the population is simply ignored and it is assumed that all persons or subjects have the same preference structure. In the models considered here the preference of an object over another object is explicitly modelled as depending on subject-specific covariates, therefore allowing for heterogeneity in the population. Since by construction the models contain a large number of parameters we propose to use penalized estimation procedures to obtain estimates of the parameters. The used regularized estimation approach penalizes the differences between the parameters corresponding to single covariates. It enforces variable selection and allows to find clusters of objects with respect to covariates. We consider simple binary but also ordinal paired comparisons models. The method is applied to data from a pre-election study from Germany.
Collapse
Affiliation(s)
| | - Gerhard Tutz
- Department of Statistics, Ludwig-Maximilians-Universität München, Germany
| |
Collapse
|
21
|
George NI, Lu TP, Chang CW. Cost-sensitive Performance Metric for Comparing Multiple Ordinal Classifiers. ARTIFICIAL INTELLIGENCE RESEARCH 2017; 5:135-143. [PMID: 28070226 DOI: 10.5430/air.v5n1p135] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The surge of interest in personalized and precision medicine during recent years has increased the application of ordinal classification problems in biomedical science. Currently, accuracy, Kendall's τb , and average mean absolute error are three commonly used metrics for evaluating the effectiveness of an ordinal classifier. Although there are benefits to each, no single metric considers the benefits of predictive accuracy with the tradeoffs of misclassification cost. In addition, decision analysis that considers pairwise analysis of the metrics is not trivial due to inconsistent findings. A new cost-sensitive metric is proposed to find the optimal tradeoff between the two most critical performance measures of a classification task - accuracy and cost. The proposed method accounts for an inherent ordinal data structure, total misclassification cost of a classifier, and imbalanced class distribution. The strengths of the new methodology are demonstrated through analyses of three real cancer datasets and four simulation studies. The new cost-sensitive metric proved better performance in its ability to identify the best ordinal classifier for a given analysis. The performance metric devised in this study provides a comprehensive tool for comparative analysis of multiple (and competing) ordinal classifiers. Consideration of the tradeoff between accuracy and misclassification cost in decisions regarding ordinal classification problems is imperative in real-world application. The work presented here is a precursor to the possibility of incorporating the proposed metric into a prediction modeling algorithm for ordinal data as a means of integrating misclassification cost in final model selection.
Collapse
Affiliation(s)
- Nysia I George
- Division of Bioinformatics and Biostatistics, National Center for Toxicological Research, FDA, Jefferson, AR 72079, USA
| | - Tzu-Pin Lu
- YongLin Biomedical Engineering Center, National Taiwan University, Taipei, Taiwan
| | - Ching-Wei Chang
- Division of Bioinformatics and Biostatistics, National Center for Toxicological Research, FDA, Jefferson, AR 72079, USA
| |
Collapse
|
22
|
Ma J, Stingo FC, Hobbs BP. Bayesian predictive modeling for genomic based personalized treatment selection. Biometrics 2015; 72:575-83. [PMID: 26575856 DOI: 10.1111/biom.12448] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2015] [Revised: 08/01/2015] [Accepted: 10/01/2015] [Indexed: 01/15/2023]
Abstract
Efforts to personalize medicine in oncology have been limited by reductive characterizations of the intrinsically complex underlying biological phenomena. Future advances in personalized medicine will rely on molecular signatures that derive from synthesis of multifarious interdependent molecular quantities requiring robust quantitative methods. However, highly parameterized statistical models when applied in these settings often require a prohibitively large database and are sensitive to proper characterizations of the treatment-by-covariate interactions, which in practice are difficult to specify and may be limited by generalized linear models. In this article, we present a Bayesian predictive framework that enables the integration of a high-dimensional set of genomic features with clinical responses and treatment histories of historical patients, providing a probabilistic basis for using the clinical and molecular information to personalize therapy for future patients. Our work represents one of the first attempts to define personalized treatment assignment rules based on large-scale genomic data. We use actual gene expression data acquired from The Cancer Genome Atlas in the settings of leukemia and glioma to explore the statistical properties of our proposed Bayesian approach for personalizing treatment selection. The method is shown to yield considerable improvements in predictive accuracy when compared to penalized regression approaches.
Collapse
Affiliation(s)
- Junsheng Ma
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, Texas 77030, U.S.A
| | - Francesco C Stingo
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, Texas 77030, U.S.A
| | - Brian P Hobbs
- Department of Biostatistics, The University of Texas MD Anderson Cancer Center, Houston, Texas 77030, U.S.A
| |
Collapse
|
23
|
Ferber K, Archer KJ. Modeling discrete survival time using genomic feature data. Cancer Inform 2015; 14:37-43. [PMID: 25861216 PMCID: PMC4360712 DOI: 10.4137/cin.s17275] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2014] [Revised: 12/18/2014] [Accepted: 12/25/2014] [Indexed: 01/31/2023] Open
Abstract
Researchers have recently shown that penalized models perform well when applied to high-throughput genomic data. Previous researchers introduced the generalized monotone incremental forward stagewise (GMIFS) method for fitting overparameterized logistic regression models. The GMIFS method was subsequently extended by others for fitting several different logit link ordinal response models to high-throughput genomic data. In this study, we further extended the GMIFS method for ordinal response modeling using a complementary log-log link, which allows one to model discrete survival data. We applied our extension to a publicly available microarray gene expression dataset (GSE53733) with a discrete survival outcome. The dataset included 70 primary glioblastoma samples from patients of the German Glioma Network with long-, intermediate-, and short-term overall survival. We tested the performance of our method by examining the prediction accuracy of the fitted model. The method has been implemented as an addition to the ordinalgmifs package in the R programming environment.
Collapse
Affiliation(s)
- Kyle Ferber
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - Kellie J Archer
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| |
Collapse
|
24
|
Archer KJ, Hou J, Zhou Q, Ferber K, Layne JG, Gentry AE. ordinalgmifs: An R Package for Ordinal Regression in High-dimensional Data Settings. Cancer Inform 2014; 13:187-95. [PMID: 25574124 PMCID: PMC4266195 DOI: 10.4137/cin.s20806] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2014] [Revised: 11/06/2014] [Accepted: 11/08/2014] [Indexed: 11/24/2022] Open
Abstract
High-throughput genomic assays are performed using tissue samples with the goal of classifying the samples as normal < pre-malignant < malignant or by stage of cancer using a small set of molecular features. In such cases, molecular features monotonically associated with the ordinal response may be important to disease development; that is, an increase in the phenotypic level (stage of cancer) may be mechanistically linked through a monotonic association with gene expression or methylation levels. Though traditional ordinal response modeling methods exist, they assume independence among the predictor variables and require the number of samples (n) to exceed the number of covariates (P) included in the model. In this paper, we describe our ordinalgmifs R package, available from the Comprehensive R Archive Network, which can fit a variety of ordinal response models when the number of predictors (P) exceeds the sample size (n). R code illustrating usage is also provided.
Collapse
Affiliation(s)
- Kellie J Archer
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - Jiayi Hou
- Clinical and Translational Research Institute, University of California San Diego, San Diego, CA
| | - Qing Zhou
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - Kyle Ferber
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| | - John G Layne
- Center for the Study of Biological Complexity, Virginia Commonwealth University, Richmond, VA, USA
| | - Amanda E Gentry
- Department of Biostatistics, Virginia Commonwealth University, Richmond, VA, USA
| |
Collapse
|
25
|
Shen Y, Liao KP, Cai T. Sparse kernel machine regression for ordinal outcomes. Biometrics 2014; 71:63-70. [PMID: 25196727 DOI: 10.1111/biom.12223] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2013] [Revised: 07/01/2014] [Accepted: 07/01/2014] [Indexed: 12/24/2022]
Abstract
Ordinal outcomes arise frequently in clinical studies when each subject is assigned to a category and the categories have a natural order. Classification rules for ordinal outcomes may be developed with commonly used regression models such as the full continuation ratio (CR) model (fCR), which allows the covariate effects to differ across all continuation ratios, and the CR model with a proportional odds structure (pCR), which assumes the covariate effects to be constant across all continuation ratios. For settings where the covariate effects differ between some continuation ratios but not all, fitting either fCR or pCR may lead to suboptimal prediction performance. In addition, these standard models do not allow for nonlinear covariate effects. In this article, we propose a sparse CR kernel machine (KM) regression method for ordinal outcomes where we use the KM framework to incorporate nonlinearity and impose sparsity on the overall differences between the covariate effects of continuation ratios to control for overfitting. In addition, we provide data driven rule to select an optimal kernel to maximize the prediction accuracy. Simulation results show that our proposed procedures perform well under both linear and nonlinear settings, especially when the true underlying model is in-between fCR and pCR models. We apply our procedures to develop a prediction model for levels of anti-CCP among rheumatoid arthritis patients and demonstrate the advantage of our method over other commonly used methods.
Collapse
Affiliation(s)
- Yuanyuan Shen
- Department of Biostatistics, Harvard School of Public Health, Boston, Massachusetts 02115, U.S.A
| | - Katherine P Liao
- Brigham and Women's Hospital, Boston, Massachusetts 02115, U.S.A
| | - Tianxi Cai
- Department of Biostatistics, Harvard School of Public Health, Boston, Massachusetts 02115, U.S.A
| |
Collapse
|
26
|
Tran T, Luo W, Phung D, Harvey R, Berk M, Kennedy RL, Venkatesh S. Risk stratification using data from electronic medical records better predicts suicide risks than clinician assessments. BMC Psychiatry 2014; 14:76. [PMID: 24628849 PMCID: PMC3984680 DOI: 10.1186/1471-244x-14-76] [Citation(s) in RCA: 73] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/19/2014] [Accepted: 03/10/2014] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND To date, our ability to accurately identify patients at high risk from suicidal behaviour, and thus to target interventions, has been fairly limited. This study examined a large pool of factors that are potentially associated with suicide risk from the comprehensive electronic medical record (EMR) and to derive a predictive model for 1-6 month risk. METHODS 7,399 patients undergoing suicide risk assessment were followed up for 180 days. The dataset was divided into a derivation and validation cohorts of 4,911 and 2,488 respectively. Clinicians used an 18-point checklist of known risk factors to divide patients into low, medium, or high risk. Their predictive ability was compared with a risk stratification model derived from the EMR data. The model was based on the continuation-ratio ordinal regression method coupled with lasso (which stands for least absolute shrinkage and selection operator). RESULTS In the year prior to suicide assessment, 66.8% of patients attended the emergency department (ED) and 41.8% had at least one hospital admission. Administrative and demographic data, along with information on prior self-harm episodes, as well as mental and physical health diagnoses were predictive of high-risk suicidal behaviour. Clinicians using the 18-point checklist were relatively poor in predicting patients at high-risk in 3 months (AUC 0.58, 95% CIs: 0.50 - 0.66). The model derived EMR was superior (AUC 0.79, 95% CIs: 0.72 - 0.84). At specificity of 0.72 (95% CIs: 0.70-0.73) the EMR model had sensitivity of 0.70 (95% CIs: 0.56-0.83). CONCLUSION Predictive models applied to data from the EMR could improve risk stratification of patients presenting with potential suicidal behaviour. The predictive factors include known risks for suicide, but also other information relating to general health and health service utilisation.
Collapse
Affiliation(s)
- Truyen Tran
- Centre for Pattern Recognition and Data Analytics, Deakin University, Geelong 3220, Australia,Department of Computing, Curtin University, Bentley, Australia
| | - Wei Luo
- Centre for Pattern Recognition and Data Analytics, Deakin University, Geelong 3220, Australia
| | - Dinh Phung
- Centre for Pattern Recognition and Data Analytics, Deakin University, Geelong 3220, Australia
| | - Richard Harvey
- Mental Health Services, Barwon Health, Geelong, Australia,School of Medicine, Deakin University, Geelong, Australia
| | - Michael Berk
- School of Medicine, Deakin University, Geelong, Australia,Barwon Health, Geelong, Australia,Mental Health Research Institute, University of Melbourne, Parkville, Australia,Orygen Youth Health Research Centre, Parkville, Australia
| | - Richard Lee Kennedy
- School of Medicine, Deakin University, Geelong, Australia,Barwon Health, Geelong, Australia
| | - Svetha Venkatesh
- Centre for Pattern Recognition and Data Analytics, Deakin University, Geelong 3220, Australia.
| |
Collapse
|
27
|
Affiliation(s)
| | - Gerhard Tutz
- Department of Statistics; Ludwig-Maximilians-University Munich; Akademiestrasse 1 D-80799 Munich Germany
| |
Collapse
|
28
|
Chen CK. The classification of cancer stage microarray data. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2012; 108:1070-1077. [PMID: 22925656 DOI: 10.1016/j.cmpb.2012.07.001] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2012] [Revised: 06/20/2012] [Accepted: 07/17/2012] [Indexed: 06/01/2023]
Abstract
Correctly diagnosing the cancer stage is most important for selecting an appropriate cancer treatment option for a patient. Recent advances in microarray technology allow the cancer stage to be predicted using gene expression patterns. The cancer stage is in ordinal scale. In this paper, we employ strict ordinal regressions including cumulative logit model in traditional statistics with data dimensionality reduction, and distribution free approaches of large margin rank boundaries implemented by the support vector machine, as well as an ensemble ranking scheme to model the cancer stage using gene expression microarray data. Predictive genes included in models are selected by univariate feature ranking, and recursive feature elimination. We perform cross-validation experiments to assess and compare classification accuracies of ordinal and non-ordinal algorithms on five cancer stage microarray datasets. We conclude that a strict ordinal classifier trained by a validated approach can predict the cancer stage more accurately than traditional non-ordinal classifiers without considering the order of cancer stages.
Collapse
Affiliation(s)
- Chi-Kan Chen
- Department of Applied Mathematics, National Chung Hosing University, Taiwan.
| |
Collapse
|