Dear Group,
since I got no reply from the last post, I try to explain the problem properly.
regarding to FLICA, I have three different modalities (VBM (FSL) and two modalities from FS) I'm trying to calculate the variance explained by each component including ALL modalities. I found this for each single modality
"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')"
--> but it just gives me the variance explained for each modality in a specific component (which is fine, too), but how can I calculate the variance explained by the whole component (network) including all modalities? (Like in Douaud et al, 2014, PNAS)
One approach might be a file called "weights" in matlab (1 x "number of components" - matrice). I am wondering if I can take these values for my desired calculation.. but I am really not sure, if this is the right way doing it.
I'm happy for any input. Thanks a lot :)
Matthias
Thanks for your help
Matthias
|