The stability matters in clinical prediction models because it makes the model to be interpretable and generalizable. It is paramount for high dimensional data, which employ sparse models with feature selection ability. We propose a new method to stabilize sparse support vector machines using intrinsic graph structure of the electronic medical records. The graph structure is exploited using the Jaccard similarity among features. Our method employs a convex function to penalize the pairwise l ∞ -norm of connected feature coefficients in the graph. We apply the alternating direction method of multipliers to solve the proposed formulation. Our experiments are conducted on a synthetic and three real-world hospital datasets. We show that our proposed method is more stable than the state-of-the-art feature selection and classification techniques in terms of three stability measures namely, Jaccard similarity measure, Spearman's rank correlation coefficient and Kuncheva index. We further show that our method has resulted in better classification performance compared to the baselines.