attach(prostate)
library(Design,T)
library(Hmisc,T)
"%in%"<-function(a,b) match(a,b nomatch=0) > 0
subset<-prostate$status %in% c("dead - heart or vascular","dead -cerebrovascular","dead-prostatic ca")
attach(prostate[subset,])
levels(ekg)[levels(ekg) %in% c("old MI", "recent MI")]<-"MI"
pf.lin<-as.integer(pf)-1
levels(pf)<-levels(pf)[c(1:3,3)]
levels(pf)<-abbreviate(levels(pf))
w<-prostate.transcan[subset, ]
sz<-impute(w,sz)
sg<-impute(w,sg)
age<-impute(w,age)
wt<-impute(w,wt)
ekg<-impute(w,ekg)


cvd<-status %in% c("dead - heart or vascular","dead - cerebrovascular")
ekg.norm<-ekg %in% c("normal","benign")
levels(ekg)<-abbreviate(levels(ekg))
sz.t<-w[,"sz"]
sg.t<-w[,"sg"]
ap.t<-w[,"ap"]
sbp.t<-w[,"sbp"]
dbp.t<-w[,"dbp"]
age.t<-w[,"age"]
wt.t<-w[,"wt"]
hg.t<-w[,"hg"]
ekg.t<-w[,"ekg"]
pf.t<-w[,"pf"]