# NOT RUN {
xx <- rep(1:10,7)
yy <- rep(1:7, rep(10,7))
zz <- rnorm(70)
p.res.2x(xx,yy,zz, restricted = 2, main = "i.i.d. N(0,1) random residuals")
# }
# NOT RUN {
example(lm.influence, echo = FALSE)
op <- mult.fig(2, marP=c(-1,-1,-1,0), main="p.res.2x(*,*, residuals(lm.SR))")$old.par
with(LifeCycleSavings,
{ p.res.2x(pop15, ddpi, residuals(lm.SR), scol=c("red", "blue"))
p.res.2x(pop75, dpi, residuals(lm.SR), scol=2:1)
})
## with formula interface:
p.res.2x(~ pop15 + ddpi, lm.SR, scol=c("red", "blue"))
p.res.2x(~ pop75 + dpi, lm.SR, scol=2:1)
par(op) # revert par() settings above
# }
Run the code above in your browser using DataLab