chap6 <- read.table("hemophilia.data",header=T)
n <- nrow(chap6)
attach(chap6)
myHIV <- rep("HIVneg",n)
myHIV[hiv==2] <- "HIVpos"
mySev <- rep("high",n)
mySev[severity==2] <- "kmoderate"
mySev[severity==3] <- "low"
mySev[severity==4] <- "unknown"
mySev[severity==5] <- "none"
x <- data.frame(hiv=myHIV,sev=mySev,year=chap6$year,age=chap6$age,py=chap6$py,deaths=chap6$deaths)
detach("chap6")
attach(x)