%% load data % Dependencies used: % colormap: cbrewer: https://de.mathworks.com/matlabcentral/fileexchange/34087-cbrewer-colorbrewer-schemes-for-matlab % brainspace: download here: https://brainspace.readthedocs.io/en/latest/pages/install.html % surfstat: https://mica-mni.github.io/surfstat/ % scientific colourmaps: https://www.fabiocrameri.ch/colourmaps/ % raincloudplots: https://github.com/RainCloudPlots/RainCloudPlots % Prepare the MPC matrices according to % https://github.com/MICA-MNI/micaopen/tree/master/MPC % Prepare the resting-state matrices like this: % dummy code of my approach: for load_fc = 1 HCP400_fc1 = zeros(length(ID),1200,400); HCP400_fc2 = HCP400_fc1; HCP400_fc3 = HCP400_fc1; HCP400_fc4 = HCP400_fc1; matric_HCP400 = zeros(1206,400,400); for i = 1:length(ID) i try mycifti1 = ft_read_cifti([ID '*rfMRI_REST1_RL_Atlas_MSMAll_hp2000_clean.dtseries.nii']) mycifti2 = ft_read_cifti([ID '*rfMRI_REST1_LR_Atlas_MSMAll_hp2000_clean.dtseries.nii']) mycifti3 = ft_read_cifti([ID '*rfMRI_REST2_RL_Atlas_MSMAll_hp2000_clean.dtseries.nii']) mycifti4 = ft_read_cifti([ID '*rfMRI_REST2_LR_Atlas_MSMAll_hp2000_clean.dtseries.nii']) ts1 = mycifti1.dtseries(mycifti1.brainstructure <= 2,:); ts2 = mycifti2.dtseries(mycifti2.brainstructure <= 2,:); ts3 = mycifti3.dtseries(mycifti3.brainstructure <= 2,:); ts4 = mycifti4.dtseries(mycifti4.brainstructure <= 2,:); HCP400_fc1(i,:,:) = labelmean(ts1',HCP400.cdata,'ignorewarning'); HCP400_fc2(i,:,:) = labelmean(ts2',HCP400.cdata,'ignorewarning'); HCP400_fc3(i,:,:) = labelmean(ts3',HCP400.cdata,'ignorewarning'); HCP400_fc4(i,:,:) = labelmean(ts4',HCP400.cdata,'ignorewarning'); fc400(i,:,:) = (corr(squeeze(HCP400_fc1(i,:,:))) + corr(squeeze(HCP400_fc2(i,:,:))) + corr(squeeze(HCP400_fc3(i,:,:))) + ... corr(squeeze(HCP400_fc4(i,:,:))))./4; %Fischer z transform fc400z(i,:,:) = 0.5*log((1+ fc400(i,:,:))./(1- fc400(i,:,:))); catch end end end % define the IDs to keep keep_mpc = find(squeeze(mean(MPC(:,1,1:400),3))>0); keep_fc = find(squeeze(mean(fc400z(:,1,1:400),3))>0); keep = intersect(keep_mpc,keep_fc); hr_fc(eye(400)==1)=0; % zeros in the diagonal hr_mpc(eye(400)==1)=0; % zeros in the diagonal % make heritability map: % use the pedigree structure of HCP in combination with solar-eclipse % scripts: https://www.nitrc.org/projects/se_linux