line <-
function(x, y=NULL)
{
	xy <- xy.coords(x, y)
	ok <- complete.cases(xy$x,xy$y)
	n <- length(ok)
	if(n <= 1) stop("insufficient observations")
	z <- .C("tukeyline",
		as.double(xy$x[ok]),
		as.double(xy$y[ok]),
		double(n),
		double(n),
		n,
		double(2),
		DUP=FALSE)
	value <- list(call=sys.call(), coefficients=z[[6]], residuals=z[[3]],
		fitted.values=xy$y-z[[3]])
	class(value) <- "tukeyline"
	value
}
coefficients.tukeyline <- coefficients.lm
residuals.tukeyline <- residuals.lm
fitted.values.tukeyline <- fitted.values.lm
print.tukeyline <- print.lm
