In this paper, the antenna selection problem is investigated for the relaying system with the base station having massive multiple antennas, where the relay station works in the full duplex mode and the two-way protocol. The antenna selection scheme is optimized by maximizing the minimum of the two signal to noise ratios(SNRs) with each corresponding to the uplink and the downlink traffics. With the optimized antenna selection scheme, the approximate probability density function(PDF) of the received SNR for each receiver is derived. In what follows, both the overall outage probability and the bit error rate are obtained in the analytical expression for the whole massive multi-input multi-output(MIMO) system. Moreover, the asymptotic overall outage probability is derived in an analytical expression with respect to the growing number of the antennas on the massive MIMO base station. Numerical simulations verify the derived analytical results. In this paper, the antenna selection problem is investigated for the relaying system with the base station having massive multiple antennas, where the relay station works in the full duplex mode and the two-way protocol. The antenna selection scheme is optimized by maximizing the minimum of the two signal to noise ratios (SNRs) with each corresponding to the uplink and the downlink traffics. With the optimized antenna selection scheme, the approximate probability density function (PDF) of the received SNR for each receiver is derived. both the overall outage probability and the bit error rate are obtained in the analytical scale for the whole massive multi-input multi-output (MIMO) system. Moreover, the asymptotic overall outage probability is derived in an analytical expression with respect to the growing number of the antennas on the massive MIMO base station. Numerical simulations verify the derived analytical results.