Two-step monotone missing data is a special case of missing data that attracted the attention of researchers since it not only appears frequently in applications but also since under the hypothesis of multivariate normality a unique and explicit solution to the likelihood equations exists. This importance prompted Yamada et al. [Kurtosis tests for multivariate normality with monotone incomplete data. Test. 2015;24:532–557. doi: 10.1007/s11749-014-0423-1] and Kurita and Seo [Multivariate normality test based on kurtosis with two-step monotone missing data. J Multivar Anal. 2022;188:104–824. doi: 10.1016/j.jmva.2021.104824] to propose four procedures in total for testing multivariate normality with two-step monotone missing data. The aim of this paper is twofold. On the one, to compare, for the first time, the performance of the four aforementioned multivariate normality tests via a Monte Carlo study. The second aim of the paper is to compare the performance of the four existing testing procedures with the respective procedures derived as a combination of 12 popular methods for handling incomplete data and 6 multivariate normality tests for complete data. We show that in several cases the tests proposed by Yamada et al. [Kurtosis tests for multivariate normality with monotone incomplete data. Test. 2015;24:532–557. doi: 10.1007/s11749-014-0423-1] and Kurita and Seo [Multivariate normality test based on kurtosis with two-step monotone missing data. J Multivar Anal. 2022;188:104–824. doi: 10.1016/j.jmva.2021.104824] are unable to achieve the nominal significance level. At the same time, no single procedure was found to be the most powerful in all situations of multivariate alternatives considered. As a consequence, the use of specific combinations of methods of handling missing data and testing multivariate normality with complete data is recommended.
Read full abstract