先检验下样本的分布情况，
如果样本分布服从正态分布的话，那么就用正态核来估计吧，hopt= 1.06*sigma*n^(-1/5)，当然也可以选择 ?density 里的那些选择 nrd, bcv, ucv 等等
sinon, 那就可以用CV(h): validation croisée 来估计hopt; 可以安装程序包 sm，用hcv 就OK了， 还可以用 Plug-in方法

楼主你好，请问这个问题你是怎么解决的

tvReg包里面可以根据数据进行计算

在R语言中，运用density()命令计算核密度估计，其中选项：bw给定窗宽，缺省时为拇指法则；kernel指定核函数，可以选择"gaussian" 、 "rectangular" 、 "triangular" 、 "epanechnikov" 、 "biweight" 、 "cosine" 或 "optcosine" 之一，默认为 "gaussian"。bw = "nrd0"为Silverman(1986)中的拇指法则，"nrd"为Scott(1992)给出的更常见的拇指法则，使用因子1.06。 "ucv"和"bcv"分别实现无偏和有偏交叉验证。"SJ"是 Sheather and Jones(1991)的方法，使用插入法选择带宽。