Hi, How can I define a new psi function in the rlm comand? In particular I would like to implement the case for \rho(u)=0.5*(u^2) and \psi(u)=u in order to assume normally distributed errors. Any help? Thanks -- ======================================================= Angelo Secchi PGP Key ID:EA280337