llh = rbind(llh,
c( res$LLF, sum(Rank[j,] * (length(active)- Rank[j,] + m)) ) )
phi[active,,,] = rbind(phi[active,,,], res$phi)
llh = rbind(llh,
c( res$LLF, sum(Rank[j,] * (length(active)- Rank[j,] + m)) ) )
phi[active,,,] = rbind(phi[active,,,], res$phi)