XL<-data.frame(TBCplapyr)
p<-ncol(XL)
N<-nrow(XL)
censor<-matrix(XL$event)
especiales<-matrix(data=0,nrow(XL),3)
especiales[,1]<-matrix(XL$id)
especiales[,2]<-matrix(XL$Tcal)
especiales[,3]<-matrix(XL$event)
niveles<-levels(factor(especiales[,1]))
for(i in 1:N){
for(j in 1:nrow(matrix(niveles))){
if (as.character(especiales[i,1])==niveles[j]) especiales[i,1]<-j}}
StudyPeriod<-matrix(data=0,nrow(matrix(niveles)),1)
start<-matrix(data=0,nrow(matrix(niveles)),1)
k<-0
for(j in 1:N){if (XL$event[j]==0){k<-k+1;StudyPeriod[k,1]<-XL$Tcal[j]}}
units<-matrix(1:nrow(matrix(niveles)),nrow(matrix(niveles)),1)
Plot.Data.Events(XL,units,start,StudyPeriod,censor,
especiales,"black","blue")
Plot.Data.Events(XL,units,start,StudyPeriod,censor,
especiales,"red","black")
Run the code above in your browser using DataLab