 is often the case
that sampling from the null parameter space is impossible because it is a set of
measure 0. While it may have measure 0, it is not empty, and will contain a finite
set of points. Thus you should provide the argument <code>theta_null_points</code> which is
a matrix where the rows contain the finite set (sometimes 1) of points
\(\theta\) such that \(\tau(\theta) = \psi_0\). There is also an argument called
<code>p_value_limits</code> that can be used to improve performance of confidence int