split original function in 2 functions: 1 for creating ITS datasets and 1 for comparing OLS and MLM estimation