Biomedical Engineering Reference
In-Depth Information
h[K+1,K+1]=-sum(tal02*colSums((z^2)*Tnt*exp(tbeta2*z)*R))#beta2
h[K+2,K+2]=-sum(tal02*colSums((z^2)*Tt*exp(tbeta3*z)*R))#beta3
h[1:K,K+1]=h[K+1,1:K]=-colSums(z*Tnt*exp(tbeta2*z)*R)#al02*beta2
h[1:K,K+2]=h[K+2,1:K]=-colSums(z*Tt*exp(tbeta3*z)*R)#al02*beta3
h
}
#q(s)functionforE-step1###########################
q_s=function(ti,dij,de,d,interval,tal01,tal02,tal12)
{
Ji=which.max(dij)
#Ji=numberofintervalincludingt[i]
Ex=q=matrix(0,nrow=1,ncol=length(interval[,1]))
#qk(x,ti)=qvalueofintervalkfortimet[i](1xk)
Del=tal01+tal02-tal12#Delta
if(de==1)
{
for(kin(1:Ji))
{
if(k!=Ji)qtemp=exp(log(tal01[k])-tal12[Ji]*(ti-interval[Ji,1])+
((tal01[k]+tal02[k])*interval[k,1])-tal12[k]*interval[k,2])#
if(Ji>=3&k!=1&k<Ji)
qtemp=qtemp*exp(-sum(((tal01+tal02)[1:(k-1)])*
((interval[,2]-interval[,1])[1:(k-1)])))
if(Ji>=4&k<=Ji-2)qtemp=qtemp*exp(-sum(((tal12)[(k+1):(Ji-1)])*
((interval[,2]-interval[,1])[(k+1):(Ji-1)])))
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)
Search WWH ::
Custom Search