Generate many, many samples (say, B in all) from a distribution with value \(\hat{\theta}\); and for each sample construct an estimate of \(\theta\), say \(\theta_j^*\), j = 1, 2, ..., B. The distribution of \(\Delta=\hat{\theta}-\theta_{0}\) is then approximated by that of \(\theta^*-\hat{\theta}\), the quantiles of which are used to form an approximate confidence interval.