Biomedical Engineering Reference
In-Depth Information
{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)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])+
Search WWH ::
Custom Search