Abstract:
|
Two developments in fMRI magnitude time series modeling, namely, the incorporation of temporal dependence and the Ricean distribution, have been separated by a distributional “mismatch”: such time series modeling is largely based upon Gaussian-distributional-based extensions to the general linear model, which precludes its use under Ricean modeling. We bridge this gap by extending independent AR(p) errors to the latent, Gaussian-distributed real and imaginary components from which the Ricean-distributed magnitudes are computed by augmenting the observed magnitude data by missing phase data in an EM algorithm framework. We use the EM algorithm for parameter estimation and extend it to compute approximate standard errors and test statistics for activation and AR order detection. We compare the performance of this “AR(p) Ricean model” to the standard Gaussian AR(p) model for simulated and experimental fMRI data, focusing on cases where the signal-to-noise ratio is low.
|