A consistent and numerically efficient variable selection method for sparse Poisson regression with applications to learning and signal recovery