debug test.EMGLLF
[valse.git] / R / indicesSelection.R
... / ...
CommitLineData
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
9indicesSelection = 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}