function make_dice_coeff_calculation

project_dir = pwd;

cd([project_dir '/supplementary_material/dice_coeff_numk_does_not_matter/dice_coefficient'])

func_file_1 = 'R_frontal_pref_data_hierarchy_F26_P15_8_branches.func.gii';
func_file_2 = 'R_frontal_minS_data_hierarchy_F40_P15_8_branches.func.gii';
func_file_3 = 'R_frontal_medianS_data_hierarchy_F67_P15_8_branches.func.gii';

data_1 = double(nonzeros(readimgfile(func_file_1)));
data_2 = double(nonzeros(readimgfile(func_file_2)));
data_3 = double(nonzeros(readimgfile(func_file_3)));

% Calculate dice 
jac_1_2             = jaccard(data_1,data_2);
similarity_1_2      = 2 * jac_1_2 ./ (1 + jac_1_2);
average_dice_1_2    = mean(similarity_1_2);

jac_1_3             = jaccard(data_1,data_3);
similarity_1_3      = 2 * jac_1_3 ./ (1 + jac_1_3);
average_dice_1_3    = mean(similarity_1_3);

jac_2_3             = jaccard(data_2,data_3);
similarity_2_3      = 2 * jac_2_3 ./ (1 + jac_2_3);
average_dice_2_3    = mean(similarity_2_3);

cd([project_dir '/supplementary_material/dice_coeff_numk_does_not_matter/dice_coefficient'])
save dice_coeff

cd(project_dir)

end