R/indicesSelection.R

   1 #' Construct the set of relevant indices -> ED: je crois que cette fonction n'est pas utile
   2 #'
   3 #' @param phi regression matrix, of size p*m
   4 #' @param thresh threshold to say a cofficient is equal to zero
   5 #'
   6 #' @return a list with A, a matrix with relevant indices (size = p*m) and B, a
   7 #'                  matrix with irrelevant indices (size = p*m)
   8 #' @export
   9 indicesSelection = function(phi, thresh = 1e-6)
  10 {
  11     dim_phi = dim(phi)
  12     p = dim_phi[1]
  13     m = dim_phi[2]
  14
  15     A = matrix(0, p, m)
  16     B = matrix(0, p, m)
  17
  18     for(j in 1:p)
  19     {
  20         cpt1 = 0
  21         cpt2 = 0
  22         for(mm in 1:m)
  23         {
  24             if(max(phi[j,mm,]) > thresh)
  25             {
  26                 cpt1 = cpt1 + 1
  27                 A[j,cpt] = mm
  28             } else
  29             {
  30                 cpt2 = cpt2+1
  31                 B[j, cpt2] = mm
  32             }
  33         }
  34     }
  35     return (list(A=A,B=B))
  36 }