File: stananim.q

package info (click to toggle)
sm 2.0.14-1
  • links: PTS
  • area: main
  • in suites: sarge
  • size: 992 kB
  • ctags: 26
  • sloc: fortran: 133; sh: 28; makefile: 13
file content (21 lines) | stat: -rw-r--r-- 609 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
provide.data(stanford)

x            <- Age
y            <- Log.time
status       <- Status
st.code      <- 1

hseq         <- seq(5, 10, by = 1)

model        <- sm.survival(x, y, status, h = hseq[1])
estimate.old <- model$estimate
for (i in (2:length(hseq))) {
  model <- sm.survival(x, y, status, h = hseq[i], 
  				display = "none")
  estimate.new <- model$estimate
  lines(model$eval.points, estimate.old, col = 0)
  lines(model$eval.points, estimate.new)
  text(x[status == st.code], y[status == st.code], "x")
  text(x[status != st.code], y[status != st.code], "o")
  estimate.old <- estimate.new
  }