- for(r in 1:k){
- p1 = p-1
- for(j in 1:p1){
- for(mm in 1:m){
- 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)){
+ for (r in 1:k)
+ {
+ for (j in 1:p)
+ {
+ for (mm in 1:m)
+ {
+ S[j,mm,r] = -rho[mm,mm,r]*ps2[j,mm,r] +
+ (if(j>1) sum(phi[1:(j-1),mm,r] * Gram2[j,1:(j-1),r]) else 0) +
+ (if(j<p) sum(phi[(j+1):p,mm,r] * Gram2[j,(j+1):p,r]) else 0)
+ if (abs(S[j,mm,r]) <= n*lambda*(pi[r]^gamma))