Print

Print


Hi pals,

I have a question about the way partial correlations were estimated for the
HCP megatrawl using FSLnets. Going through the 'ridgep' case in
nets_netmats.m, it looks like the covariance matrix for each subject is
scaled by the square root of the mean of the variances squared - yes, no,
maybe so?

from line 88 of nets_netmats.m:

grot = cov(<component timeseries>);
grot = grot/sqrt(mean(diag(grot).^2));

I'm just trying to figure out why I need to square the variances before I
take the mean. Can someone give me a quick explanation or point me towards
a good reference?

Thanks in advance for your help!
mb