R中ks中kde的带宽

时间:2018-12-06 17:06:31

标签: r kde

我正在使用R中的ks软件包中的kde函数进行三维内核密度估计:https://cran.r-project.org/web/packages/ks/ks.pdf。该函数接受作为带宽的参数“ H”-但是,我不明白给出该参数的格式。我对数学的理解建议带宽应该是单个值(实际上是内核函数的输出乘以该值)或三个值(即每个维度上的带宽不同)。但是,函数实际接受的值为3x3矩阵。使用同一个程序包中的Hucv或Hucv.diag之类的带宽选择函数,我可以生成kde函数随后愉快地接受的值,例如:

            [,1]        [,2]        [,3]
[1,]  6679417.38  72656228.3  69791.9709
[2,] 72656228.26 790328476.1 759145.3598
[3,]    69791.97    759145.4    816.4863

           [,1]    [,2]     [,3]
[1,] 3532545052 0.00000  0.00000
[2,]          0 4.58034  0.00000
[3,]          0 0.00000 33.42598

但是,我得到警告,这些带宽选择功能在诸如我的数据(有很多重复值)的情况下表现不佳-实际上,从输出来看,这些带宽显然在一维上太大而另外两个太小。但是由于我不了解这个矩阵的确切含义以及它在kde中的用法,因此无法有效地对其进行修改。包装的小插图和手册没有对此进行解释。

我的问题是:有人知道为什么这个R函数接受3x3矩阵作为其带宽参数吗?应该如何修改此矩​​阵以更改内核的输出?

0 个答案:

没有答案