- 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:p)
- {
- b = A2[j,2:dim(A2)[2],lambdaIndex]
- b = b[b!=0]
- if (length(b) > 0)
- phi[A2[j,1,lambdaIndex],b,,lambdaIndex] = 0.
- c = A1[j,2:dim(A1)[2],lambdaIndex]
- dimension = dimension + sum(c!=0)
- }