estimatevars.m

The missing values are left out when estimating variances.

N = size(x, dim);
x(find(missing)) = probdist(0, 0);
xval = sum(x.e .^ 2 + x.var, dim);
N = N - sum(missing, dim);