You need evaluate deparse(substitute(x)) _before_ you change x. so have
xlab; ylab
early in the body of your function.
However, NEWS in R-devel says
o scatter.smooth() and loess.smooth() now handle missing values in their inputs.
so why reinvent that wheel?
On Thu, 10 Feb 2005, William Briggs wrote:
I rewrote scatter.smooth to handle missing values, but I have a question about a move I had to make. Here's the code:
Mscatter.smooth<-function (x, y, span = 2/3, degree = 1, family = c("symmetric",
"gaussian"), xlab = deparse(substitute(x)), ylab = deparse(substitute(y)),
ylim = range(y, prediction$y), evaluation = 50, ...)
{
if (inherits(x, "formula")) {
if (length(x) < 3)
stop("need response in formula")
thiscall <- match.call()
thiscall$x <- x[[3]]
thiscall$y <- x[[2]]
return(invisible(eval(thiscall, sys.parent())))
}
##################
plot(x, y, xlab = xlab, ylab = ylab, ...)
i<-complete.cases(x,y)
x<-x[i]
y<-y[i]
##################
prediction <- loess.smooth(x, y, span, degree, family, evaluation)
lines(prediction)
invisible()
}
The changes I made are between the '###########'. The idea was only to pass complete cases on to loess.smooth and thus avoice the NA error you'd usually get by calling scatter.smooth with missing values.
I had to move the plot above the loess.smooth, whereas in the orginal scatter.smooth is was just below to take adavantage of prediction$y to set the ylim in the plotting function.
But if I try to plot after creating the complete cases, the names of the x and y variables disappear and instead the plotting labels on the x and y axis look like, for example, 'c(1,2,1.2,5.2, ...)', that is, the vector of values being plotted.
Why does that happen?
Because it is documented that way!
-- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595
______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
