Previous work has been demonstrated that tracking features describing the dynamic and time-varying patterns in brain monitoring signals provide additional predictive information beyond that derived from static features based on snapshot measurements. To achieve more accurate predictions of outcomes of patients with traumatic brain injury (TBI), we proposed a statistical framework to extract dynamic features from brain monitoring signals based on the framework of Gaussian processes (GPs). GPs provide an explicit probabilistic, nonparametric Bayesian approach to metric regression problems. This not only provides probabilistic predictions, but also gives the ability to cope with missing data and infer model parameters such as those that control the function’s shape, noise level and dynamics of the signal. Through experimental evaluation, we have demonstrated that dynamic features extracted from GPs provide additional predictive information in addition to the features based on the pressure reactivity index (PRx). Significant improvements in patient outcome prediction were achieved by combining GP-based and PRx-based dynamic features. In particular, compared with the a baseline PRx-based model, the combined model achieved over 30 % improvement in prediction accuracy and sensitivity and over 20 % improvement in specificity and the area under the receiver operating characteristic curve.