PROTECT(S = allocArray(REALSXP, dimPhiS));
PROTECT(affec = allocVector(INTSXP, n));
double *pPhi=REAL(phi), *pRho=REAL(rho), *pPi=REAL(pi), *pLLF=REAL(LLF), *pS=REAL(S);
PROTECT(S = allocArray(REALSXP, dimPhiS));
PROTECT(affec = allocVector(INTSXP, n));
double *pPhi=REAL(phi), *pRho=REAL(rho), *pPi=REAL(pi), *pLLF=REAL(LLF), *pS=REAL(S);