File: debug1.s

package info (click to toggle)
survival 2.29-1
  • links: PTS
  • area: main
  • in suites: etch, etch-m68k
  • size: 3,204 kB
  • ctags: 1,077
  • sloc: asm: 8,713; ansic: 6,928; sh: 22; makefile: 2
file content (14 lines) | stat: -rw-r--r-- 465 bytes parent folder | download | duplicates (3) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
#
# Simplest printout to debug frailty models
#
fit1 <- survreg(Surv(time, status) ~ frailty(group, theta=.1, sparse=F), 
		data=leukemia)

fit2 <-survreg(Surv(time, status) ~ frailty(group, theta=.1, sparse=F),
	       data=leukemia, debug=2, init=c(fit1$coef, log(fit1$scale)),
	       iter=0)

fit3 <-survreg(Surv(time, status) ~ frailty(group, theta=.1, sparse=T),
	       data=leukemia, debug=2, 
	       init=c(fit1$coef[c(2,3,1)], log(fit1$scale)), iter=0)