问题标签 [kernel-density]

尝试绘制我的数据时，我收到以下错误消息：

这个可重现的代码有效：

但是使用我的真实数据不会：

然而，这两个数据集的结构是相同的（虚拟数据 [dat] 和真实数据 [temp]）：

它们都存储在数据帧中，我知道的唯一区别是长度，它以andtemp$V1为界。-11

两个数据集的输出kdexy <- kde2d()不同。在示例数据中，这些Z部分填充了非常小的数字；在真实数据集中，每个点都填充了“NaN”。

在对包使用多重插补后，我试图产生类似于densityplot()from 的东西。这是一个可重现的示例：lattice packageggplot2mice

产生：

我想对输出有更多的控制（我也用它作为 ggplot 的学习练习）。所以，对于bmi变量，我尝试了这个：

产生这个：

所以它有几个问题：

1. 颜色不对。看来我控制颜色的尝试完全错误/被忽略了
2. 有不需要的水平线和垂直线
3. 我希望图例显示 Imputed 和 Observed 但我的代码给出了错误invalid argument to unary operator

此外，要完成一条线完成的工作似乎需要做很多工作densityplot(impute)- 所以我想知道我是否可能完全以错误的方式去做这件事？

编辑：我应该添加第四个问题，正如@ROLO 所指出的：

.4. 地块的范围似乎不正确。

我正在尝试在 ggplot2 中的 R 中创建一个简单的密度图。这是我的代码，效果很好。

问题是我无法将 x 轴调整为负数。

效果很好，但是

根本不起作用！它给了我这个错误：

if (zero_range(range)) { 中的错误：需要 TRUE/FALSE 的缺失值

请帮忙！

我有一个包含 224900 个观察值和 10 个变量的数据集，这些变量是不同泰勒级数反向转换到原始数据值的结果。我希望覆盖这 10 个变量中每一个变量的密度图，以显示泰勒级数反变换对数据估计的稳健性水平。而不是只有 10 条线，我认为应用颜色会很好，这样每个密度图都会贡献 10% 的灰度。如果数据仅与其中一个图相关，则灰色为 10%，两个图的暗度为 20% 时的两倍，直到所有密度图重叠的地方，即 100%。

我曾经melt得到一个 2249000 行长的数据框。共有三列，第一列是人员 ID，第二列是分组变量（variable），第三列是每日 kJ 摄入量值（value）。

我使用下面的代码来覆盖密度图，ggplot2但它对组使用不同的颜色。如何更改此代码以获得我的灰度？我希望所有 10 个组具有相同的颜色和颜色密度；该图的目的只是使用灰度直观地显示密度图上的重叠量。

一些测试数据可供希望提供帮助的人使用，使用 5 组而不是 10 组：

我认为答案可能与修改scale_colour_grey和/或有关，scale_manual​​但我对自己的理解还不够。

可能重复：
将密度曲线拟合到 R 中的直方图

x是一个 NAs 自由数值向量。

我跑：

有人建议我在调用 hist 时设置 prob =TRUE。如果你也能解释这一点，那就太好了。谢谢你。

我有一个以天为单位的年龄列表，我希望在密度图上以年为单位显示它们。

我做了这两种方法 - 将 x 轴上的标签更改为年，并将数据除以 365。这些方法给了我不同的密度估计：

第一个图是使用以下方法生成的：

y 轴上显示的密度范围为 0 到 0.0002

但是，当我这样做时（将年龄除以 365 得到年份 - 而不仅仅是像上面那样更改 x 标签）：

该图看起来相同，但密度范围从 0 到 0.08 我很难理解发生了什么 - 为什么两个图之间的密度不同？

我想将密度曲线添加到直方图和累积直方图，如下所示：

这是我能做到的：

我正在覆盖几个密度图以在 ggplot2 中创建一个图形。假设我有一系列年份的每个 12 月的数据点（在本例中为 2004-2012 年），我想绘制每个 12 月 + 年的密度函数并将它们叠加。我想突出显示特定月份的一条密度线，方法是在所有其他密度线都是实心的地方用虚线表示。我在下面有一个可重现的示例：

目前，我能够制作这个。不幸的是，右边的图例显示的是 08 年 12 月的蓝色实线，而不是虚线。如何使图例对应于仅 08 年 12 月的虚线密度线（例如，所有其他图例键都是实线）？我曾尝试使用 legend.key 但这仅适用于轮廓和背景，而不适用于图例中的实际对象。

我正在使用均值偏移，此过程计算数据集中每个点的收敛位置。我还可以计算两个不同点收敛的坐标之间的欧几里得距离，但我必须给出一个阈值，也就是说，如果（距离 < 阈值）那么这些点属于同一个簇，我可以合并它们。

如何找到用作阈值的正确值？
（我可以使用每个值，这取决于结果，但我需要最佳值）

我正在尝试获取 R 中股票价格对数的密度估计值。我知道我可以使用plot(density(x)). 但是，我实际上想要函数的值。

我正在尝试实现内核密度估计公式。这是我到目前为止所拥有的：

任何帮助将非常感激。一直在这几天！