next up previous
Next: Data set (D) Up: Simulated data Previous: Data set (B)

5.1.3 Data set (C)

In data set (C), different parameters for the HRF convolutions resulted in small differences in the time courses. As a consequence, each of the three spatial maps has a slightly different associated time-course in each subject. As such, there is no longer a single 'true' associated time course across subjects. The induced variation in the temporal domain, however, is small (time courses have $r>0.75$ temporal correlation) and the different time courses are well approximated by their dominant Eigenvector, i.e. by the best rank-1 approximation. The data permits two different representations: firstly, the signal content can either be approximated as a linear combination of 3 processes (where in the temporal domain the rank-1 approximation to the 3 slightly different time courses is used), or can fully be expressed as a linear combination of 9 processes with large co-linearity in $\mbox{\protect\boldmath$A$}$ and 3 multiple versions for each of the true spatial sources in $\mbox{\protect\boldmath$B$}$.

Figure 4[*] shows the estimated set of source processes for tensor-PICA and PARAFAC. The tensor-PICA decomposition, due to the independence assumption in the spatial domain, represents the data via a set of 3 source processes. Compared to the tensor-PICA results, the PARAFAC spatial estimates exhibit some cross-talk, e.g. the first spatial map is visibly confounded by map 2 and map 3. Also, in the temporal- and subject domains, PARAFAC finds less accurate estimates of the true source processes. The two approaches differ most significantly in the way in which true spatial maps correlate with each of the $R=14$ estimated maps: while the spatial tensor-PICA decomposition always results in only one source which correlates strongly with the true spatial map, the PARAFAC decomposition shows that, especially for sources 2 and 3, multiple PARAFAC estimates correlate with the true maps. As such, PARAFAC does not represent the signal of interest via 3 different source processes but equally does not find the representation by 9 sources: almost all of the estimated correlated maps show significant amount of cross-talk. Convergence in this case is particularly slow, with 47 times the number of floating point operations compared to tensor-PICA.

Figure 4: Tensor-PICA and PARAFAC decomposition results for data set (C) with $R=14$. Each of the 'true' spatial maps no longer has a single 'true' signal time course associated and estimated time courses (blue) are shown together with the dominant Eigenvector of all possible time courses (red).
\includegraphics[width=\figwidth]{fig4}
(i) tensor-PICA (ii) PARAFAC


next up previous
Next: Data set (D) Up: Simulated data Previous: Data set (B)
Christian Beckmann 2004-12-14