Most of the optimization is
Icens have a one dimensional root-finding component. Since the quantities involved are generally restricted to a subset of [0, 1] we use bisection to find the roots. This function is borrowed from package
We search from
pvec in the direction
ndir to obtain the new value of
pvec that maximizes the likelihood.
It is used in
ModifiedEMICMmac and not to be called by the user.
Alain Vandal and Robert Gentleman
Any book on optimization.