- S[j,mm,r] = -rho[mm,mm,r]*ps2[j,mm,r] + phi[1:j-1,mm,r]%*%t(Gram2[j,1:j-1,r]) + phi[j+1:p,mm,r]%*%t(Gram2[j,j+1:p,r])
- if(abs(S(j,mm,r)) <= n*lambda*(Pi[r]^gamma)){
+ j1 = j-1
+ j2 = j+1
+ v1 = c(1:j1)
+ v2 = c(j2:p)
+ S[j,mm,r] = -rho[mm,mm,r]*ps2[j,mm,r] + phi[v1,mm,r]%*%(Gram2[j,v1,r]) + phi[v2,mm,r]%*%(Gram2[j,v2,r]) #erreur indice
+ if(abs(S[j,mm,r]) <= n*lambda*(Pi[r]^gamma)){