Causal mechanism of relationship between the clinical outcome (efficacy or safety endpoints) and putative biomarkers, clinical baseline, and related predictors is usually unknown and must be deduced empirically from experimental data. Such relationships enable the development of tailored therapeutics and implementation of a precision medicine strategy in clinical trials to help stratify patients in terms of disease progression, clinical response, treatment differentiation, and so on. These relationships often require complex modeling to develop the prognostic and predictive signatures. For the purpose of easier interpretation and implementation in clinical practice, defining a multivariate biomarker signature in terms of thresholds (cutoffs/cut points) on individual biomarkers is preferable. In this paper, we propose some methods for developing such signatures in the context of continuous, binary and time‐to‐event endpoints. Results from simulations and case study illustration are also provided. Copyright © 2017 John Wiley & Sons, Ltd.