The R Project SVN R

Rev

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