par(mfrow=c(1,1)) plot(sspot,ty="l") readline() fitted = sspot-fit2$resid lines(fitted,col="red")