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

}
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)*R))#beta1
h[2,2]=-sum(al02*colSums((z^2)*Tnt*exp(tbeta2*z)*R))#beta2
h[3,3]=-sum(al02*colSums((z^2)*Tt*exp(tbeta3*z)*R))#beta3
h
}
#N-Rfor"theta"#####################################################
gr_theta=function(ttheta)
sum(R-logR)/(ttheta)^2-N*(1/(ttheta)^2-log(ttheta)/(ttheta)^2-
digamma(1/ttheta)/(ttheta)^2)
he_theta=function(ttheta)
-2*sum(R-logR)/(ttheta)^3-N*(-3/(ttheta)^3+2*log(ttheta)/(ttheta)^3+
2*digamma(1/ttheta)/(ttheta)^3+trigamma(1/ttheta)/(ttheta)^4)
#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)*R))#beta1
h[1:K,K+1]=h[K+1,1:K]=-colSums(z*Tnt*exp(tbeta1*z)*R)#al01*beta1
h
}
#SEfor"al02,theta2,theta3"#######################################
he_al02beta=function(param)
{
tal02=param[1:K];tbeta2=param[K+1];tbeta3=param[K+2];
h=matrix(0,ncol=length(param),nrow=length(param))
h[1:K,1:K]=-diag(colSums(dij*d)/tal02^2)#al02