plotbyfactor(x, y, f, data, col, pch, add, lg, ...)
data(iris) plotbyfactor(petal.wid,petal.len,species,data=iris)
locfit