pf_stats.m

％ ％ set some variables in the workspace to control behaviour: ％ ％ graphicsMode 0 no graphics, ％

function [Mean, Var] = pf_stats(pf)

	z = [pf.x.s]';
	prob = [pf.x.pi]';

	mn = sum( z .* [prob prob] );
	s2 = sum( (z - ones(pf.n,1)*mn)'.^2)';
	var = sum( s2 .* prob );

	if nargout == 0,
		fprintf('mean = %f %f, var = %f\n', mn, var);
	else
		Mean = mn;
		Var = var;
	end