Interval coefficients have been introduced in OR and CP to specify uncertain data in order to provide reliable solutions to convex models. The output is generally a solution set, guaranteed to contain all solutions possible under any realization of the data. This set can be too large to be meaningful. Furthermore, each solution has equal uncertainty weight, thus does not reflect any possible degree of knowledge about the data. To overcome these problems we propose to extend the notion of interval coefficient by introducing a second dimension to each interval bound. Each bound is now specified by its data value and its degree of knowledge. This is formalized using the cumulative distribution function of the data set. We define the formal framework of constraint reasoning over this cdf-intervals. The main contribution of this paper concerns the formal definition of a new interval arithmetic and its implementation. Promising results on problem instances demonstrate the approach.