for (j in 1:length(a))
phi[a[j],,,lambdaIndex] = res$phi[j,,]
rho[,,,lambdaIndex] = res$rho
pi[,lambdaIndex] = res$pi
dimension = 0
for (j in 1:length(a))
phi[a[j],,,lambdaIndex] = res$phi[j,,]
rho[,,,lambdaIndex] = res$rho
pi[,lambdaIndex] = res$pi
dimension = 0