fn_parcorrModelFC.m

function[rho_all, p_all, rho_intra, p_intra, rho_inter, p_inter] = ...
    fn_parcorrModelFC(FC, modelFC, SC, cfg)

LH_parc = cfg.LH_parc; 
RH_parc = cfg.RH_parc; 
corrType = cfg.corrType;

% Extract and vectorize SC elements
[SC_vec,SC_intra_vec,SC_inter_vec] = fn_interIntraConnVec(SC, LH_parc, RH_parc);

% Extract and vectorize FC elements
[FC_vec,FC_intra_vec,FC_inter_vec] = fn_interIntraConnVec(FC, LH_parc, RH_parc);

% Extract and vectorize model FC elements
[modelFC_vec,modelFC_intra_vec,modelFC_inter_vec] = fn_interIntraConnVec(modelFC, LH_parc, RH_parc);

% Correlate between actual and model FC, controlling for SC
[rho_all, p_all] = ...
    partialcorr(modelFC_vec,FC_vec,SC_vec,'type',corrType,'rows','pairwise');

[rho_intra, p_intra] = ...
    partialcorr(modelFC_intra_vec,FC_intra_vec,SC_intra_vec,'type',corrType,'rows','pairwise');

[rho_inter, p_inter] = ...
    partialcorr(modelFC_inter_vec,FC_inter_vec,SC_inter_vec,'type',corrType,'rows','pairwise');

end