Biomedical Engineering Reference
In-Depth Information
sigma2=sum(Rsq)/N
if(w==nMax){flag=1;error_flag=1}else#error
{
w=w+1#thenumberofEMiteration
conv[w,]=c(al01,al02,beta1,beta2,beta3,sigma2)
convLik[w,1]=logFlik(conv[w,])
if(max(abs(c(al01,al02,beta1,beta2,beta3)-
c(pal01,pal02,pbeta1,pbeta2,pbeta3)))<tol)flag=1
}
#if(as.integer(w/10)==w/10)
cat("EMestimation....iteration#",w,convLik[w,1],
round(conv[w,(K2+1):(K2+4)],4),"\n")
#forchecking
}
#ENDofEM
#Result
if(error_flag==1|w==nMax)print("Notcovergence")else
{
est=c(al01,al02,beta1,beta2,beta3,sigma2)
var_al01beta=
-diag(qr.solve(he_al01beta(c(al01,beta1)),tol=1e-17))
var_al02beta=
-diag(qr.solve(he_al02beta(c(al02,beta2,beta3)),tol=1e-17))
var_sigma2=-qr.solve(he_sigma2(sigma2),tol=1e-17)
se=sqrt(c(var_al01beta[1:K],var_al02beta[1:K],var_al01beta[K+1],
var_al02beta[(K+1):(K+2)],var_sigma2))
lower=est-(1.96*(se))
upper=est+(1.96*(se))
z_value=est/se
p_z=2*(1-pnorm(abs(z_value),0,1))
result=round(cbind(est,se,lower,upper,z_value,p_z),3)
rownames(result)=c(paste(
'
hazard.tumor(
'
,int1[-length(int1)],
'
,
'
,
int1[-1],
'
]
'
,sep=
''
),
Search WWH ::
Custom Search