Biomedical Engineering Reference
In-Depth Information
he_al02beta=function(param)
{
tal02=param[1:K];tbeta2=param[K+1];tbeta3=param[K+2];
h=matrix(0,ncol=length(param),nrow=length(param))
h[1:K,1:K]=-diag(colSums(dij*d)/tal02^2)#al02
h[K+1,K+1]=-sum(tal02*colSums((z^2)*Tnt*exp(tbeta2*z)*expR))#beta2
h[K+2,K+2]=-sum(tal02*colSums((z^2)*Tt*exp(tbeta3*z)*expR))#beta3
h[1:K,K+1]=h[K+1,1:K]=-colSums(z*Tnt*exp(tbeta2*z)*expR)#al02*beta2
h[1:K,K+2]=h[K+2,1:K]=-colSums(z*Tt*exp(tbeta3*z)*expR)#al02*beta3
h
}
#SEfor"sigma2"###########################################
he_sigma2=function(sigma2)
N/(2*sigma2^2)-sum(Rsq)/(sigma2*sigma2^2)
#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)])))
Search WWH ::
Custom Search