Internal function: Detrend Y by column-wise centering
detrend(Y, is_Y_detrended)
Y
Data matrix
Detrended data matrix