We develop a new technique to integrate realistic uncertainties in probabilistic models of past sea-level change. The new framework performs better than past methods (in precision, accuracy, bias, and model fit) because it enables the incorporation of previously unused data and exploits correlations in the data. This method has the potential to assess the validity of past estimates of extreme sea-level rise and highstands providing better context in which to place current sea-level change.