Full_Name: Bin Nan Version: 1.3.0 OS: Win32 Submission from: (NULL) (141.211.15.110) The function resid(fit, "pearson") seems not giving the same Pearson's residuals for logistic regression as what Splus does. I found the problem when I fitted glm with family = binomial using Crowder's seed data.