Non-negative low-rank approximations for multi-dimensional arrays on statistical manifold