Computational issues with fitting joint location/dispersion models in unreplicated 2k factorials

Maximum likelihood estimation for a joint location/dispersion model has been found occasionally to experience convergence problems when applied to experiments of the 2^k factorial series. We explore these problems and identify models for which the likelihood diverges or is multimodal. We derive the conditions under which this occurs and provide simple ways to check for problems both before and during computation.