car (version 2.0-11)

invResPlot: Inverse Response Plots to Transform the Response

Description

For a lm model, draws an inverse.response plot with the response $Y$ on the vertical axis and the fitted values $\hat{Y}$ on the horizontal axis. Uses nls to estimate $\lambda$ in the function $\hat{Y}=b_0+b_1Y^{\lambda}$. Adds the fitted curve to the plot. invResPlot is an alias for inverseResponsePlot.

Usage

inverseResponsePlot(model, lambda=c(-1,0,1), xlab=NULL, ...)

## S3 method for class 'lm':
inverseResponsePlot(model, lambda=c(-1,0,1), xlab=NULL,
   labels=names(residuals(model)), ...)

invResPlot(model, ...)

Arguments

model
A lm regression object
lambda
A vector of values for lambda. A plot will be produced with curves corresponding to these lambdas and to the least squares estimate of lambda
xlab
The horizontal axis label. If NULL, it is constructed by the function.
labels
Case labels if labeling is turned on; see invTranPlot and showLabels for arguments.
...
Other arguments passed to invTranPlot and then to plot.

Value

  • As a side effect, a plot is produced with the response on the horizontal axis and fitted values on the vertical axis. Several lines are added to be plot as the ols estimates of the regression of $\hat{Y}$ on $Y^{\lambda}$, interpreting $\lambda$ = 0 to be natural logarithms. Numeric output is a list with elements
  • lambdaEstimate of transformation parameter for the response
  • RSSThe residual sum of squares at the minimum

References

Fox, J. and Weisberg, S. (2011) An R Companion to Applied Regression, Second Edition, Sage. Weisberg, S. (2005) Applied Linear Regression, Third Edition, Wiley, Chapter 7.

See Also

invTranPlot, powerTransform, showLabels

Examples

Run this code
m2 <- lm(rate ~ log(len) + log(ADT) + slim + shld + log(sigs1), Highway1)
invResPlot(m2)

Run the code above in your browser using DataCamp Workspace