This paper is concerned with the computational estimation of the error of numerical solutions of potentially degenerate reaction-diffusion equations. The underlying motivation is a desire to compute accurate estimates as opposed to deriving inaccurate analytic upper bounds. The authors outline, analyze and test an approach to obtain computational error estimates based on the introduction of the residual error of the numerical solution and in which the effects of the accumulation of errors are estimated computationally.