We propose a probabilistic prediction method, based on local polynomial regressions, which complements the point forecasts with robust estimates of the corresponding forecast uncertainty. The reliability, practicability and generality of the method is demonstrated by applying it to astronomical, physiological, economic, and geophysical time series.