R/plot.mra.R

In zhangh12/iva: Instrumental Variable Analysis in Case-Control Association Studies

plot.mra <- function(x, ...){

  resid <- resid(x)
  fitted <- fitted(x)

  plot(fitted, resid, type = 'p',
       xlab = 'Fitted values', ylab = 'Residuals',
       main = 'Residuals vs Fitted')
  abline(h = 0, lty = 'dotted', col = 'grey')
  lines(lowess(fitted, resid), col = 'red')

}