Statistical Analysis of Dependent Current Status Data with Application to Tumorigenicity Experiments - Interval-Censored-Time-to-Event Data: Methods and Applications - page 133

if(k<=(Ji-1))
{qtemp1=qtemp*(1/Del[k])*exp(-Del[k]*interval[k,1])
qtemp2=qtemp*(1/Del[k])*exp(-Del[k]*interval[k,2])
q[k]=qtemp1-qtemp2
Ex[k]=qtemp1*interval[k,1]-qtemp2*
interval[k,2]+(1/Del[k])*(qtemp1-qtemp2)
}
if(k==Ji&Ji==1)
{qtemp1=exp(log(tal01[k])+(tal01[k]+
tal02[k])*interval[k,1]-tal12[k]*ti-Del[Ji]*
interval[Ji,1])*(1/Del[Ji])#
qtemp2=exp(log(tal01[k])+((tal01[k]+tal02[k])*
interval[k,1])-tal12[k]*ti-Del[Ji]*ti)*(1/Del[Ji])#
q[k]=qtemp1-qtemp2
Ex[k]=qtemp1*interval[Ji,1]-qtemp2*ti+(1/Del[Ji])*(qtemp1-qtemp2)
}
if(k==Ji&Ji!=1)
{qtemp1=exp(log(tal01[k])-sum(((tal01+tal02)[1:(Ji-1)])*
((interval[,2]-interval[,1])[1:(Ji-1)]))+#
((tal01[k]+tal02[k])*interval[k,1])-
tal12[k]*ti-Del[Ji]*interval[Ji,1])*(1/Del[Ji])
qtemp2=exp(log(tal01[k])-sum(((tal01+tal02)[1:(Ji-1)])*
((interval[,2]-interval[,1])[1:(Ji-1)]))+#
(tal01[k]+tal02[k])*interval[k,1]-tal12[k]*
ti-Del[Ji]*ti)*(1/Del[Ji])
q[k]=qtemp1-qtemp2
Ex[k]=qtemp1*interval[Ji,1]-qtemp2*ti+
(1/Del[Ji])*(qtemp1-qtemp2)
}
}
}
for(iin1:length(q)){if(is.nan(q[i])==TRUE|
is.na(q[i])==TRUE)q[i]=0;if(is.infinite(q[i])==TRUE)