Biomedical Engineering Reference
In-Depth Information
if(q[i]>0)q[i]=1e300elseq[i]=-1e300;}
if(de==0|sum(q)==0)p=rep(0,length(q))
elsep=(q/sum(q))[1,]#pj(ti)(1xJ)
Ejx=Ex/q#Ej(x|ti)(1xJ)
NT=de*p#E(NT)(1xJ)
Ttot=Tnt=matrix(0,nrow=1,ncol=length(interval[,1]))
Ttot=(1-cumsum(dij))*(interval[,2]-interval[,1])+
dij*(ti-interval[,1])#E(Ttot)(1xJ)
#E(Tt)(1xJ)
Tt=cumsum(p)*(1-cumsum(dij))*(interval[,2]-interval[,1])+
cumsum(p)*cumsum(dij)*dij*(ti-interval[,1])-p*(Ejx-interval[,1])[1,]
for(iin1:length(Tt))if(is.nan(Tt[i])==TRUE)Tt[i]=0
#E(Tnt)(1xJ)
Tnt=Ttot-Tt
result=list(p=p,Tnt=Tnt,Tt=Tt)
return(result)
}
############################################################
t=data$time#deathorsacrifiedtimeofindividuali
de=data$tumor#delta:indicatorwithtumoronset
d=data$death#d:indicatorwithdeath
z=data$trt#covariate
#Intervalandnumberofinterval(K)######################
int1=unique(c(0,int,max(t)));K=length(int1)-1
int=sort(rep(int1,2));int=int[2:(length(int)-1)]
############################################################
K2=2*K;N=length(t)
interval=matrix(int,nrow=K,ncol=2,byrow=TRUE)
dij=matrix(0,nrow=length(t),ncol=length(interval[,1]))#nxK
for(iin1:length(t))
for(jin1:length(interval[,1]))
if(interval[j,1]<t[i]&t[i]<=interval[j,2])
{dij[i,j]=1;break;}
 
Search WWH ::




Custom Search