[en] Data-driven methodologies based on the restoring force method have been developed over the past few decades for building predictive reduced-order models (ROMs) of nonlinear dynamical systems. These methodologies involve fitting a polynomial expansion of the restoring force in the dominant state variables to observed states of the system. ROMs obtained in this way are usually prone to errors and uncertainties due to the approximate nature of the polynomial expansion and experimental limitations. We develop in this article a stochastic methodology that endows these errors and uncertainties with a probabilistic structure in order to obtain a quantitative description of the proximity between the ROM and the system that it purports to represent. Specifically, we propose an entropy maximization procedure for constructing a multi-variate probability distribution for the coefficients of power-series expansions of restoring forces. An illustration in stochastic aeroelastic stability analysis is provided to demonstrate the proposed framework.