Discovery of optimal factors in binary data via a novel method of matrix decomposition