We illustrate how the Hill relation and the notion of quasi-stationary distribution can be used to analyse the error introduced by many algorithms that have been proposed in the literature, in particular in molecular dynamics, to compute mean reaction times between metastable states for Markov processes. We present in particular how this analysis gives rigorous foundations to methods using splitting algorithms to sample the reactive trajectories.