[CuriouslyC]

Gaussian processes are great for this. Another easy option is to bootstrap the estimator, though for large models that isn't feasible.