-
Notifications
You must be signed in to change notification settings - Fork 1
/
prostate.Rd
70 lines (61 loc) · 1.91 KB
/
prostate.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/data.R
\docType{data}
\name{prostate}
\alias{prostate}
\title{Prostate Cancer Data}
\format{
A data frame with 97 observations on the following 10 variables.
\describe{
\item{lcavol}{log cancer volume}
\item{lweight}{log prostate weight}
\item{age}{in years}
\item{lbph}{log of the amount of benign prostatic hyperplasia}
\item{svi}{seminal vesicle invasion}
\item{lcp}{log of capsular penetration}
\item{gleason}{a numeric vector}
\item{pgg45}{percent of Gleason score 4 or 5}
\item{lpsa}{response}
\item{train}{a logical vector}
}
}
\source{
Stamey, T., Kabalin, J., McNeal, J., Johnstone, I., Freiha, F.,
Redwine, E. and Yang, N (1989) Prostate specific antigen in the diagnosis
and treatment of adenocarcinoma of the prostate II. Radical prostatectomy
treated patients, \emph{Journal of Urology}, \bold{16}: 1076--1083.
}
\description{
Data to examine the correlation between the level of prostate-specific
antigen and a number of clinical measures in men who were about to receive a
radical prostatectomy.
}
\details{
This data set came originally from the (now defunct) ElemStatLearn package.
The last column indicates which 67 observations were used as the "training
set" and which 30 as the test set, as described on page 48 in the book.
}
\note{
There was an error in this dataset in earlier versions of the package,
as indicated in a footnote on page 3 of the second edition of the book. As
of version 2012.04-0 this was corrected.
}
\examples{
data(prostate)
str( prostate )
cor( prostate[,1:8] )
prostate <- prostate[, -10]
prostate.mod <- lm(lpsa ~ ., data=prostate)
vif(prostate.mod)
py <- prostate[, "lpsa"]
pX <- data.matrix(prostate[, 1:8])
pridge <- ridge(py, pX, df=8:1)
pridge
# univariate ridge trace plots
traceplot(pridge)
traceplot(pridge, X="df")
# bivariate ridge trace plots
plot(pridge)
pairs(pridge)
}
\keyword{datasets}