trans.plot: Diagnostic plot to find power transform to remove...

Usage Arguments Value See Also Examples

Usage

1
2
trans.plot(fit, predictors=all.predictors(fit)[1:2],
  terms=predict(fit, type = "terms"), xlab, ylab, ...)

Arguments

fit

linear model (lm) fit

predictors

predictors to compare. They must be in the model, as well as theirinteraction in the specified order.

terms

estimates of model terms

xlab,ylab

x and y axis labels

...

Plot parameters

Value

Power = 1 - slope estimate. SE comes from simple linear regression and is at best advisory. Consider tranformation $y^p$ = exp(p*log(y)), or transform log(y) if p=0. In practice, examine nearest half integer.

See Also

effect.plot.

Examples

1
2
3
4
5
6
7
8
9
## Not run: 
# Bekk data from Hoaglin, Mosteller and Tukey (1991)
bekk.fit <- lm(smooth~lab*mat, bekk)
effect.plot(bekk.fit)trans.plot(bekk.fit)
bekk.lfit <- lm(log(smooth)~lab*mat, bekk)
effect.plot(bekk.lfit)
trans.plot(bekk.lfit)

## End(Not run)

byandell/pda documentation built on May 13, 2019, 9:27 a.m.