loess.control(surface=c("interpolate", "direct"), statistics= c("approximate", "exact"), trace.hat=c("wait.to.decide", "approximate", "exact"), cell=0.2, iterations=4, method=c("loess", "model.frame"), ...)
loess(NOx ~ C * E, data = ethanol, span = 1/2, drop.square = "C", parametric = "C", control = loess.control(surface = "direct")) # The following has the same result: loess(NOx ~ C * E, data = ethanol, span = 1/2, drop.square = "C", parametric = "C", surface = "direct")