The R Project SVN R

Rev

Rev 85953 | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
42333 ripley 1 
% File src/library/stats/man/plot.ppr.Rd
68948 ripley 2 
% Part of the R package, https://www.R-project.org
70615 maechler 3 
% Copyright 1995-2016 R Core Team
42333 ripley 4 
% Distributed under GPL 2 or later
5 
 
27464 ripley 6 
% file modreg/man/plot.ppr.Rd
7 
% copyright (C) 1995-8 B. D. Ripley
8 
%
85953 hornik 9 
 
27464 ripley 10 
\name{plot.ppr}
56186 murdoch 11 
\alias{plot.ppr}
70615 maechler 12 
\title{Plot Ridge Functions for Projection Pursuit Regression Fit}
27464 ripley 13 
\description{
70615 maechler 14 
  Plot the ridge functions for a projection pursuit regression
15 
  (\code{\link{ppr}}) fit.
27464 ripley 16 
}
17 
\usage{
70615 maechler 18 
\method{plot}{ppr}(x, ask, type = "o", cex = 1/2,
19 
     main = quote(bquote(
20 
         "term"[.(i)]*":" ~~ hat(beta[.(i)]) == .(bet.i))),
21 
     xlab = quote(bquote(bold(alpha)[.(i)]^T * bold(x))),
22 
     ylab = "", \dots)
27464 ripley 23 
}
24 
\arguments{
70615 maechler 25 
  \item{x}{an \R object of class \code{"ppr"} as produced by a call to
26 
    \code{ppr}.}
27464 ripley 27 
  \item{ask}{
70615 maechler 28 
    the graphics parameter \code{ask}: see \code{\link{par}} for details.
27464 ripley 29 
    If set to \code{TRUE} will ask between the plot of each cross-section.
30 
  }
70615 maechler 31 
  \item{type}{the type of line (see \code{\link{plot.default}}) to draw.}
32 
  \item{cex}{plot symbol expansion factor (\emph{relative} to
33 
    \code{\link{par}("cex")}).}
34 
  \item{main, xlab, ylab}{axis annotations, see also
35 
    \code{\link{title}}.  Can be an expression (depending on \code{i} and
85953 hornik 36 
    \code{bet.i}), as by default which will be \I{\code{eval()}uated}.}
70615 maechler 37 
  \item{\dots}{further graphical parameters, passed to
38 
    \code{\link{plot}()}.}
27464 ripley 39 
}
40 
\value{
41 
  None
42 
}
43 
\section{Side Effects}{
44 
  A series of plots are drawn on the current graphical device, one for
45 
  each term in the fit.
46 
}
47 
\seealso{
48 
  \code{\link{ppr}}, \code{\link{par}}
49 
}
50 
\examples{
41508 ripley 51 
require(graphics)
52 
 
70646 maechler 53 
rock1 <- within(rock, { area1 <- area/10000; peri1 <- peri/10000 })
74363 maechler 54 
par(mfrow = c(3,2)) # maybe: , pty = "s"
27464 ripley 55 
rock.ppr <- ppr(log(perm) ~ area1 + peri1 + shape,
70646 maechler 56 
                data = rock1, nterms = 2, max.terms = 5)
61160 ripley 57 
plot(rock.ppr, main = "ppr(log(perm)~ ., nterms=2, max.terms=5)")
58 
plot(update(rock.ppr, bass = 5), main = "update(..., bass = 5)")
59 
plot(update(rock.ppr, sm.method = "gcv", gcvpen = 2),
27464 ripley 60 
     main = "update(..., sm.method=\"gcv\", gcvpen=2)")
61 
}
62 
\keyword{hplot}