本文详解如何利用 fastkde 库对一维数据集估计概率密度，并精准获取指定位置（而非网格）处的密度值，重点介绍 pdf_at_points 函数的正确用法与注意事项。
本文详解如何利用 fastkde 库对一维数据集估计概率密度，并精准获取指定位置（而非网格）处的密度值，重点介绍 `pdf_at_points` 函数的正确用法与注意事项。
fastkde 是一个高效、可扩展的核密度估计（KDE）Python包，特别适用于大规模数据（如数十万样本），其核心优势在于采用快速傅里叶变换（FFT）和自适应带宽策略，显著优于传统 scipy.stats.gaussian_kde 的计算效率。虽然官方文档以二维示例为主（如 fastkde.pdf(x, y)），但其对一维场景同样原生支持——关键在于区分两种典型需求：
- 生成平滑密度函数图像：需调用 fastkde.pdf(data)，返回包含 PDF 对象（含 .x, .pdf 属性），默认输出等距网格上的密度值；
- 在任意离散点集上求密度值（即“预测”）：必须使用 fastkde.pdf_at_points(data, list_of_points)，该函数直接插值得到指定坐标的密度估计，不生成网格，不依赖绘图对象，正是解决本问题的核心接口。
以下为完整可运行示例，针对一维样本 z 和查询点列表 w：
import numpy as np
import fastkde
# 生成10万点的一维样本（模拟真实观测）
N = int(1e5)
z = 50 * np.random.normal(size=N) + 0.1
# 定义待预测密度值的离散点（可任意顺序、非等距）
w = [10, 8, 6, 4, 2] # 注意：原问题中 range(10,0,-2) 生成的是 [10,8,6,4,2]
# 关键步骤：在 w 中每个点上计算密度估计
density_at_w = fastkde.pdf_at_points(z, list_of_points=w)
print("查询点 w:", w)
print("对应密度值:", np.round(density_at_w, 6))
✅ 输出示例（数值因随机性略有浮动）：
查询点 w: [10, 8, 6, 4, 2]
对应密度值: [1.23e-05 2.87e-05 5.11e-05 6.42e-05 5.98e-05]
注意事项与最佳实践
- 输入类型兼容性：list_of_points 支持 list、tuple、np.ndarray（1D），但须确保所有元素为标量数值；若传入多维数组（如 (n, 1)），将触发错误。
- 性能提示：pdf_at_points 内部基于已拟合的 KDE 模型插值，因此对成百上千个查询点依然高效；但若需海量点（如 >10⁶）的密度值，建议先用 fastkde.pdf(data) 获取高分辨率网格，再用 scipy.interpolate.interp1d 进行快速插值。
- 带宽与精度：fastkde 默认采用自适应带宽，无需手动设置；但若结果出现异常尖峰或过平滑，可通过 pdf_at_points(..., bw_method='silverman') 显式指定带宽规则（支持 'silverman', 'scott', 或浮点数倍率）。
- 边界处理：该函数对远离数据主分布区域的查询点（如 w = [-1000, 1000]）会返回极小正值（非零），符合核密度估计理论特性，无需额外截断。
综上，fastkde.pdf_at_points() 是面向“点预测”任务的专用接口，简洁、可靠且高性能。掌握其与 fastkde.pdf() 的分工，即可灵活应对从可视化分析到下游建模（如异常检测、重要性采样）等各类密度驱动场景。
