我希望有人可以在这里提供帮助。我一直在使用R中的ks
软件包来开发我的模拟二维样本的多变量核密度估计值,$ W =(W_1,W_2)$,组件是iid。与之相关的代码是我使用的代码:
init.bw = Hpi.diag(W,nstage='2',pilot='amse',pre='scale') #Initial plugin bandwidth
cv.bw = Hscv.diag(W,nstage=2,pre='scale',Hstart=init.bw) #Cross - validation bandwidth estimator
hat.fw.mult = kde(W,H=cv.bw,eval.points=W)
问题在于,当我试图绘制这个估算器时,我已经遵循了这个估算,它看起来像是:
plot(hat.fw.mult, display="slice", cont=c(25,50,75), abs.cont,
approx.cont=FALSE, xlab, ylab, zlab="Density function", add=FALSE,
drawpoints=FALSE, drawlabels=TRUE, theta=-30, phi=40, d=4,
ptcol="blue", col="black")
问题是我遇到以下错误:
Error in plotkde.2d(fhat, ...) : Need a grid of density estimates
我不知道如何解决这个问题..任何帮助都将不胜感激!