R/ryp.R

In asuR: Advanced statistics using R

ryp <- function(mymodel, id = c("all", "none"), ...){
### ================================= checking and argument matching
  if(class(mymodel)[1] != "lm"){stop("ERROR: the model you want to inspect should be fitted with lm")}
  id <- match.arg(id)
### ================================= calculations
  ri <- rstudent(mymodel)
  yi <- predict(mymodel)
### ================================= plot
  plot1 <- xyplot(ri ~ yi, ...,
                  ylab = "studentized residuals",
                  xlab = "predicted values",
                  panel = function(x,y,...){
                    panel.xyplot(x,y,...)
                    panel.abline(a=0)
                    panel.loess(x,y,...)
                  }
                  )
  print(plot1)
### ================================= IDENTIFICATION
  identifyControl(panel.matrix=trellis.currentLayout(),
                  original.row.names=dimnames(model.matrix(mymodel))[[1]],
                  id=id
                  )
}