fre = [1 2 3; 3 4 3];
orien = "r";
refM = [1 1 1];
-refV = [0.889522663062 0.593015108708 0.593015108708];
+refV = [0.200138037385 0.334558519179 0.125031231272];
[v, m] = variancef( x, fre, orien, ones(meanf(x,fre,orien)) );
assert_checkalmostequal([v m], [refV refM]);
orien = "c";
refM = [1; 1];
-refV = [0.288922678414; 0.019966608736];
+refV = [0.492838922640; 0.060105711640];
[v, m] = variancef( x, fre, orien, ones(meanf(x,fre,orien)) );
assert_checkalmostequal([v m], [refV refM]);