faca por partes... conserte o codigo a seguir onde necessario (nao testei) set.seed(1) x <- rnorm(100) y <- rnorm(100) plot(x, y, type='n') args <- append(as.list(par('usr')), 'grey70') names(args) <- c('xleft', 'xright', 'ybottom', 'ytop', 'col') do.call(rect, args) points(x, y) --b