have an issue here:
I have 8 different datasets, combined in a list EE_list.
TFP_EE <- list() # create list for estimation results Omegas <