log(y + alpha) ~ x1 + x2 + ....
logtrans(object, ..., alpha = seq(0.5, 6, by = 0.25) - min(y),
plotit = <<see below>>, interp = <<see below>>,
xlab = "alpha", ylab = "log Likelihood")
y ~ x1 + x2 + .... The function is generic.
object is a formula, this argument may specify a data frame
as for
lm.
T
if a non-null device is currently active, else
F.)
T if plotting
is to be done and the number of real points is less than 100.)
plot.
plot.
x (for alpha) and
y (for the marginal log-likelihood
values).
Venables & Ripley, Chapter 6.
logtrans(Days ~ Age*Sex*Eth*Lrn, data = quine,
alpha = seq(0.75, 6.5, len=20), singular.ok = T)