以下是R命令,这些命令突出显示了使用无偏交叉验证选择高斯内核密度估计器的带宽的两个函数的输出之间的差异:
> library(MASS);
>
> x <- rnorm(10000)
>
> stats::bw.ucv(x)
[1] 0.09771488
>
> MASS::ucv(x)
[1] 0.3908595
>
> sessionInfo()
R version 3.3.1 (2016-06-21)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: CentOS release 6.6 (Final)
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] MASS_7.3-45
在这里,MASS::ucv
的估算值大于stats:bw.ucv
的估算值。
有趣的是,在R stats软件包的“作者”部分的“内核密度估计的带宽选择器”下,说明如下:
B。 D. Ripley,摘自MASS软件包的早期版本。