I downloaded group average functional correlation
file: HCP_S900_820_rfMRI_MSMAll_groupPCA_d4500ROW_zcorr.dconn.nii
Some diagonal elements of the square matrix (91282x91282) are infinite
(Please see below).
I want to use this matrix in ananalysis; however, I am not sure how to
understand or deal with infinite diagonal values.
I appreciate any insight
Thanks,
Sang
======================
In [1]: import nibabel
In [2]: asdf =
nibabel.load('HCP_S900_820_rfMRI_MSMAll_groupPCA_d4500ROW_zcorr.dconn.nii')
In [3]: img = asdf.get_data()
In [4]: img.shape
Out[4]: (1, 1, 1, 1, 91282, 91282)
In [5]: S = img[0,0,0,0,:,:]
In [6]: S
Out[6]:
memmap([[ 8.66434002e+00, 1.96847185e-01, 1.66294336e-01, ...,
1.01449557e-01, 7.45474100e-02, 1.15624115e-01],
[ 1.96847185e-01, inf, 3.36383432e-01, ...,
-5.70017472e-03, -5.49946353e-02, 3.72834280e-02],
[ 1.66294336e-01, 3.36383432e-01, inf, ...,
-4.45242636e-02, -6.07097335e-02, -1.51601573e-02],
...,
[ 1.01449557e-01, -5.70017472e-03, -4.45242636e-02, ...,
inf, 1.91883039e+00, 9.20160294e-01],
[ 7.45474100e-02, -5.49946353e-02, -6.07097335e-02, ...,
1.91883111e+00, 8.31776619e+00, 8.82132888e-01],
[ 1.15624115e-01, 3.72833721e-02, -1.51601573e-02, ...,
9.20160294e-01, 8.82132888e-01, 8.66434002e+00]],
dtype=float32)
In [7]: S.diagonal()
Out[7]:
memmap([ 8.66434002, inf, inf, ..., inf,
8.31776619, 8.66434002], dtype=float32)
_______________________________________________
HCP-Users mailing list
[email protected]
http://lists.humanconnectome.org/mailman/listinfo/hcp-users
