The R Project SVN R

Rev

Rev 70615 | Rev 74363 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 70615 Rev 70646
Line 47... Line 47...
47 
  \code{\link{ppr}}, \code{\link{par}}
 47 
  \code{\link{ppr}}, \code{\link{par}}
48 
}
 48 
}
49 
\examples{
 49 
\examples{
50 
require(graphics)
 50 
require(graphics)
51 
 
 51 
 
52 
with(rock, {
 - 
 
53 
area1 <- area/10000; peri1 <- peri/10000
 52 
rock1 <- within(rock, { area1 <- area/10000; peri1 <- peri/10000 })
54 
par(mfrow = c(3,2)) # maybe: , pty = "s")
 53 
par(mfrow = c(3,2)) # maybe: , pty = "s")
55 
rock.ppr <- ppr(log(perm) ~ area1 + peri1 + shape,
 54 
rock.ppr <- ppr(log(perm) ~ area1 + peri1 + shape,
56 
                data = rock, nterms = 2, max.terms = 5)
 55 
                data = rock1, nterms = 2, max.terms = 5)
57 
plot(rock.ppr, main = "ppr(log(perm)~ ., nterms=2, max.terms=5)")
 56 
plot(rock.ppr, main = "ppr(log(perm)~ ., nterms=2, max.terms=5)")
58 
plot(update(rock.ppr, bass = 5), main = "update(..., bass = 5)")
 57 
plot(update(rock.ppr, bass = 5), main = "update(..., bass = 5)")
59 
plot(update(rock.ppr, sm.method = "gcv", gcvpen = 2),
 58 
plot(update(rock.ppr, sm.method = "gcv", gcvpen = 2),
60 
     main = "update(..., sm.method=\"gcv\", gcvpen=2)")
 59 
     main = "update(..., sm.method=\"gcv\", gcvpen=2)")
61 
})
 - 
 
62 
}
 60 
}
63 
\keyword{hplot}
 61 
\keyword{hplot}