fof-visually-guided/elutils/+stats/withinstderr.m

function [mu,se]=withinstderr(M)

col_mean=nanmean(M,2);
row_mean=nanmean(M,1);
grand_mean=mean(M(:));
nM=bsxfun(@minus,M,col_mean)+grand_mean;
mu = nanmean(nM,1);
se=stats.nanstderr(nM,1);