105
cum.hazd < — matrix(NA,l,n);
cum.hazdc < — matrix(NA,l,n);
item < — matrix(NA,l,n)
itemc < — matrix(NA,l,n)
item.den < — matrix (0,1,n)
surv.fnc.s < — matrix (NA,l,n)
surv.fnc.r < — matrix (NA,l,n)
for (i in l:n) { #11 i represents t
if (deltaln[i]!=0){ #13
for (k in i:n){#14
item.den[i] < — item.den[i] + exp(betaa.fitem[i] < — deltaln[i]∕(item.den[i]);
} else { #13 & 15
item[i] < — 0 } #15
if (delta2n[i]!=0){ #13
for (k in i:n){#14
item.denc[i] < — item.denc[i] + exp(betac.f% * %cow[,k]) } #14
itemc[i] < — delta2n[i]∕(item.denc[i])
} else { #13 & 15
itemc[i] < — 0} }#11
for (j in l:n) {#16
cum.hazd∣j] < — cumsum(item)[j] ;