lyman.glm.estimate_residual_autocorrelation

lyman.glm.estimate_residual_autocorrelation(Y, X, tukey_m=None)

Fit OLS model and estimate residual autocorrelation with regularization.

Parameters

Yn_tp x n_vox array

Array of time series data for multiple voxels.

Xn_tp x n_ev array

Design matrix for the model.

tukey_m: int or None

Size of tukey taper window or None to use default rule.

Returns

acftukey_m x n_vox array

Regularized autocorrelation function estimate for each voxel.