rent <- read.table(hh("datasets/rent.dat"), ## Weisberg's file alr162
col.names=c("rnt.alf","rnt.till",
"cow.dens","prop.past","lime"))
rent$lime <- factor(rent$lime, labels=c("no.lime","lime"))
rent$alf.till <- rent$rnt.alf / rent$rnt.till
rent.lm <- lm(rnt.alf ~ rnt.till + cow.dens + lime, data=rent)
rent$resid.rent <- resid(rent.lm)
xysplom(resid.rent ~ rnt.till + cow.dens | lime, data=rent,
layout=c(2,2))
xysplom(resid.rent ~ rnt.till + cow.dens | lime, data=rent,
layout=c(2,2),
xlab="", ylab="",
x.label="", y.label="",
group.label.side="",
par.strip.text=list(cex=1.2),
panel=panel.cartesian,
axis3.line=2.4,
scales=list(
relation="same",
alternating=FALSE, labels=FALSE, ticks=FALSE),
between=list(x=1, y=3))
xysplom(resid.rent ~ rnt.till + cow.dens | lime, data=rent,
layout=c(2,2),
xlab="", ylab="",
x.label="", y.label="",
group.label.side="",
par.strip.text=list(cex=1.2),
panel=panel.cartesian,
axis3.line=3.6,
scales=list(
relation="same",
alternating=FALSE, labels=FALSE, ticks=FALSE),
rescale=list(x=FALSE, y=FALSE),
between=list(x=1, y=3))
xysplom(resid.rent ~ rnt.till + cow.dens | lime, data=rent,
layout=c(2,2),
xlab="", ylab="",
x.label="", y.label="",
group.label.side="",
par.strip.text=list(cex=1.2),
panel=panel.cartesian,
axis3.line=3.6,
scales=list(
relation="free",
alternating=FALSE, labels=FALSE, ticks=FALSE),
between=list(x=1, y=3))
xysplom(resid.rent ~ rnt.till + cow.dens | lime, data=rent,
layout=c(2,2),
xlab="", ylab="",
y.label="resid",
group.label.side="top",
par.strip.text=list(cex=1.2),
panel=panel.cartesian,
axis3.line=3.6,
scales=list(alternating=FALSE, labels=FALSE, ticks=FALSE),
rescale=list(x=FALSE, y=FALSE),
between=list(x=1, y=5))
Run the code above in your browser using DataCamp Workspace