Massive multiple-input multiple-output (MIMO) and full-duplex (FD) are promising candidates for achieving the spectral efficiency to meet the needs of 5G communications. One essential key to realizing practical FD massive MIMO systems is how to effectively mitigate the self-interference (SI). Conventionally, however, the performance comparison of different SI methods by reflecting the actual channel characteristics was insufficient in the literature. Accordingly, this paper presents a performance analysis of SI cancellation (SIC) methods in FD massive MIMO systems. Analytical and numerical results confirm that, in an imperfect channel-estimation case, the ergodic rates performance of the spatial suppression in the uplink outperforms those of the SI subtraction, due to the correlation between the precoder and the estimation error of the SI channel. In addition, we discuss which method performs better under different given system constraints such as uplink and downlink sum rates, the total transmit power, and the power scaling law.