- }
- EMGLLf = EMGLLF(phiInit[a,,],rhoInit,piInit,gamInit,mini,maxi,gamma,0,X[,a],Y,tau)
-
- phiLambda = EMGLLf$phi
- rhoLambda = EMGLLf$rho
- piLambda = EMGLLf$Pi
-
- for(j in 1:length(a)){
- phi[a[j],,,lambdaIndex] = phiLambda[j,,]
- }
- rho[,,,lambdaIndex] = rhoLambda
- Pi[,lambdaIndex] = piLambda
-
+
+ res = EMGLLF(phiInit[a,,],rhoInit,piInit,gamInit,mini,maxi,gamma,0.,X[,a],Y,tau)
+
+ for (j in 1:length(a))
+ phi[a[j],,,lambdaIndex] = res$phi[j,,]
+ rho[,,,lambdaIndex] = res$rho
+ pi[,lambdaIndex] = res$pi
+