Setup light cone geometry
This function makes weights of a mixture model more sparse using gradient based penalty methods.
1 2  sparsify_weights(weight.matrix.proposed, weight.matrix.current = NULL,
penalty = "entropy", lambda = 0)

weight.matrix.proposed 
N \times K weight matrix 
weight.matrix.current 
N \times K weight matrix 
penalty 
type of penalty: 
lambda 
penalization parameter: larger

compute_mixture_penalty
,
mixed_LICORS
1 2 3 4 5 6 7  WW < matrix(c(rexp(10, 1/10), runif(10)), ncol = 5, byrow = FALSE)
WW < normalize(WW)
WW_sparse < sparsify_weights(WW, lambda = 0.1)
WW_more_sparse < sparsify_weights(WW, lambda = 0.5)
compute_mixture_penalty(WW)
compute_mixture_penalty(WW_sparse)
compute_mixture_penalty(WW_more_sparse)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.