HpM implement(1)

If you need to index this paper, please show the track adress! All rights reserved by qianjiang!

function InitialMode=Course_Sampling(xi, htheta, n,KernelType)
% find location with max density estimate
xi = sort(xi)          % ----------------- Added 23:50,1/3/2006
k=1:10;
j=k*round(n*1.0/11);
switch(KernelType)
case 'Epanechnikov', density_xi = DensityEstimate(@Epanechnikov,xi,xi(j),htheta);;  
case 'Biweight', density_xi = DensityEstimate(@Biweight,xi,xi(j),htheta);
otherwise, error(sprintf('Missing Density Estimate for %s function.',wfun));
end
disp('density=');
disp(density_xi);
[value,j0] = max(density_xi)
if length(j0) > 1
    j0 = j0(1);
end
InitialMode=xi(j(j0));    % ----------------- Modified by qianjiang at 23:50,1/3/2006, the oringinal one is
                          % ----------------- InitialMode=xi(j0), meanwhile
                          % ----------------- there was no order xi = sort(xi).