Termplot: Produce a Term Plot Lattice

Description

Termplot produces a lattice plot of termplots. Terms are not plotted individually, rather the terms in which a variable appears are summed and plotted.
Usage

## S3 method for class 'default':
Termplot(object,
        ...,
        variables=NULL,
        col.term = 2,
        lty.term = 1,
        lwd.term = 1.5,
        se = TRUE,
        col.se = "orange",
        lty.se = 2,
        lwd.se = 1,
        col.res = "gray",
        residuals = c("deviance","none","pearson","working"),
        cex = 1,
        pch = 1,
        jitter.resid=FALSE,
        smooth = TRUE,
        col.smth = "darkred",
        lty.smth = 2,
        lwd.smth = 1,
        span.smth = 2/3,
        aspect="fill",
        xlab=NULL,
        ylab=NULL,
        main=paste(deparse(object$call),collapse=""),
        models=c("rows","columns"),
        xrot = 0,
        layout=NULL)
object{an model fit object, or a list of model objects}
  ...{further model objects.}
  variables{a character vector giving the names of 
    independent variables; note that
    the combined effect of all terms containing the respective
    values will be plotted; if empty, the effect of each independent
    variable will be plotted. Currently, higher-order terms
    will be ignored.}
  col.term,lty.term,lwd.term{parameters for the graphical representation
    of the terms with the same meaning as in termplot.}
  se{a logical value; should standard error curves added to the plot?}
  col.se,lty.se,lwd.se{graphical parameters for the
    depiction of the standard error curves, see termplot.}
  residuals{a character string, to select the type of
    residuals added to the plot.}
  col.res,cex,pch{graphical parameters for the depiction
        the residuals as in termplot.}
  jitter.resid{a logical vector of at most length 2.
      Determines whether residuals should be jittered
      along the x-axis (first element) and along the y-axis.
      If this argument has length 1, its setting applies to
      both axes.}
  smooth{a logical value; should a LOWESS smooth added to the plot?}
  span.smth{a numerical value, the span of the smoother.}
  col.smth,lty.smth,lwd.smth{graphical parameters for the smoother curve.}
  aspect{aspect ratio of the plot(s), see xyplot.}
  xlab{label of the x axis, see xyplot.}
  ylab{label of the y axis, see xyplot.}
  main{main heading, see xyplot.}
  models{character; should models arranged in rows or columns?}
  xrot{angle by which labels along the x-axis are rotated.}
  layout{layout specification, see xyplot.}
A trellis object.
library(lattice)
library(grid)
if(interactive())
  old.prompt <- devAskNewPage(TRUE)
lm0 <- lm(sr ~ pop15 + pop75,              data = LifeCycleSavings)
lm1 <- lm(sr ~                 dpi + ddpi, data = LifeCycleSavings)
lm2 <- lm(sr ~ pop15 + pop75 + dpi + ddpi, data = LifeCycleSavings)
berkeley <- Aggregate(Table(Admit,Freq)~.,data=UCBAdmissions)
berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial")
berk1 <- glm(cbind(Admitted,Rejected)~Gender,data=berkeley,family="binomial")
berk2 <- glm(cbind(Admitted,Rejected)~Gender+Dept,data=berkeley,family="binomial")
Termplot(lm2)
Termplot(berk2)
Termplot(lm0,lm1,lm2)
Termplot(berk0,berk1,berk2)
Termplot(By(~Gender,glm(cbind(Admitted,Rejected)~Dept,family="binomial"),
                    data=berkeley))
Termplot(By(~Dept,glm(cbind(Admitted,Rejected)~Gender,family="binomial"),
                    data=berkeley))
require(splines)
xyz <- data.frame(
  x = 1:100,
  z = factor(rep(LETTERS[1:4],25))
)
xyz <- within(xyz,
  y <- rnorm(100,sin(x/10)+x/50+as.numeric(z))
)
yxz.lin <- glm(y ~ x + z, data=xyz)
yxz.bs <- glm(y ~ bs(x,6) + z, data=xyz)
yxz.ns <- glm(y ~ ns(x,6) + z, data=xyz)
yxz.poly <- glm(y ~ poly(x,6) + z, data=xyz)
yxz.sincos <- glm(y ~ sin(x/10) + cos(x/10) + x + z, data=xyz)
# Terms containing
# the same variable are not plotted
# individually but their combined effect is plotted
#
Termplot(yxz.lin,yxz.bs,yxz.ns,yxz.poly,yxz.sincos,models="columns",
  span.smth=1/3)
Termplot(yxz.lin,yxz.bs,yxz.ns,yxz.poly,yxz.sincos,variables="x",
  span.smth=1/3)

if(interactive())
  devAskNewPage(old.prompt)
hplot
Description

Usage

Arguments
