2830

Lomb-Scargle your way to RSFC parameter estimation in AFNI-FATCAT
Paul A Taylor1, Gang Chen2, Daniel R Glen2, Richard C Reynolds1, and Robert W Cox1

1NIMH, NIH, Bethesda, MD, United States, 2NIH, Bethesda, MD, United States

Synopsis

We propose a new tool in AFNI-FATCAT to estimate the above RSFC parameters even when time series are censored, using the Lomb-Scargle (L-S) periodogram. The L-S approach for estimating RSFC parameters is useful and generalizable for FMRI data, where censoring is nearly always performed during processing. The method shows minimal bias of parameter estimation, and also allows for the estimation of confidence intervals for the parameters.

Introduction

Estimating frequency spectra from time series is central to many aspects of FMRI analysis. Several resting state functional connectivity (RSFC) parameters have been defined in terms of the low frequency fluctuation (LFF). For example, ALFF (the amplitude of LFFs [1]) is the sum of LFF spectral amplitudes, and fALFF (the fractional ALFF [2]) is the ratio of ALFF to the sum over all spectral amplitudes. Additionally, RSFA (resting state BOLD fluctuation amplitude [3]) is the square-root of the LFF power spectral density (PSD), and fRSFA (fractional RSFA) is the RSFA normalized by the square-root of the unfiltered PSD’s sum. Traditionally, the Fast Fourier Transform (FFT) is applied to FMRI time series to estimate the above spectral parameters. However, censoring time points is also a common step in FMRI processing, removing large motion or outliers; this effectively produces a time series with a variable sampling rate, so the FFT can no longer be applied. Replacing points with mean time series values would bias the estimate toward increased low frequency estimates, and other interpolations can also bias the RSFC estimates with assumptions. Therefore, we propose a new tool in AFNI-FATCAT to estimate the above RSFC parameters even when time series are censored, using the Lomb-Scargle (L-S) periodogram.

Methods

The Lomb-Scargle (L-S) periodogram was developed to examine periodicity in non-uniformly sampled time series [4-5]. Here, we use the efficient form of L-S calculation from [6], whose results in the presence of no censoring exactly match those from FFT, so the L-S approach can be applied generally. We assume the time series have been de-meaned, with either no or non-uniform (i.e., approximately random) censoring. Amplitudes are calculated up to a maximum frequency (pseudo-Nyquist limit), with an appropriate value shown to be $$$f_{max} = 1/(2*TR)$$$, where TR is sampling rate of acquisition [6]; with no censoring, this upper frequency precisely matches with the FFT case (the Nyquist limit). In traditional (uncensored) calculations, the ALFF and RSFA parameters from a (zero mean) time series $$$x_n$$$ of $$$N$$$ time points can be expressed as a sum of spectral amplitudes $$$|X_k|$$$ with harmonics in the interval $$$k\in[A, B]$$$: $$RSFA=\sqrt{\frac{1}{(N-1)N}\sum_{k=A}^{B}|X_k|^2}\,,~~~~~~~~ALFF=\frac{1}{\sqrt{N}}\sum_{k=A}^{B}|X_k|\,.$$ When censoring, the forms and scaling of the RSFC quantities must be expressed so the parameters are unbiased with respect to the the number of censored time points. Here, we derive the quantities starting from Parseval’s theorem, a standard relation between $N$ time points $x_n$ (here, with zero mean) and their $L = N/2$ independent spectral amplitudes $|X_k|$:

$$\sum_{n=0}^{N-1}|x_n|^2=\frac{1}{N}\sum_{k=0}^{N-1}|X_k|^2\,.$$

For real-valued, zero-mean time series, symmetries in the spectrum lead to:

$$\sum_{n=0}^{N-1}|x_n|^2=\frac{1}{L}\sum_{k=1}^{L-1}|X_k|^2\,.$$

These relations can be expressed with reference to the definition of variance:

$$\sigma^2=\frac{1}{N-1}\sum_{n=0}^{N-1}|x_n|^2=\frac{1}{(N-1)L}\sum_{k=1}^{L-1}|X_k|^2\,.$$

RSFA is the standard deviation of the LFF-bandpassed time series $$$x’_n$$$ with frequency harmonics in the interval [A, B]: $$RSFA^2=\frac{1}{N-1}\sum_{n=0}^{N-1}|x’_n|^2=\frac{1}{(N-1)L}\sum_{k=A}^{B}|X_k|^2\,.$$

The RSFA value is essentially the $L_2$ norm of amplitudes, scaled by an additional factor of $$$\sqrt{N-1}$$$. ALFF, above, is the $$$L_1$$$ norm. The L-S periodogram approximately maintains this relation in the presence of censored data that has $M<N$ time points. $L$ frequencies are still calculated, estimating the PSD in the absence of censoring. From this L-S spectrum, both the number of censored time points $M$ and spectral magnitudes $$$L$$$ must be considered, for both RSFA and ALFF. We demonstrate that the following are robust estimators of ALFF, fALFF, RSFA and fRSFA:

$$RSFA=\sqrt{\frac{1}{(M-1)L}\sum_{k=A}^{B}|X_k|^2}\,,~~~~~fRSFA=\sqrt{\frac{\sum_{k=A}^{B}|X_k|^2}{\sum_{k=1}^{L-1}|X_k|^2}}\,,$$

$$ALFF=\frac{1}{\sqrt{(M-1)L}}\sum_{k=A}^{B}|X_k|\,,~~~~~fALFF=\frac{\sum_{k=A}^{B}|X_k|}{\sum_{k=1}^{L-1}|X_k|}\,.$$

We examine the estimators using RS-FMRI data (196 volumes) from the publicly available FATCAT Demo data set [8], which was processed with standard steps in AFNI [9]. Time series were taken from within the brain mask. Random censoring was performed with $$$C=5,10,20,40,60$$$ points, each with 500 repetitions. The mean and standard deviation values were calculated, and the distributions were compared.

Results

Distributions for each set and RSFC parameter are shown in Fig. 1. The distributions are unimodal and strongly peaked, peaked at approximately the same values. As expected, the distribution widens with increased censoring.

The mean and standard deviations for each level of censoring are shown in Fig. 2 for each RSFC parameter, along with the uncensored value for reference. Even at very high levels of censoring (up to nearly one third of the points), there does not appear to be significant bias in parameter estimate.

Conclusions

The L-S approach for estimating RSFC parameters is useful and generalizable for FMRI data, where censoring is nearly always performed during processing. The method shows minimal bias of parameter estimation, and also allows for the estimation of confidence intervals for the parameters.

Acknowledgements

No acknowledgement found.

References

1. Zang YF, He Y, Zhu CZ, Cao QJ, Sui MQ, Liang M, Tian LX, Jiang TZ, Wang YF. 2007. Altered baseline brain activity in children with ADHD revealed by resting-state functional MRI. Brain Dev 29:83–91.

2. Zuo, X.N., Di Martino, A., Kelly, C., Shehzad, Z.E., Gee, D.G., Klein, D.F., Castellanos, F.X., Biswal, B.B., Milham, M.P., 2010. The oscillating brain: complex and reliable. Neuroimage 49:1432–1445.

3. Kannurpatti SS, Biswal BB. 2008. Detection and scaling of task induced fMRI-BOLD response using resting state fluctuations. Neuroimage 40:1567–1574.

4. Lomb, N. R. 1976, Ap&SS, 39, 447

5. Scargle, J. D. 1982, ApJ, 263, 835

6. Press, W. H., & Rybicki, G. B. 1989, ApJ, 338, 277

7. Eyer, L., & Bartholdi, P. 1999, A&AS, 135, 1

8. Taylor PA, Gohel S, Di X, Walter M, Biswal BB. 2012. Functional covariance networks: obtaining resting-state networks from intersubject variability. Brain Connect. 2(4):203-17.

9.

Cox RW. 1996. Comput Biomed Res 29:162-173.

Figures

Distributions of RSFC parameters from a RS-FMRI time series with various amounts of (random) censoring: number of points censored = 5, 10, 20, 40, 60. In each case, 500 repetitions were made.

The mean and standard deviation of RSFC parameters from the examples in Fig. 1. Results for each number of (randomly) censored points are shown. For reference, the RSFC parameter value for the uncensored time series is shown in bold.

Proc. Intl. Soc. Mag. Reson. Med. 26 (2018)
2830