We propose a method to overcome the usual limitation of data reduction pipelines and of the Optical Interferometry FITS format: the absence of correlated and non Gaussian errors. To do so, we use a bootstrap method that allows to sample the probability density function of the interferometric observables. We have applied the method to the accurate assessment of uncertainties on the stellar diameters of under-resolved stars.