Incorporate ComBat into ModelArray

[smeisler]

Hi,

I think it would be great if ComBat harmonization could be incorporated into ModelArray to account for inter-site variance. Please let me know if I can help in any way.

Thanks,
Steven