The target of this paper is to demonstrate the benefits of using tolerance regions statistics in risk analysis. In particular, adopting the expected beta content tolerance regions as an alternative approach for choosing the optimal order of a response polynomial it is possible to improve results in reference class forecasting methodology. Reference class forecasting tries to predict the result of a planned action based on actual outcomes in a reference class of similar actions to that being forecast. Scientists/analysts do not usually work with a best fitting polynomial according to a prediction criterion. The present paper proposes an algorithm, which selects the best response polynomial, as far as a future prediction is concerned for reference class forecasting. The computational approach adopted is discussed with the help of an example of a relevant application.