res1.std<-
  function(a, b = 3, d = 0.95) {
    c <- round((residuals(a) * sqrt(a$df.residual/deviance(a))), b)
    outlier <- as.character(ifelse(abs(c) > qnorm(d + (1 - d)/2), "***",
        " "))
    list(rs = c, o = outlier)
  }