Ill
denomc < — matrix(NA,n,n)
denomc.term < — rep(0,n)
Ioglik < — О
der.l < — matrix(O,m,l)
der.2 < — matrix(O,m,m)
Ioglikc < — О
der.lc < — matrix(O,m,l)
der.2c < — matrix(O,m,m)
for (j in l:n){ #64 j-time
Ioglik.i < — O
loglikc.i < — O
der.l.i < — matrix(O,m,l)
der.2.i < — matrix(O,m,m)
der.lc.i < — matrix(O,m,l)
der.2c.i < — matrix(O,m,m)
zye < — matrix(O,m,l)
zzye < — matrix(O,m,m)
zyec < — matrix(O,m,l)
zzyec < — matrix(O,m,m)
for (k in l:n){ #65.1
denom.term[j] < — denom.term∣j] + weight[k,j]