For a lm
model, draws an inverse.response plot with the response Y on the
vertical axis and the fitted values Yhat
on the horizontal axis. Uses nls
to
estimate lambda in the function
Yhat = b0 + b1(Y)^(lambda).
Adds the fitted curve to the plot.
invResPlot
is an alias for inverseResponsePlot
.
1 2 3 4 5 6 7  inverseResponsePlot(model, lambda=c(1,0,1), robust=FALSE, xlab=NULL, ...)
## S3 method for class 'lm'
inverseResponsePlot(model, lambda=c(1,0,1), robust=FALSE,
xlab=NULL, labels=names(residuals(model)), ...)
invResPlot(model, ...)

model 
A 
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 
labels 
Case labels if labeling is turned on; see

robust 
If TRUE, then estimation uses Huber Mestimates with the median absolute deviation to estimate scale and k= 1.345. The default is FALSE. 
... 
Other arguments passed to 
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 Yhat on Y^(lambda), interpreting lambda = 0 to be natural logarithms.
Numeric output is a list with elements
lambda 
Estimate of transformation parameter for the response 
RSS 
The residual sum of squares at the minimum if robust=FALSE. If robust = TRUE, the value of Huber objective function is returned. 
Sanford Weisberg, sandy@umn.edu
Fox, J. and Weisberg, S. (2011) An R Companion to Applied Regression, Second Edition, Sage.
Pendergast, L, and Sheather, S. (in press). On sensitivity of response plot estimation of a robust estimation approach. Scandinavian Journal of Statistics.
Weisberg, S. (2014) Applied Linear Regression, Fourth Edition, Wiley, Chapter 7.
invTranPlot
, powerTransform
,
showLabels
1 2 
Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
All documentation is copyright its authors; we didn't write any of that.