Deep Evidential Regression
Abstract
Deterministic neural networks (NNs) are increasingly being deployed in safety critical domains, where calibrated, robust and efficient measures of uncertainty are crucial. While it is possible to train regression networks to output the parameters of a probability distribution by maximizing a Gaussian likelihood function, the resulting model remains oblivious to the underlying confidence of its predictions. In this paper, we propose a novel method for training deterministic NNs to not only estimate the desired target but also the associated evidence in support of that target. We accomplish this by placing evidential priors over our original Gaussian likelihood function and training our NN to infer the hyperparameters of our evidential distribution. We impose priors during training such that the model is penalized when its predicted evidence is not aligned with the correct output. Thus the model estimates not only the probabilistic mean and variance of our target but also the underlying uncertainty associated with each of those parameters. We observe that our evidential regression method learns wellcalibrated measures of uncertainty on various benchmarks, scales to complex computer vision tasks, and is robust to adversarial input perturbations.
1 Introduction
Recent advances in deep supervised learning have yielded super human level performance and precision (Liu et al., 2015; Gebru et al., 2017). While these models empirically generalize well when placed into new test enviornments, they are often easily fooled by adversarial perturbations (Goodfellow et al., 2014), and have difficulty understanding when their predictions should not be trusted. Today, regression based neural networks are being deployed in safety critical domains of computer vision (Godard et al., 2017; Alahi et al., 2016) as well as in robotics and control (Bojarski et al., 2016) where the ability to infer model uncertainty is crucial for eventual widescale adoption. Furthermore, precise uncertainty estimates are useful both for human interpretation of confidence and anomaly detection, and also for propagating these estimates to other autonomous components of a larger, connected system.
Existing approaches to uncertainty estimation are roughly split into two main categories: (1) learning aleatoric uncertainty (uncertainty in the data) and (2) epistemic uncertainty (uncertainty in the prediction). While representations for aleatoric uncertainty can be learned directly from data, approaches for estimating epistemic uncertainty primarily focus on placing probabilistic priors over all weights and sampling many times to obtain a measure of variance. In practice, many challenges arise with this approach, such as the computational expense of sampling during inference, how to pick an appropriate weight prior, or even how to learn such a representation given your prior.
We approach the problem of uncertainty estimation in regression from an evidential state of mind, where the model can acquire evidence during learning as it sees training examples. Every training example adds support to a learned higherorder, evidential distribution. Sampling from this evidential distribution yields instances of lowerorder, likelihood functions from which the data was drawn (cf. Fig. 1). We demonstrate that, by placing priors over our likelihood function (instead of all weights), we can learn a grounded representation of epistemic and aleatoric uncertainty that can be computed without sampling during inference.
In summary, this work makes the following contributions:

A novel and scalable method for learning representations of epistemic and aleatoric uncertainty, specifically on regression problems, by placing evidential priors over our likelihood function;

Evaluation of learned epistemic uncertainty on benchmark regression tasks and comparison against other stateoftheart uncertainty estimation techniques for neural networks;

Robustness evaluation against out of distribution and adversarially perturbed test data.
2 Modelling uncertainties from data
2.1 Preliminaries
Consider the following supervised optimization problem: given a dataset, , of paired training examples, , we aim to learn a function , parameterized by a set of weights, , which approximately solves the following optimization problem:
(1) 
where describes a loss function. In this work, we consider deterministic regression problems, which commonly optimize the sum of squared errors, . In doing so, the model is encouraged to learn the average correct answer for a given input, but does not explicitly model any underlying noise or uncertainty in the data when making its estimation.
2.2 Maximum likelihood estimation
Alternatively, we can approach our optimization problem from a maximum likelihood perspective, where we learn model parameters that maximize the likelihood of observing the particular set of training datapoints. In the context of deterministic regression, if we assume our targets, , were drawn i.i.d. from a Gaussian distribution with mean and variance parameters , then the likelihood of observing a single target, , can be expressed as
(2) 
In maximum likelihood estimation, we aim to learn a model to infer the that maximize the likelihood of observing our targets, . Equivalently, instead of maximizing the likelihood function, in practice it is common to instead minimize the negative log likelihood by setting
(3) 
In learning the parameters , this likelihood function allows us to successfully model the uncertainty of our data, also known as the aleatoric uncertainty. However, our model remains oblivious to the predictive model uncertainty. This metric, known as epistemic uncertainty, corresponds to the model’s uncertainty in its own output prediction (Kendall and Gal, 2017). In this paper, we present a novel approach for estimating the evidence in support of network predictions by directly learning both the inferred aleatoric uncertainty as well as the underlying epistemic uncertainty over its predictions. We achieve this by placing higherorder prior distributions over the learned parameters governing the distribution from which our observations are drawn.
3 Evidential uncertainty for regression
3.1 Problem setup
We consider the problem where our observed targets, , are drawn i.i.d. from a Gaussian distribution with unknown mean and variance , which we seek to probabilistically estimate. We model this by placing a conjugate prior distribution on . If we assume our observations are drawn from a Gaussian, this leads to placing a Gaussian prior on our unknown mean and an InverseGamma prior on our unknown variance:
(4)  
where , , , .
From a variational Bayesian perspective, our aim is to estimate a posterior distribution of the parameters . To obtain an approximation for the true posterior, we assume that the estimated distribution can be factorized into independent factors such that . In this case, the true distribution takes the form of a Normal InverseGamma (N.I.G.) distribution:
(5) 
The mean of this distribution can be interpreted as being estimated from observations with sample mean while its variance was estimated from observations with sample mean and sum of squared deviations . We denote the total evidence as the sum of all inferred observations counts, .
Thus, we can interpret the estimated posterior as an evidential, higherorder probability distribution on top of the unknown lowerorder likelihood distribution from which observations are drawn. Drawing a single sample from our evidential posterior yields representing a single instance of our likelihood function, namely . Thus, the parameters of the posterior, specifically (), determine not only the location but also the dispersion concentrations, or uncertainty, associated with our inferred likelihood function.
For example, in Fig. 2A we visualize the posterior of different evidential N.I.G. distributions with varying model parameters. We illustrate that by increasing the evidential parameters (i.e. ) of this distribution, the p.d.f. becomes more tightly concentrated about its inferred likelihood function. Considering a single parameter realization of this higherorder distribution, cf. Fig. 2B, we can subsequently sample many lowerorder realizations of our likelihood function, as shown in Fig. 2C.
In this work, we use neural networks to learn a higherorder evidential distribution that directly captures prediction uncertainties and evaluate our method on various regression tasks. This approach presents several distinct advantages. First, this method enables simultaneous learning of the desired regression task, along with uncertainty estimation, built in, due to our evidential priors. Second, since the evidential prior is a higherorder N.I.G. distribution, the maximum likelihood Gaussian can be computed analytically from the expected values of the parameters, without the need for sampling. Third, by explicitly modeling the evidence, we effectively capture the epistemic or model uncertainty associated with the network’s prediction. This can be done by simply evaluating the variance of our inferred evidential distribution.
3.2 Learning the evidential distribution
Having formalized the problem of using an evidential distribution to capture model uncertainty, we next describe our approach for actually learning this distribution. Given a set of observations, variational inference methods aim to approximate a posterior distribution over unobserved variables or parameters by maximizing the evidence lower bound (ELBO) (Kingma and Welling, 2013). Similarly, here we seek to estimate the posterior distribution governed by the higherorder distribution parameters to maximize the likelihood of our observations. Applying the principle of variational inference, we have:
(6) 
where . Similar to the principle of variational inference, in the remainder of this section we will discuss how we learn evidential distributions for regression by maximizing the loglikelihood of model evidence and minimizing the distance to an uncertainty prior. As we will see, maximizing the loglikelihood allows our model to fit the data, while the regularization provides an “uncertainty” penalty so the model can express when it does not know the answer.
We define the “model evidence” as the likelihood of an observation, , given the evidential distribution parameters , as . We apply Bayes’ theorem and marginalize over the likelihood parameters to obtain an equation for the model evidence:
(7) 
The model evidence is not, in general, straightforward to evaluate since computing it involves integrating out the dependence on model parameters. However, by placing a N.I.G. prior on our Gaussian likelihood function an analytical solution does exist.
(8) 
For computational reasons it is common to instead minimize the negative loglikelihood of the model evidence (). For a complete derivation please refer to the appendix.
(9) 
Alternatively, we can also derive the negative log likelihood of model evidence from the sumofsquare deviations to compute :
(10)  
(11) 
In our experiments, using resulted in greater training stability and increased performance, as opposed to the loss. Therefore, the remainder of this paper present results using .
3.3 Expressing “I don’t know”
In the previous subsection, we outlined a loss function for training a NN to output parameters of a N.I.G. distribution which maximize the loglikelihood of our data. In this subsection we describe how we regularize training against a prior where the model does not have any evidence (i.e. maximum uncertainty). In variational inference this is done by minimizing the KLdivergence between the inferred posterior, , and a prior, , cf. Eq. 6. In the evidential setting, our prior is also a Normal InverseGamma distribution, but with zero evidence (or infinite uncertainty). Therefore, during training we aim to minimize our evidence (or maximize our uncertainty) everywhere except where we have training data, as enforced by the negative loglikelihood loss term.
Unfortunately, the KLdivergence between an arbitrary N.I.G. distribution and another with infinitely low evidence is not well defined (Soch and Allefeld, 2016). To address this, we formulate a custom evidence regularizer, , based on the error of the th prediction,
where represents the Lp norm of .
This regularization loss imposes a penalty whenever there is an error in the prediction that scales with the total evidence of our inferred posterior. Conversely, large amounts of predicted evidence will not be penalized as long as the prediction is close to the target observation.
The combined loss function employed during training consists of the two loss terms for maximizing model evidence and regularizing evidence,
(12) 
3.4 Evaluating aleatoric and epistemic uncertainty
The aleatoric uncertainty, also referred to as statistical or data uncertainty, is representative of unknowns that differ each time we run the same experiment. We evaluate the aleatoric uncertainty from . The epistemic, also known as the model uncertainty, describes the estimated uncertainty in the learned model and is defined as , based on the N.I.G. definition.
4 Experiments
4.1 Predictive accuracy and uncertainty benchmarking
We first qualitatively compare the performance of our approach against a set of benchmarks on a onedimensional toy regression dataset. The training set consists of training examples drawn from , where in the region , whereas the test data is unbounded. Not only determinisitic or maximum likelihood regression, but also techniques using empirical variance of the networks’ predictions such as MCdropout, modelensembles, and BayesbyBackprop underestimate the uncertainty outside the training distribution. In contrast, evidential regression estimates uncertainty appropriately and grows the uncertainty estimate with increasing distance from the training data (Figure 3).
Additionally, we compare our approach to stateoftheart methods for predictive uncertainty estimation using NNs on common real world datasets used in (HernándezLobato and Adams, 2015; Lakshminarayanan et al., 2017; Gal and Ghahramani, 2016). We evaluate our proposed evidential regression method against modelensembles and BBB based on root mean squared error (RMSE), and negative loglikelihood (NLL). We do not provide results for MCdropout since it consistently performed inferior to the other baselines. The results in Table 1 indicate that although the loss function for evidential regression is more complex than competing approaches, it is the top performer in RMSE and NLL in 8 out of 9 datasets.
RMSE  NLL  

Dataset  Ensembles  BBB  Evidential  Ensembles  BBB  Evidential 
Boston  0.09 4.3e4  0.09 3.7e4  0.09 1.0e6  0.89 6.5e2  0.67 1.5e2  0.87 2.2e2 
Concrete  0.07 4.4e3  0.06 3.3e6  0.06 7.0e7  1.29 4.1e2  1.32 4.3e3  1.31 1.9e2 
Energy  0.10 2.3e4  0.10 1.6e5  0.10 9.0e7  0.61 8.9e2  0.60 2.0e2  0.75 1.4e2 
Kin8nm  0.07 3.5e4  0.17 3.5e4  0.08 3.8e3  0.78 1.4e2  0.32 6.3e3  1.17 2.6e2 
Naval  0.01 1.0e7  0.04 1.2e2  0.01 3.4e4  2.55 3.3e2  1.83 2.4e1  3.17 2.1e3 
Power  0.06 4.0e7  0.06 2.3e6  0.06 5.3e6  1.29 6.9e2  1.33 2.5e3  1.40 6.2e3 
Protein  0.17 1.0e6  0.17 8.0e4  0.17 1.6e6  0.27 6.7e2  0.32 5.9e2  0.29 1.1e2 
Wine  0.10 3.0e4  0.10 2.9e4  0.10 3.8e5  0.46 2.5e1  0.89 2.4e3  0.85 6.9e3 
Yacht  0.07 1.3e3  0.07 3.4e3  0.06 6.2e5  1.16 6.3e2  0.74 5.8e2  1.28 9.4e3 
4.2 Depth estimation
After establishing benchmark comparison results, in this subsection we demonstrate the scalability of our evidential learning by extending to the complex, highdimensional task of depth estimation. Monocular endtoend depth estimation is a central problem in computer vision which aims to learn a representation of depth directly from an RGB image of the scene. This is a challenging learning task since the output target is very highdimensional, , where are the height and width of the input image respectively. For every pixel in the image we regress over the desired depth and simultaneously want to estimate the uncertainty associated to that individual pixel estimate.
Our training data consists of over 27k RGBtodepth pairs of indoor scenes (e.g. kitchen, bedroom, etc.) from the NYU Depth v2 dataset (Nathan Silberman and Fergus, 2012). We train a UNet style NN (Ronneberger et al., 2015) for inference. Spatial dropout (Tompson et al., 2015) (with ) is used for the dropout baseline. The final layer of our model outputs a single activation map in the case of deterministic regression, dropout, ensembling and BBB. However, for our evidential model, we infer four outputs, each corresponding to respectively.
# Parameters  Inference Speed  RMSE  NLL  

Absolute  Relative  Seconds  Relative  
Evidential (Ours)  7,846,776  1  0.013  1  0.02 0.04  1.05 0.35 
Spatial Dropout  7,846,657  0.99  0.093  7.21  0.03 0.03  1.22 0.46 
Ensembles  39,233,285  4.99  0.071  5.49  0.03 0.03  0.99 0.28 
BBB  15,693,314  1.99  0.102  7.84     
Table 2 summarizes the size and speed of all models. Evidential models contain significantly fewer trainable parameters than ensembles (where the number of parameters scales linearly with the size of the ensemble). BBB maintains a trainable mean and variance for every weight in the network, so its size is roughly larger as well. The number of trainable parameters for evidential regression is closest to that of dropout, which has fewer as it contains a slightly smaller final output layer. Since evidential regression models do not require sampling in order to estimate their uncertainty, their forwardpass inference times are also significantly more efficient when compared to the baselines. Finally, we demonstrate that we achieve comparable predictive accuracy (through RMSE and NLL) to the other models. Note that the output size of the depth estimation problem presented significant learning challenges for the BBB baseline, and it was unable to converge during training. As a result, for the remainder of this analysis we compare against only spatial dropout and ensembles.
We evaluate these models in terms of both their accuracy and their predictive uncertainty on previously unseen test set examples. Fig. 4AC visualizes the predicted depth, absolute error from ground truth, and predictive uncertainty across three randomly picked test images. Ideally, a strong predictive uncertainty would capture any errors in the prediction (i.e., roughly correspond to where the model is making errors). We note that, compared to dropout and ensembling approaches, evidential uncertainty modeling captures the depth errors while providing clear and localized predictions of confidence, cf. Fig. 4. In general, dropout drastically underestimates the amount of uncertainty present, while ensembling occasionally overestimates the uncertainty, cf. Fig. 4A,C.
To evaluate how calibrated the predictive uncertainty was to the groundtruth errors, we fit receiver operating characteristic (ROC) curves to normalized estimates of error and uncertainty. Thus, we test the network’s ability to detect how likely it is to make an error at a given pixel using its predictive uncertainty. ROC curves take into account sensitivity and specificity of the uncertainties towards error predictions and are stronger if they contain greater area under their curve (AUC). Fig. 4D demonstrates that our evidential model provides uncertainty estimates which are the most attuned to where the model is making the errors.
4.3 Robustness to adversarial samples
A key use case of uncertainty estimation is to understand when a model is faced with test examples that fall outside of its training distribution or when the model’s output cannot be trusted. In the previous subsection, we showed that our evidential uncertainties were well calibrated with the model’s errors. In this subsection, we evaluate the uncertainty response for the depth estimation task under the extreme case where our model is presented with adversarially perturbed inputs.
We compute adversarial perturbations to our test set using the fast gradient sign method (Goodfellow et al., 2014), with increasing scales, , of noise. Fig. 5A confirms that the absolute error of all methods increasing as adversarial noise is added. We also observe a positive effect noise on our predictive uncertainty estimates in Fig. 5B. An additional desirable property of evidential uncertainty modeling is that it presents a higher overall uncertainty when presented with adversarial inputs compared to dropout and ensembling methods. Furthermore, we observe this strong overall uncertainty estimation despite the model losing calibration accuracy from the adversarial examples (Fig. 5C).
The robustness of evidential uncertainty against adversarial perturbations is visualized in greater detail in Fig. 6, which illustrates the predicted depth, error, and estimated pixelwise uncertainty as we perturb the input image with greater amounts of noise (lefttoright). Note that the predictive uncertainty not only steadily increases as we increase the noise, but the spatial concentrations of uncertainty throughout the image maintain tight correspondence with the error.
5 Discussion and Related work
Uncertainty estimation has a long history in neural networks, from modeling probability distribution parameters over outputs (Bishop, 1994) to Bayesian deep learning (Kendall and Gal, 2017). Our work builds on this foundation and presents a scalable representation for inferring the parameters of an evidential uncertainty distribution while simultaneously learning regression tasks via MLE.
In Bayesian deep learning, priors are placed over network weights and estimated using variational inference (Kingma et al., 2015). Dropout (Gal and Ghahramani, 2016; Molchanov et al., 2017) and BayesbyBackprop (Blundell et al., 2015) rely on multiple sampling iterations to estimate a predictive variance. Ensembles (Lakshminarayanan et al., 2017) provide a tangential approach where sampling occurs over multiple trained instances of the model. In contrast, we place uncertainty priors directly over our likelihood output function and thus only a single forward pass to evaluate both prediction and uncertainty. Additionally, our approach of uncertainty estimation proved to be better calibrated and capable of predicting where the model fails.
A large topic of research in Bayesian inference focuses on placing prior distributions over hierarchical models to estimate uncertainty (Gelman and others, 2006; Gelman et al., 2008). Our methodology falls under the class of evidential deep learning which leverages the Theory of Evidence to model prior distributions over neural network predictions and interpret uncertainty. Prior works in this field (Sensoy et al., 2018; Malinin and Gales, 2018) have focused exclusively on modeling uncertainty in the classification domain with Dirichlet prior distributions. Our work extends this field into the broad range of regression learning tasks and demonstrates generalizability to outofdistribution test samples.
6 Conclusion
In this paper, we develop a novel method for training deterministic NNs that both estimates a desired target and evaluates the evidence in support of the target to generate robust metrics of model uncertainty. We formalize this in terms of learning evidential distributions, and achieve stable training by penalizing our model for prediction errors that scale with the available evidence. Our approach for evidential regression is validated on a benchmark regression task. We further demonstrate that this method robustly scales to a key task in computer vision, depth estimation, and that the predictive uncertainty increases with increasing outofdistribution adversarial perturbation. This framework for evidential representation learning provides a means to achieve the precise uncertainty metrics required for robust neural network deployment in safetycritical domains.
References
 Social lstm: human trajectory prediction in crowded spaces. In Proceedings of the IEEE conference on computer vision and pattern recognition, pp. 961–971. Cited by: §1.
 Mixture density networks. In Tech. Rep. NCRG/94/004, Neural Computing Research Group, Cited by: §5.
 Weight uncertainty in neural networks. arXiv preprint arXiv:1505.05424. Cited by: §5.
 End to end learning for selfdriving cars. arXiv preprint arXiv:1604.07316. Cited by: §1.
 Dropout as a bayesian approximation: representing model uncertainty in deep learning. In international conference on machine learning, pp. 1050–1059. Cited by: §4.1, §5.
 Using deep learning and google street view to estimate the demographic makeup of neighborhoods across the united states. Proceedings of the National Academy of Sciences 114 (50), pp. 13108–13113. Cited by: §1.
 A weakly informative default prior distribution for logistic and other regression models. The Annals of Applied Statistics 2 (4), pp. 1360–1383. Cited by: §5.
 Prior distributions for variance parameters in hierarchical models (comment on article by browne and draper). Bayesian analysis 1 (3), pp. 515–534. Cited by: §5.
 Unsupervised monocular depth estimation with leftright consistency. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pp. 270–279. Cited by: §1.
 Explaining and harnessing adversarial examples. arXiv preprint arXiv:1412.6572. Cited by: §1, §4.3.
 Practical variational inference for neural networks. In Advances in neural information processing systems, pp. 2348–2356. Cited by: Table 1.
 Probabilistic backpropagation for scalable learning of bayesian neural networks. In International Conference on Machine Learning, pp. 1861–1869. Cited by: §4.1.
 What uncertainties do we need in bayesian deep learning for computer vision?. In Advances in neural information processing systems, pp. 5574–5584. Cited by: §2.2, §5.
 Autoencoding variational bayes. arXiv preprint arXiv:1312.6114. Cited by: §3.2.
 Variational dropout and the local reparameterization trick. In Advances in Neural Information Processing Systems, pp. 2575–2583. Cited by: §5.
 Simple and scalable predictive uncertainty estimation using deep ensembles. In Advances in Neural Information Processing Systems, pp. 6402–6413. Cited by: §4.1, Table 1, §5.
 Targeting ultimate accuracy: face recognition via deep embedding. arXiv preprint arXiv:1506.07310. Cited by: §1.
 Predictive uncertainty estimation via prior networks. In Advances in Neural Information Processing Systems, pp. 7047–7058. Cited by: §5.
 Variational dropout sparsifies deep neural networks. In Proceedings of the 34th International Conference on Machine LearningVolume 70, pp. 2498–2507. Cited by: §5.
 Indoor segmentation and support inference from rgbd images. In ECCV, Cited by: §4.2.
 Unet: convolutional networks for biomedical image segmentation. In International Conference on Medical image computing and computerassisted intervention, pp. 234–241. Cited by: §4.2.
 Evidential deep learning to quantify classification uncertainty. In Advances in Neural Information Processing Systems, pp. 3179–3189. Cited by: §5.
 Kullbackleibler divergence for the normalgamma distribution. arXiv preprint arXiv:1611.01437. Cited by: §3.3.
 Efficient object localization using convolutional networks. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pp. 648–656. Cited by: §4.2.
7 Appendix
7.1 Model evidence derivations
For convenience, define be the precision of a Gaussian distribution.
7.1.1 Type II Maximum Likelihood Loss
(13)  
(14)  
(15)  
(16) 
For computational reasons it is common to instead minimize the negative logarithm of the model evidence.
(17) 
7.1.2 Sum of Squares Loss
(18)  
(19)  
(20)  
(21)  
(22) 
(23) 