Biomedical Engineering Reference
In-Depth Information
############################################################
#initialvalue:constant#################################
beta1=init[1];beta2=init[1];beta3=init[1];sigma2=init[3]
al01=init[2]*1:K;al02=init[2]*1:K
R=rep(0,N);expR=rep(1,N);Rsq=rep(0,N);
#ConditionforEMinteration##############################
nMax=1000#maximumEMiteration
#tol#convergencecriteria
#ConditionforE-step2###################################
GH_n=9;
Xj=c(3.190993,2.266581,1.468553,0.723551,0.000000,
-0.723551,-1.468553,-2.266581,-3.190993);
Wj=c(0.000040,0.004944,0.088475,0.432652,0.720235,
0.432652,0.088475,0.004944,0.000040);
L=f0=f1=f2=f3=rep(NA,GH_n)
#Memorystorage
p=Tnt=Tt=matrix(NA,ncol=K,nrow=N)
conv=matrix(NA,nrow=nMax,ncol=K2+4)
convLik=matrix(NA,nrow=nMax,ncol=1)
#E-M######################################################
w=0;flag=0;#w:countofEMiteration,flag:stopindicator
error_flag=0;#error_flag:errorindicator
#cat("LLbeta1beta2beta3sigma2","\n")
while(w<nMax&flag==0)
{
#####previousdatasaving#################################
pal01=al01;pal02=al02;pbeta1=beta1;pbeta2=beta2;pbeta3=beta3;
psigma2=sigma2
#####Step:E-step1:E[Nt]E[Nnt]process#################
for(iin1:length(t))
{
y=q_s(t[i],dij[i,],de[i],d[i],interval,al01*exp(beta1*z[i])*
expR[i],al02*exp(beta2*z[i])*expR[i],al02*exp(beta3*z[i])*expR[i])
Search WWH ::
Custom Search