ISSN : 1229-067X
When using a multilevel model to analyze experimental data, the results of hypothesis testing can vary depending on the random effects structure. Previous studies have shown that the likelihood ratio test (LRT) for selecting the random effect structure does not match the shape of the null distribution with the theoretical chi-squared distribution due to the parameter space constraints of the variance. In this study, we investigate the shape, the type 1 error rate, and convergence failure rate of null distribution of LRT to select random effects structures through Monte Carlo simulations. Before proceeding with this study, the validity of the methodology was secured by checking the shape of null distribution of LRT for testing fixed effects. We conducted simulations by manipulating the experimental design, number of participants, and number of items, and found that the mean, variance, 95th and 99th percentiles of null distribution of LRT for selecting random effects structures were smaller than the theoretical chi-squared distribution under all conditions. When LRT was performed using the theoretical chi-squared distribution, the type 1 error rate was smaller than .05 significance level. The convergence failure occurred frequently in models with complex experimental designs and correlation structure in the structure of the random effects. The caveats of choosing the model’s random effects structure were discussed for empirical researchers when using multilevel models to test for experimental condition effects.
As many real datasets do not satisfy the normality assumption, the need for simulation studies related to nonnormality has arisen. To conduct such simulation studies, it is important to generate appropriate non-normal data. Although many non-normal data generation methods have been developed, there is a lack of studies that discuss the entire process of non-normal data generation. Therefore, this study offers guidelines that researchers can consider for conducting simulation studies by outlining the process of non-normal data generation with specified steps. Specifically, the non-normal data generation process will be discussed based on the power method, which is the most widely used technique for generating non-normal data, and the program and functions that researchers can utilize will also be examined. Furthermore, by providing a real example of non-normal data generation based on a two-factor CFA model along with the program code, we expect that this study will enhance understanding of the data generation process while providing tools that researchers can use.