Hi Shahrzad,
It used to be the case that this was automatically calculated, however this was unfortunately lost on the last versions of FLICA.
Once you've got the k number of outputs (or modalities, in our PNAS paper we had m=3), each with n number of components, and you are in the relevant FLICA directory, you can simply do this for each k in [1,m]:
fslmaths niftiOut_mik.nii.gz -sqr /tmp/tmpk.nii.gz && fslstats -t
/tmp/tmpk.nii.gz -m >/tmp/energyk.txt
And in matlab:
load /tmp/energyk.txt
energyk = energyk * diag(1./sum(energyk))
Mk=[(0:n-1)', (energyk*100)]
dlmwrite('variancek.txt', Mk, 'delimiter','\t','precision','%.2f')
You'll see in variancek.txt that the second column sums up to 100% for that given k modality.
Hope this helps,Gwenaëlle
Thanks dear Gwenaëlle,
Yes, you are right, I see it now. and thanks for the components.
I have another question:
How does one estimate how much of the variance in the data is being explained by each component?
From the outputs flica gives, I guess I should be using the weights of each component over sum of weight of all components multiplied by 100. Am I right?
Thanks a lot again.
Shahrzad
Hi Shahrzad,
The simple answers are:
I haven't added the first age-related map, as it is fairly unspecific and whole-brain, but I'm happy to add it on the webpage if you feel this would be helpful!
Cheers,
Gwenaëlle
Dear FSL group,
- In Flica program, is it possible to choose the number of independent components automatically?
My other question is with regard to the spatial maps of the components. Specifically, the recent Douaud et al. (2014, PNAS) paper presents 2 main components that are correlating significantly with age.
Are these spatial maps freely available ? (like the spatial maps provided by smith et al for resting state fmri components)
Shahrzad