getneff.m

来自「Non-parametric density estimation」· M 代码 · 共 11 行

function Neff = getNeff(dens)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% getNeff(P)
%  returns an estimate of Neff, the 'effective' number of points in the NPDE
%   (= sum(1/weights^2) )
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Copyright (C) 2003 Alexander Ihler; distributable under GPL -- see README.txt

  Neff = 1/sum(getWeights(dens).^2);