正文
fdata =7,]
fdata$age2 fdata$group 7) + 1*(fdata$flc.grp >9), levels=1:3, labels=c("FLC < 3.38", "3.38 - 4.71", "FLC > 4.71"))
sfit1 plot(sfit1, mark.time=F, col=c(1,2,4), lty=1, lwd=2, xscale=365.25, xlab="Years from Sample", ylab="Survival")
text(c(11.1, 10.5, 7.5)*365.25, c(.88, .57, .4),c("FLC < 3.38", "3.38 - 4.71", "FLC > 4.71"), col=c(1,2,4))
temp dd # Select subjects
set.seed(1978)
select for (j in 1:dd[2]) {
for (k in 1:dd[3]) {
n for (i in 1:2) {
indx select[i,j,k] temp[i,j,k])))
}
indx select[3,j,k] }
}
data2 sfit2 plot(sfit2,col=c(1,2,4), lty=1, lwd=2,xscale=365.25, xlab="Years from Sample", ylab="Survival")
lines(sfit1, col=c(1,2,4), lty=2, lwd=1,xscale=365.25)
legend(730, .4, levels(fdata$group), lty=1, col=c(1,2,4),bty='n', lwd=2)