Biomedical Engineering Reference
In-Depth Information
source("CSD.r")
CSD(data=bladder,int=c(12,18),frailty=
'
gamma
'
,
init=c(0,0.1,1),tol=0.001,plots=T)
CSD(data=bladder,int=c(12,18),frailty=
'
normal
'
,
init=c(0,0.1,1),tol=0.001,plots=T)
#####<Example1>Rawdata:Lung###############################
#time
#tumor=1(notumor),0(tumor)
#death=1(death),0(sacrified)
#trt
lung=read.table(
'
lung.txt
'
,col.names=
c(
'
id
'
,
'
time
'
,
'
tumor
'
,
'
trt
'
,
'
death
'
))
lung$tumor=1-lung$tumor
###############################################
source("CSD.r")
CSD(data=lung,int=c(12,18),frailty=
'
gamma
'
,init=c(0,0.1,1),
tol=0.001,plots=T)
CSD(data=lung,int=c(12,18),frailty=
'
normal
'
,init=c(0,0.1,1),
tol=0.001,plots=T)
model_n<-function(data,int,init,tol,plots)
{
#############################################################
#Fulllog-likelihood#######################################
logFlik<-function(param)
{
tal01=param[1:K];tal02=param[(K+1):K2]
tbeta1=param[K2+1];tbeta2=param[K2+2];tbeta3=param[K2+3];tsigma2
=param[K2+4]
logL=sum(colSums(p)*log(tal01))+sum(colSums(dij*d)*log(tal02))-
sum(tal01*colSums(Tnt*exp(tbeta1*z)*expR)+
tal02*colSums(Tnt*exp(tbeta2*z)*expR)+
tal02*colSums(Tt*exp(tbeta3*z)*expR))+
sum((d+de)*R)+sum((de*tbeta1+(1-de)*d*tbeta2+de*d*tbeta3)*z)-
Search WWH ::
Custom Search