X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=R%2Fselectionindice.R;h=97014b1ff419e228b7f232157e3cc283c6c72683;hb=6bd2e869a17f3980d52820643c1c1b5f3725738e;hp=6782b2eef9f3fa391ad56f7d7e85446a601ecba6;hpb=493a35bfea6d1210c94ced8fbfe3e572f0389ea5;p=valse.git diff --git a/R/selectionindice.R b/R/selectionindice.R index 6782b2e..97014b1 100644 --- a/R/selectionindice.R +++ b/R/selectionindice.R @@ -1,24 +1,28 @@ -selectionindice = function(phi, seuil){ - dim_phi = dim(phi) - pp = dim_phi[1] - m = dim_phi[2] - - A = matrix(0, pp, m) - B = matrix(0, pp, m) - - for(j in 1:pp){ - cpt1 = 0 - cpt2 = 0 - for(mm in 1:m){ - if(max(phi[j,mm,])> seuil){ - cpt1 = cpt1 + 1 - A[j,cpt] = mm - } - else{ - cpt2 = cpt2+1 - B[j, cpt2] = mm - } - } - } - return(list(A,B)) -} \ No newline at end of file +selectionindice = function(phi, seuil) +{ + dim_phi = dim(phi) + pp = dim_phi[1] + m = dim_phi[2] + + A = matrix(0, pp, m) + B = matrix(0, pp, m) + + for(j in 1:pp) + { + cpt1 = 0 + cpt2 = 0 + for(mm in 1:m) + { + if(max(phi[j,mm,]) > seuil) + { + cpt1 = cpt1 + 1 + A[j,cpt] = mm + } else + { + cpt2 = cpt2+1 + B[j, cpt2] = mm + } + } + } + return (list(A,B)) +}