The R Project SVN R

Rev

Rev 14707 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 14707 Rev 15702
Line 11... Line 11...
11
         terms=NULL, se=FALSE, xlabs=NULL, ylabs=NULL, main = NULL,
11
         terms=NULL, se=FALSE, xlabs=NULL, ylabs=NULL, main = NULL,
12
         col.term = 2, lwd.term = 1.5,
12
         col.term = 2, lwd.term = 1.5,
13
         col.se = "orange", lty.se = 2, lwd.se = 2,
13
         col.se = "orange", lty.se = 2, lwd.se = 2,
14
         col.res= "gray", cex = 1, pch = par("pch"),
14
         col.res= "gray", cex = 1, pch = par("pch"),
15
         ask = interactive() && nb.fig < n.tms && .Device != "postscript",
15
         ask = interactive() && nb.fig < n.tms && .Device != "postscript",
16
         ...)
16
         \dots)
17
}
17
}
18
\arguments{
18
\arguments{
19
  \item{model}{fitted model object}
19
  \item{model}{fitted model object}
20
  \item{data}{data frame in which the variables in \code{model} can be found}
20
  \item{data}{data frame in which the variables in \code{model} can be found}
21
  \item{partial.resid}{logical; should partial residuals be plotted?}
21
  \item{partial.resid}{logical; should partial residuals be plotted?}
Line 62... Line 62...
62
z <- factor(rep(1:4,25))
62
z <- factor(rep(1:4,25))
63
y <- rnorm(100,sin(x/10)+as.numeric(z))
63
y <- rnorm(100,sin(x/10)+as.numeric(z))
64
model <- glm(y ~ ns(x,6) + z)
64
model <- glm(y ~ ns(x,6) + z)
65
 
65
 
66
par(mfrow=c(2,2)) ## 2 x 2 plots for same model :
66
par(mfrow=c(2,2)) ## 2 x 2 plots for same model :
67
termplot(model, main = paste("termplot( ", deparse(model$call)," ..)"))
67
termplot(model, main = paste("termplot( ", deparse(model$call)," ...)"))
68
termplot(model, rug=TRUE)
68
termplot(model, rug=TRUE)
69
termplot(model, partial=TRUE, rug= TRUE,
69
termplot(model, partial=TRUE, rug= TRUE,
70
         main="termplot(.., partial = TRUE, rug = TRUE)")
70
         main="termplot(..., partial = TRUE, rug = TRUE)")
71
termplot(model, partial=TRUE, se = TRUE, main = TRUE)
71
termplot(model, partial=TRUE, se = TRUE, main = TRUE)
72
if(rs) detach("package:splines")
72
if(rs) detach("package:splines")
73
}
73
}
74
\keyword{hplot}
74
\keyword{hplot}
75
\keyword{regression}
75
\keyword{regression}