# An altered log transformation
warp.lm1 <- lm(log(breaks + 1) ~ wool*tension, data = warpbreaks)
rg1 <- update(ref.grid(warp.lm1),
tran = list(linkinv = function(eta) exp(eta) - 1,
mu.eta = function(eta) exp(eta)),
inv.lbl = "pred.breaks")
summary(rg1, type = "response")
lsm.options(ref.grid = list(level = .90),
contrast = list(infer = c(TRUE,FALSE)),
estble.tol = 1e-6)
# Sets default confidence level to .90 for objects created by ref.grid
# AS WELL AS lsmeans called with a model object (since it creates a
# reference grid). In addition, when we call 'contrast', 'pairs', etc.,
# confidence intervals rather than tests are displayed by default.
lsm.options(disable.pbkrtest = TRUE)
# This forces use of asymptotic methods for lmerMod objects.
# Set to FALSE or NULL to re-enable using pbkrtest.
print(lsm.options()) # see the current settings
Run the code above in your browser using DataLab