ndwmean(array,weights,directions) calculates the n-dimensional weighted mean of an array. The size of the weighting array and the optional input "directions" determine the dimensions in which the weighted mean is applied (see examples). Supports NaN values.

Examples: x = rand(3,3,3,2)
ndwmean(x,rand(1,2)) % 4th dimension
ndwmean(x,rand(3,2),[2 4]) % 2nd and 4th dimensions
ndwmean(x,rand(1,3,1,2)) % 2nd and 4th dimensions
ndwmean(x,rand(1,3)) % 3rd dimension
ndwmean(x,rand(3,3)) % 2nd and 3rd dimensions