quartic.m

来自「计量工具箱」· M 代码 · 共 25 行

function heights=quartic(data,center,bw);
% PURPOSE:
%     Quartic kernel function
% 
% USAGE:
%     heights=triangular(data,center,bw)
% 
% INPUTS:
%     data: Data used for height calculation
%     center: The location for height calculation
%     bw:  The bandwidth
% 
% OUTPUTS:
%     heights: Scalar height for kernel density
% 
% COMMENTS:
%     For kernel density estimation
% 
% Author: Kevin Sheppard
% [e-mail]
% Revision: 2    Date: 12/31/2001
u=(data-center)/bw;
heights=(15/16)*((1-u.^2).^2).*(abs(u)<1);