Hello, Is the "std.dev" component of ls.diag( lsfit(x,y) ) the sample standard deviation of the residuals of the fit? I have ls.diag(lsfit(xx,yy))$std.dev different from sd(lsfit(xx,yy)$residuals) where xx and yy are vectors of 5 elements. Regards,