The name is absent

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)