ceeds p_target. The algorithm
begins by sampling uniformly from the unit simplices defining the parameter space, but
alternatives can be specified in <code>theta_sampler</code>. By default
gradient ascent (<code>ga = TRUE</code>) is performed during the p-value maximization
procedure, and <code>ga_gfactor</code> and <code>ga_lrate</code> control options for the gradient
ascent. At each iteration, the gradient of the multinomial probability at the current maximum
theta is computed, and a step is taken to <co