如何获取stats :: density中特定值的密度估计？
戴维·帕萨雷蒂（Davide Passaretti）
假设我有如下数据：
val <- .65
set.seed(1)
distr <- replicate(1000, jitter(.5, amount = .2))
d <- density(distr)
由于stats::density使用的是特定值bw，因此不会在间隔中包含所有可能的值（因为它们是无限的）：
d$x[ d$x > .64 & d$x < .66 ]
[1] 0.6400439 0.6411318 0.6422197 0.6433076 0.6443955 0.6454834 0.6465713 0.6476592 0.6487471
[10] 0.6498350 0.6509229 0.6520108 0.6530987 0.6541866 0.6552745 0.6563624 0.6574503 0.6585382
[19] 0.6596261
我想找到一种提供val给密度函数的方法，以便它返回其d$y估计值（然后将其用于为密度图的区域着色）。
我无法猜测这个问题有多愚蠢，但我找不到快速解决方案。
我想通过d$y对应于的两个值d$x更接近的线性插值来获得它val。有没有更快的方法？
红外光谱
这说明了approxfun：
> Af <- approxfun(d$x, d$y)
> Af(val)
[1] 2.348879
> plot(d(
+
> plot(d)
> points(val,Af(val) )
> png();plot(d); points(val,Af(val) ); dev.off()
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
编辑于
相关文章
Related 相关文章
- 1
如何整合核密度估计
- 2
如何在R中重用我的内核密度估计函数？
- 3
如何使用scikit归一化内核密度估计？
- 4
如何制作这些顺序的直方图/密度估计图
- 5
如何将内核密度估计作为scikit中的一维聚类方法来学习？
- 6
R中的密度估计
- 7
如何估计2D GMM及其梯度的核密度估计？
- 8
在density2d中计算密度估计值？
- 9
在numpy中创建密度估计
- 10
R中的高斯核密度估计
- 11
R中的核密度估计器
- 12
在numpy中创建密度估计
- 13
从 sklearn 核密度估计中采样
- 14
如何使用sklearn在2D图像/阵列上获得内核密度估计器？
- 15
将内核密度估计值提取到R中的新采样点
- 16
从 R 中的 3D 核密度估计和提取点值
- 17
Python中的多变量内核密度估计
- 18
Python中的多处理空间内核密度估计
- 19
Matlab中的knn（k最近邻）密度估计源
- 20
python中的加权高斯核密度估计
- 21
Matlab中的knn（k最近邻）密度估计源
- 22
如何从SciPy中的scipy.stats.distribution对象获取分发名称？
- 23
如何从冻结的Spicer.stats分布中获取参数自变量？
- 24
如何从SciPy中的scipy.stats.distribution对象获取分发名称？
- 25
sm.density.compare（）：在一个图中显示多个密度估计
- 26
如何将带有内核密度估计线的图表从 d3.js 版本 3 转换为 d3.js 版本 4
- 27
数据流的密度估计
- 28
使用Seaborn为python中的matplotlib进行内核密度估计的下限
- 29
一张d3.js图表中的多个内核密度估计
我来说两句