Biomedical Engineering Reference
In-Depth Information
paste(
'
hazard.death(
'
,int1[-length(int1)],
'
,
'
,int1[-1],
'
]
'
,sep=
''
),
"beta1","beta2","beta3","theta")
colnames(result)=c("Estimate","SE","95%CI(L",",U)","Z","P")
cat("Normalfrailtymodel","\n")
print(result)
cat("Fulllog-likelihood=",logFlik(est),"n=",N,"\n")
#print("relativeriskoftumorondeath");
print((1+theta)*exp(beta1-beta2))
}
#EndofProgram
###############
if(plots==TRUE)
{
par(mar=c(4,4,4,4)+0.1)
par(mfrow=c(2,2))
plot(conv[1:(w-1),(K2+1)],type='l',
'
l
'
,
ylab=
'
Estimateofbeta1
'
,
xlab=
'
Thenumberofiterations
'
)
plot(conv[1:(w-1),(K2+2)],type='l',
'
l
'
,
ylab=
'
Estimateofbeta2
'
,
xlab=
'
Thenumberofiterations
'
)
plot(conv[1:(w-1),(K2+3)],type='l',
'
l
'
,
ylab=
'
Estimateofbeta3
'
,
xlab=
'
Thenumberofiterations
'
)
plot(conv[1:(w-1),(K2+4)],type='l',
'
l
'
,
ylab=
'
Estimateoftheta
'
,
xlab=
'
Thenumberofiterations
'
)
}
}#EndofFunction
model_g<-function(data,int,init,tol,plots)
{
##############################################################
#Fulllog-likelihood#######################################
Search WWH ::
Custom Search