Biomedical Engineering Reference
In-Depth Information
logFlik<-function(param)
{
tal01=param[1:K];tal02=param[(K+1):K2]
tbeta1=param[K2+1];tbeta2=param[K2+2];tbeta3=
param[K2+3];ttheta=param[K2+4]
logL=sum(colSums(p)*log(tal01))+sum(colSums(dij*d)*log(tal02))-
sum(tal01*colSums(Tnt*exp(tbeta1*z)*R)+tal02*
colSums(Tnt*exp(tbeta2*z)*R)+tal02*colSums(Tt*exp(tbeta3*z)*R))+
sum((d+de)*logR)+sum((de*tbeta1+(1-de)*d*tbeta2+de*d*tbeta3)*z)+
sum((1/ttheta-1)*logR-R/ttheta)-N*(log(ttheta)/
ttheta+lgamma(1/ttheta))logL
}
#Marginallog-likelihood##########################################
logMlik<-function(param)
{
tal01=param[1:K];tal02=param[(K+1):K2]
tbeta1=param[K2+1];tbeta2=param[K2+2];tbeta3
=param[K2+3];ttheta=param[K2+4]
marL=sum(colSums(p)*log(tal01)+colSums(dij*d)*log(tal02))-
sum((1/ttheta+de+d)*log(1+ttheta*(rowSums((Tnt*exp(tbeta1*z))%*%al01)+
rowSums((Tnt*exp(tbeta2*z))%*%al02)+rowSums((Tt*exp(tbeta3*z))%*%al02))))
+sum(de*d*(1+ttheta))+sum((de*tbeta1+(1-de)*d*tbeta2+de*d*tbeta3)*z)
marL
}
#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)*R))#beta1
g[2]=sum((1-de)*d*z)-sum(al02*colSums(z*Tnt*exp(tbeta2*z)*R))#beta2
g[3]=sum(de*d*z)-sum(al02*colSums(z*Tt*exp(tbeta3*z)*R))#beta3
g
Search WWH ::
Custom Search