Plotting a 3D surface in R

I have a lm of the form lm(log(z) ~ x*y). I can make a 3D scatterplot of the data using s3d <- scatterplot3d(x,y,log(z)). How do I plot the surface of my regression function? I know I can do ...