"controlBMAgamma0" <-
function(maxIter = .Machine$integer.max, eps = sqrt(.Machine$double.eps),
nEsteps = 1, transformation = function(x) x^(1/3),
inverseTransformation = function(x) x^3,
start = list(varCoefs = NULL, weights = NULL))
{
list(maxIter = maxIter, eps = eps, nEsteps = nEsteps,
transformation = transformation,
inverseTransformation = inverseTransformation,
start = start)
}