[valse.git] / R / indicesSelection.R

#' Construct the set of relevant indices -> ED: je crois que cette fonction n'est pas utile
#'
#' @param phi regression matrix, of size p*m
#' @param thresh threshold to say a cofficient is equal to zero
#'
#' @return a list with A, a matrix with relevant indices (size = p*m) and B, a 
#'					matrix with irrelevant indices (size = p*m)
#' @export
indicesSelection = function(phi, thresh = 1e-6)
{
	dim_phi = dim(phi)
	p = dim_phi[1]
	m = dim_phi[2]

	A = matrix(0, p, m)
	B = matrix(0, p, m)

	for(j in 1:p)
	{
		cpt1 = 0
		cpt2 = 0
		for(mm in 1:m)
		{
			if(max(phi[j,mm,]) > thresh)
			{
				cpt1 = cpt1 + 1
				A[j,cpt] = mm
			} else
			{
				cpt2 = cpt2+1
				B[j, cpt2] = mm
			}
		}
	}
	return (list(A=A,B=B))
}
Commit	Line	Data
	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	}