Statistics for UN torus link invariants in the large N limit from the matrix model approach