IlO
return(surv.fnc.f) } #10
В.6 Calculate the likelihood function as shown in
Equation 3.12
likel.h < — function(covv,deltain,delta2n,betaa.f,betac.f,rww,rwwc,PP,PPc){#61
betaa.f < — matrix(betaa.f,l,m)
betac.f < — matrix(betac.f,l,m)
zbetaa < — rep(NA,n)
zbetacc < — rep (N A,n)
for (z in l:n){
zbetaa[z] < — betaa.f% * %covv[,z]
zbetacc[z] < — betac.f% * %covv[,z] }
weight < — matrix(NA,n,n)
weightc < — matrix(NA,n,n)
for ( i in l:n){ #62
for ( j in l:n) { #63
weight[i,j] < — (exp(zbetaa[i]))*rww[i,j];
weightc[i,j] < — (exp(zbetacc[i]))*rwwc[i,j] } } #62
denom < — matrix(NA,n,n)
denom.term < — rep(0,n)