Statistical Analysis of Dependent Current Status Data with Application to Tumorigenicity Experiments - Interval-Censored-Time-to-Event Data: Methods and Applications - page 131

(N/2)*log(2*pi*tsigma2)-sum(Rsq)/(2*tsigma2)
logL
}
#N-Rfor"beta1beta2beta3"###############################
gr_beta=function(param)
{
tbeta1=param[1];tbeta2=param[2];tbeta3=param[3];
g=matrix(0,ncol=length(param),nrow=1)
g[1]=sum(de*z)-sum(al01*colSums(z*Tnt*exp(tbeta1*z)*expR))#beta1
g[2]=sum((1-de)*d*z)-sum(al02*colSums(z*Tnt*exp(tbeta2*z)*expR))#beta2
g[3]=sum(de*d*z)-sum(al02*colSums(z*Tt*exp(tbeta3*z)*expR))#beta3
g
}
he_beta=function(param)
{
tbeta1=param[1];tbeta2=param[2];tbeta3=param[3];
h=matrix(0,ncol=length(param),nrow=length(param))
h[1,1]=-sum(al01*colSums((z^2)*Tnt*exp(tbeta1*z)*expR))#beta1
h[2,2]=-sum(al02*colSums((z^2)*Tnt*exp(tbeta2*z)*expR))#beta2
h[3,3]=-sum(al02*colSums((z^2)*Tt*exp(tbeta3*z)*expR))#beta3
h
}
#SEfor"al01,theta1"######################################
he_al01beta=function(param)
{
tal01=param[1:K];tbeta1=param[K+1];
h=matrix(0,ncol=length(param),nrow=length(param))
h[1:K,1:K]=-diag(colSums(p)/tal01^2)#al01
h[K+1,K+1]=-sum(tal01*colSums((z^2)*Tnt*exp(tbeta1*z)*expR))#beta1
h[1:K,K+1]=h[K+1,1:K]=-colSums(z*Tnt*exp(tbeta1*z)*expR)#al01*beta1
h
}
#SEfor"al02,theta2,theta3"#############################