p.res.2fact {sfsmisc} | R Documentation |
Plots a numeric (“residual like”) variable against two factor covariates, using boxplots.
p.res.2fact(x, y, z, restricted, notch = FALSE, xlab = NULL, ylab = NULL, main = NULL)
x,y |
two factors or numeric vectors giving the levels of factors. |
z |
numeric vector of same length as |
restricted |
positive value which truncates the size. The corresponding symbols are marked by stars. |
notch |
logical indicating if the boxplots should be notched, see
|
xlab,ylab |
axis labels, see |
main |
main title passed to |
if values are restricted, this make use of the auxiliar
function u.boxplot.x
.
Lorenz Gygax logyg@wild.unizh.ch and Martin Maechler, Jan.95;
starting from p.res.2x()
.
p.res.2x
, boxplot
,
plot.lm
, TA.plot
.
I <- 8; J <- 3; K <- 20 xx <- factor(rep(rep(1:I, rep(K,I)),J)) yy <- factor(rep(1:J, rep(I*K,J))) zz <- rt(I*J*K, df=5) #-- Student t with 5 d.f. p.res.2fact(xx,yy,zz, restr= 4, main= "i.i.d. t <- 5 random |.| <= 4") mtext("p.res.2fact(xx,yy,zz, restr= 4, ..)", line=2, adj=1, outer=TRUE, cex=1) ## Real data data(warpbreaks) (fm1 <- lm(breaks ~ wool*tension, data = warpbreaks)) with(warpbreaks, p.res.2fact(wool, tension, residuals(fm1)))