In this paper, outage probability (OP) of a joint transmit and receive antenna selection (JTRAS) scheme is analyzed in multiple-input multiple-output non orthogonal multiple access based downlink energy harvesting (EH) relaying networks. In this dual-hop and amplify-and-forward relaying based network, since the first and second hops are types of single-user and multi-user systems, respectively, the optimal JTRAS and suboptimal majority-based JTRAS schemes are employed in the first and second hops. The theoretical OP analysis is carried out over Nakagami-m fading channels in the cases of perfect and imperfect successive interference cancellation. Finally, Monte Carlo simulations are performed to substantiate the accuracy of the theoretical analysis. It is shown that the optimal power splitting ratios at the EH relay are different for users and the users with good channel conditions have minimum optimal ratios.