Biomedical Engineering Reference
In-Depth Information
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
Search WWH ::
Custom Search