Wow! great response Alfonso --

There are all sorts of issues I ran into trying to implement this on my own (memory issues nonwithstanding) so I have been using existing implementations. I've managed to estimate MDL for each individual subject's denoised data using the icatb_estimate_dimension.m function in the GIFT toolbox.  However I'm not too confident in the results since the MDL estimate hasn't been performed using the group level data but rather individual subjects (which all turned out to be 56...).  I've attached a plot of the MDL/AIC estimate with shaded error bars across subjects.  Does this look right to you?

Thank you so much for your continued support!!