samer@4: function x=meanlastdims(x,n)
samer@4: % meanlastdims - compute average over last n dimensions of multidim array
samer@4: %
samer@4: % Usage: y=meanlastdims(x,n)
samer@4: %
samer@4: % If size1(x) = [m(1) m(2) ... m(k)], then size1(y) = [m(1) ... m(k-n)]
samer@4: % MEANLASTDIMS ignores any trailing singleton dimensions including 2nd
samer@4: % dimension of a column vector.
samer@4: 
samer@4: sz=size1(x);
samer@4: for d=length(sz)-(0:n-1);
samer@4: 	x=sum(x,d)/sz(d);
samer@4: end